Problem with tick marks in lines.survfit (package survival)
Sorry, yes I think that would work, for the example I gave, but I am trying to superimpose a curve onto an existing plot. The option of using "plot" to draw both curves is not available, since the data are (or were) in two different columns. For that reason I was trying to use "lines" for the second curve. I used the example given with the package so that I didn't have to publish my own data (which are not mine to publish, if you see what I mean). I have now given up and created a new data frame with OS and PFS stacked on top of each other, but this is a very time-consuming workaround. Do you (or anyone) know how to make "lines.survfit" draw tick marks? Rachel
On 19/09/05, Robert Baer <rbaer at atsu.edu> wrote:
How about: plot(fit, mark.time=TRUE, xscale=365.24,xlab='Years', ylab='Survival',lwd=c(2,1)) Rob ----- Original Message ----- From: "Rachel Pearce" <pearce.rachel at gmail.com> To: <r-help at stat.math.ethz.ch> Sent: Monday, September 19, 2005 4:59 AM Subject: [R] Problem with tick marks in lines.survfit (package survival)
I have attempted to follow posting guidelines but I have failed to find
out
what I am doing wrong here. I am trying to use lines.survfit to plot a second curve onto a survival curve produced by plot.survfit. In my case this is to be a progression
free
survival curve superimposed upon an overall survival curve, but I will illustrate my problem using the example given in the help for
lines.survfit.
I would like to have tick marks for censored data points on both curves,
so
I would like to do something like:
> fit <- survfit(Surv(time, status) ~ sex, pbc,subset=1:312) plot(fit[2], mark.time=TRUE, xscale=365.24,
+ xlab='Years', ylab='Survival')
lines(fit[1], lwd=2, xscale=365.24, mark.time=TRUE)
but when I do this, no tick marks appear on the second curve. Can anyone see what I am doing wrong here? Rachel [[alternative HTML version deleted]]
______________________________________________ R-help at stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!