在不使用 Lumley 调查包的情况下,我非常感谢在 R 中指定概率权重的任何帮助。我正在使用 Imai 等人的中介包在 R 中进行中介分析,该软件包目前不支持 svyglm。
我目前正在运行的代码是:
olsmediator_basic<-lm(poledu ~ gateway_strict_alt + gender_n + spline1 + spline2 + spline3,
data = unifiedanalysis, weights = designweight).
但是,我不确定这是否正确地加权数据。原因是这段代码产生的标准错误与我在 Stata 中得到的不同。我正在运行的Stata代码是:
reg poledu gateway_strict_alt gender_n spline1 spline2 spline3 [pweight=designweight]).
我想知道 R 中的权重选项是否可能不适用于逆概率权重,但我无法从文档、本论坛或其他地方确定这一点。如果我遗漏了什么,我真的很抱歉 - 我对 R 以及这个论坛都是新手。
预先感谢您的帮助。