分散分析（ANOVA）は測定値の誤差の分析

測定値には誤差を伴います。その誤差を分析する方法が、分散分析（Analysis of Variance; ANOVA）です。

一元配置分散分析（one-way ANOVA）の例として、ある細胞培養用の40℃に設定した振とう機の温度測定を４回繰り返し行いました。

この４回の繰り返し測定を３日間にわたって実施した（初日は実験１、２日目は実験２、３日目は実験３）ところ、以下の結果が得られました。

さて、これらのデータから何が言えるでしょうか。

先ず平均値を見てみると、３回にわたる実験で、それぞれの平均値は概ね４０〜４１度で誤差（ばらつき）が小さいことがわかるでしょう。

しかし、個々の値を見てみると、結構誤差（ばらつき）が大きいことわかります。

特に実験２などは、１回目が４５．５２度で、３回目が３６．２９度とかなり異なっています。

さて、この誤差（ばらつき）を、客観的にどう説明したらよいでしょうか。

実はこれは、

観測値＝真の値＋誤差

という簡単な一次のモデル式で説明することができます。

真の値というのは神のみぞ知る値で、これを仮に振とう機の設定温度40℃としましょう。

しかしながら、実際の測定温度は、温度計の性能、温度計の読み取り、測定環境、当日の室温、測定者の手技など、様々な誤差要因に左右されますので、決して唯一無二の同じ40℃にはなり得ないわけです。

しかしそうも言っておられず、温度があまりにもばらつくと、細胞の培養環境が変わるため、その後実施する細胞実験の精度に悪影響を及ぼす可能性があります。

ではどうすればいいのでしょうか。

先ずは誤差の原因を分析することです。

すると、誤差として、１回の実験のなかで発生する誤差（４回の繰り返しの誤差）と、各実験（実験１〜３）の間の誤差、が考えられます。

すなわち、

観測値＝真の値＋実験間の誤差＋実験内の誤差

と誤差を分解して考えることができるわけです。

もし実験間の誤差が大きければ、各実験日の室温を一定に保とうとか、実験日ごとに測定者を変えない、などの次の対策が考えられます。

また、実験内の誤差が大きければ、例えば測定を朝、昼、夜のように間隔を空けるのではなくある時間内にまとめて４回測定しようとか、電源をオンにしてから１時間以上経って安定してから測定しよう、などの次の対策が考えられます。

測定値は真の値から必ず誤差を伴います。

先ず測定値の誤差の原因を分析してみましょう。

実験間と実験内、それぞれの誤差の要因について考察してみましょう。

分散分析（ANOVA）：全体の誤差平方和を分解

では、実験間誤差と実験内誤差をどのように計算したらよいでしょうか。

まず、データを図のように分解します。

つまり、各実験（各行）の平均値と全部の平均値との差により実験間誤差を推定します。

そして個々の値と各実験（各行）の平均値との差により実験内誤差を推定します。

全体の誤差は、個々の値と全部の平均値との差により推定します。

すると、以下の式が成立します。

全体の誤差平方和＝実験間の誤差平方和＋実験内誤差平方和

そして、帰無仮説：実験間に差がない　とする仮説検定を行います。

以下のようにExcelで簡単に計算し検定結果を算出することができます。

実験間の誤差平方和は、実験間誤差の２乗和で、SUMSQという関数で計算できます。

実験内の誤差平方和は、実験内誤差の２乗和で、これもSUMSQという関数で計算できます。

全体の誤差平方和は、全データの全平均値との差の２乗の和で、これはDEVSQという関数で計算できます。

Excel関数のSUMSQは単なる２乗の和であるのに対し、DEVSQは全平均値を引いた上での２乗の和を計算する関数です。便利なので覚えておきましょう。

そして、平方和を自由度で割って平均平方、すなわち分散を計算し、これらの比がすなわちF値です。

このF値をF分布表に照らして有意かどうか判定することもできますが、FDISTという関数を使えば簡単にp値を計算できます。

結果としてp=0.907となり、帰無仮説：実験間に差がない、が棄却されなかったことになります。

棄却されなかったので、実験間に差があったと積極的に主張することはできないということです。

検定結果にばかり注目がいきがちですが、分散分析では平方和と平均平方（分散）に着目することも重要です。

平方和を見ると、全体が80.31に対して実験間が1.73と、非常に小さいことがわかります。

平均平方（分散）を見ると、実験内が8.73に対して実験間0.86と、非常に小さいことがわかります。

つまり実験内誤差の方が圧倒的に大きいということです。

この結果から、実験日ごとの条件を厳密にコントロールするよりは、毎回の実験の条件をしっかりコントロールすることの方が大切である、という次のアクションに結びつけることができます。

これがいわゆる一元配置分散分析で、実験間、実験内に限らず、一般に繰り返しを伴う群間の誤差の構造を調べるに適した方法です。

先ずは全体の誤差を、実験間誤差と実験内誤差に分解しましょう。

一元配置分散分析は、群間に差がない、という帰無仮説に対する仮説検定です。

検定の結果だけでなく、平方和や分散にも注目し、誤差の構造について考察しましょう。

ChatGPTに聞いてみた

プロンプト

あなたは最高のデータサイエンティストです。
一元配置分散分析について解説してください。
#制約条件：
サルでもわかる平易な日本語で
文字数は300字程度

ChatGPTの回答

一元配置分散分析（ANOVA）は、3つ以上のグループの平均値が、偶然かどうかを調べる統計手法です。例えば、3つの異なる肥料を使った植物の成長を比較する場合に使います。ポイントは、各グループのデータがどれだけ散らばっているか（分散）を見て、グループ間で平均値に違いがあるかを調べることです。全てのグループの平均値が同じなら、肥料による効果はないと考えられます。しかし、一つでも平均値が異なるグループがあれば、少なくとも一つの肥料には効果があると言えます。簡単に言うと、一元配置分散分析は、「グループ間で、本当に差があるのか？」を数学的に判断する方法です。これにより、偶然の結果なのか、実際に違いがあるのかを明らかにできます。