[R] only need the p-value

Boks, M.P.M. M.P.M.Boks at umcutrecht.nl
Sat Sep 30 15:09:15 CEST 2006


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)



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!,


More information about the R-help mailing list