Re: [R] Multiple comparison on lme model with 2 fixed factors

2009-08-18 Thread Billeke
Sorry to come back again, but with this and other (real) data, I'm not able to obtain those comparisons. I'm using the following command and obtain the following result (lmer model). glht(lmelac.ex, linfct = mcp("ath1:tim20 - ath1:tim40 = 0"="Tukey"))->comp Erreur dans mcp2matrix(model, linfct =

Re: [R] Multiple comparison on lme model with 2 fixed factors

2009-08-18 Thread Billeke
Hi again, Reproducing part of my data (which indeed are here well standardized), I obtain similar results (with regard to the R output) than previously. These are "part" of an experimentation where people with different athletic (ath) training exercised during 120 minutes. Plasma lactate was one

Re: [R] Multiple comparison on lme model with 2 fixed factors

2009-08-18 Thread Dieter Menne
Billeke wrote: > > Reproducing part of my data (which indeed are here well standardized), I > obtain similar results (with regard to the R output) than previously. > > These are "part" of an experimentation where people with different > athletic (ath) training exercised during 120 minutes. Pla

Re: [R] Multiple comparison on lme model with 2 fixed factors

2009-08-18 Thread Dieter Menne
Billeke wrote: > > > mult2 <- glht(lmelac.ex, linfct = mcp("ath:time" = "Tukey")) > > I don't have the time to check, but got a mail for Torsten (author) with the following idea: glht(lmeglu, linfct = c("Ath1:tim2 - Ath1:tim3 = 0", ...)) He suggested to build the matrix explicitly instead

Re: [R] Multiple comparison on lme model with 2 fixed factors

2009-08-17 Thread Dieter Menne
Billeke wrote: > > I'm a bit lost while performing multiple comparisons on a lme model of > that type: > > lmeglu=lme(glucose~Ath*tim,random=~1|Vol,na.action=na.omit,data=data) > multc = glht(lmeglu, linfct = mcp(Ath = "Tukey", tim = "Tukey")) > > This works fine for identifying the effect of