如何解决sim.survdata
我最近使用“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 举报，一经查实，本站将立刻删除。
