Как рассчитать хи-квадрат в 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