[R] How to use restricted cubic spline in survfit.cph function in survival package?

2014-03-06 Thread Zhiyuan Sun
Hello, I used PBC data set included in the survival package to fit a cox model. The model included a restricted cubic spline transformation on age. Then I tried using survfit function to predict a survival curve using the first row of the original data. I got an error message. The R code and

Re: [R] How to use restricted cubic spline in survfit.cph function in survival package?

2014-03-06 Thread Andrews, Chris
Try cph in rms (where rcs is defined). library(rms) fit-cph(Surv(time,status==2)~rcs(age,4)+sex, data=pbc, y=TRUE, x=TRUE) id1-pbc[1,] surv.id1-survfit(fit,newdata=id1) plot(surv.id1) summary(surv.id1) -Original Message- From: Zhiyuan Sun [mailto:sam.d@gmail.com] Sent: Wednesday,