1

我使用倾向得分匹配将病例与对照进行匹配(matchIt 包,比例 1:2)。现在我想使用权重(权重已知)获得两组年龄变量的独立 t 检验统计量。

我之前使用过相同的代码并且它有效。它适用于 20 个 obs 的样本数据集;确实适用于 10 个 obs 或整个数据集(856 个 obs)。我现在不知道出了什么问题。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