Displaying 1 result from an estimated 1 matches for "height2008_mean".
2013 Jul 25
1
lme (weights) and glht
...ies=="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(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="T...