问题描述
我最近使用“coxed”包中的“sim.survdata”函数来模拟生存数据。我用一个协变量和系数 = 1 进行了 2k 次模拟。每次我拟合 Cox 回归并记录估计的 beta 和 95% CI。然而，均值（betahat）不等于 1，并且 CI 不会以 95% 的次数覆盖 1。请参阅我在下面使用的代码。我想知道原因，欢迎任何输入。提前致谢。
bhat=NULL
lb=NULL
ub=NULL
for(i in 1:2000){
set.seed(20110312+i)
simdata <- sim.survdata(N=2633,T=200,xvars=1,censor=.85,num.data.frames = 1,beta = 1)
simdata$data$event=as.numeric(simdata$data$Failed)
fit=coxph(Surv(y,event==1)~X,data=simdata$data)
bhat=c(bhat,summary(fit)$coefficients[1])
lb=c(lb,log(summary(fit)$conf.int[3]))
ub=c(ub,log(summary(fit)$conf.int[4]))
}
> mean(bhat)
[1] 0.9408025
> mean(lb<=1 & ub>=1)
[1] 0.904
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）