在 R 中执行最近邻匹配时,是否可以查看哪些案例与哪些控件匹配的身份?
When performing nearest neighbour matching in R, is it possible to view the identity of which cases matched with which controls?
我首先在 R Studio 中使用此处找到的小型练习数据集(584 个观测值,5 个变量)(https://scholarworks.umass.edu/cgi/viewcontent.cgi?article=1330&context=pare)
进行尝试
使用此代码,我可以使用最近邻匹配来查找匹配案例和控件 (1:1) 的平均差异,其中 stw 是我的分组变量,tot、min 和 dis 是匹配变量:
m.out = matchit(stw ~ tot + min + dis,
data = mydata, method = "nearest",
ratio = 1)
我想知道的是如何查看哪些案例与哪些控件匹配(即向我显示确切的身份)?
我也想在匹配 5 个最近的邻居(即比率 = 5)并查看这些确切身份的情况下执行此操作。为此需要额外的代码吗?
非常感谢
匹配中的单位标识存储在 matchit
输出对象的 match.matrix
组件中。这是一个矩阵,其行数与处理的单元数一样多,列数为 ratio
。每行对应一个处理过的单元,行中的值对应于与该处理过的单元匹配的控制单元的标识。下面是一个例子:
> head(m$match.matrix)
1 2
NSW1 "PSID305" "PSID369"
NSW2 "PSID198" "PSID273"
NSW3 "PSID333" "PSID282"
NSW4 "PSID400" "PSID92"
NSW5 "PSID426" "PSID262"
NSW6 "PSID391" "PSID261"
处理单元 NSW1
与控制单元 PSID305
和 PSID369
匹配。如果原始数据集没有行名称,则值对应数据集中出现相应单位的行。
我首先在 R Studio 中使用此处找到的小型练习数据集(584 个观测值,5 个变量)(https://scholarworks.umass.edu/cgi/viewcontent.cgi?article=1330&context=pare)
进行尝试使用此代码,我可以使用最近邻匹配来查找匹配案例和控件 (1:1) 的平均差异,其中 stw 是我的分组变量,tot、min 和 dis 是匹配变量:
m.out = matchit(stw ~ tot + min + dis,
data = mydata, method = "nearest",
ratio = 1)
我想知道的是如何查看哪些案例与哪些控件匹配(即向我显示确切的身份)?
我也想在匹配 5 个最近的邻居(即比率 = 5)并查看这些确切身份的情况下执行此操作。为此需要额外的代码吗?
非常感谢
匹配中的单位标识存储在 matchit
输出对象的 match.matrix
组件中。这是一个矩阵,其行数与处理的单元数一样多,列数为 ratio
。每行对应一个处理过的单元,行中的值对应于与该处理过的单元匹配的控制单元的标识。下面是一个例子:
> head(m$match.matrix)
1 2
NSW1 "PSID305" "PSID369"
NSW2 "PSID198" "PSID273"
NSW3 "PSID333" "PSID282"
NSW4 "PSID400" "PSID92"
NSW5 "PSID426" "PSID262"
NSW6 "PSID391" "PSID261"
处理单元 NSW1
与控制单元 PSID305
和 PSID369
匹配。如果原始数据集没有行名称,则值对应数据集中出现相应单位的行。