要因計画の2値応答の分析の例

食品科学者が、食品の腐敗に影響する要因を研究するとします。2水準要因実験を使用して、食品の腐敗率に影響を及ぼす複数の要因を評価します。

2水準要因計画を分析して、防腐剤の種類、真空封止圧力、汚染度、冷却温度がどのように果物の腐敗に影響を及ぼすかを判定します。800個の果物容器の標本では、腐敗の有無にかかわらず、応答は2値です。

  1. 標本データを開く、食品の損傷.MTW.
  2. 統計 > 実験計画法(DOE) > 要因計画 > 2値応答の分析を選択します。
  3. 事象名(V)で、事象を入力します。
  4. 事象数損傷を入力します。
  5. 試行回数コンテナを入力します。
  6. をクリックします。
  7. モデルに含める項の次数2を選択します。
  8. 各ダイアログボックスでOKをクリックします。

結果を解釈する

逸脱度表で、3つの主効果の項である防腐剤、真空プレス、汚染度のp値は有意です。p値は有意水準0.05より小さいため、食品科学者はこれらの要因を統計的に有意だと結論付けます。二元交互作用のどれも有意ではありません。食品科学者はモデルの縮約を検討できます。

逸脱R2値は、モデルが応答における総逸脱の97.95%を説明していることを示しており、モデルが良好にデータに適合することを示しています。

ほとんどの分散拡大係数(VIF)は小さく、モデルに含まれる項は相関していないことを示しています。

効果のパレート図を使用すると、重要な効果を視覚的に識別して、さまざまな効果の相対的重要度を比較することができます。この結果では、防腐剤の種類(A)、真空封止圧力(B)、汚染度(C)の3つの主効果は統計的に有意です(α = 0.05)。また、防腐剤の種類(A)は最も長く伸びているため、これが最大の効果を持つこともわかります。冷却温度相互作用(AD)による防腐剤の効果の伸びは最も短いため、効果は最小です。

方法

リンク関数Logit
使用中の行16

応答情報

変数計数事象名
損傷事象506Event
  非事象7482 
コンテナ合計7988 

コード化係数

効果係数係数の標準誤差VIF
定数  -2.73700.0479 
防腐剤0.44970.22490.04771.03
真空プレス0.25740.12870.04771.06
汚染レベル0.29540.14770.04781.06
冷却温度-0.1107-0.05540.04781.07
防腐剤*真空プレス-0.0233-0.01170.04731.05
防腐剤*汚染レベル0.07220.03610.04741.06
防腐剤*冷却温度0.00670.00340.04721.05
真空プレス*汚染レベル-0.0430-0.02150.04691.04
真空プレス*冷却温度-0.0115-0.00580.04651.02
汚染レベル*冷却温度0.15730.07860.04671.02

連続予測変数のオッズ比

変更ユニットオッズ比95%信頼区間
真空プレス10.0*(*, *)
汚染レベル22.5*(*, *)
冷却温度5.0*(*, *)
交互作用項に含まれる予測変数のオッズ比は、交互作用項の他の予測変数の値に依存するため、計算されません。

カテゴリ予測変数のオッズ比

水準A水準Bオッズ比95%信頼区間
防腐剤     
  任意の水準任意の水準*(*, *)
水準Bと比べた水準Aのオッズ比
交互作用項に含まれる予測変数のオッズ比は、交互作用項の他の予測変数の値に依存するため、計算されません。

モデル要約

逸脱
(deviance)
R二乗
逸脱
(deviance)
R二乗 (調整済み)
AICAICc(修正済み
赤池情報量基準)
BIC(ベイズ
情報量基準)
97.95%76.75%105.98171.98114.48

適合度検定

検定自由度カイ二乗p値
逸脱 (deviance)50.970.965
ピアソン50.970.965
Hosmer-Lemeshow60.101.000

分散分析

要因自由度調整偏差調整平均カイ二乗p値
モデル1046.21304.621346.210.000
  防腐剤122.683522.683522.680.000
  真空プレス17.33137.33137.330.007
  汚染レベル19.62099.62099.620.002
  冷却温度11.34411.34411.340.246
  防腐剤*真空プレス10.06080.06080.060.805
  防腐剤*汚染レベル10.57800.57800.580.447
  防腐剤*冷却温度10.00510.00510.010.943
  真空プレス*汚染レベル10.21060.21060.210.646
  真空プレス*冷却温度10.01530.01530.020.902
  汚染レベル*冷却温度12.84752.84752.850.092
誤差50.96740.1935   
合計1547.1804     

非コード化単位の回帰式

P(Event)=exp(Y')/(1 + exp(Y'))
Y'=-2.721 + 0.188 防腐剤 + 0.0172 真空プレス - 0.00249 汚染レベル - 0.0286 冷却温度 - 0.00117 防腐剤*真空プレス
+ 0.00160 防腐剤*汚染レベル + 0.00067 防腐剤*冷却温度 - 0.000096 真空プレス*汚染レベル - 0.000115 真空プレス*冷却温度
+ 0.000699 汚染レベル*冷却温度