将Stata代码翻译为R:调查权重

时间:2017-10-13 20:31:39

标签: r stata survey weighting

我试图在R中加权调查数据。我使用Stata代码作为参考。

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

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

我将这段代码翻译成R代码:

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

但我不确定这是否是正确的翻译,如果我在R中正确应用这一点。任何见解和建议将不胜感激!

1 个答案:

答案 0 :(得分:0)

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

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

对于具有重复权重,抽样权重svrepdesign和重复权重wgtp和“连续差异”权重的调查,就是wgtp1-wgtp80