search for: g_h08_lme

...t 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)~as.factor(Diversity), data=ME, random=~1|Plot/SubPlot, weights=varPower...