Goodness of fit

From Wikipedia, the free encyclopedia

The goodness of fit of a statistical model describes how well it fits a set of observations. Measures of goodness of fit typically summarize the discrepancy between observed values and the values expected under the model in question. Such measures can be used in statistical hypothesis testing, e.g. to test for normality of residuals, to test whether two samples are drawn from identical distributions (see Kolmogorov-Smirnov test), or whether outcome frequencies follow a specified distribution (see Pearson's chi-square test). In the analysis of variance, one of the components into which the variance is partitioned may be a lack-of-fit sum of squares.

Contents 1 Fit of distributions 2 Regression analysis 2.1 Example 3 Categorical data 3.1 Example 1 3.2 Binomial case 4 Other measures of fit 5 See also 6 References

[edit] Fit of distributions

In assessing whether a given distribution is suited to a data-set, the following tests and their underlying measures of fit can be used:

• Kolmogorov–Smirnov test;
• Cramér–von-Mises criterion;
• Anderson-Darling test.

[edit] Regression analysis

In regression analysis, the following topics relate to goodness of fit:

• Coefficient of determination (The R squared measure of goodness of fit);
• lack-of-fit sum of squares.

[edit] Example

One way in which a measure of goodness of fit statistic can be constructed, in the case where the variance of the measurement error is known, is to construct a weighted sum of squared errors:

where σ² is the known variance of the observation.^[1] This definition is only useful when one has estimates for the error on the measurements, but it leads to a situation where a chi-square distribution can be used to test goodness of fit, provided that the errors can be assumed to have a normal distribution.

The reduced chi-squared statistic is simply the chi-squared divided by the number of degrees of freedom: ^{[1] [2] [3] [4]}

where ν is the number of degrees of freedom, usually given by N − n − 1, where N is the number of observations, and n is the number of fitted parameters, assuming that the mean value is an additional fitted parameter. The advantage of the reduced chi-squared is that it already normalizes for the number of data points and model complexity.

As a rule of thumb, a large indicates a poor model fit. However indicates that the model is 'over-fitting' the data (either the model is improperly fitting noise, or the error variance has been over-estimated). A indicates that the fit has not fully captured the data (or that the error variance has been under-estimated). In principle a value of indicates that the extent of the match between observations and estimates is in accord with the error variance.

[edit] Categorical data

The following are examples that arise in the context of categorical data.

[edit] Example 1

Pearson's chi-square test uses a measure of goodness of fit which is the sum of differences between observed and expected outcome frequencies (that is, counts of observations), each squared and divided by the expectation:

where:

O_i = an observed frequency (ie count) for the i^th bin

E_i = an expected (theoretical) frequency for the i^th bin, asserted by the null hypothesis.

The resulting value can be compared to the chi-square distribution to determine the goodness of fit. In order to determine the degrees of freedom of the chi-squared distribution, one takes the total number of observed frequencies and subtracts one. For example, if there are eight different frequencies, one would compare to a chi-squared with seven degrees of freedom.

[edit] Binomial case

A binomial experiment is a sequence of independent trials in which the trials can result in one of two outcomes, success or failure. There are n trials each with probability of success, denoted by p. Provided that np_i ≫ 1 for every i (where i = 1, 2, ..., k), then

This has approximately a chi-squared distribution with k − 1 df. The fact that df = k − 1 is a consequence of the restriction . We know there are k observed cell counts, however, once any k − 1 are known, the remaining one is uniquely determined. Basically, one can say, there are only k − 1 freely determined cell counts, thus df = k − 1.

[edit] Other measures of fit

The likelihood ratio test statistic is a measure of the goodness of fit of a model, judged by whether an expanded form of the model provides a substantially improved fit.

[edit] See also

• Overfitting

[edit] References

1. ^ ^{a b} Charlie Laub and Tonya L. Kuhl: Chi-Square Data Fitting. University California, Davis.
2. ^ John Robert Taylor: An introduction to error analysis, page 268. University Science Books, 1997.
3. ^ Kirkman, T.W.: Chi-Square Curve Fitting.
4. ^ David M. Glover, William J. Jenkins, and Scott C. Doney: Least Squares and regression techniques, goodness of fit and tests, non-linear least squares techniques. Woods Hole Oceanographic Institute, 2008.