匹配数据R的加权t检验
我使用倾向评分匹配来匹配病例和对照组(匹配包,比率1:2)。现在我想使用权重(权重已知)得到两组中年龄变量的独立t检验统计量 我以前也用过同样的代码,而且很有效。它适用于20个obs的样本数据集;不适用于10个obs或整个数据集(856个obs)。我现在想不出怎么了。var age的类型为整数。 以下是我使用的代码(所有受试者都在同一数据集中“匹配的_数据”,w=权重变量): 以下是警告信息:匹配数据R的加权t检验,r,statistics,matching,weighted,R,Statistics,Matching,Weighted,我使用倾向评分匹配来匹配病例和对照组(匹配包,比率1:2)。现在我想使用权重(权重已知)得到两组中年龄变量的独立t检验统计量 我以前也用过同样的代码,而且很有效。它适用于20个obs的样本数据集;不适用于10个obs或整个数据集(856个obs)。我现在想不出怎么了。var age的类型为整数。 以下是我使用的代码(所有受试者都在同一数据集中“匹配的_数据”,w=权重变量): 以下是警告信息: Warning message: In wtd.t.test(matched_data$age[mat
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
提前谢谢 只有当x和y的长度不相同时,才会出现警告。这不是错误,只是警告。警告是为了让你知道有些东西可能是错的,但不一定是错的。谢谢大家,我最终上传并清理了数据集,效果更好。相同的变量,相同的观察次数。
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