Lilith, >> No the big mystery is the Tukey test. I just can't find the mistake, it >> keeps telling me, that >> there are " less than two groups" >> ... >> ### Tukey test ## >> summary(glht(PAM.lme, linfct = mcp(Provenancef = "Tukey"))) >> >> Error message: >> Fehler in glht.matrix(model = list(modelStruct = list(reStruct = >> list(Code .808654423456211, : >> ncol(linfct) is not equal to length(coef(model))
You need to look, in particular, at what's in your variable "Provenancef." ## PAMdata$Provenancef levels(PAMdata$Provenancef) And if the call to glht() is returning an object you need to inspect obj$linfct. It contains the contrast matrix (for Tukey contrasts) and needs to match the number of coefficients in your model. Try the following for further clues: glht(PAM.lme, linfct = mcp(Provenancef = "Tukey"))$linfct coef(PAMaov) Regards, Mark. -- View this message in context: http://r.789695.n4.nabble.com/Tukey-Test-lme-error-less-than-two-groups-tp3069789p3071678.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.
