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 > [EMAIL PROTECTED] > > ______________________________________________ > 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