1

傾向スコア マッチングを使用して、ケースとコントロールを一致させました (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

前もって感謝します!

4

0 に答える 0