为什么我通过调用boot.ci来获取错误消息?我想获得AUC的置信区间。我也不太了解d< - data [indices,]是否有用并且应该在那里。
欢迎提出所有意见和建议。
期待您的回复,
library(survivalROC)
library(boot)
data(mayo)
attach(mayo)
bs <- function(data, indices) {
# d <- data[indices,] # allows boot to select sample
surv.res = survivalROC(Stime = time, status = censor, marker =mayoscore5, predict.time = 100,span = 0.25*NROW(mayo)^(-0.20))
return(surv.res$AUC)
}
boot.res <- boot(data=mayo, statistic=bs, R=100)
boot.ci(boot.res)
1] "All values of t are equal to 0.975215095269041 \n Cannot calculate confidence intervals"
NULL