Lumley 調査パッケージを使用せずに、R で確率の重みを指定する方法を教えていただければ幸いです。現在svyglmをサポートしていない今井らのメディエーションパッケージを使用して、Rでメディエーション分析を行っています。
私が現在実行しているコードは次のとおりです。
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とこのフォーラムに不慣れです。
よろしくお願いいたします。