Well, I am very proud of myself. Here is something to help you.
The dataset is attached, with "anonymized" column names. The code to reproduce the error is below:
library(rms) d <- datadist(x) options(datadist="d") f <- cph(Surv(start,stop,status)~P1*P2 + P3, data=x, x=T, y=T, surv=T) plot(Predict(f)) # Works fine plot(Predict(f, time=3)) # Produces the error Best regards, -- Michal J. Figurski, PhD HUP, Pathology & Laboratory Medicine Biomarker Research Laboratory 3400 Spruce St. 7 Maloney Philadelphia, PA 19104 tel. (215) 662-3413 On 2010-05-27 10:44, Terry Therneau wrote:
but 'plot(Predict(f, time=3))' throws an error:"Error in summary.survfit(g, print.it = FALSE, times = times) : object 'n.enter' not found"Congratulations! You just discovered a new, never-before-encountered bug in the survival package. I should have a repair posted to Rforge by the end of the day. Terry Therneau (author of the package)
______________________________________________ [email protected] 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.
