Displaying 1 result from an estimated 1 matches for "ace_pse".
Did you mean:
ace_p
2013 Jul 25
1
lme (weights) and glht
...(summary(g_H08_lme))
--> using lme with weights I changed the order of factor() and
introduced as.factor in the model
H08_lme<-lme(log(Height2008_mean)~as.factor(Diversity), data=ME,
random=~1|Plot/SubPlot, weights=varPower(form=~Diversity),
na.action=na.omit, subset=ME$Species=="Ace_pse", method="ML")
summary(H08_lme)
anova(H08_lme)
ME$Diversity=factor(ME$Diversity)
g_H08_lme<-glht(H08_lme, linfct=mcp(Diversity="Tukey"))
Error in mcp2matrix(model, linfct = linfct) :
Variable(s) ‘Diversity’ have been specified in ‘linfct’ but cannot
be found in ‘mo...