search for: ace_pse

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...