Regression validation

In statistics, regression validation is the process of deciding whether the numerical results quantifying hypothesized relationships between variables, obtained from regression analysis, are acceptable as descriptions of the data. The validation process can involve analyzing the goodness of fit of the regression, analyzing whether the regression residuals are random, and checking whether the model's predictive performance deteriorates substantially when applied to data that were not used in model estimation.

Goodness of fit

Main article: Goodness of fit

One measure of goodness of fit is the coefficient of determination, often denoted, R². In ordinary least squares with an intercept, it ranges between 0 and 1. However, an R² close to 1 does not guarantee that the model fits the data well. For example, if the functional form of the model does not match the data, R² can be high despite a poor model fit. Anscombe's quartet consists of four example data sets with similarly high R² values, but data that sometimes clearly does not fit the regression line. Instead, the data sets include outliers, high-leverage points, or non-linearities.

One problem with the R² as a measure of model validity is that it can always be increased by adding more variables into the model, except in the unlikely event that the additional variables are exactly uncorrelated with the dependent variable in the data sample being used. This problem can be avoided by doing an F-test of the statistical significance of the increase in the R², or by instead using the adjusted R2.

Analysis of residuals

Main article: residual analysis

The residuals from a fitted model are the differences between the responses observed at each combination of values of the explanatory variables and the corresponding prediction of the response computed using the regression function. Mathematically, the definition of the residual for the i^th observation in the data set is written

${\displaystyle e_{i}=y_{i}-f(x_{i};{\hat {\beta }}),}$

An illustrative plot of a fit to data (green curve in top panel, data in red) plus a plot of residuals: red points in bottom plot. Dashed curve in bottom panel is a straight line fit to the residuals. If the functional form is correct then there should be little or no trend to the residuals - as seen here.

with y_i denoting the i^th response in the data set and x_i the vector of explanatory variables, each set at the corresponding values found in the i^th observation in the data set.

If the model fit to the data were correct, the residuals would approximate the random errors that make the relationship between the explanatory variables and the response variable a statistical relationship. Therefore, if the residuals appear to behave randomly, it suggests that the model fits the data well. On the other hand, if non-random structure is evident in the residuals, it is a clear sign that the model fits the data poorly. The next section details the types of plots to use to test different aspects of a model and gives the correct interpretations of different results that could be observed for each type of plot.

Graphical analysis of residuals

See also: Statistical graphics

A basic, though not quantitatively precise, way to check for problems that render a model inadequate is to conduct a visual examination of the residuals (the mispredictions of the data used in quantifying the model) to look for obvious deviations from randomness. If a visual examination suggests, for example, the possible presence of heteroscedasticity (a relationship between the variance of the model errors and the size of an independent variable's observations), then statistical tests can be performed to confirm or reject this hunch; if it is confirmed, different modeling procedures are called for.

Different types of plots of the residuals from a fitted model provide information on the adequacy of different aspects of the model.

1. sufficiency of the functional part of the model: scatter plots of residuals versus predictors
2. non-constant variation across the data: scatter plots of residuals versus predictors; for data collected over time, also plots of residuals against time
3. drift in the errors (data collected over time): run charts of the response and errors versus time
4. independence of errors: lag plot
5. normality of errors: histogram and normal probability plot

Graphical methods have an advantage over numerical methods for model validation because they readily illustrate a broad range of complex aspects of the relationship between the model and the data.

Quantitative analysis of residuals

Main article: Regression diagnostic

Numerical methods also play an important role in model validation. For example, the lack-of-fit test for assessing the correctness of the functional part of the model can aid in interpreting a borderline residual plot. One common situation when numerical validation methods take precedence over graphical methods is when the number of parameters being estimated is relatively close to the size of the data set. In this situation residual plots are often difficult to interpret due to constraints on the residuals imposed by the estimation of the unknown parameters. One area in which this typically happens is in optimization applications using designed experiments. Logistic regression with binary data is another area in which graphical residual analysis can be difficult.

Serial correlation of the residuals can indicate model misspecification, and can be checked for with the Durbin–Watson statistic. The problem of heteroskedasticity can be checked for in any of several ways.

Out-of-sample evaluation

Main article: Cross-validation

Cross-validation is the process of assessing how the results of a statistical analysis will generalize to an independent data set. If the model has been estimated over some, but not all, of the available data, then the model using the estimated parameters can be used to predict the held-back data. If, for example, the out-of-sample mean squared error, also known as the mean squared prediction error, is substantially higher than the in-sample mean square error, this is a sign of deficiency in the model.

A development in medical statistics is the use of out-of-sample cross validation techniques in meta-analysis. It forms the basis of the validation statistic, Vn, which is used to test the statistical validity of meta-analysis summary estimates. Essentially it measures a type of normalized prediction error and its distribution is a linear combination of χ² variables of degree 1. ^[1]

See also

• All models are wrong
• Model selection
• Prediction error
• Prediction interval
• Resampling (statistics)
• Statistical conclusion validity
• Statistical model specification
• Statistical model validation
• Validity (statistics)
• Coefficient of determination
• Lack-of-fit sum of squares
• Reduced chi-squared

References

1. Willis BH, Riley RD (2017). "Measuring the statistical validity of summary meta-analysis and meta-regression results for use in clinical practice". Statistics in Medicine. 36 (21): 3283–3301. doi:10.1002/sim.7372. PMC   5575530 . PMID   28620945.

Further reading

• Arboretti Giancristofaro, R.; Salmaso, L. (2003), "Model performance analysis and model validation in logistic regression", Statistica , 63: 375–396
• Kmenta, Jan (1986), Elements of Econometrics (Second ed.), Macmillan, pp. 593–600; republished in 1997 by University of Michigan Press

External links

• How can I tell if a model fits my data? (NIST)
• NIST/SEMATECH e-Handbook of Statistical Methods
• Model Diagnostics (Eberly College of Science)

 This article incorporates public domain material from the National Institute of Standards and Technology