Example of Fit Poisson Model

A quality engineer is concerned about two types of defects in molded resin parts: discoloration and clumping. Discolored streaks in the final product can result from contamination in hoses and from abrasions to resin pellets. Clumping can occur when the process is run at higher temperatures and faster rates of transfer. The engineer identifies three possible predictor variables for the responses (defects). The engineer records the number of each type of defect in hour long sessions, while varying the predictor levels.

The engineer wants to study how several predictors affect discoloration defects in resin parts. Because the response variable describes the number of times that an event occurs in a finite observation space, the engineer fits a Poisson model.

  1. Enter the sample data, ResinDefects.MTW.
  2. Choose Stat > Regression > Poisson Regression > Fit Poisson Model.
  3. In Response, enter 'Discoloration Defects'.
  4. In Continuous predictors, enter 'Hours Since Cleanse' Temperature.
  5. In Categorical predictors, enter 'Size of Screw'.
  6. Click Graphs.
  7. In Residuals for plots, select Standardized.
  8. Under Residuals plots, select Four in one.
  9. Click OK in each dialog box.

Interpret the results

The plot of the standardized deviance residuals versus the fitted values shows a distinct curve. In the plot of the residuals versus order, the residuals in the middle tend to be higher than the residuals at the beginning and end of the data set. For these data, both patterns are because of a missing interaction term between the size of the screw and the temperature. The pattern is visible on the residuals versus order plot because the engineer did not collect the data in random order. The engineer refits the model with the interaction between temperature and the size of the screw to model the defects more accurately.

Method

Link functionNatural log
Categorical predictor coding(1, 0)
Rows used36

Regression Equation

Discoloration Defects=exp(Y')
Size of
Screw
largeY'=4.398 + 0.01798 Hours Since Cleanse - 0.001974 Temperature
       
smallY'=4.244 + 0.01798 Hours Since Cleanse - 0.001974 Temperature

Coefficients

TermCoefSE CoefZ-ValueP-ValueVIF
Constant4.39820.062870.020.000 
Hours Since Cleanse0.017980.008262.180.0291.00
Temperature-0.0019740.000318-6.200.0001.00
Size of Screw         
  small-0.15460.0427-3.620.0001.00

Model Summary

Deviance
R-Sq
Deviance
R-Sq(adj)
AICAICcBIC
64.20%60.80%253.29254.58259.62

Goodness-of-Fit Tests

TestDFEstimateMeanChi-SquareP-Value
Deviance3231.607220.9877331.610.486
Pearson3231.267130.9771031.270.503

Analysis of Variance



Wald Test
SourceDFChi-SquareP-Value
Regression356.290.000
  Hours Since Cleanse14.740.029
  Temperature138.460.000
  Size of Screw113.090.000

Fits and Diagnostics for Unusual Observations

ObsDiscoloration
Defects
FitResidStd Resid
3343.0058.18-2.09-2.18R
R  Large residual
  1. Press Ctrl+E, or click the Edit Last Dialog button on the Standard toolbar.
  2. Click Model.
  3. In Predictors, select Temperature and 'Size of Screw'.
  4. Next to Interactions through order, choose 2 and click Add.
  5. Click OK in each dialog box.

For the model with the interaction, the AIC is approximately 236, which is lower than the model without the interaction. The AIC criterion indicates that the model with the interaction is better than the model without the interaction. The curvature in the residuals versus fits plot is gone. The engineer notices that some coefficients have VIF values that are > 5. In this case, an analysis with standardized continuous predictors to reduce the effect of collinearity gives the same conclusions about the statistical significance of the terms in the model. (For more information, go to Multicollinearity in regression.) The engineer decides to interpret this model rather than the model without the interaction.

Method

Link functionNatural log
Categorical predictor coding(1, 0)
Rows used36

Regression Equation

Discoloration Defects=exp(Y')
Size of
Screw
largeY'=4.576 + 0.01798 Hours Since Cleanse - 0.003285 Temperature
       
smallY'=4.032 + 0.01798 Hours Since Cleanse - 0.000481 Temperature

Coefficients

TermCoefSE CoefZ-ValueP-ValueVIF
Constant4.57600.073662.150.000 
Hours Since Cleanse0.017980.008262.180.0291.00
Temperature-0.0032850.000441-7.460.0001.92
Size of Screw         
  small-0.54440.0990-5.500.0005.37
Temperature*Size of Screw         
  small0.0028040.0006404.380.0006.64

Model Summary

Deviance
R-Sq
Deviance
R-Sq(adj)
AICAICcBIC
85.99%81.46%236.05238.05243.97

Goodness-of-Fit Tests

TestDFEstimateMeanChi-SquareP-Value
Deviance3112.365980.3989012.370.999
Pearson3112.316110.3972912.320.999

Analysis of Variance



Wald Test
SourceDFChi-SquareP-Value
Regression478.770.000
  Hours Since Cleanse14.740.029
  Temperature155.600.000
  Size of Screw130.210.000
  Temperature*Size of Screw119.170.000