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 <laomen...@gmail.com> > Yes. > According to your suggestion,I modified my code.It works well. > > Thanks you very much. > > > My best. > > > > 2011/7/26 Peter Ehlers <ehl...@ucalgary.ca> > >> 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, linfct = mcp(f_GROUP="Tukey") ) >>> Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected >>> >>> >>> I can't find out the reason for Error. >>> >> >> I think glht() is looking for f_GROUP in data_ori. >> You have defined f_GROUP in your global environment but >> have not added it to your dataframe. >> >> Peter Ehlers >> >> >>> >>> Any help from you are welcome. >>> >>> >>> Many thanks! >>> >>> [[alternative HTML version deleted]] >>> >>> ______________________________**________________ >>> R-help@r-project.org mailing list >>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help> >>> PLEASE do read the posting guide http://www.R-project.org/** >>> posting-guide.html <http://www.R-project.org/posting-guide.html> >>> and provide commented, minimal, self-contained, reproducible code. >>> >> >> > [[alternative HTML version deleted]] ______________________________________________ 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.