r - 将Stata代码转换为R：调查权重

Question

我正在尝试对 R 中的调查数据进行加权。我使用 Stata 代码作为参考。

该程序上用于加权我的数据的代码是：

svyset [iw=wgtp], vce(sdr), sdrweight(wgtp1-wgtp80) mse

我将该代码翻译成这样的 R 代码：

svydesign(ids = ~1, data = house.data, weights = (~wgtp1-wgtp80))

但是，我不确定这是否是正确的翻译，以及我是否在 R 中正确应用它。任何见解和建议将不胜感激！

Answer 1

使用当前版本的“调查”，您想要

svrepdesign(data=house.data, weights=~wgtp, repweights="wgtp[1-9]+",type="successive-difference")

这svrepdesign适用于具有重复权重、抽样权重wgtp和重复权重wgtp1-wgtp80以及“连续差异”权重的调查。