簡介 - 引數生存模型與生存包的基本擬合和繪圖

survival 是 R 中最常用的生存分析包。使用內建的 lung 資料集，我們可以通過將回歸模型與 survreg() 函式擬合，建立一條帶有 survfit() 的曲線，並通過呼叫繪製預測的生存曲線來開始生存分析。具有新資料的此包的 predict 方法。

在下面的示例中，我們繪製了 2 條預測曲線，並在 2 組新資料之間改變了 sex，以顯示其效果：

require(survival)
s <- with(lung,Surv(time,status))

sWei  <- survreg(s ~ as.factor(sex)+age+ph.ecog+wt.loss+ph.karno,dist='weibull',data=lung)

fitKM <- survfit(s ~ sex,data=lung)
plot(fitKM)

lines(predict(sWei, newdata = list(sex      = 1, 
                                   age      = 1, 
                                   ph.ecog  = 1, 
                                   ph.karno = 90,
                                   wt.loss  = 2),
                                 type = "quantile",
                                 p    = seq(.01, .99, by = .01)),
                                 seq(.99, .01, by        =-.01),
                                 col = "blue")

lines(predict(sWei, newdata = list(sex      = 2, 
                                   age      = 1, 
                                   ph.ecog  = 1, 
                                   ph.karno = 90,
                                   wt.loss  = 2),
                                 type = "quantile",
                                 p    = seq(.01, .99, by = .01)),
                                 seq(.99, .01, by        =-.01),
                                 col = "red")

http://i.stack.imgur.com/G0SCt.jpg