加权 Kaplan Meier 曲线上的 R 比例与具有 svymean 函数的比例

问题描述

我正在使用调查设计对象 (RORpps) 绘制 6m 至 17 岁儿童 ROR 疫苗接种覆盖率 (VC) 的 Kaplan Meier 曲线。这是我的代码：

fit <- svykm(Surv(dROR2,ROR2) ~ 1,design = RORpps,se =  T)
svyjskm(fit,ci = F,cumhaz = T,ylims = c(0,0.8),xlims = c(0,13),ystrataname = "",timeby = 1,showpercent = T,ylab = "VC rate (%)",xlab = "Age",surv.scale = "percent",pval = T,table = F,pval.testname = T)

图中的百分比表示使用 svymean 函数时的覆盖率为 71.1%

svymean(~ ROR2,RORpps,na.rm = TRUE)

我得到了 70.1%。为什么会有这种差异？不应该是一样的吗？预先感谢您的帮助。

解决方法

我是 jskm 软件包的开发者。

svymean计算出的比例没有考虑时间(dPOR2)，因此与kaplan-meier估计不同。