Чи-квадрат добротности пригонки в R

У меня есть вектор наблюдаемых значений, а также вектор значений, вычисленных с помощью модели:

actual <- c(1411,439,214,100,62,38,29,64)
expected <- c(1425.3,399.5,201.6,116.9,72.2,46.3,30.4,64.8)

Теперь я использую критерий пригодности Chi-squared, чтобы увидеть, насколько хорошо работает моя модель. Я написал следующее:

chisq.test(expected,actual)

но это не сработает. Можете ли вы мне помочь?

Ответы

Ответ 1

X ^ 2 = 10,2 при 7 градусах свободы даст вам p ~ 0,18.

> 1-pchisq(10.2, df = 7)
[1] 0.1775201

Вы должны передать ожидаемые значения в аргументе p. Убедитесь, что ваши значения суммируются до 1.

> chisq.test(actual, p = expected/sum(expected))

    Chi-squared test for given probabilities

data:  actual 
X-squared = 10.2581, df = 7, p-value = 0.1744

Это о том, что делает тест X ^ 2. Вы даете функции модель (expected) и спрашиваете: насколько вероятно, что мои данные observed пришли из популяции, которая "сгенерирована" expected?