分析の手順Rで分散分析を行う際には,以下の3つのステップを経る。対象となるデータを変数(オブジェクト)に格納する統計分析関数に,そのデータを渡す結果を解釈する前述のt検定において,2の実行と同時に結果が画面に表示された。Rにおける分散分析では,結果を解釈するステップが加わる。詳しいことは後述する。サンプルデータこのページでは,調理中のドーナツが吸収する脂肪分の量を分析する。4種類の脂肪質を使い,それぞれ6個のドーナツを調理して,その脂肪吸収量を測定した。脂肪質を効果(4水準)として,一元配置の分散分析を実施してみよう。一元配置の分散分析も,前述の検定と同様の手順で行う。ファイルのデータを読み取って変数に格納し,分析関数を呼び出せばよい。ところが,データファイルの作成にあたり,注意しなければならないことがある。1行につき1個の測定値とするその行には,観測値のほか,処置群(効果)を識別する記号を含める詳しくは,サンプルのデータファイルを参照してほしい。グループを識別する記号は何でもよい。この例では「A,B,C,D」を用いたが,「1,2,3,4」あるいは「FAT1,FAT2,FAT3,FAT4」でも問題ない。データの読み込みサンプルファイルは,ヘッダ付きのカンマ区切り(CSV)形式で保存されている。するが,変数名が長くて入力に不便なので,別の変数に代入しておく。分かりやすく,入力しやすい変数名にしておくとよい。下の例では,短い1文字の変数に置き換えた。分散分析の実行Rで分散分析を行うには,aov()という関数を使用する。観測値に,どのような要因が関与しているかを与える必要がある。unbalanced分散分析は,無事に完了し,いくつかの数値が表示されている。残念ながら,これだけでは分散分析表は表示されない。以下のように入力する。それについて,詳しい結果を表示するこの2段階の入力が面倒ならば,以下のように全体をsummary()で囲んでもよい。上と同じ内容が表示されるなお,同じような働きをするコマンドとして,summary()がある。分析が成功すると,以下のように表示される。これは,分散分析表である。必要な部分を日本語で書き直して整形すると,以下のようになる。-------------------------------------------------------要因Fについて,P値が0.05を下回っていることから,4水準の平均値は有意差が認められる。obihiro.ac.jp---./index.html。