使用＆＃39; survcomp＆＃39;未能获得D指数的95％CI和＆＃39; boot＆＃39;包

Question

我使用了＆＃39; survcomp＆＃39;计算预测模型的D指数。代码如下：

source("http://bioconductor.org/biocLite.R")
biocLite("survcomp")
library(survcomp)
set.seed(12345)
age <- rnorm(100, 50, 10)
stime <- rexp(100)
cens <- runif(100,.5,2)
sevent <- as.numeric(stime <= cens)
stime <- pmin(stime, cens)
D.index(x=age, surv.time=stime, surv.event=sevent)

当我使用＆＃39; boot＆＃39;包计算模型的D指数的95％CI。

D指数CI

的计算

library(boot)
df <- data.frame(age, size, stime, sevent, cens)
boot_Dindex <- function(x, surv.time, surv.event, data, indices){
  d <- data[indices,]
  Dindex <- D.index(x=x, surv.time=surv.time, surv.event=surv.event)
  return(Dindex$d.index)
} 
set.seed(1234)
results <- boot(data=df,x=df$age, surv.time=df$stime, surv.event=df$sevent, statistic = boot_Dindex, R= 1000)  
print(results)

普通非参数引导程序 呼叫： boot（data = small_intestine_2004，statistic = boot_Dindex，R = 1000，     x =诺模图，生存时间=时间，幸存者= CSS）

Bootstrap统计：

 original  bias    std. error

t1 * 0.9096271 0 0

我不知道这个问题的原因是什么，没有偏见和标准。可以计算错误并将其重新设为0。