Learn and talk about Kolmogorov–Smirnov test, Non-parametric statistics, Normality tests, Statistical methods, Statistical tests

Illustration of the Kolmogorov-Smirnoff statistic. Red line is CDF, blue line is an ECDF, and the black arrow is the K-S statistic.

In statistics, the Kolmogorov–Smirnov test (K–S test) is a nonparametric test for the equality of continuous, one-dimensional probability distributions that can be used to compare a sample with a reference probability distribution (one-sample K–S test), or to compare two samples (two-sample K–S test). The Kolmogorov–Smirnov statistic quantifies a distance between the empirical distribution function of the sample and the cumulative distribution function of the reference distribution, or between the empirical distribution functions of two samples. The null distribution of this statistic is calculated under the null hypothesis that the samples are drawn from the same distribution (in the two-sample case) or that the sample is drawn from the reference distribution (in the one-sample case). In each case, the distributions considered under the null hypothesis are continuous distributions but are otherwise unrestricted.

The two-sample KS test is one of the most useful and general nonparametric methods for comparing two samples, as it is sensitive to differences in both location and shape of the empirical cumulative distribution functions of the two samples.

The Kolmogorov–Smirnov test can be modified to serve as a goodness of fit test. In the special case of testing for normality of the distribution, samples are standardized and compared with a standard normal distribution. This is equivalent to setting the mean and variance of the reference distribution equal to the sample estimates, and it is known that using these to define the specific reference distribution changes the null distribution of the test statistic: see below. Various studies have found that, even in this corrected form, the test is less powerful for testing normality than the Shapiro–Wilk test or Anderson–Darling test.^[1]

1 Kolmogorov–Smirnov statistic
2 Kolmogorov distribution
- 2.1 Test with estimated parameters
3 Two-sample Kolmogorov–Smirnov test
4 Setting confidence limits for the shape of a distribution function
5 The Kolmogorov–Smirnov statistic in more than one dimension
6 See also
7 Footnotes
8 References
9 External links

Kolmogorov–Smirnov statistic [edit]

The empirical distribution function F_n for n iid observations X_i is defined as

$F_n(x)={1 \over n}\sum_{i=1}^n I_{X_i\leq x}$

where $I_{X_i\leq x}$ is the indicator function, equal to 1 if X_i ≤ x and equal to 0 otherwise.

The Kolmogorov–Smirnov statistic for a given cumulative distribution function F(x) is

$D_n=\sup_x |F_n(x)-F(x)|$

where sup _x is the supremum of the set of distances. By the Glivenko–Cantelli theorem, if the sample comes from distribution F(x), then D_n converges to 0 almost surely. Kolmogorov strengthened this result, by effectively providing the rate of this convergence (see below). The Donsker's theorem provides yet a stronger result.

In practice, the statistic requires a relatively large number of data points to properly reject the null hypothesis.

Kolmogorov distribution [edit]

The Kolmogorov distribution is the distribution of the random variable

$K=\sup_{t\in[0,1]}|B(t)|$

where B(t) is the Brownian bridge. The cumulative distribution function of K is given by^[2]

$\operatorname{Pr}(K\leq x)=1-2\sum_{k=1}^\infty (-1)^{k-1} e^{-2k^2 x^2}=\frac{\sqrt{2\pi}}{x}\sum_{k=1}^\infty e^{-(2k-1)^2\pi^2/(8x^2)}.$

Both the form of the Kolmogorov–Smirnov test statistic and its asymptotic distribution under the null hypothesis were published by Andrey Kolmogorov,^[3] while a table of the distribution was published by Nikolai Vasilyevich Smirnov.^[4] Recurrence relations for the distribution of the test statistic in finite samples are available.^[3]

Under null hypothesis that the sample comes from the hypothesized distribution F(x),

$\sqrt{n}D_n\xrightarrow{n\to\infty}\sup_t |B(F(t))|$

in distribution, where B(t) is the Brownian bridge.

If F is continuous then under the null hypothesis $\sqrt{n}D_n$ converges to the Kolmogorov distribution, which does not depend on F. This result may also be known as the Kolmogorov theorem; see Kolmogorov's theorem for disambiguation.

The goodness-of-fit test or the Kolmogorov–Smirnov test is constructed by using the critical values of the Kolmogorov distribution. The null hypothesis is rejected at level $\alpha$ if

$\sqrt{n}D_n>K_\alpha,\,$

where K_α is found from

$\operatorname{Pr}(K\leq K_\alpha)=1-\alpha.\,$

The asymptotic power of this test is 1.

Test with estimated parameters [edit]

If either the form or the parameters of F(x) are determined from the data X_i the critical values determined in this way are invalid. In such cases, Monte Carlo or other methods may be required, but tables have been prepared for some cases. Details for the required modifications to the test statistic and for the critical values for the normal distribution and the exponential distribution have been published,^[5] and later publications also include the Gumbel distribution.^[6] The Lilliefors test represents a special case of this for the normal distribution. If kolmogorov test data is not fit for assumption (that is data is not normally distributed). So the logarithm may help to overcome the assumption.

Two-sample Kolmogorov–Smirnov test [edit]

Illustration of the Two-Sample Kolmogorov-Smirnoff statistic. Red and blue lines each correspond to an empirical distribution function, and the black arrow is the two-sample K-S statistic.

The Kolmogorov–Smirnov test may also be used to test whether two underlying one-dimensional probability distributions differ. In this case, the Kolmogorov–Smirnov statistic is

$D_{n,n'}=\sup_x |F_{1,n}(x)-F_{2,n'}(x)|,$

where $F_{1,n}$ and $F_{2,n'}$ are the empirical distribution functions of the first and the second sample respectively.

The null hypothesis is rejected at level $\alpha$ if

$\sqrt{\frac{n + n'}{n n'}}D_{n,n'}>K_\alpha.$ ^[7]

The value of $D_{n,n'}$ is given in the table below for each level of $\alpha$ ^[7]

$\alpha$	0.10	0.05	0.025	0.01	0.005	0.001
$D_{n,n'}$	1.22	1.36	1.48	1.63	1.73	1.95

Note that the two-sample test checks whether the two data samples come from the same distribution. This does not specify what that common distribution is (e.g. normal or not normal). Again, tables of critical values have been published.^[5]^[7]

Setting confidence limits for the shape of a distribution function [edit]

While the Kolmogorov–Smirnov test is usually used to test whether a given F(x) is the underlying probability distribution of F_n(x), the procedure may be inverted to give confidence limits on F(x) itself. If one chooses a critical value of the test statistic D_α such that P(D_n > D_α) = α, then a band of width ±D_α around F_n(x) will entirely contain F(x) with probability 1 − α.

The Kolmogorov–Smirnov statistic in more than one dimension [edit]

A distribution free multivariate Kolmogorov-Smirnov goodness of fit test has been proposed by Justel, Peña and Zamar (1997). The test uses a statistic which is built using Rosenblatt's transformation and an algorithm is developed to compute it in the bivariate case. An approximate test, that can be easily computed in any dimension, is also presented.

The Kolmogorov–Smirnov test statistic needs to be modified if a similar test is to be applied to multivariate data. This is not straightforward because the maximum difference between two joint cumulative distribution functions is not generally the same as the maximum difference of any of the complementary distribution functions. Thus the maximum difference will differ depending on which of $\Pr(x < X \and y < Y)$ or $\Pr(X < x \and Y > y)$ or any of the other two possible arrangements is used. One might require that the result of the test used should not depend on which choice is made.

One approach to generalizing the Kolmogorov–Smirnov statistic to higher dimensions which meets the above concern is to compare the cdfs of the two samples with all possible orderings, and take the largest of the set of resulting K-S statistics. In d dimensions, there are 2^d−1 such orderings. One such variation is due to Peacock^[8] and another to Fasano & Franceschini^[9] (see Lopes et al. for a comparison and computational details).^[10] Critical values for the test statistic can be obtained by simulations, but depend on the dependence structure in the joint distribution.

Footnotes [edit]

^ Stephens, M. A. (1974). "EDF Statistics for Goodness of Fit and Some Comparisons". Journal of the American Statistical Association (American Statistical Association) 69 (347): 730–737. doi:10.2307/2286009. JSTOR 2286009.
^ Marsaglia G, Tsang WW, Wang J (2003). "Evaluating Kolmogorov’s Distribution". Journal of Statistical Software 8 (18): 1–4.
^ ^a ^b Kolmogorov A (1933). "Sulla determinazione empirica di una legge di distribuzione". G. Ist. Ital. Attuari 4: 83.
^ Smirnov NV (1948). "Tables for estimating the goodness of fit of empirical distributions". Annals of Mathematical Statistics 19: 279.
^ ^a ^b Pearson E.S. and Hartley, H.O., ed. (1972). Biometrika Tables for Statisticians 2. Cambridge University Press. pp. 117–123, Tables 54, 55. ISBN 0-521-06937-8.
^ Galen R. Shorack and Jon A. Wellner (1986). Empirical Processes with Applications to Statistics. Wiley. p. 239. ISBN 047186725X.
^ ^a ^b ^c Table of critical values for the two-sample test
^ Peacock J.A. (1983). "Two-dimensional goodness-of-fit testing in astronomy". Monthly Notices of the Royal Astronomical Society 202: 615–627. Bibcode:1983MNRAS.202..615P.
^ Fasano, G., Franceschini, A. (1987). "A multidimensional version of the Kolmogorov–Smirnov test". Monthly Notices of the Royal Astronomical Society (ISSN 0035-8711) 225: 155–170.
^ Lopes, R.H.C., Reid, I., Hobson, P.R. (April 23–27, 2007). "The two-dimensional Kolmogorov-Smirnov test". XI International Workshop on Advanced Computing and Analysis Techniques in Physics Research. Amsterdam, the Netherlands.

References [edit]

Justel, A., Peña, D. and Zamar, R. (1997) A multivariate Kolmogorov-Smirnov test of goodness of fit, Statistics & Probability Letters, 35(3), 251-259.
Eadie, W.T.; D. Drijard, F.E. James, M. Roos and B. Sadoulet (1971). Statistical Methods in Experimental Physics. Amsterdam: North-Holland. pp. 269–271. ISBN 0-444-10117-9.
Stuart, Alan; Ord, Keith; Arnold, Steven [F.] (1999). Classical Inference and the Linear Model. Kendall's Advanced Theory of Statistics 2A (Sixth ed.). London: Arnold. pp. 25.37–25.43. ISBN 0-340-66230-1. MR 1687411. Unknown parameter |publisher2= ignored (help); Unknown parameter |location2= ignored (help)
Corder, G.W., Foreman, D.I. (2009).Nonparametric Statistics for Non-Statisticians: A Step-by-Step Approach Wiley, ISBN 978-0-470-45461-9
Stephens, M.A. (1979) Test of fit for the logistic distribution based on the empirical distribution function, Biometrika, 66(3), 591-5.

External links [edit]

Hazewinkel, Michiel, ed. (2001), "Kolmogorov-Smirnov test", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
Short introduction
KS test explanation
JavaScript implementation of one- and two-sided tests
Online calculator with the K-S test
Open-source C++ code to compute the Kolmogorov distribution and perform the K-S test
Paper on Evaluating Kolmogorov’s Distribution; contains C implementation. This is the method used in Matlab.

Statistics

Descriptive statistics

Continuous data

Location	Mean (Arithmetic, Geometric, Harmonic) Median Mode

Dispersion	Range Standard deviation Coefficient of variation Percentile Interquartile range

Shape	Variance Skewness Kurtosis Moments L-moments

Count data

Index of dispersion

Summary tables

Dependence

Statistical graphics

Data collection

Designing studies	Effect size Standard error Statistical power Sample size determination

Survey methodology	Sampling Stratified sampling Opinion poll Questionnaire

Controlled experiment	Design of experiments Randomized experiment Random assignment Replication Blocking Factorial experiment Optimal design

Uncontrolled studies	Natural experiment Quasi-experiment Observational study

Statistical inference

Statistical theory	Sampling distribution Order statistics Sufficiency Completeness Exponential family Permutation test (Randomization test) Empirical distribution Bootstrap U statistic Efficiency Asymptotics Robustness

Frequentist inference	Unbiased estimator (Mean unbiased minimum variance, Median unbiased) Biased estimators (Maximum likelihood, Method of moments, Minimum distance, Density estimation) Confidence interval Testing hypotheses Power Parametric tests (Likelihood-ratio, Wald, Score)

Specific tests	Z (normal) Student's t-test F Chi-squared Signed-rank (1-sample, 2-sample, 1-way anova) Shapiro–Wilk Kolmogorov–Smirnov

Bayesian inference	Bayesian probability Prior Posterior Credible interval Bayes factor Bayesian estimator Maximum posterior estimator

Correlation and regression analysis

Correlation	Pearson product–moment correlation Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust

Generalized linear model	Exponential families Logistic (Bernoulli) Binomial Poisson

Partition of variance	Analysis of variance (ANOVA) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical, multivariate, time-series, or survival analysis

Categorical data

Multivariate statistics

Time series analysis

General	Decomposition Trend Stationarity Seasonal adjustment

Time domain	ACF PACF XCF ARMA model ARIMA model Vector autoregression

Frequency domain	Spectral density estimation

Survival analysis

Applications

Biostatistics	Bioinformatics Clinical trials & studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process & Quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Original courtesy of Wikipedia: http://en.wikipedia.org/wiki/Kolmogorov–Smirnov_test — Please support Wikipedia.
A portion of the proceeds from advertising on Digplanet goes to supporting Wikipedia.

Kolmogorov–Smirnov test

Contents