Example of Analyze Binary Response for Factorial Design

A food scientist is studying factors that affect food spoilage. The scientist uses a 2-level factorial experiment to assess several factors that could impact the rate of food spoiling.

The scientist analyzes a 2-level factorial design to determine how preservative type, vacuum packaging pressure, contamination level, and cooling temperature affect the spoilage of fruit. The response is binary—whether spoilage is detected or not—in a sample of 500 containers of fruit.

  1. Open the sample data, FoodSpoilage.MTW.
  2. Choose Stat > DOE > Factorial > Analyze Binary Response.
  3. In Event name, enter Event.
  4. In Number of events, enter Spoilage.
  5. In Number of trials, enter Containers.
  6. Click Terms.
  7. Under Include terms in the model up through order, choose 2.
  8. Click OK in each dialog box.

Interpret the results

In the Deviance table, the p-values for three of the main effect terms—Preservative, VacuumPress, and ContaminationLevel—are significant. Because the p-values are less than the significance level of 0.05, the scientist concludes that these factors are statistically significant. None of the two-way interactions are significant. The scientist can consider reducing the model.

The Deviance R2 value shows that the model explains 97.95% of the total deviance in the response, which indicates that the model fits the data well.

Most of the VIFs are small, which indicates that the terms in the model are not correlated.

The Pareto plot of the effects allow you to visually identify the important effects and compare the relative magnitude of the various effects. In these results, three main effects are statistically significant (α = 0.05) - preservative type (A), vacuum seal pressure (B), and contamination level (C). In addition, you can see that the largest effect is preservative type (A) because it extends the farthest. The effect for the preservative by cooling temperature interaction (AD) is the smallest because it extends the least.

Method

Link functionLogit
Rows used16

Response Information

VariableValueCountEvent Name
SpoilageEvent506Event
  Non-event7482 
ContainersTotal7988 

Coded Coefficients

TermEffectCoefSE CoefVIF
Constant  -2.73700.0479 
Preservative0.44970.22490.04771.03
VacuumPress0.25740.12870.04771.06
ContaminationLevel0.29540.14770.04781.06
CoolTemp-0.1107-0.05540.04781.07
Preservative*VacuumPress-0.0233-0.01170.04731.05
Preservative*ContaminationLevel0.07220.03610.04741.06
Preservative*CoolTemp0.00670.00340.04721.05
VacuumPress*ContaminationLevel-0.0430-0.02150.04691.04
VacuumPress*CoolTemp-0.0115-0.00580.04651.02
ContaminationLevel*CoolTemp0.15730.07860.04671.02

Odds Ratios for Continuous Predictors

Unit of
Change
Odds Ratio95% CI
VacuumPress10.0*(*, *)
ContaminationLevel22.5*(*, *)
CoolTemp5.0*(*, *)
Odds ratios are not calculated for predictors that are included in interaction terms because
     these ratios depend on values of the other predictors in the interaction terms.

Odds Ratios for Categorical Predictors

Level ALevel BOdds Ratio95% CI
Preservative     
  Any levelAny level*(*, *)
Odds ratio for level A relative to level B
Odds ratios are not calculated for predictors that are included in interaction terms because
     these ratios depend on values of the other predictors in the interaction terms.

Model Summary

Deviance
R-Sq
Deviance
R-Sq(adj)
AICAICcBIC
97.95%76.75%105.98171.98114.48

Goodness-of-Fit Tests

TestDFChi-SquareP-Value
Deviance50.970.965
Pearson50.970.965
Hosmer-Lemeshow60.101.000

Analysis of Variance

SourceDFAdj DevAdj MeanChi-SquareP-Value
Model1046.21304.621346.210.000
  Preservative122.683522.683522.680.000
  VacuumPress17.33137.33137.330.007
  ContaminationLevel19.62099.62099.620.002
  CoolTemp11.34411.34411.340.246
  Preservative*VacuumPress10.06080.06080.060.805
  Preservative*ContaminationLevel10.57800.57800.580.447
  Preservative*CoolTemp10.00510.00510.010.943
  VacuumPress*ContaminationLevel10.21060.21060.210.646
  VacuumPress*CoolTemp10.01530.01530.020.902
  ContaminationLevel*CoolTemp12.84752.84752.850.092
Error50.96740.1935   
Total1547.1804     

Regression Equation in Uncoded Units

P(Event)=exp(Y')/(1 + exp(Y'))
Y'=-2.721 + 0.188 Preservative + 0.0172 VacuumPress - 0.00249 ContaminationLevel
- 0.0286 CoolTemp - 0.00117 Preservative*VacuumPress
+ 0.00160 Preservative*ContaminationLevel + 0.00067 Preservative*CoolTemp
- 0.000096 VacuumPress*ContaminationLevel - 0.000115 VacuumPress*CoolTemp
+ 0.000699 ContaminationLevel*CoolTemp