13

観測値のベクトルと、モデルで計算された値のベクトルがあります。

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)

現在、カイ2乗適合度検定を使用して、モデルのパフォーマンスを確認しています。私は次のように書いた:

chisq.test(expected,actual) 

しかし、それは機能しません。これを手伝ってくれませんか。

4

1 に答える 1

15

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ますか?

于 2012-07-10T08:25:11.347 に答える