Interpret the key results for Fit Binary Logistic Model

Complete the following steps to interpret a binary logistic model. Key output includes the p-value, the coefficients, R2, and the goodness-of-fit tests.

Step 1: Determine whether the association between the response and the term is statistically significant

To determine whether the association between the response and each term in the model is statistically significant, compare the p-value for the term to your significance level to assess the null hypothesis. The null hypothesis is that there is no association between the term and the response. Usually, a significance level (denoted as α or alpha) of 0.05 works well. A significance level of 0.05 indicates a 5% risk of concluding that an association exists when there is no actual association.
P-value ≤ α: The association is statistically significant
If the p-value is less than or equal to the significance level, you can conclude that there is a statistically significant association between the response variable and the term.
P-value > α: The association is not statistically significant
If the p-value is greater than the significance level, you cannot conclude that there is a statistically significant association between the response variable and the term. You may want to refit the model without the term.
If there are multiple predictors without a statistically significant association with the response, you can reduce the model by removing terms one at a time. For more information on removing terms from the model, go to Model reduction.
If a model term is statistically significant, the interpretation depends on the type of term. The interpretations are as follows:
  • If a continuous predictor is significant, you can conclude that the coefficient for the predictor is different from zero.
  • If a categorical predictor is significant, you can conclude that not all of the levels of the factor have the same probability.
  • If an interaction term is significant, you can conclude that the relationship between a predictor and the probability of the event depends on the other predictors in the term.
  • If a polynomial term is significant, you can conclude that the relationship between a predictor and the probability of the event depends on the magnitude of the predictor.

Analysis of Variance



Wald Test
SourceDFChi-SquareP-Value
Regression17.830.005
  Dose (mg)17.830.005

Coefficients

TermCoefSE CoefZ-ValueP-ValueVIF
Constant-5.251.99-2.640.008 
Dose (mg)3.631.302.800.0051.00
Key Results: P-Value, Coefficients

In these results, the dosage is statistically significant at the significance level of 0.05. You can conclude that changes in the dosage are associated with changes in the probability that the event occurs.

Assess the coefficient to determine whether a change in a predictor variable makes the event more likely or less likely. The relationship between the coefficient and the probability depends on several aspects of the analysis, including the link function. Generally, positive coefficients indicate that the event becomes more likely as the predictor increases. Negative coefficients indicate that the event becomes less likely as the predictor increases. For more information, go to Coefficients and regression equation for Fit Binary Logistic Model.

The coefficient for Dose is 3.63, which suggests that higher dosages are associated with higher probabilities that the event will occur.

If an interaction term is statistically significant, the relationship between a predictor and the response differs by the level of the other predictor. In this case, you should not interpret the main effects without considering the interaction effect. To obtain a better understanding of the main effects, interaction effects, and curvature in your model, go to Factorial Plots and Response Optimizer.

Step 2: Understand the effects of the predictors

Use the odds ratio to understand the effect of a predictor. The interpretation of the odds ratio depends on whether the predictor is categorical or continuous. Minitab calculates odds ratios when the model uses the logit link function.
Odds Ratios for Continuous Predictors

Odds ratios that are greater than 1 indicate that the event is more likely to occur as the predictor increases. Odds ratios that are less than 1 indicate that the event is less likely to occur as the predictor increases.

Odds Ratios for Continuous Predictors

Unit of
Change
Odds Ratio95% CI
Dose (mg)0.56.1279(1.7218, 21.8087)
Key Result: Odds Ratio

In these results, the model uses the dosage level of a medicine to predict the presence or absence of bacteria in adults. In this example, the absence of bacteria is the Event. Each pill contains a 0.5 mg dose, so the researchers use a unit change of 0.5 mg. The odds ratio is approximately 6. For each additional pill that an adult takes, the odds that a patient does not have the bacteria increase by about 6 times.

Odds Ratios for Categorical Predictors

For categorical predictors, the odds ratio compares the odds of the event occurring at 2 different levels of the predictor. Minitab sets up the comparison by listing the levels in 2 columns, Level A and Level B. Level B is the reference level for the factor. Odds ratios that are greater than 1 indicate that the event is more likely at level A. Odds ratios that are less than 1 indicate that the event is less likely at level A. For information on coding categorical predictors, go to Coding schemes for categorical predictors.

Binary Logistic Regression: Cancellation versus Month

Odds Ratios for Categorical Predictors

Level ALevel BOdds Ratio95% CI
Month     
  211.1250(0.0600, 21.0834)
  313.3750(0.2897, 39.3165)
  417.7143(0.7461, 79.7592)
  512.2500(0.1107, 45.7172)
  616.0000(0.5322, 67.6397)
  323.0000(0.2547, 35.3325)
  426.8571(0.6556, 71.7169)
  522.0000(0.0976, 41.0019)
  625.3333(0.4679, 60.7946)
  432.2857(0.4103, 12.7323)
  530.6667(0.0514, 8.6389)
  631.7778(0.2842, 11.1200)
  540.2917(0.0252, 3.3719)
  640.7778(0.1464, 4.1326)
  652.6667(0.2124, 33.4861)
Odds ratio for level A relative to level B
Key Result: Odds Ratio

In these results, the categorical predictor is the month from the start of a hotel's busy season. The response is whether or not a guest cancels a reservation. In this example, a cancellation is the Event. The largest odds ratio is approximately 7.71, when level A is month 4 and level B is month 1. This indicates that the odds that a guest cancels a reservation in month 4 is approximately 8 times higher than the odds that a guest cancels a reservation in month 1.

For more information, go to Odds Ratios for Fit Binary Logistic Model.

Step 3: Determine how well the model fits your data

To determine how well the model fits your data, examine the statistics in the Model Summary table.

Many of the model summary and goodness-of-fit statistics are affected by how the data are arranged in the worksheet and whether there is one trial per row or multiple trials per row. The Hosmer-Lemeshow test is unaffected by the data format and is comparable between formats. For more information, go to How data formats affect goodness-of-fit in binary logistic regression.

Deviance R-sq

The higher the deviance R2, the better the model fits your data. Deviance R2 is always between 0% and 100%.

Deviance R2 always increases when you add additional predictors to a model. For example, the best 5-predictor model will always have an R2 that is at least as high as the best 4-predictor model. Therefore, deviance R2 is most useful when you compare models of the same size.

For binary logistic regression, the format of the data affects the deviance R2 value. The deviance R2 is usually higher for data in Event/Trial format. Deviance R2 values are comparable only between models that use the same data format.

Goodness-of-fit statistics are just one measure of how well the model fits the data. Even when a model has a desirable value, you should check the residual plots and goodness-of-fit tests to assess how well a model fits the data.

Deviance R-sq (adj)

Use adjusted deviance R2 to compare models that have different numbers of predictors. Deviance R2 always increases when you add a predictor to the model. The adjusted deviance R2 value incorporates the number of predictors in the model to help you choose the correct model.

AIC, AICc, and BIC

Use AIC, AICc, and BIC to compare different models. For each statistic, smaller values are desirable. However, the model with the smallest value for a set of predictors does not necessarily fit the data well. Also use goodness-of-fit tests and residual plots to assess how well a model fits the data.

Area Under ROC Curve

The area under the ROC curve values range from 0.5 to 1. When the binary model can perfectly separate the classes, then the area under the curve is 1. When the binary model cannot separate the classes better than a random assignment, then the area under the curve is 0.5.

Model Summary

Deviance
R-Sq
Deviance
R-Sq(adj)
AICAICcBICArea Under
ROC Curve
96.04%91.81%10.6314.6310.220.9398
Key Results: Deviance R-Sq, Deviance R-Sq (adj), AIC, Area Under ROC Curve

In these results, the model explains 96.04% of the total deviance in the response variable. For these data, the Deviance R2 value indicates the model provides a good fit to the data. The area under the ROC curve is 0.9398. This value indicates that the model classifies much of the data correctly. If additional models are fit with different predictors, use the adjusted Deviance R2 value, the AIC value, the AICc value, the BIC value, and the area under the ROC curve to compare how well the models fit the data.

Step 4: Determine whether the model does not fit the data

Use the goodness-of-fit tests to determine whether the predicted probabilities deviate from the observed probabilities in a way that the binomial distribution does not predict. If the p-value for the goodness-of-fit test is lower than your chosen significance level, the predicted probabilities deviate from the observed probabilities in a way that the binomial distribution does not predict. This list provides common reasons for the deviation:
  • Incorrect link function
  • Omitted higher-order term for variables in the model
  • Omitted predictor that is not in the model
  • Overdispersion

If the deviation is statistically significant, you can try a different link function or change the terms in the model.

For binary logistic regression, the format of the data affects the p-value because it changes the number of trials per row.

  • Deviance: The p-value for the deviance test tends to be lower for data that are in the Binary Response/Frequency format compared to data in the Event/Trial format. For data in Binary Response/Frequency format, the Hosmer-Lemeshow results are more trustworthy.
  • Pearson: The approximation to the chi-square distribution that the Pearson test uses is inaccurate when the expected number of events per row in the data is small. Thus, the Pearson goodness-of-fit test is inaccurate when the data are in Binary Response/Frequency format.
  • Hosmer-Lemeshow: The Hosmer-Lemeshow test does not depend on the number of trials per row in the data as the other goodness-of-fit tests do. When the data have few trials per row, the Hosmer-Lemeshow test is a more trustworthy indicator of how well the model fits the data.

Response Information

VariableValueCountEvent Name
EventEvent160Event
  Non-event340 
TrialTotal500 

Goodness-of-Fit Tests

TestDFChi-SquareP-Value
Deviance23.780.151
Pearson23.760.152
Hosmer-Lemeshow33.760.288
Key Results for Event/Trial Format: Response Information, Deviance Test, Pearson Test, Hosmer-Lemeshow Test

In these results, the Response Information table shows Event and Trial in the Variable column. These labels indicate that the data are in Event/Trial format. All of the goodness-of-fit tests have p-values higher than the usual significance level of 0.05. The tests do not provide evidence that the predicted probabilities deviate from the observed probabilities in a way that the binomial distribution does not predict.

Response Information

VariableValueCount
YEvent160(Event)
  Non-event340 
  Total500 

Goodness-of-Fit Tests

TestDFChi-SquareP-Value
Deviance497552.030.044
Pearson497504.420.399
Hosmer-Lemeshow33.760.288
Key Results for Binary Response/Frequency Format: Response Information, Deviance Test, Pearson Test, Hosmer-Lemeshow Test

In these results for the same data, the Response Information table shows Y in the variable column. This label indicates that the data are in Binary Response/Frequency format. The deviance test has a p-value less than the usual significance level of 0.05, but the Hosmer-Lemeshow test is the most trustworthy test. The Hosmer-Lemeshow test does not provide evidence that the predicted probabilities deviate from the observed probabilities in a way that the binomial distribution does not predict.