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 <thern...@mayo.edu> wrote: > > > On 11/03/2013 05:00 AM, r-help-requ...@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]] ______________________________________________ R-help@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.