下面是一个简单的例子,假设你有一个包含两个分类变量的数据集,你想要检验它们之间的关联性:
# 创建一个示例数据框
data <- data.frame(
Category1 = c("A", "B", "A", "B", "A"),
Category2 = c("X", "Y", "X", "Y", "Z")
)
# 执行卡方检验
result <- chisq.test(table(data$Category1, data$Category2))
# 打印卡方检验的结果
print(result)
在这个例子中,table(data$Category1, data$Category2)用于创建一个列联表,然后chisq.test()函数用于执行卡方检验。
卡方检验的输出将包括卡方统计量、自由度和p值等信息。p值小于显著性水平(通常为0.05),则可以拒绝原假设,表示两个分类变量之间存在显著关联。
请确保你的数据满足卡方检验的假设,特别是每个单元格的预期频数不应太小。
转载请注明出处:http://www.zyzy.cn/article/detail/6429/R语言