deaR user,
I found an unexpected behaviour of the rms::survplot.survfit function, that
is giving me inverted labels for the patient-at-risk rows. The problem is
that, for some reason, the survival::summary.survfit function changes the
order of two of the suvfit object's strata when called (on my dataset) with
the times= option (is this another unexpected behaviour?). survplot takes
the labels from the fit passed as argument and the numbers at risk from a
summary(fit, time=stime) object, thus inverting the two labels.
I suppose the problem is caused by some feature of my dataset but I didn't
manage to reproduce the problem on simulated data (although I didn't try
hard) I'm thinking of debugging also summary.survfit to understand why this
is happening. Before doing doing that, however, I'd like to know if some of
you have some hints that could help me either to solve the problem or to
understand it.
Thank you
Marco.
[[alternative HTML version deleted]]
______________________________________________
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.
