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