使用 Cochran–Mantel–Haenszel 检验评估分类变量的独立性 - R语言

使用 Cochran–Mantel–Haenszel 检验评估分类变量的独立性 - R语言

Cochran–Mantel–Haenszel (CMH) 检验是一种用于评估两个分类变量之间独立性的统计方法。它广泛应用于医学研究和流行病学领域，通常用于分析二元变量在控制其他相关变量后的关联性。本文将介绍如何使用R语言中的 mantelhaen.test 函数执行CMH检验，并提供相应的源代码。

首先，我们需要确保已经安装并加载了适当的R软件包，以便能够使用 mantelhaen.test 函数。在R环境中，使用以下命令来安装和加载 “vcd” 软件包：

install.packages("vcd")
library(vcd)

完成安装和加载后，我们可以使用 mantelhaen.test 函数来执行 CMH 检验。此函数的基本语法如下：

mantelhaen.test(x, y, stratify = NULL)

其中，x 和 y 是我们要比较的两个分类变量。stratify 参数是一个可选的分层变量，用于控制其他相关变量的影响。如果没有分层变量，可以将 stratify 参数设置为 NULL。

下面我们通过一个示例来演示如何使用 mantelhaen.test 函数执行 CMH 检验。假设我们有一个名为 data 的