[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@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Problem with tick marks in lines.survfit (package survival)
How about: plot(fit, mark.time=TRUE, xscale=365.24,xlab='Years', ylab='Survival',lwd=c(2,1)) Rob - Original Message - From: Rachel Pearce [EMAIL PROTECTED] To: r-help@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@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] 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 [EMAIL PROTECTED] 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 [EMAIL PROTECTED] To: r-help@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@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html