Расширенные Линии Участка Выживания в R

Question

Я получил график выживания из следующего кода:

    s = Surv(outcome.[,1], outcome.[,2])
    survplot= (survfit(s ~ person.list[,1]))
    plot(survplot, mark.time = FALSE)

person.list - это просто список из 15 человек.

Когда я строю график, все линии на моем графике заканчиваются в разные моменты времени.Есть ли способ расширить все строки, чтобы они заканчивались в определенный момент времени?(т.е. исход. [, 1] - переменная времени до события, и я хотел бы, чтобы линии выживания на графике расширились до 5 (лет))

Спасибо, Мэтт

Answer 1

Это код, который я отправил на похожий вопрос в rhelp некоторое время назад: http://finzi.psych.upenn.edu/Rhelp10/2010-September/253817.html

  ?survfit   # to get a working example since you did not provide one
  lsurv2 <- survfit(Surv(time, status) ~ x, aml, type='fleming')
  plot(lsurv2, lty=2:3, xmax=300)  # drats, no effect of xmax
  str(lsurv2)      # so see the structure of the survfit object
  lsurv2$time[21] <- 300         #add a time value
  lsurv2$n.censor[21] <- 1       # mark as censoring time
  lsurv2$strata[2] <- 11         # add to count of group 2

  plot(lsurv2, lty=2:3, xmax=300)  # horizontal line to 300 for group 2

И это был более поздний ответ Терно (предположительно, лучше, чем у меня): http://finzi.psych.upenn.edu/Rhelp10/2010-September/253879.html

plot(surv, mark.time=F, fun='event', xlim=c(0, 54))
  for (i in 1:length(surv$strata)) { #number of curves
     temp <- surv[i]
     lines(c(max(temp$time), 54), 1- rep(min(temp$surv),2))
     }

Answer 2

Это не ответ о том, как делать то, что вы просите, а скорее объяснение того, почему вы не должны делать то, что вы просите.

Строки останавливаются там, где заканчиваются данные. По истечении этого времени у вас нет информации для оценки выживаемости (это делается в традиционном анализе выживаемости Каплана-Мейера, как вы его установили). Таким образом, оценка Каплана-Мейера не очень хорошо определена за пределами этого времени, и поэтому расширение этой кривой не имеет какого-либо конкретного значения. Хотя графически вы можете просто нарисовать горизонтальную линию на том же уровне, что и последнее значение выживания, на самом деле это не имеет смысла.