Hi, I am new to R and I am trying to perform a post hoc tukey test using the multcomp package's ghlt for a lme model. My first attempt at doing so gave me an output, HOWEVER I have tried to do this again, it keeps coming up with the error: (Error in contrMat(table(mf[[nm]]), type = types[pm]) : less than two groups)
My model is looking at effect of incubation temperature (3 groups) on PHA with hen as a my random factor. Temperature treatment has been made factorial (Tempfact). my code is as follows: library(multcomp) > PHA <- lme(DiffPHA48 ~ Tempfact, random= ~1|Hen, na.action=na.omit) > summary(glht(PHA, linfct=mcp(Tempfact = "Tukey"))) Error in contrMat(table(mf[[nm]]), type = types[pm]) : less than two groups I am confused why it would work the first time and now everyother time it comes up with this error. Can anyone please help me?? Thanks! -- View this message in context: http://r.789695.n4.nabble.com/Trouble-with-performing-post-hoc-analysis-Tukey-for-lme-model-using-ghlt-tp3849982p3849982.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.