分散分析は表 1 のようになる。
| 要因 | 平方和 | 自由度 | 平均平方 | $F$ 値 | 有意確率 |
|---|---|---|---|---|---|
| 回帰 | 6.67164 | 2 | 3.33582 | 823.47652 | $\lt$ 0.00001 |
| 残差 | 0.02836 | 7 | 0.00405 | ||
| 全体 | 6.70000 | 9 | 0.74444 |
R を用いると以下のようになる。
> ## 回帰の分散分析
> MSr <- Sr/p # 回帰で説明される平均平方
> F <- MSr/MSe # 回帰の分散分析の F 値
> p.F <- pf(F, p, df, lower.tail=FALSE) # P 値
> c(F=F, df1=p, df2=df, p=p.F)
F df1 df2 p
8.234765e+02 2.000000e+00 7.000000e+00 4.931875e-09
直前のページへ戻る E-mail to Shigenobu AOKI