0

在不使用 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 以及这个论坛都是新手。

预先感谢您的帮助。

4

2 回答 2

2

R 文档指定函数的参数weightslm观测值的方差成反比。这是分析权重的定义,或者aweights在 Stata 中。

看看ipw逆概率加权的包。

于 2012-06-15T23:58:27.497 回答
-1

为了更正先前的答案-我查阅了有关权重的手册,并找到了以下有关lm中权重的描述

非 NULL 权重可用于表示不同的观测值具有不同的方差(权重中的值与方差成反比);或者等价地,当权重的元素是正整数 w_i 时,每个响应 y_i 是 w_i 个单位权重观测值的平均值(包括 w_i 个观测值等于 y_i 并且数据已经汇总的情况)。

这些实际上是频率权重(stata 中的 fweights)。它们将观察值乘以权重向量定义的n次。另一方面,概率权重是指观察组包含在总体中的概率。这样做会调整观测值对系数的影响,但不会调整对标准误差的影响,因为它们不会改变样本中表示的观测值数量。

于 2017-02-23T20:26:42.113 回答