r语言 卡方检验可视化

在R语言中，可以使用`chisq.test()`函数进行卡方检验，并使用可视化工具来展示结果。下面是一个示例代码：

# 创建一个数据框，包含两个分类变量
data <- data.frame(
  category1 = c("A", "B", "A", "B", "A"),
  category2 = c("X", "Y", "X", "X", "Y")
)

# 进行卡方检验
result <- chisq.test(table(data$category1, data$category2))

# 打印卡方检验结果
print(result)

# 可视化卡方检验结果
library(ggplot2)

# 提取卡方检验结果的预期频数
expected <- as.data.frame(result$expected)

# 添加预期频数和实际频数的列名称
colnames(expected) <- colnames(table(data$category1, data$category2))
expected$category1 <- rownames(expected)

# 将数据转换为长格式以方便绘图
expected_long <- tidyr::pivot_longer(expected, -category1, names_to = "category2", values_to = "frequency")

# 绘制柱状图
ggplot(expected_long, aes(x = category1, y = frequency, fill = category2)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(x = "Category 1", y = "Frequency", fill = "Category 2") +
  theme_minimal()

这段代码首先创建了一个包含两个分类变量的数据框，并使用`chisq.test()`函数进行卡方检验。然后，通过`ggplot2`库将卡方检验结果的预期频数可视化为柱状图。你可以根据自己的数据进行相应的修改和调整。