I’ve obtained a survival plot from the following code:
s = Surv(outcome.[,1], outcome.[,2])
survplot= (survfit(s ~ person.list[,1]))
plot(survplot, mark.time = FALSE)
person.list is just a list of 15 people.
When I plot this, the lines on my plot all end at different time points. Is there a way to extend all the lines to make them end at a certain time point? (i.e outcome.[,1] is a time to event variable and I would like the survival lines on the plot to extend out to say 5(years) )
Thanks,
Matt
This isn’t an answer of how to do what you ask, but rather an explanation of why you should not do what you ask.
The lines stop where the data stops. Beyond that time, you have no information in order to make an estimate of the survival (this is in a traditional Kaplan-Meier survival analysis, as you have set it up). Therefore, the Kaplan-Meier estimate is not well defined beyond that time, and so extending that curve does not have any particular meaning. While graphically you could just draw a horizontal line at the same level as the last survival value, this is not really meaningful.