「一元配置実験のデータをエクセルで解析したい」
「エクセルのアドイン機能を使って手軽にやりたい」
「自力で関数を使って計算できるようになりたい」
このような悩みをお持ちの方に向けた記事です。
一元配置実験は一つの因子を変化させて特性値への影響を調べるための実験手法で、実験計画法のなかでも最も基本的なものの一つです。
この記事では、エクセル関数を使った一元配置実験のデータ解析の方法、アドイン機能を使った手軽な解析方法の2種類について詳しく解説しています。
ぜひ最後まで読んで参考にしていただければ幸いです。
一元配置実験とは？
一元配置実験とは一つの因子を変化させて特性値への影響を調べるための実験手法です。
例えば次の実験データのように、温度を一定間隔で変えて繰り返しデータを取得するような実験計画のことを指します。
特性値への影響の度合いとしては、分散を用いて水準どうしを比較し、水準間の分散が水準内の分散よりも十分に大きい場合に、水準による効果が大きいと判断します。
分散の比はF検定を用いることで、定量的な基準で判定することができます。
一元配置実験の基本的な考え方や手順については、以下の記事で詳しく解説していますので、合わせてご覧ください。
アドイン機能を用いたデータ解析の方法
エクセルを使った一元配置実験のデータ解析のやり方としては、アドイン機能を使って手軽に行う方法と、関数を使って自力で解く方法の2つがあります。
まずは、簡単にできるアドイン機能を使った方法から紹介します。
手順
アドイン機能を有効にする
まず、Excelの設定をしましょう。
「ファイル」→「その他」→「オプション」を選択します。
「アドイン」→「設定」を指定し、「分析ツール」にチェックを入れて有効化します。
データをマトリックス状に入力する
次に、一元配置実験で得られたデータをマトリックス状に入力します。
今回の例では、行に水準、列に試行回数を取りましたが、行と列はデータ分析の設定の際に自分で選べるので、行に試行回数、列に水準を取っても問題ありません。
データ範囲と有意水準を設定する
次に、アドイン機能を用いて一元配置実験のデータ解析を行います。
「データ」タブの「データ分析」をクリックし、「分散分析：一元配置」を選択します。
入力範囲に特性値のデータ範囲を設定します。データ方向は水準の配置の向きを表すので、今回の場合は「行」を選択します。
また、F検定の有意水準として$α$を設定します。
デフォルトでは、0.05（5%）が設定されています。
最後に結果の出力先を設定したら、「OK」を押して完了です。
あっという間で便利な機能だね
結果の見かた
各項目の意味
それでは、得られた結果が何を意味するのか、見ていきましょう。
まず、「概要」で示す表は見て分かる通り、基本統計量が算出されています。
「データの個数」は繰り返しの試行回数、「分散」は不偏分散を意味します。
ここまでは、わざわざアドイン機能を使わなくても簡単に求められますね。
次に、本題の「分散分析表」を見てみましょう。
「変動要因」の「グループ間」と「グループ内」というのは、それぞれ水準間と水準内ということを意味しています。
冒頭にも触れた通り、分散分析は水準間の分散と水準内の分散の比によって有意性を検定するものですので、要因を分割して計算する必要があります。
「合計」というのは水準間変動と水準内変動を合わせたもので、全変動を意味します。
「変動」は偏差平方和（平均値と各データの差分の平方和）のことを表します。
水準間変動、水準内変動、全変動の定義と求め方は次のようになります。
水準間変動$S_{A}$：（全体の平均と水準の平均との差）の二乗和
水準内変動$S_{e}$：（各データと水準内平均との差）の二乗和
全変動$S_{T}$：（各データと全体平均との差）の二乗和
全変動$S_{T}$＝水準間変動$S_{A}$＋水準内変動$S_{e}$
また、自由度、分散、観測された分散比は、次のように求められます。
観測された分散比がF検定での検定統計量となります。
結果の判定
アドイン機能では便利なことに、F検定の検定結果まで計算してくれます。
その判定指標となるのが「P-値」と「F境界値」です。
P値は先ほど求めた検定統計量（観測された分散比）に該当する確率のことであり、分散比が大きいほど値が小さくなります。
つまり、今回の場合、水準間変動が大きいとは言えない確率が0.025（=2.5%）しかなく、裏を返すと水準間の影響が十分に大きいと判断できるのです。
判定の基準としては一般的に5%や1%を用いるので、これより値が小さければ、水準の違いがあると考えて良いでしょう。
なお、「F境界値」というのは、最初に設定した有意水準に該当するF値であり、今回の場合は5%に相当する値が出力されています。
分散比がこれよりも大きい値なので、F境界値と分散比の大小関係からも有意性を判断することができます。
大小の比較だけだから簡単
エクセル関数を使ったデータ解析の方法
次にエクセル関数を用いて自力で計算をする場合の手順を説明します。
アドイン機能を用いれば簡単に答えを出せますが、その反面どのような計算から導き出されたものか理解しないまま使ってしまう恐れがあるので、応用が利きません。
計算自体は複雑なものではないので、一度は自力で計算してみることをおススメします。
手順
「概要」の表で出力される項目は、いずれも基本的な関数ばかりなので詳細は割愛します。
次に「分散分析表」ですが、まずは全変動と水準間変動の偏差平方和を求めるところから始めます。
全変動に関しては便利な関数があり、DEVSQ関数を用いて全データ範囲を指定すれば、偏差平方和を一発で求めることができます。
水準間変動については関数がないので、自力で求めるしかありません。
数式が長くなって一見複雑に見えますが、先ほど説明した定義の式を一つずつ分解して足し合わせただけなので、特に難しいものではありません。
水準内変動については、全変動から水準間変動を引けば求められます。
自由度、分散、観測された分散比は、先ほどにも説明した通り関数を使わずに単純な計算で求められます。
次にP-値ですが、これはF分布の上側（右側）確率を返すF.DIST.RT関数を用いれば、あとは分散比と自由度を入れるだけで求められます。
また、F境界値はF分布の上側逆関数を返すF.INV.RT関数を用いれば、有意水準と自由度を入れて求められます。
今回の場合、有意水準を0.05としています。
これで、アドイン機能と同じ表を関数計算で作ることができました。
信頼区間の求め方
最後に信頼区間の求め方を紹介しておきます。
実験データから得られた平均値はあくまで標本平均であり、母集団の平均値（母平均）の推定値に過ぎません。
そのため、母集団の平均値には推定の「幅」があり、例えば信頼度95%で○○～△△の範囲に収まるといった表現をします。
一元配置実験では、次の計算式から信頼区間を求めることができます。
水準内変動の自由度と信頼度に対応するｔ値、水準内変動の不偏分散から計算でき、エクセルの関数を用いて計算すると次のようになります。
T.INV関数はｔ分布の逆関数を返す関数で、区間の信頼度と自由度から求められます。
残念ながら、信頼区間はアドイン機能では計算してくれませんので、自力で導出するしかありません。
ただ、計算自体は複雑なものではないので、計算式だけ覚えておけば特に問題ないことと思います。
関数の使い方を知っておけば応用もバッチリだね
まとめ
自力で計算するのに使える主なエクセル関数を整理しておきます。
- COUNT関数：データの個数
- SUM関数：総和
- VAR.S関数：不偏分散
- DEVSQ関数：偏差平方和
- F.DIST.RT関数：F分布の上側確率
- F.INV.RT関数：F分布の上側逆関数
- T.INV関数：ｔ分布の逆関数
- SQRT関数：平方根
最後までご覧いただきありがとうございました。
数式の苦手な方でも安心の入門編。
回帰分析、検定、実験計画法まで幅広く、エクセルを用いた実践方法も習得したい方に。
乱塊法、分割法もお任せ！ QC検定1級の教材としてもおススメ。
コメント