[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)
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
More information about the R-help
mailing list