Weibull 모형 계산
공변량이 있는 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)