[R] only need the p-value

[Boks, M.P.M.]

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