从匹配的同类群组中提取数据并运行mcnemar / chi sq测试

时间:2019-04-25 23:30:12

标签: r statistics extract chi-squared

如何创建一个数据框,其中包含匹配子集的信息,但包含不包含为协变量但包含在原始数据集中的变量?

数据集包含: y(曝光变量) x1,x2,x3(协变量) z(结果变量)

暴露变量(y)与基线协变量的对照的倾向匹配: mod <-y〜x1 + x2 + x3,数据,方法,游标卡尺

我对y和结果变量(z)之间的关系感兴趣,这两个变量都是二进制(0/1)。 我可以获取匹配矩阵的信息,但是由于z不在倾向匹配中,因此如何创建一个仅包含来自匹配样本的主题的子集。之后,我希望仅提取匹配样本的z数据,并对y和z进行卡方检验或mcnemar测试。

我可以通过执行以下操作来查看匹配示例的详细信息: mod.data <-match.data(mod) 但它不包含感兴趣的变量(z),但z在原始数据集中。 我也通过执行mod $ match.matrix知道匹配主题的ID 我尝试使用似乎有效的data [mod,],但它没有给出2x2表,而仅给出1x2,其中y ==0。但是,我需要它报告y == 1,所以我可以执行卡方测试/麦克内马尔测试。

我是R和倾向得分匹配的初学者,因此,任何建议/帮助都将不胜感激。预先非常感谢。

0 个答案:

没有答案