R help - Mar 2003 - help--Cox ph model

Dear r-users,

I want to use the Cox's ph model to analyze survival data set.  How can I
extract the model coefs. and Wald test p-value or Score?

For example: I use the data set melanom in iSwR package.
> library(survival)
> data(melanom)
> attach(melanom)
> cox.model <- coxph(Surv(days,status==1)~sex)
> summary(cox.model)
Call:
coxph(formula = Surv(days, status == 1) ~ sex)

  n= 205 

     coef exp(coef) se(coef)    z     p
sex 0.662      1.94    0.265 2.50 0.013

    exp(coef) exp(-coef) lower .95 upper .95
sex      1.94      0.516      1.15      3.26

Rsquare= 0.03   (max possible= 0.937 )
Likelihood ratio test= 6.15  on 1 df,   p=0.0131
Wald test            = 6.24  on 1 df,   p=0.0125
Score (logrank) test = 6.47  on 1 df,   p=0.0110


Thanks.

Guangchun

"Song, Guangchun" <Guangchun.Song at stjude.org> writes:
> Dear r-users,
> 
> I want to use the Cox's ph model to analyze survival data set.  How can
I
> extract the model coefs. and Wald test p-value or Score?
Unlike most other modelling code, summary.coxph() just prints its
results rather than store them in an object that can be queried. So
basically, you need to copy code from summary.coxph. It's not all that
hard. 

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907