This article includes a list ofgeneral references, butit lacks sufficient correspondinginline citations. Please help toimprove this article byintroducing more precise citations.(September 2016) (Learn how and when to remove this message) |
| Part of a series on |
| Regression analysis |
|---|
| Models |
| Estimation |
| Background |
Instatistics andoptimization,errors andresiduals are two closely related and easily confused measures of thedeviation of anobserved value of anelement of astatistical sample from its "true value" (not necessarily observable). Theerror of anobservation is the deviation of the observed value from the true value of a quantity of interest (for example, apopulation mean). Theresidual is the difference between the observed value and theestimated value of the quantity of interest (for example, asample mean). The distinction is most important inregression analysis, where the concepts are sometimes called theregression errors andregression residuals and where they lead to the concept ofstudentized residuals.Ineconometrics, "errors" are also calleddisturbances.[1][2][3]
Suppose there is a series of observations from aunivariate distribution and we want to estimate themean of that distribution (the so-calledlocation model). In this case, the errors are the deviations of the observations from the population mean, while the residuals are the deviations of the observations from the sample mean.
Astatistical error (ordisturbance) is the amount by which an observation differs from itsexpected value, the latter being based on the wholepopulation from which the statistical unit was chosen randomly. For example, if the mean height in a population of 21-year-old men is 1.75 meters, and one randomly chosen man is 1.80 meters tall, then the "error" is 0.05 meters; if the randomly chosen man is 1.70 meters tall, then the "error" is −0.05 meters. The expected value, being themean of the entire population, is typically unobservable, and hence the statistical error cannot be observed either.
Aresidual (or fitting deviation), on the other hand, is an observableestimate of the unobservable statistical error. Consider the previous example with men's heights and suppose we have a random sample ofn people. Thesample mean could serve as a good estimator of thepopulation mean. Then we have:
Note that, because of the definition of the sample mean, the sum of the residuals within a random sample is necessarily zero, and thus the residuals are necessarilynotindependent. The statistical errors, on the other hand, are independent, and their sum within the random sample isalmost surely not zero.
One can standardize statistical errors (especially of anormal distribution) in az-score (or "standard score"), and standardize residuals in at-statistic, or more generallystudentized residuals.
If we assume a normally distributed population with mean μ andstandard deviation σ, and choose individuals independently, then we have
and thesample mean
is a random variable distributed such that:
Thestatistical errors are then
withexpected values of zero,[4] whereas theresiduals are
The sum of squares of thestatistical errors, divided byσ2, has achi-squared distribution withndegrees of freedom:
However, this quantity is not observable as the population mean is unknown. The sum of squares of theresiduals, on the other hand, is observable. The quotient of that sum by σ2 has a chi-squared distribution with onlyn − 1 degrees of freedom:
This difference betweenn andn − 1 degrees of freedom results inBessel's correction for the estimation ofsample variance of a population with unknown mean and unknown variance. No correction is necessary if the population mean is known.
It is remarkable that thesum of squares of the residuals and the sample mean can be shown to be independent of each other, using, e.g.Basu's theorem. That fact, and the normal and chi-squared distributions given above form the basis of calculations involving the t-statistic:
where represents the errors, represents the sample standard deviation for a sample of sizen, and unknownσ, and the denominator term accounts for the standard deviation of the errors according to:[5]
Theprobability distributions of the numerator and the denominator separately depend on the value of the unobservable population standard deviationσ, butσ appears in both the numerator and the denominator and cancels. That is fortunate because it means that even though we do not know σ, we know the probability distribution of this quotient: it has aStudent's t-distribution withn − 1 degrees of freedom. We can therefore use this quotient to find aconfidence interval for μ. This t-statistic can be interpreted as "the number of standard errors away from the regression line."[6]
Inregression analysis, the distinction betweenerrors andresiduals is subtle and important, and leads to the concept ofstudentized residuals. Given an unobservable function that relates the independent variable to the dependent variable – say, a line – the deviations of the dependent variable observations from this function are the unobservable errors. If one runs a regression on some data, then the deviations of the dependent variable observations from thefitted function are the residuals. If the linear model is applicable, a scatterplot of residuals plotted against the independent variable should be random about zero with no trend to the residuals.[5] If the data exhibit a trend, the regression model is likely incorrect; for example, the true function may be a quadratic or higher order polynomial. If they are random, or have no trend, but "fan out" - they exhibit a phenomenon calledheteroscedasticity. If all of the residuals are equal, or do not fan out, they exhibithomoscedasticity.
However, a terminological difference arises in the expressionmean squared error (MSE). The mean squared error of a regression is a number computed from the sum of squares of the computedresiduals, and not of the unobservableerrors. If that sum of squares is divided byn, the number of observations, the result is the mean of the squared residuals. Since this is abiased estimate of the variance of the unobserved errors, the bias is removed by dividing the sum of the squared residuals bydf =n − p − 1, instead ofn, wheredf is the number ofdegrees of freedom (n minus the number of parameters (excluding the intercept) p being estimated - 1). This forms an unbiased estimate of the variance of the unobserved errors, and is called the mean squared error.[7]
Another method to calculate the mean square of error when analyzing the variance of linear regression using a technique like that used inANOVA (they are the same because ANOVA is a type of regression), the sum of squares of the residuals (aka sum of squares of the error) is divided by the degrees of freedom (where the degrees of freedom equaln − p − 1, wherep is the number of parameters estimated in the model (one for each variable in the regression equation, not including the intercept)). One can then also calculate the mean square of the model by dividing the sum of squares of the model minus the degrees of freedom, which is just the number of parameters. Then the F value can be calculated by dividing the mean square of the model by the mean square of the error, and we can then determine significance (which is why you want the mean squares to begin with.).[8]
However, because of the behavior of the process of regression, thedistributions of residuals at different data points (of the input variable) may varyeven if the errors themselves are identically distributed. Concretely, in alinear regression where the errors are identically distributed, the variability of residuals of inputs in the middle of the domain will behigher than the variability of residuals at the ends of the domain:[9] linear regressions fit endpoints better than the middle. This is also reflected in theinfluence functions of various data points on theregression coefficients: endpoints have more influence.
Thus to compare residuals at different inputs, one needs to adjust the residuals by the expected variability ofresiduals, which is calledstudentizing. This is particularly important in the case of detectingoutliers, where the case in question is somehow different from the others in a dataset. For example, a large residual may be expected in the middle of the domain, but considered an outlier at the end of the domain.
The use of the term "error" as discussed in the sections above is in the sense of a deviation of a value from a hypothetical unobserved value. At least two other uses also occur in statistics, both referring to observableprediction errors:
Themean squared error (MSE) refers to the amount by which the values predicted by an estimator differ from the quantities being estimated (typically outside the sample from which the model was estimated).Theroot mean square error (RMSE) is the square root of MSE.Thesum of squares of errors (SSE) is the MSE multiplied by the sample size.
Sum of squares of residuals (SSR) is the sum of the squares of the deviations of the actual values from the predicted values, within the sample used for estimation. This is the basis for theleast squares estimate, where the regression coefficients are chosen such that the SSR is minimal (i.e. its derivative is zero).
Likewise, thesum of absolute errors (SAE) is the sum of the absolute values of the residuals, which is minimized in theleast absolute deviations approach to regression.
Themean error (ME) is the bias.Themean residual (MR) is always zero for least-squares estimators.