โมเดล Weibull สำหรับประมาณเส้นโค้งการรอดชีวิตแบบราบเรียบ

Survival Analysis in R

Heidi Seibold

Statistician at LMU Munich

ทำไมต้องใช้โมเดล Weibull?

โมเดล Weibull:

wb <- survreg(Surv(time, event) ~ 1, data)

โมเดล Weibull:

wb <- survreg(Surv(time, event) ~ 1, data)

การประมาณค่า Kaplan-Meier:

km <- survfit(Surv(time, event) ~ 1, data)

wb <- survreg(Surv(time, cens) ~ 1, data = GBSG2)

ผู้ป่วย 90 เปอร์เซ็นต์รอดชีวิตเกินจุดเวลา:

predict(wb, type = "quantile", p = 1 - 0.9, newdata = data.frame(1))

       1 
384.9947

p = 1 - 0.9 เพราะฟังก์ชันการแจกแจงคือ 1 ลบด้วยฟังก์ชันการรอดชีวิต

wb <- survreg(Surv(time, cens) ~ 1, data = GBSG2)

เส้นโค้งการรอดชีวิต:

surv <- seq(.99, .01, by = -.01)

t <- predict(wb, type = "quantile", p = 1 - surv, newdata = data.frame(1))
head(data.frame(time = t, surv = surv))

       time surv
 1  60.6560 0.99
 2 105.0392 0.98
 3 145.0723 0.97
 4 182.6430 0.96
 5 218.5715 0.95
 6 253.3125 0.94

Survival Analysis in R