Hi all,
I've noticed that a number of people have had this problem, but can't seem
to find any resolutions.
I am running a PGLS analysis, with body size as the dependent variable, and
two categorical traits as the explanatory variables:
> data.pgls <- pgls(MeanBodyMass ~ Habitat_aq * TrophicLevel, species.comp,
> lambda='ML', delta='ML')
>
The pgls runs without any problems, and summary(data.pgls) gives proper
results.
However, when I run anova(data.pgls), I get the following error:
> Error in terms.formula(formula, data = data) :
> invalid model formula in ExtractVars
My trophic level categorical explanatory variable is multistate (not
binary), so the anova results are much preferred to the summary results
(which splits multistate variables into a number of binary ones).
Has anyone run into this problem before and has found a fix or workaround?
Thanks once again in advance,
Will
--
William Gearty
PhD Student, Paleobiology
Department of Geological and Environmental Sciences
Stanford School of Earth, Energy & Environmental Sciences
people.stanford.edu/wgearty
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
