Displaying 1 result from an estimated 1 matches for "height2005_mean".
2013 Jul 25
1
lme (weights) and glht
...rrors when using
weights=varPower().
The glht error makes sense as glht needs factor levels and the model
works fine without weights=. Does anyone know a solution so I do not
have to change the lme model?
Thanks
Sibylle
--> works fine
ME$Diversity=factor(ME$Diversity)
H08_lme<-lme(log(Height2005_mean)~Diversity, data=ME, random=~1|Plot/
SubPlot, na.action=na.omit, subset=ME$Species=="Pse_men", method="ML")
summary(H08_lme)
anova(H08_lme)
g_H08_lme<-glht(H08_lme, linfct=mcp(Diversity="Tukey"))
print(summary(g_H08_lme))
--> using lme with weights I changed th...