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 文档指定函数的参数weights与lm观测值的方差成反比。这是分析权重的定义，或者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 回答