教程:如何手动计算卡方检验统计量和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}