Displaying 1 result from an estimated 1 matches for "pse_men".
Did you mean:
ps_mem
2013 Jul 25
1
lme (weights) and glht
...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 the order of factor() and
introduced as.factor in the model
H08_lme<-lme(log(Height2008_mean)~a...