因此,我正在使用生存软件包3.1-7和R版本3.6.1。我使用bpcp包中的leuk数据集。
我想自己计算标准误差,并用R来检查结果(我也用Kaplan-Meier估计器也将其除掉,一切都很好),结果却有所不同。这意味着,当我们计算置信区间时,它们也有很大不同(虽然相差不大,但仍然如此)。我检查了文档,并使用与“ fleming-harrington”类型相同的公式(而不是“ fh2”)。
fh <- survfit(freireich~1,type='fleming-harrington',conf.type='plain')
summary(fh)
S_6_HF <- exp(-fh$n.event[1]/fh$n.risk[1])
sd6hf <- S_6_HF*sqrt(fh$n.event[1]/(fh$n.risk[1]*(fh$n.risk[1]-fh$n.event[1])))
当我的标准误差(sd6hf)为0.0772时,survfit函数为第一行提供的标准误差为0.0715。考虑到我的生存率(S_6_HF)与Survfit的生存率相同,那怎么可能呢? survfit在计算标准误差时还会做其他事情吗?
非常感谢您。