応答曲面計画を分析の分散分析表

分散分析表のすべての統計量の定義と解釈について解説します。

自由度(DF)

総自由度は、データに含まれる情報量のことです。この情報から、未知の母集団のパラメータの値を分析し推定します。総自由度は、サンプルに含まれる観測値の数によって決定されます。項の自由度は、その項が使う情報量を示します。サンプルサイズを大きくすると、母集団に関して提供される情報が増え、総自由度が高くなります。モデルに含める項の数を増やすと情報量が増え、パラメータ推定値の変動性を推定するのに使える自由度が低くなります。

2つの条件が一致すると、Minitabは誤差の自由度を分割します。1つ目の条件は、現在のモデルには含まれていない、データと適合できる項があることです。たとえば、3つ以上の異なる値を持つ連続予測変数がある場合、その予測変数に対して2次項を推定できます。モデルが2次項を含まない場合、データと適合できる項はモデルに含まれていないので、この条件は満たされていることになります。

2つ目の条件は、データに反復が含まれていることです。反復とは、各予測変数の値が同じ観測値のことを言います。たとえば、圧力5、温度25の観測値が3つある場合、それら3つの観測値は反復となります。

2つの条件が一致すると、誤差の自由度の2つの部分は不適合かつ純誤差となります。不適合の自由度は、モデル形式が適切かどうかの検定を可能にします。不適合検定では、不適合に対する自由度が使用されます。純粋誤差の自由度が大きいほど、不適合検定の検定力は高くなります。

調整平方和

調整平方和は、モデル内の異なる部分の変動の測度です。モデル内の予測変数の次数は、調整平方和の計算に影響を及ぼしません。分散分析表では、調整平方和は、異なる要因による変動を説明する成分に分けられます。

調整平方和モデル
モデルの調整平方和は、全体平方和と誤差平方和の差です。モデルに含まれる項のすべての逐次平方和の合計です。
項のグループの調整平方和
項のグループの調整平方和は、項のグループにより説明される応答データの変動量を定量化します。
調整平方和項
項の調整平方和は、他の項だけを持つモデルと比較した場合のモデル平方和の増加を表します。項によって説明される応答データの変動量を定量化します。
調整誤差平方和
調整誤差平方和は残差の平方和です。モデルで説明づけられないデータの変動を定量化します。
純粋誤差の調整平方和
調整純粋誤差平方和は誤差平方和の一部です。純粋誤差平方和は、純粋誤差の自由度があるときに存在します。詳細は、自由度(DF)セクションを参照してください。因子とブロックに同じ値を持つ観測値のデータの変動を定量化します。
調整全体平方和
全体調整平方和は、モデルの平方和と誤差の平方和の合計です。データの変動全体を定量化します。

解釈

Minitabは調整平方和を使用して分散分析表のp値を計算します。また、平方和を使用してR2の統計量も計算します。通常は、平方和ではなく、p値とR2統計量を解釈します。

調整平均平方

調整平均平方は、項やモデルによってどれだけの変動を説明できるかを測定するものです。このとき、その他のすべての項は、入力された順序にかかわらずモデル内に含まれると仮定します。調整平方和と異なり、調整平均平方では、自由度が考慮されます。

調整平均平方誤差(MSEまたはs2)は適合値からの分散です。

解釈

Minitabは調整平均平方を使用して分散分析表のp値を計算します。また、調整平均平方を使用して調整済みR2の統計量も計算します。通常は、調整平均平方ではなく、p値と調整済みR2の統計量を解釈します。

逐次平方和

逐次平方和は、モデル内の異なる部品の変動の測度です。調整平方和と異なり、逐次平方和は項がモデルに追加された順序に依存します。

逐次平方和モデル
モデルの逐次平方和は、モデルに含まれる全体平方和と誤差平方和の差です。モデルに含まれる項のすべての逐次平方和の合計です。
項のグループの逐次平方和
モデルに含まれる項のグループに対する逐次平方和は、グループ内のすべての項に対する逐次平方和の合計です。
逐次平方和項
項の逐次平方和は、分散分析表のモデル上に他の項だけを持つモデルと比較した場合のモデル平方和の増加を表します。
逐次誤差平方和
逐次誤差平方和は残差の平方和です。予測変数では説明できないデータの変動を定量化します。
逐次平方和純粋誤差
逐次純粋誤差平方和は誤差平方和の一部です。純粋誤差平方和は、純粋誤差の自由度があるときに存在します。詳細は、自由度(DF)セクションを参照してください。因子とブロックに同じ値を持つ観測値のデータの変動を定量化します。
全体逐次平方和
全体逐次平方和は、モデルの平方和と誤差の平方和の合計です。データの変動全体を定量化します。

解釈

計画分析時にp値を計算するためには逐次平方和は使用しませんが、回帰モデルの適合または一般線形モデルの適合を使用するときに逐次平方和を使用できます。通常、調整平方和を基に、p値とR2統計量を解釈します。

寄与度

寄与度には、各ソースが応答の総分散に寄与する割合が表示されます。

解釈

割合が高いほど、ソースがより応答変数の変動の原因になっていることを示しています。応答曲面モデルの寄与率はR2と同じです。

F値

F値は分散分析表の各検定に表示されます。

モデルのF値
F値は、モデルの項がブロックと因子項を含む応答と関係があるかどうかを判断するために使用する検定統計量です。
ブロックのF値
F値は、ブロック内のさまざまな条件が応答と関連付けられているかを判断する検定統計量です。
因子項の種類のF値
F値は、項のグループが応答に関連付けられているかを判断する検定統計量です。項のグループの例には、線形効果、二乗項、2元交互作用があります。
個々の項のF値
F値は項が応答に関連付けられているかを判断する検定統計量です。
不適合度検定におけるF値
F値は、モデルが実験に因子が含まれる欠損している項かを判断する検定統計量です。ステップワイズ手順でモデルからブロックが削除される場合、不適合検定には、これらの項も含まれます。

解釈

F値を使用してMinitabで計算されるp値に基づいて、検定の統計的有意性に関する決定を下すことができます。p値は帰無仮説を棄却するための証拠を測定する確率です。確率が低いほど、帰無仮説を棄却する強力な証拠となります。F値が十分に大きい場合、統計的有意性を示します。

F値から帰無仮説を棄却するかどうかを判断するには、F値を棄却限界値と比較します。Minitabで棄却限界値を計算することも、ほとんどの統計に関する書籍に掲載されているF分布表で棄却値を見つけることもできます。Minitabを使用した棄却値の計算に関する詳細は、逆累積分布関数(ICDF)の使用に進み、「ICDFを使用して棄却値を計算」をクリックします。

p値~モデル

p値は帰無仮説を棄却するための証拠を測定する確率です。確率が低いほど、帰無仮説を棄却する強力な証拠となります。

解釈

モデルによって応答の変動を説明できるかどうかを判断するには、モデルのp値と有意水準を比較して帰無仮説を評価します。モデルの帰無仮説は、モデルでは応答の変動は説明できないという仮定です。通常は、有意水準(αまたはアルファとも呼ばれる)として0.05が適切です。0.05の有意水準は、実際にはモデルによって応答の変動は説明できないにも関わらず、説明できると結論付ける可能性が5%であることを示しています。

p値 ≤ α:モデルにより応答での変動が説明されます
p値が有意水準以下の場合は、そのモデルにより応答での変動が説明されると結論付けます。
p値 > α:応答での変動はモデルによって説明されると結論付けるだけの十分な証拠はありません
p値が有意水準より大きい場合、そのモデルにより応答での変動が説明されると結論することはできません。新しいモデルを適合することができます。

p値~ブロック

p値は帰無仮説を棄却するための証拠を測定する確率です。確率が低いほど、帰無仮説を棄却する強力な証拠となります。

ブロックは、異なる条件下で実験が実行された場合に起こりうる差を説明します。たとえば、ある技師が溶接を分析する実験を計画し、すべてのデータを一日では収集できないとします。溶接の質は、相対湿度などの技師では制御できない、日々変わる複数の不確定要素に影響されます。これらの制御できない変数を説明するため、各日で行われた実験を個別のブロックにグループ化します。ブロックは、制御できない変数の効果とエンジニアが分析したい因子の効果が混同されないよう、制御できない変数からの変動性を説明します。Minitabで実行をブロックに割り当てる方法についての詳細は、ブロックとはを参照してください。

解釈

実行間で異なる条件により応答が変化するかを判断するには、ブロックのp値と有意水準を比較して帰無仮説を評価します。この帰無仮説は、異なる条件によって応答は変化しないという仮定です。

通常は、有意水準(αまたはアルファとも呼ばれる)として0.05が適切です。0.05の有意水準は、実際には実行における異なる条件は応答に影響しないにもかかわらず、影響すると結論付ける可能性が5%であることを示しています。

p値 ≤ α: 条件が異なることで応答が変化します
p値が有意水準以下の場合は、条件が異なることによって応答が変化すると結論付けます。
p値 > α: 条件が異なることで応答が変化すると結論付けるだけの十分な証拠はありません
p値が有意水準より大きい場合は、条件が異なることによって応答が変化すると結論付けることはできません。ブロックなしでモデルを適合してみてください。

p値~因子、交互作用、項グループ

p値は帰無仮説を棄却するための証拠を測定する確率です。確率が低いほど、帰無仮説を棄却する強力な証拠となります。

解釈

モデル項が統計的に有意な場合、解釈は項のタイプによって異なります。解釈は以下のとおりです。
  • カテゴリ因子が有意な場合は、すべての水準平均が等しいとは限らないと結論付けることができます。
  • 交互作用項が有意な場合は、因子と応答の間の関係がその項の他の因子に依存すると結論付けることができます。
  • 2次項が有意な場合、応答曲面が曲線を備えていると結論付けられます。

項グループの検定

項のグループが統計的に有意な場合、グループ内の少なくとも1つの項が応答に対して効果を持つと結論付けることができます。モデルに残す項を統計的有意性によって決定する場合、通常は一度に項のグループ全体を取り除くことはしません。個々の項の統計的有意性は、モデルに含まれる項によって変わることがあるためです。

分散分析

要因自由度調整平方和調整平均平方F値p値
モデル10447.76644.77717.610.003
  線形4428.937107.23442.180.000
    材料1181.151181.15171.250.000
    射出圧力1112.648112.64844.310.001
    射出温度173.72573.72529.000.003
    冷却温度161.41261.41224.150.004
  2元交互作用618.8283.1381.230.418
    材料*射出圧力10.3420.3420.130.729
    材料*射出温度10.7780.7780.310.604
    材料*冷却温度14.5654.5651.800.238
    射出圧力*射出温度10.0020.0020.000.978
    射出圧力*冷却温度10.0390.0390.020.906
    射出温度*冷却温度113.10113.1015.150.072
誤差512.7122.542   
合計15460.478     

このモデルにおいては、二元交互作用は水準0.05において統計的に有意ではありません。

分散分析

要因自由度調整平方和調整平均平方F値p値
モデル5442.0488.40847.950.000
  線形4428.94107.23458.160.000
    材料1181.15181.15198.240.000
    射出圧力1112.65112.64861.090.000
    射出温度173.7373.72539.980.000
    冷却温度161.4161.41233.310.000
  2元交互作用113.1013.1017.110.024
    射出温度*冷却温度113.1013.1017.110.024
誤差1018.441.844   
合計15460.48     

最大のp値の二元交互作用から初めて、モデルから項を1つずつ除外して縮約する場合、最後の二元交互作用は水準0.05において統計的に有意です。

p値~不適合度

p値は帰無仮説を棄却するための証拠を測定する確率です。確率が低いほど、帰無仮説を棄却する強力な証拠となります。

解釈

モデルによって応答と予測の関係が正しく指定されるかどうかを判断するには、不適合度検定のp値と有意水準を比較して帰無仮説を評価します。不適合度検定の帰無仮説は、モデルによって応答と予測の関係が正しく指定されるという仮定です。通常は、有意水準(αまたはアルファとも呼ばれる)として0.05が適切です。0.05の有意水準は、実際にはモデルによって応答と予測の関係が正しく指定されるのにも関わらず、正しく指定されないと結論付ける可能性が5%であることを示しています。
p値 ≤ α:その不適合は統計的に有意です
p値が有意水準以下の場合は、そのモデルでは関係が正しく指定されないと結論付けます。モデルを改善するには、項を追加するか、またはデータを変換する必要があります。
p値 > α:その不適合は統計的に有意ではありません

p値が有意水準より大きい場合は、検定で何の不適合も検出されません。