使用Lumley調查軟件包在R * without *中指定概率權重

我真的很感謝使用Lumley調查軟件包在R 沒有中指定概率權重的任何幫助。我正在使用Imai等調解軟件包進行R中的調解分析，該軟件包目前不支持svyglm。使用Lumley調查軟件包在R * without *中指定概率權重

我目前正在運行的代碼是：

olsmediator_basic<-lm(poledu ~ gateway_strict_alt + gender_n + spline1 + spline2 + spline3, 
    data = unifiedanalysis, weights = designweight).

不過，我不能確定這是否正確加權的數據。原因是這段代碼產生的標準錯誤與我在Stata中獲得的錯誤不同。我運行的是塔塔代碼：

reg poledu gateway_strict_alt gender_n spline1 spline2 spline3 [pweight=designweight]).

我在想，如果R中的權重的選擇可能不是逆概率權，但我無法從文檔，這個論壇或其他地方確定這一點。如果我錯過了一些東西，我真的很抱歉 - 我對R以及本論壇都是新手。

非常感謝您的幫助。

來源

2012-06-14 sabaya

如果您提供了一個可重複的示例，這會更容易。 'designweight'中有什麼值？調查概率或逆概率。如果你想用逆概率加權，你必須把這些作爲權重。 – mnel

感謝您回覆mnel並對此抱歉。設計權值是逆概率權重。任何我應該做的不同？ – sabaya

這裏是一個指南，以提問在Stackoverflow的好問題的鏈接：http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example通過遵循其中的建議，你可能會增加您有機會收到有用的回覆。一個可重複的例子將會有所幫助，比如在運行這個可重現的例子時，包括R和Stata的輸出。 –

R documentation指定lm函數的weights參數與觀察值的方差成反比。這是Stata中分析權重的定義或aweights。

查看用於反向概率加權的ipw包。

來源

2012-06-15 23:58:27

-1

要糾正以前的答案 - 我擡頭一看說明書上的權重，並找到了權重下面的描述中流明

非空的權重可以用來表示不同的意見有不同的方差（權重值與方差成反比）;或等價地，當權重的元素是正整數w_i時，每個響應y_i是w_i單位重量觀測值的均值（包括觀測值等於y_i且數據已被彙總的情況）。

這些實際上是頻率權重（stata中的fweights）。它們乘以由權重向量定義的觀察次數的觀察值n。另一方面，概率權重指的是觀察組包含在人羣中的概率。這樣做會調整觀察對係數的影響，但不在標準誤差上，因爲它們不會更改樣本中表示的觀察值的數量。

來源

2017-02-23 20:26:42 EconomySizeAl

我知道在我回答時已有近5年的歷史了，但是一旦我偶然發現整個辦公室進入加權討論階段 – EconomySizeAl

這是不正確的。「lm」和「glm」實際上實現了分析權重（Stata中的「aweight」），也稱爲逆方差權重（與您給出的描述相匹配）。參見例如[這個鏈接]（http://bc.bojanorama.pl/2015/09/linear-models-with-weighted-observations/）。 –

使用Lumley調查軟件包在R * without *中指定概率權重

回答

相關問題