由 Weibull 模型取得存活曲線的分位數
現在我們可以用 Weibull 模型來估計 GBSG2 資料中乳癌患者的存活情形。
predict() 函式搭配 type = "quantile" 可以計算分配函數的分位數。到目前為止,你已經學到存活函數等於 1 減去分配函數($S = 1 - D$),因此我們可以用 predict() 函式輕鬆計算存活函數的分位數。
本練習已為你載入 survival 套件與 GBSG2 資料。
本練習屬於課程
R 的生存分析
練習說明
- 估計乳癌患者的 Weibull 模型。
- 取得「存活時間長於該時間點的機率為 70%」時對應的時間點(使用
predict()並設定type = "quantile")。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Weibull model
wb <- ___(Surv(___, ___) ~ 1, data = GBSG2)
# 70 Percent of patients survive beyond time point...
predict(wb, type = ___, p = ___, newdata = data.frame(1))