Rational ArithmeticMathematicaFunctions to Evaluate the Two-Sided One Sample K-S Cumulative Sampling Distribution

Abstract

One of the most widely used goodness-of-fit tests is the two-sided one sample KolmogorovSmirnov (K-S) test which has been implemented by many computer statistical software packages. To calculate a two-sided p value (evaluate the cumulative sampling distribution), these packages use various methods including recursion formulae, limiting distributions, and approximations of unknown accuracy developed over thirty years ago. Based on an extensive literature search for the two-sided one sample K-S test, this paper identifies an exact formula for sample sizes up to 31, six recursion formulae, and one matrix formula that can be used to calculate a p value. To ensure accurate calculation by avoiding catastrophic cancelation and eliminating rounding error, each of these formulae is implemented in rational arithmetic. For the six recursion formulae and the matrix formula, computational experience for sample sizes up to 500 shows that computational times are increasing functions of both the sample size and the number of digits in the numerator and denominator integers of the rational number test statistic. The computational times of the seven formulae vary immensely but the Durbin recursion formula is almost always the fastest. Linear search is used to calculate the inverse of the cumulative sampling distribution (find the confidence interval half-width) and tables of calculated half-widths are presented for sample sizes up to 500. Using calculated half-widths as input, computational times for the fastest formula, the Durbin recursion formula, are given for sample sizes up to two thousand.

A Novel Method for Fast Change-Point Detection on Simulated Time Series and Electrocardiogram Data

“…KS test is one of the most useful and general non-parametric methods for comparing two samples, because it is sensitive to those differences in both location and shape of the empirical cumulative distribution functions (e.c.d.f) of two samples [31] , [32] . Suppose is a time series, we observe, where are discrete and centred i.i.d.…”

Computing the Cumulative Distribution Function and Quantiles of the limit of the Two-sided Kolmogorov-Smirnov Statistic

“…Drew Glen & Leemis [14] generated the collection of polynomial splines for n <= 30. Brown and Harvey [15,16,17] implemented several algorithms in both rational arithmetic and arbitrary precision arithmetic. Simard and L'Ecuyer [18] analyzed all the known algorithms for numerical stability and sped.…”