我使用倾向得分匹配将病例与对照进行匹配(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
提前致谢!