計算 Weibull 模型
對於含有共變數的 Weibull 模型，我們可以像沒有共變數的情況一樣計算存活曲線。唯一需要做的是在 predict() 函式裡，為想要的存活曲線指定共變數數值。這可以用引數 newdata 來完成。
本練習屬於課程
R 的生存分析
練習說明
- 以
GBSG2資料集與共變數horTh建立一個 Weibull 模型，用來分析荷爾蒙治療對患者存活的影響。 - 計算有接受荷爾蒙治療患者的存活曲線。
- 使用
str()檢視存活曲線。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Weibull model
wbmod <- survreg(Surv(___, ___) ~ ___, data = GBSG2)
coef(wbmod)
# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
t_yes <- predict(___, type = "quantile", p = ___,
newdata = data.frame(horTh = ___))
# Take a look at survival curve
str(t_yes)