From: Ritwik Sinha <ritwik.sinha_at_gmail.com>

Date: Sat 30 Sep 2006 - 13:38:48 GMT

Date: Sat 30 Sep 2006 - 13:38:48 GMT

This is how you go about doing this.

summary(results)$coefficients[1,5]

You will have to check this for you code. But the idea is that summary(results) is a list (?) and one of its components is called "coefficients", which is a matrix. So the problem is just to extract one element of this matrix.

I am not well versed with coxph so there may be some minor details I am missing, but that is the general idea (same as with lm, glm etc.).

Ritwik.

On 9/30/06, Boks, M.P.M. <M.P.M.Boks@umcutrecht.nl> wrote:

*>
**>
*

> 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.
**>
*

-- Ritwik Sinha Graduate Student Epidemiology and Biostatistics Case Western Reserve University ritwik.sinha_at_gmail.com | +12163682366 | http://darwin.cwru.edu/~rsinha ______________________________________________ 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:44:42 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:17 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.
*