如何在 R 编程中计算卡方？

Question

这是我的名为 mydata 的数据框。

                   Denomination     
Attendance    Protestant Catholic Jewish
Regular          182      213    203
Irregular        154      138    110

我想通过卡方检验来计算这个假设。所以我这样设定我的假设：

H0：面额与出勤率没有关系

H1：面额和出勤率有关系。

我尝试通过计算实际值和期望值并比较它们来计算卡方。但是没有得到输出。

我试过的代码：

> rowSums(mydata,na.rm = FALSE,dims = 1L)
> colSums(mydata,na.rm = FALSE,dims = 1L)
> sum(mydata)
> e = rowSums(mydata) * colSums(mydata)/ 1000
> chisq.test(mydata) = sum((mydata-e)^2 / e)

但我没有得到结果。请推荐？

Answer 1

这里有几点需要注意：

1. chisq.test 只能有数字列（或整数）。您显示数据的方式似乎 Attendance 是它自己的列，而不是 rowname。
2. 您无需计算rowSums 或colSums 即可使用测试。

这是我运行的代码：

dfs <- data.frame(Attendance = c("Regular", "Irregular"), Protestant = c(182, 154), Catholic = c(213, 138), Jewish = c(203,110))
rownames(dfs) <- dfs[,1]
dfs$Attendance <- NULL
chisq.test(dfs)

    Pearson's Chi-squared test

data:  dfs
X-squared = 7.8782, df = 2, p-value = 0.01947