傾向スコア マッチングを使用して、ケースとコントロールを一致させました (matchIt パッケージ、比率 1:2)。ここで、重みを使用して両方のグループの年齢変数の独立した t 検定統計を取得したいと考えています (重みは既知です)。
以前に同じコードを使用しましたが、うまくいきました。これは、20 個の obs のサンプル データセットで動作します。10 個の観測点またはデータセット全体 (856 個の観測点) に対して機能しません。私は今何が間違っているのか理解できません。var age の型は整数です。これが私が使用するコードです(すべての被験者は同じデータセット「matched_data」、w = weights varにあります):
wtd.t.test {weights}
wtd.t.test(matched_data$age[matched_data$group==1],
matched_data$age[matched_data$group==0],
weight=matched_data$w[matched_data$group==1],
weighty=matched_data$w[matched_data$group==0], samedata=FALSE)
警告メッセージは次のとおりです。
Warning message:
In wtd.t.test(matched_data$age[matched_data$group == 1], matched_data$age[matched_data$group == :
Treating data for x and y separately
> sessionInfo()
R version 2.15.2 (2012-10-26)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
サンプル データセット:
group age w
0 25 1.911564626
1 30 1
0 72 1.911564626
1 72 1
1 72 1
0 58 0.955782313
0 83 0.955782313
0 61 0.955782313
0 71 0.955782313
0 37 0.955782313
1 77 1
1 78 1
1 65 1
0 89 0.955782313
1 50 1
1 89 1
1 63 1
0 81 0.955782313
1 65 1
0 91 0.955782313
前もって感謝します!