Wittner, Ben, Ph.D. wrote:
Does anyone know how to make a plot like Fig. 4 of Paik, et al., New England
Journal of Medicine, Dec. 30, 2004?
Given survival data and a covariate, they plot a curve giving "Rate of Distant
Recurrence at 10 Yr (% of patients)" on the y-axis versus the covariate on the
x-axis. They also plot curves giving a 95% confidence interval.
Thanks very much.
-Ben
Such a plot is easy to do with the rms package if using a Cox or
accelerated failure time model, e.g.
require(rms)
dd <- datadist(mydata); options(datadist='dd')
f <- cph(Surv(rtime, event) ~ rcs(covariate,4) + sex + ..., x=TRUE,
y=TRUE) # restricted cubic spline with 4 knots
plot(Predict(f, covariate, sex, time=10)) # separate curves for male
and female; omit sex to make one curve; add age=50 to predict for a 50
year old
--
Frank E Harrell Jr Professor and Chairman School of Medicine
Department of Biostatistics Vanderbilt University
______________________________________________
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.
