の続きです。今回もANOVAです。explanatoly variables(説明変数)をType(最近、先行きの２値をとるカテゴリカル変数)にresponse variable(反応変数)はValue(短観の値)です。
まずは、前回と同じくグラフで見てみます。
黒い点が最近、赤い点が先行きです。黒の水平線が最近の平均値、赤の水平線が先行きの平均値です。
tapply関数でそれぞれの平均値を求めてみます。
先行きの平均値は、最近の半分以下ですね。
では、SSY、the sum of the squares of the differences between the y values and the overall meanを計算します。
SSYは68298です。
SSE、the error sum of squaresを計算します。最近の平均値、先行きの平均値はさきほど計算していますから、これを利用しましょう。
SSEは65572です。
SSA、the treatment sum of squaresはSSY - SSEで計算できます。
SSAは2726です。SSA, SSE, SSYがわかりましたので、ANOVA表を作ってみます。
D3のセルのところ、SSEのMean Squareを計算します。
E2のところがF ratioです。これは、2726/166です。
これでANOVA表が完成しました。このF ratioの16.38が大きいのか小さいのか、qf関数で判定します。
3.865よりも大きなF ratioだとSSAには意味があるということですから最近の平均と先行きの平均には有意な違いがありますね。
aov関数とsummary関数でも同じようにANOVA分析をします。
p値が5.24e-05と0.05以下ですから、最近と先行きでは、有意な差があります。先行きのほうが悲観的なのですね。
今回は以上です。
