为何重叠加权后love图的SMD与svyCreateTableOne函数绘制基线表的不一致

Question

用R语言对三组数据进行重叠加权，重叠加权后用bal.tab函数及love.plot函数去看均衡性显示有一个SMD>0.1,但是用svydesign及svyCreateTableOne函数绘制基线表得出加权后的SMD值都是<0.1，应该以哪个为标准的？

#####重叠加权#########
PS.fmu<-treatment ~ age + BMI + figo.stages + Pathological.type+ lymph.node.metastasis
SumSt<-SumStat(ps.formula=PS.fmu,data=bc,weight=c("overlap"),method="glm") 
IPW<-SumSt[["ps.weights"]][,3] 

bal_tab <- bal.tab(
  treatment ~ age + BMI + figo.stages + Pathological.type+ lymph.node.metastasis,
  data = bc,
  weights = IPW,
  method = "weighting",
  stats=c("m","v"),continuous="std",binary="std",thresholds=c(m=0.1,v =2),pairwise=TRUE,abs=FALSE,un=TRUE, s.d.denom = "pooled"
)
###SMD可视化
love.plot(bal_tab,threshold=0.1,agg.fun="max",binary="std")

###得出加权后的基线表
Svy.design<-svydesign(ids=~0,weights=bc$weight,data=bc) 
Tabipw.stas<-svyCreateTableOne(vars=allVars,strata="treatment",data=Svy.design,test=TRUE,smd=TRUE,addOverall=TRUE) 
Tab.ipw<-print(Tabipw.stas, nonnormal = c("age"), smd = TRUE,
               formatOptions = list(big.mark = ","))

Answer 1

阿里嘎多学长整理AIGC生成，因移动端显示问题导致当前答案未能完全显示，请使用PC端查看更加详细的解答过程

为何重叠加权后love图的SMD与svyCreateTableOne函数绘制基线表的不一致

你遇到的情况是，在重叠加权后，使用 love.plot 函数显示的 SMD（Standardized Mean Difference）大于 0.1，但是在使用 svyCreateTableOne 函数绘制基线表时，结果不同。这个问题可能是由于两种方法对样本权重的处理方式不同所致。

love.plot 函数使用的是简单的加权平均，而 svyCreateTableOne 函数使用的是urveyighted 加权平均。前者对所有 sampling unit 使用相同的加权，而后者使用 weighted 加权 averaging， weight 是根据 survey design 的结果确定的。

可能的解决方法是，在使用 svyCreateTableOne 函数时，指定使用之前计算的 weighted 加权平均结果，例如：

library(RSurvey)
# ...

# 重叠加权
weight <- ...

# 使用 svyCreateTableOne 函数
svyCreateTableOne(~variable, design = svydesign(weight = weight))

这样可以确保 svyCreateTableOne 函数使用的是之前计算的 weighted 加权平均结果，而不是重新计算的结果。