Я хотел бы узнать о том, как получить доверительный интервал из предсказания значения и стандартной ошибки из модели выживания Вейбулла в R "выживание".
R doc для предиката.survreg содержит пример, показывающий график, который дает не только fit прогноза модели выживания Вейбулла, но также fit+2*se.fit и fit-2*se.fit. Тогда какому доверительному интервалу соответствует диапазон от fit-2*se.fit до fit+2*se.fit? И для данного доверительного интервала, как выбрать n для fit +/- n*se.fit, чтобы получить этот доверительный интервал? Ниже приведен код:
lfit <- survreg(Surv(time, status) ~ ph.ecog, data=lung)
pct <- 1:98/100 # The 100th percentile of predicted survival is at +infinity
ptime <- predict(lfit, newdata=data.frame(ph.ecog=2), type='quantile', p=pct, se=TRUE)
matplot(cbind(ptime$fit, ptime$fit + 2*ptime$se.fit, ptime$fit - 2*ptime$se.fit)/30.5, 1-pct, xlab="Months", ylab="Survival", type='l', lty=c(1,2,2), col=1)
Для нормального распределения, я думаю, для 95% доверительного интервала можно использовать fit +/- 1.96*se.fit. Я не уверен, работает ли он для модели Weibull Survreg.