Re: [R] a question about glht function

2011-07-27 Thread Lao Meng
Here's another question: What's kind of difference does glht tests? In other words, If I get 4 intercepts and 4 slopes for all levels of GROUP,then glht tests for the difference of the 4 intercepts or 4 slopes,or something else? Thanks. 2011/7/27 Lao Meng > Yes. > According to your suggestio

Re: [R] a question about glht function

2011-07-26 Thread Lao Meng
Yes. According to your suggestion,I modified my code.It works well. Thanks you very much. My best. 2011/7/26 Peter Ehlers > On 2011-07-26 00:16, Lao Meng wrote: > >> Hi all: >> There's a question about glht function. >> >> My data:data_ori,which inclue CD4, GROUP,time. >> >> >> f_GROUP<-fac

Re: [R] a question about glht function

2011-07-26 Thread Achim Zeileis
On Tue, 26 Jul 2011, Peter Ehlers wrote: On 2011-07-26 01:13, Achim Zeileis wrote: On Tue, 26 Jul 2011, Lao Meng wrote: Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3)

Re: [R] a question about glht function

2011-07-26 Thread Peter Ehlers
On 2011-07-26 01:13, Achim Zeileis wrote: On Tue, 26 Jul 2011, Lao Meng wrote: Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result<- lme(sqrt(CD4) ~ f_GROUP*time ,r

Re: [R] a question about glht function

2011-07-26 Thread Achim Zeileis
On Tue, 26 Jul 2011, Lao Meng wrote: Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result <- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori) glht(resu

Re: [R] a question about glht function

2011-07-26 Thread Mark Difford
On Jul 26, 2011 Lao Meng wrote: > glht(result, linfct = mcp(f_GROUP="Tukey") ) > Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected It is almost certainly the underscore in the name ("_") that is causing the problem. Try putting the term in quotes ("f_GROUP"). Regards, Mar

Re: [R] a question about glht function

2011-07-26 Thread Peter Ehlers
On 2011-07-26 00:16, Lao Meng wrote: Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result<- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori) glht(result,

[R] a question about glht function

2011-07-26 Thread Lao Meng
Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result <- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori) glht(result, linfct = mcp(f_GROUP="Tukey") ) Error