Vars indicates the number of predictors contained in each model. For each number of predictors, by default, Minitab selects the two models with the highest R^{2 }values. The right side of the table indicates which predictors are in the models indicated by "X" symbols.
Vars | R-Sq | R-Sq (adj) | R-Sq (pred) | Mallows Cp | S | I n s o l a t i o n | E a s t | S o u t h | N o r t h | T i m e o f D a y |
---|---|---|---|---|---|---|---|---|---|---|
1 | 72.1 | 71.0 | 66.9 | 38.5 | 12.328 | X | ||||
1 | 39.4 | 37.1 | 26.3 | 112.7 | 18.154 | X | ||||
2 | 85.9 | 84.8 | 81.4 | 9.1 | 8.9321 | X | X | |||
2 | 82.0 | 80.6 | 74.2 | 17.8 | 10.076 | X | X | |||
3 | 87.4 | 85.9 | 79.0 | 7.6 | 8.5978 | X | X | X | ||
3 | 86.5 | 84.9 | 81.4 | 9.7 | 8.9110 | X | X | X | ||
4 | 89.1 | 87.3 | 80.6 | 5.8 | 8.1698 | X | X | X | X | |
4 | 88.0 | 86.0 | 79.3 | 8.2 | 8.5550 | X | X | X | X | |
5 | 89.9 | 87.7 | 78.8 | 6.0 | 8.0390 | X | X | X | X | X |
R^{2} is the percentage of variation in the response that is explained by the model. It is calculated as 1 minus the ratio of the error sum of squares (which is the variation that is not explained by model) to the total sum of squares (which is the total variation in the model).
Use R^{2} to determine how well the model fits your data. The higher the R^{2} value, the better the model fits your data. R^{2} is always between 0% and 100%.
R^{2} always increases when you add additional predictors to a model. For example, the best five-predictor model will always have an R^{2} that is at least as high as the best four-predictor model. Therefore, R^{2} is most useful when you compare models of the same size.
Small samples do not provide a precise estimate of the strength of the relationship between the response and predictors. For example, if you need R^{2} to be more precise, you should use a larger sample (typically, 40 or more).
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 to verify that the model meets the model assumptions.
Adjusted R^{2} is the percentage of the variation in the response that is explained by the model, adjusted for the number of predictors in the model relative to the number of observations. Adjusted R^{2} is calculated as 1 minus the ratio of the mean square error (MSE) to the mean square total (MS Total).
Use adjusted R^{2} when you want to compare models that have different numbers of predictors. R^{2} always increases when you add a predictor to the model, even when there is no real improvement to the model. The adjusted R^{2} value incorporates the number of predictors in the model to help you choose the correct model.
Model | % Potato | Cooling rate | Cooking temp | R^{2} | Adjusted R^{2} |
---|---|---|---|---|---|
1 | X | 52% | 51% | ||
2 | X | X | 63% | 62% | |
3 | X | X | X | 65% | 62% |
The first model yields an R^{2} of more than 50%. The second model adds cooling rate to the model. Adjusted R^{2} increases, which indicates that cooling rate improves the model. The third model, which adds cooking temperature, increases the R^{2} but not the adjusted R^{2}. These results indicate that cooking temperature does not improve the model. Based on these results, you consider removing cooking temperature from the model.
The prediction error sum of squares (PRESS) is a measure of the deviation between the fitted values and the observed values. PRESS is similar to the sum of squares of the residual error (SSE), which is the summation of the squared residuals. However, PRESS uses a different calculation for the residuals. The formula used to calculate PRESS is equivalent to a process of systematically removing each observation from the data set, estimating the regression equation, and determining how well the model predicts the removed observation.
Use PRESS to assess your model's predictive ability. Usually, the smaller the PRESS value, the better the model's predictive ability. Minitab uses PRESS to calculate the predicted R^{2}, which is usually more intuitive to interpret. Together, these statistics can prevent over-fitting the model. An over-fit model occurs when you add terms for effects that are not important in the population, although they may appear important in the sample data. The model becomes tailored to the sample data and therefore, may not be useful for making predictions about the population.
Predicted R^{2} is calculated with a formula that is equivalent to systematically removing each observation from the data set, estimating the regression equation, and determining how well the model predicts the removed observation. The value of predicted R^{2} ranges between 0% and 100%. (While the calculations for predicted R^{2} can produce negative values, Minitab displays zero for these cases.)
Use predicted R^{2} to determine how well your model predicts the response for new observations. Models that have larger predicted R^{2} values have better predictive ability.
A predicted R^{2} that is substantially less than R^{2} may indicate that the model is over-fit. An over-fit model occurs when you add terms for effects that are not important in the population. The model becomes tailored to the sample data and, therefore, may not be useful for making predictions about the population.
Predicted R^{2} can also be more useful than adjusted R^{2} for comparing models because it is calculated with observations that are not included in the model calculation.
For example, an analyst at a financial consulting company develops a model to predict future market conditions. The model looks promising because it has an R^{2} of 87%. However, the predicted R^{2} is 52%, which indicates that the model may be over-fit.
Mallows' Cp can help you choose between competing multiple regression models. Mallows' Cp compares the full model to models with the subsets of predictors. It helps you strike an important balance with the number of predictors in the model. A model with too many predictors can be relatively imprecise while a model with too few predictors can produce biased estimates. Using Mallows' Cp to compare regression models is only valid when you start with the same complete set of predictors.
A Mallows' Cp value that is close to the number of predictors plus the constant indicates that the model produces relatively precise and unbiased estimates.
A Mallows' Cp value that is greater than the number of predictors plus the constant indicates that the model is biased and does not fit the data well.
S represents the standard deviation of the distance between the data values and the fitted values. S is measured in the units of the response.
Use S to assess how well the model describes the response. S is measured in the units of the response variable and represents how far the data values fall from the fitted values. The lower the value of S, the better the model describes the response. However, a low S value by itself does not indicate that the model meets the model assumptions. You should check the residual plots to verify the assumptions.
For example, you work for a potato chip company that examines the factors that affect the percentage of crumbled potato chips per container. You reduce the model to the significant predictors, and S is calculated as 1.79. This result indicates that the standard deviation of the data points around the fitted values is 1.79. If you are comparing models, values that are lower than 1.79 indicate a better fit, and higher values indicate a worse fit.
The corrected Akaike’s Information Criterion (AICc) and the Bayesian Information Criterion (BIC) are measures of the relative quality of a model that account for fit and the number of terms in the model.
Use AICc and BIC to compare different models. Smaller values are desirable. However, the model with the least value for a set of predictors does not necessarily fit the data well. Also use tests and residual plots to assess how well the model fits the data.
Both AICc and BIC assess the likelihood of the model and then apply a penalty for adding terms to the model. The penalty reduces the tendency to overfit the model to the sample data. This reduction can yield a model that performs better in general.
As a general guideline, when the number of parameters is small relative to the sample size, BIC has a larger penalty for the addition of each parameter than AICc. In these cases, the model that minimizes BIC tends to be smaller than the model that minimizes AICc.
In some common cases, such as screening designs, the number of parameters is usually large relative to the sample size. In these cases, the model that minimizes AICc tends to be smaller than the model that minimizes BIC. For example, for a 13-run definitive screening design, the model that minimizes AICc will tend to be smaller than the model that minimizes BIC among the set of models with 6 or more parameters.
For more information on AICc and BIC, see Burnham and Anderson.^{1}
The condition number measures the multicollinearity among model terms. Multicollinearity is a condition that occurs when some terms in the model are correlated with other terms. When you compare models, a smaller condition number is better.
Use the condition number to compare models with different terms. A condition number of 1 indicates that the model terms are uncorrelated. Larger values indicate more multicollinearity.
While the consequences of a certain condition number depend on several conditions, values greater than 100 usually indicate the need for investigation. When terms in the model are multicollinear, the interpretation of the model is less straightforward than the interpretation of a model with uncorrelated terms. For more information, go to Multicollinearity in regression.
In these results, the condition number is 1 when the model has only 1 term. (The condition number is always 1 when the model has 1 continuous predictor.) None of the models have a condition number greater than 100, so the multicollinearity among the predictors is unlikely to have a large effect on the results.
Vars | R-Sq | R-Sq (adj) | PRESS | R-Sq (pred) | Mallows Cp | S | AICc | BIC | Cond No | I n s o l a t i o n | E a s t | S o u t h | N o r t h | T i m e o f D a y |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 72.1 | 71.0 | 4855.9 | 66.9 | 38.5 | 12.328 | 232.873 | 236.015 | 1.000 | X | ||||
1 | 39.4 | 37.1 | 10822.6 | 26.3 | 112.7 | 18.154 | 255.321 | 258.463 | 1.000 | X | ||||
2 | 85.9 | 84.8 | 2736.5 | 81.4 | 9.1 | 8.9321 | 215.798 | 219.600 | 1.807 | X | X | |||
2 | 82.0 | 80.6 | 3786.4 | 74.2 | 17.8 | 10.076 | 222.788 | 226.590 | 5.344 | X | X | |||
3 | 87.4 | 85.9 | 3089.7 | 79.0 | 7.6 | 8.5978 | 215.390 | 219.618 | 2.428 | X | X | X | ||
3 | 86.5 | 84.9 | 2725.9 | 81.4 | 9.7 | 8.9110 | 217.466 | 221.693 | 5.141 | X | X | X | ||
4 | 89.1 | 87.3 | 2847.2 | 80.6 | 5.8 | 8.1698 | 214.454 | 218.840 | 5.988 | X | X | X | X | |
4 | 88.0 | 86.0 | 3045.7 | 79.3 | 8.2 | 8.5550 | 217.127 | 221.512 | 20.427 | X | X | X | X | |
5 | 89.9 | 87.7 | 3109.9 | 78.8 | 6.0 | 8.0390 | 215.799 | 220.037 | 22.614 | X | X | X | X | X |