分散分析パッケージ
このパッケージは,集団間での平均値の違いを測定する一変量分散分析(ANOVA)を実行するための関数を提供するものである.関数ANOVAは,直交計画において任意数の固定因子を持つモデルを扱うことができる.また,要素が揃っているかどうかにかかわらずバランスデータ・アンバランスデータのどちらを扱うこともできる.結果はすべてType I平方和で与えられる.ANOVAでは,比較のための多様な事後検定を実行することもできる.
ANOVA関数
data は{{α1,β1,…,y1},{α2,β2,…,y2},…}という形式でなければならない.ここで αi,βi 等は i 番目の応答 yi と関連したカテゴリ変数 vars の値である.
引数 vars は,モデルのカテゴリ変数を表すシンボルのリストである.
引数 model は,モデルを共同で指定する主効果および交互作用のリストでなければならない.交互作用項は変数の積として与えられる.例えば,三元配置分散分析のすべての因子モデルは{α,β,γ,α β,α γ,β γ,α β γ}と書き表すことができる.ここで α,β,γ とは主効果を,α β,α γ,β γ とは二因子交互作用,α β γ とは三因子交互作用を示す.モデルはAllを使って指定されたカテゴリ変数間のすべての主効果と交互作用を表すように書くこともできる.したがって,三元配置分散分析のすべての因子モデルは{α,β,γ,All}と書くこともできる.
<<ANOVA`onewaydata = {{1, 7.0}, {1, 5.3}, {1, 5.9}, {1, 6.6}, {1, 4.9}, {2, 4.4}, {2, 6.8}, {2, 7.7}, {2, 8.3}, {2, 6.6}, {3, 8.1}, {3, 10.4}, {3, 8.0}, {3, 6.8}, {3, 9.2}, {4, 5.7}, {4, 3.9}, {4, 6.2}, {4, 5.5}, {4, 6.2}};ANOVA[onewaydata]twowaydata = {{1, 1, 10.1}, {1, 1, 10.5}, {1, 1, 11.3}, {1, 2, 13.1}, {1, 2, 14.7}, {1, 3, 14.1}, {1, 3, 12.6}, {2, 1, 10.7}, {2, 1, 15.3}, {2, 1, 17.9}, {2, 1, 18.0}, {2, 2, 28.7}, {2, 3, 16.0}, {2, 3, 9.2}, {2, 3, 12.1}};ANOVA[twowaydata, {factor1, factor2, All}, {factor1, factor2}]ANOVA[Drop[twowaydata, {-4}], {factor1, factor2, All}, {factor1, factor2}]threewaydata = {{1, 1, 1, 50}, {1, 1, 1, 50}, {1, 1, 1, 54}, {1, 1, 2, 40}, {1, 1, 2, 36}, {1, 1, 2, 40}, {1, 2, 1, 48}, {1, 2, 1, 48}, {1, 2, 1, 44}, {1, 2, 2, 14}, {1, 2, 2, 18}, {1, 2, 2, 14}, {2, 1, 1, 40}, {2, 1, 1, 36}, {2, 1, 1, 36}, {2, 1, 2, 18}, {2, 1, 2, 14}, {2, 1, 2, 18}, {2, 2, 1, 6}, {2, 2, 1, 2}, {2, 2, 1, 2}, {2, 2, 2, 20}, {2, 2, 2, 16}, {2, 2, 2, 20}};ANOVA[threewaydata, {α, β, γ, α β, α γ, β γ}, {α, β, γ}]
オプション名
|
デフォルト値
| |
| CellMeans | True | 各集合の平均値のリストを与えるかどうか |
| WorkingPrecision | MachinePrecision | 計算で使用する演算精度 |
| PostTests | {} | 実行する事後検定のリスト |
| SignificanceLevel | 0.05 | 多重比較に使用する有意水準の基準 |
ANOVAのオプション
| Bonferroni | 集合の数に基づいて修正された α を持つスチューデントの t 分布による平均値比較検定 |
| Duncan | スチューデント化された範囲の分布に基づくリベラルな範囲検定 |
| StudentNewmanKeuls | スチューデント化された範囲の分布に基づく保守的な範囲検定 |
| Tukey | スチューデント化された範囲の分布に基づく平均値比較検定 |
| Dunnett | 最初の1つの集合をコントロール群とした場合の,これに対する他群の平均値の比較検定 |
PostTestsオプションで使用できる検定
ANOVA[onewaydata, PostTests -> Tukey]ANOVA[onewaydata, PostTests -> {Tukey, Bonferroni}, SignificanceLevel -> .01, CellMeans -> False]