Josh,
Good point. And it actually accidentally solved my problem.
Previously I attached .Rdata files. After your email I exported this
dataset into a csv file and then re-loaded it to see if the problem
persists, and... the problem disappeared.
Nevertheless, if I load the original .Rdata fil
Actually, I have another problem with the same data - this time with
plotting simple KM lines. The dataset is attached.
In this excercise I was testing two ways of representing time from the
previously attached interval dataset. Here t1 is the "stop" time at the
last interval for each patient,
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(Pre
> 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 e
On 05/26/2010 10:55 AM, Michal Figurski wrote:
Dear R-helpers,
I am working with 'cph' models from 'rms' library. When I build simple
survival models, based on 'Surv(time, event)', everything is fine and I
can make nice plots using plot(Predict(f, time=3)).
However, recently I tried to be more
5 matches
Mail list logo
