У меня есть вектор наблюдаемых значений, а также вектор значений, вычисленных с помощью модели:
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)
но это не сработает. Можете ли вы мне помочь?
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?