[R] only need the p-value

From: Boks, M.P.M. <M.P.M.Boks_at_umcutrecht.nl>
Date: Sat 30 Sep 2006 - 13:09:15 GMT

 

Dear R users,  

I am calculating several cox proportional hazard models after each other (I know this is unusual, but I am just exploring the data). For the purpose of multiple testing correction I need to construct an array of these p-values. However since the output is not an array in itself, I cannot find a way to obtain the p-value only.  

> attach(tms)

>

> goal<-rep(0.7*FREQUENC[1:13],6)

> event<- Surv(TIJD,FREQUENC<goal)

> results<-coxph(event~ TYPETREA)

> summary(results)
 

Call:

coxph(formula = event ~ TYPETREA)  

  n=76 (2 observations deleted due to missing)

                     coef exp(coef) se(coef)     z     p

TYPETREAnon-guided -0.826     0.438    0.484 -1.71 0.088

 

                   exp(coef) exp(-coef) lower .95 upper .95

TYPETREAnon-guided     0.438       2.28     0.169      1.13

 

Rsquare= 0.041 (max possible= 0.829 )

Likelihood ratio test= 3.2 on 1 df, p=0.0737

Wald test = 2.91 on 1 df, p=0.088

Score (logrank) test = 3.08 on 1 df, p=0.0794  

Does anyone now how to extract the p-value?

Many thanks!,

Marco



R-help@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. Received on Sat Sep 30 23:16:37 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 30 Sep 2006 - 14:31:16 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.