Thanks for Terry to reply for my qustion.
That comes to my second question, is there a simple way to plot the Hazard
ratio from a coxph?
anyone familiar with the muhaz package?
thanks a bunch
Maggie
On Tuesday, September 27, 2005, at 06:24AM, Terry Therneau <therneau@mayo.edu>
wrote:
> The simplest way that I know to get an estimate of the
> hazard is to use smoothing splines.
>
> fit <- survfit(Surv(time, status) ~ sex, data=lung)
>
> temp1 <- smooth.spline(fit[1]$time, 1-fit[1]$surv, df=5)
> temp2 <- smooth.spline(fit[2]$time, 1-fit[2]$surv, df=5)
>
> plot( predict(temp1, deriv=1), type='l')
> lines(predict(temp2, deriv=1), col=2)
>
>
>------
> The lung cancer data set is part of Splus (it's used in several
> manual pages as an example). The plot shows the males to have a
> higher initial hazard than the females.
>
> You choice of 5 degrees of freedom for the spline was completely
>arbitrary.
>
> Terry Therneau
>
>--------------------------------------------------------------------
>This message was distributed by s-news@lists.biostat.wustl.edu. To
>unsubscribe send e-mail to s-news-request@lists.biostat.wustl.edu with
>the BODY of the message: unsubscribe s-news
>
>
|