我观察到的值的矢量,并且还与模型计算的值的矢量:
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)
现在,我使用的是适合测试的卡方善良,看我的模型执行情况。 我写了下面:
chisq.test(expected,actual)
但它不工作。 你能帮助我吗?
我观察到的值的矢量,并且还与模型计算的值的矢量:
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)
现在,我使用的是适合测试的卡方善良,看我的模型执行情况。 我写了下面:
chisq.test(expected,actual)
但它不工作。 你能帮助我吗?
X ^ 2 = 10.2 7个自由度会给你AP〜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
?