如何手动计算卡方检验统计量和p值

教程：如何手动计算卡方检验统计量和p值

步骤1：创建观察频率表和期望频率表

首先，您需要创建一个观察频率表（observed frequency table），该表包含您的观察数据，以及一个期望频率表（expected frequency table），根据假设的独立性来计算。观察频率表如下所示：

	Not Stopped	Bribe requested	Stopped/given warning
Upper class	14	6	7
Lower class	7	7	1

在R中，您可以使用以下代码来创建这些表格：

# 创建观察频率表
observed <- matrix(c(14, 6, 7, 7, 7, 1), nrow = 2, byrow = TRUE)

# 计算行总计
row_totals <- rowSums(observed)

# 计算列总计
col_totals <- colSums(observed)

# 计算总计
total <- sum(observed)

# 计算期望频率表
expected <- outer(row_totals, col_totals) / total

步骤2：计算卡方检验统计量

卡方检验统计量的计算方式如下：

[ \chi^{2}