使用Lumley调查包在没有*的情况下指定R *中的概率权重

Question

我真的很感激使用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以及这个论坛。

提前感谢您的帮助。

Answer 1

R documentation指定weights函数的lm参数与观察值的方差成反比。这是分析权重的定义，或Stata中的aweights。

查看ipw包的逆概率加权。