array chip <arrayprofile <at> yahoo.com> writes: (Real names are prefered here)
> I am trying to test some contrasts, using glht() in > multcomp package on fixed effects in a linear mixed > model fitted with lme() in nlme package. The command I > used is: > data=dat) > glht(dat.lme, linfct = mcp(type =c("b+t-2*m=0"))) > > The lme model fit is ok, but I got an error message > with glht(): > Error in eval(expr, envir, enclos) : object > "info.index" not found > Error in factor_contrasts(model) : no 'model.matrix' > method for 'model' found! > > according to help page of glht(), it should work with > linear mixed model, what is the problem here? Here an example from Torsten Hothorn (author of multcomp), both for lme and lmer : To get a better answer, please provide the (simulated?) data for your example as required by the posting guide, # multcomplme.r nlmeOK <- require("nlme") lme4OK <- require("lme4") library("multcomp") K <- rbind(c(0,1,-1,0),c(0,1,0,-1),c(0,0,1,-1)) data("ergoStool", package = "nlme") stool.lmer <- lmer(effort ~ Type + (1 | Subject), data = ergoStool) glme4 <- glht(stool.lmer,K) summary(glme4) #Linear Hypotheses: # Estimate Std. Error z value p value #1 == 0 1.6667 0.5187 3.213 0.00376 ** #2 == 0 3.2222 0.5187 6.212 < 0.001 *** #3 == 0 1.5556 0.5187 2.999 0.00761 ** stool.lme <- lme(effort ~ Type, data = ergoStool, random = ~ 1 | Subject) gnlme <- glht(stool.lme,K) summary(gnlme) ______________________________________________ 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.