Anqi Qiu wrote: > I am not sure whether there is a bug. When I tested the example given for > "glht" > in the help, I entered the following error: > > Running commands: > > amod <- aov(minutes ~ blanket, data = recovery) > rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), > alternative = "less") > > Errors are: > > Error in try(coef.(model)) : could not find function "coef." > > Error in modelparm.default(model, ...) : no 'coef' method for 'model' > found! > > Thanks!
It works for me. > library(multcomp) Loading required package: mvtnorm > amod <- aov(minutes ~ blanket, data = recovery) > rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), + alternative = "less") > summary(rht) Simultaneous Tests for General Linear Hypotheses Multiple Comparisons of Means: Dunnett Contrasts Fit: aov(formula = minutes ~ blanket, data = recovery) Linear Hypotheses: Estimate Std. Error t value p value b1 - b0 >= 0 -2.1333 1.6038 -1.330 0.2412 b2 - b0 >= 0 -7.4667 1.6038 -4.656 <0.001 *** b3 - b0 >= 0 -1.6667 0.8848 -1.884 0.0924 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Adjusted p values reported) > sessionInfo() R version 2.6.1 RC (2007-11-20 r43507) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] multcomp_0.992-6 mvtnorm_0.8-1 loaded via a namespace (and not attached): [1] tools_2.6.1 > Regards, > > Anqi Qiu > > Anqi Qiu > Center for Imaging Science > Johns Hopkins University > Email: [EMAIL PROTECTED] > tel: 410-516-8103 > > Anqi Qiu > Center for Imaging Science > Johns Hopkins University > Email: [EMAIL PROTECTED] > tel: 410-516-8103 > ______________________________________________ > 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. -- Chuck Cleland, Ph.D. NDRI, Inc. 71 West 23rd Street, 8th floor New York, NY 10010 tel: (212) 845-4495 (Tu, Th) tel: (732) 512-0171 (M, W, F) fax: (917) 438-0894 ______________________________________________ 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.