R语言频数检验：Cochran-Mantel-Haenszel检验用于检验重复测量两离散变量的独立性

R语言频数检验：Cochran-Mantel-Haenszel检验用于检验重复测量两离散变量的独立性

Cochran-Mantel-Haenszel（CMH）检验是一种常用的统计方法，用于检验两个离散变量之间的独立性，尤其适用于重复测量的情况。本文将详细介绍如何使用R语言进行CMH检验，并提供相应的源代码。

首先，我们需要确保已经安装了R语言环境，并加载所需的包。在R中，我们可以使用vcd包进行CMH检验，因此需要先安装该包并加载。

# 安装vcd包
install.packages("vcd")

# 加载vcd包
library(vcd)

接下来，我们需要准备用于分析的数据。假设我们有两个离散变量X和Y，分别包含重复测量的观察值。我们可以将这些观察值组织成一个交叉表格，其中行表示X的水平，列表示Y的水平。

# 创建交叉表格
data <- matrix(c(10, 20, 30, 40, 15, 25, 35, 45), nrow = 2, ncol = 4, byrow = TRUE)
rownames(data) <- c("X1", "X2")
colnames(data) <- c("Y1", "Y