Re: [R] Help on extract paramters from fitted models

From: Marc Schwartz <marc_schwartz_at_comcast.net>
Date: Tue, 29 Apr 2008 14:43:51 -0500

Lisa wrote:
> Hi, I have a question about how to extract paramters from a fitted model. I
> can extract coefficients and std, but from some other statistics, I dont
> know how to extract. Can anyone help?
>
> Here it is an example:
>
>

>>  coxout<-coxph(Surv(t,t.censor)~x)

>
>
>> coxout

> Call:
> coxph(formula = Surv(t, t.censor) ~ x)
>
> coef exp(coef) se(coef) z p
> x 0.349 1.42 0.257 1.36 0.17
>
> Likelihood ratio test=1.84 on 1 df, p=0.175 n= 200
>
>> coxout$coef

> x
> 0.3490500
>
>> coxout$p

> NULL
>
> Any convienent way to extract the LRT and the p value? thank you!

Use the summary method for coxph objects:

   summary(coxout)$logtest["test"]

   summary(coxout)$logtest["pvalue"]

Unfortunately, this is not documented in the built-in help for the survival package, but you can see the structure for the summary.coxph method by using:

   str(summary(coxout))

?summary.coxph actually points to the structure details of a coxph object, rather than the summary.coxph object.

HTH, Marc Schwartz



R-help_at_r-project.org 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 Tue 29 Apr 2008 - 19:47:25 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 29 Apr 2008 - 20:30:44 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.

list of date sections of archive