Hi Dr. Therneau,

Thanks very much for your kind help! Does survfit() just give me the
survival curve? What I wanted is the baseline hazard estimates, i.e.,
lambda_{0} (t). How can I obtain this estimate from coxph()? Or using
basehaz()?

Thanks,
YH










On Mon, Nov 4, 2013 at 5:29 PM, Terry Therneau <therneau@mayo.edu> wrote:

>
>
> On 11/03/2013 05:00 AM, r-help-request@r-project.org wrote:
>
>> I can easily get the other parameter estimates by using coxph() but don't
>> know how to get the baseline hazard estimates from it. Does anyone know
>> how
>> to obtain nonparametric baseline hazard estimates in the gamma frailty
>> model?
>>
>> Thanks,
>> YH
>>
>>
> I don't see what the problem is.
>
>  fit1 <- coxph(Surv(time, status) ~ age + ph.ecog + frailty(inst),
> data=lung)
>  sfit <- survfit(fit1)
>  plot(sfit)
>
>
> Please give an actual example of the problem.
>
> Terry Therneau
>

	[[alternative HTML version deleted]]

