[R] multcomp

Nair, Murlidharan T mnair at iusb.edu
Wed Mar 21 02:35:05 CET 2007


Thanks for the hint. 
Here is what I had done earlier according to the old package
 
 mult.comp<-simint(a~b, data=z, conf.level=0.99, type=c("Tukey") )
 isoforms<-as.vector(rownames(mult.comp$estimate))
 estimate<-as.vector(mult.comp$estimate)
 lower<-as.vector(mult.comp$conf.int[,1])
 upper<-as.vector(mult.comp$conf.int[,2])
 p.val.raw<-as.vector(mult.comp$p.value.raw)
 
Here is how I modified the above so that I can use the new package. Can someone comment if what I have done is correct or not?
amod<-aov(a~b, data=z)
mult.comp<-glht(amod, linfct=mcp(b="Tukey"))
estimate<-confint(mult.comp)[1,1]
lower<-confint(mult.comp)[1,2]
upper<-confint(mult.comp)[1,3]
summary(mult.comp) does give the adjusted pvalues.  I tried all the subscripts to extract it to store in a variable but could not, can any one help me here. Also,  how do I specify the conf.level to be equal to 0.99. The default is 0.95
 
Thanks for your help,
Cheers../Murli

 
________________________________

From: talepanda [mailto:talepanda at gmail.com]
Sent: Sun 3/18/2007 11:07 PM
To: Nair, Murlidharan T
Cc: r-help at stat.math.ethz.ch
Subject: Re: [R] multcomp



?glht says
>with 'print', 'summary',  'confint', 'coef' and 'vcov'  methods
>being available.

try:

example(glht)
summary(glht(amod, linfct = mcp(tension = "Tukey")))
confint(glht(amod, linfct = mcp(tension = "Tukey")))

On 3/19/07, Nair, Murlidharan T <mnair at iusb.edu> wrote:
> I used the multcomp package sometime back for doing multiple
> comparisons. I see that it has been updated and the methods like simint
> are no longer supported. When I run the program it prompts to me to use
> glht. How do I get the lower and upper conf int and the pValues using
> glht? Does anyone have an example?
>
> Thanks ../Murli
>
>
>
>
>       [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>



More information about the R-help mailing list