Re: [R] analyzing suvival data using splines (a.k.a., piecewise log-hazard-ratio models)

From: Frank E Harrell Jr <f.harrell_at_vanderbilt.edu>
Date: Thu 23 Jun 2005 - 07:05:46 EST

Wittner, Ben wrote:
> I'm looking for software that makes plots such as fig 4 (a)-(e), fig 5 anf fig 7
> of
> Gray, Robert, "Flexible Methods for Analyzing Survival Data Using Splines, with
> Applications to Breast Cancer Prognosis," 1992, J Am Stat Assoc, pp 942-51.
>
> In other words, I'm looking for software that takes survival data and a
> continuous
> covariate as input and computes a curve giving log hazard ratio (or rate of
> failure for
> a specific time) as a function of the continuous covariate, as well as curves
> giving
> +-SE for that curve.

library(Design)
d <- datadist(mydata); options(datadist='d') f <- cph(Surv(dtime,death) ~ age + sex + rcs(blood.pressure,5)) plot(f, blood.pressure=NA)

Frank

>
> A plot of this nature can also be found as figure 4 in Paik, et al., "A
> Multigene Assay to
> Predict Recurrence of Tamoxifen-Treated, Node-Negative Breast Cancer," 2004, New
> England
> Journal of Medicine, pp 2817-26.
>
> Any help would be greatly appreciated.
>
> Thanks.
>
> -Ben
>
> Ben Wittner
> Research fellow, MGH & Harvard Medical School
> wittner.ben@mgh.harvard.edu
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>

-- 
Frank E Harrell Jr   Professor and Chair           School of Medicine
                      Department of Biostatistics   Vanderbilt University

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Thu Jun 23 07:07:36 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:57 EST