PROC NPAR1WAY: Tests Based on the Empirical Distribution Function :: SAS/STAT(R) 9.2 User's Guide, Second Edition

The NPAR1WAY Procedure

Tests Based on the Empirical Distribution Function

If you specify the EDF option, PROC NPAR1WAY computes tests based on the empirical distribution function. These include the Kolmogorov-Smirnov and Cramer-von Mises tests, and also the Kuiper test for two-sample data. This section gives formulas for these test statistics. For further information about the formulas and the interpretation of EDF statistics, see Hollander and Wolfe (1999) and Gibbons and Chakraborti (1992). For details about the $\text{[math]}$ -sample analogs of the Kolmogorov-Smirnov and Cramer-von Mises statistics used by NPAR1WAY, see Kiefer (1959).

The empirical distribution function (EDF) of a sample $\text{[math]}$ , $\text{[math]}$ , is defined as

$\text{[math]}$

where $\text{[math]}$ is an indicator function. PROC NPAR1WAY uses the subsample of values within the $\text{[math]}$ th class level to generate an EDF for the class, $\text{[math]}$ . The EDF for the overall sample, pooled over classes, can also be expressed as

$\text{[math]}$

where $\text{[math]}$ is the number of observations in the $\text{[math]}$ th class level, and $\text{[math]}$ is the total number of observations.

Kolmogorov-Smirnov Test

The Kolmogorov-Smirnov statistic measures the maximum deviation of the EDF within the classes from the pooled EDF. PROC NPAR1WAY computes the Kolmogorov-Smirnov statistic as

$\text{[math]}$

The asymptotic Kolmogorov-Smirnov statistic is computed as

$\text{[math]}$

For each class level $\text{[math]}$ and overall, PROC NPAR1WAY displays the value of $\text{[math]}$ at the maximum deviation from $\text{[math]}$ and the value $\text{[math]}$ at the maximum deviation from $\text{[math]}$ . PROC NPAR1WAY also gives the observation where the maximum deviation occurs.

If there are only two class levels, PROC NPAR1WAY computes the two-sample Kolmogorov-Smirnov test statistic $\text{[math]}$ as

$\text{[math]}$

The p-value for this test is the probability that $\text{[math]}$ is greater than the observed value $\text{[math]}$ under the null hypothesis of no difference between class levels or samples. PROC NPAR1WAY computes the asymptotic p-value for $\text{[math]}$ with the approximation

$\text{[math]}$

where

$\text{[math]}$

The quality of this approximation has been studied by Hodges (1957).

If you specify the D option, or if you request exact Kolmogorov-Smirnov p-values with the KS option in the EXACT statement, PROC NPAR1WAY also computes the one-sided Kolmogorov-Smirnov statistics $\text{[math]}$ and $\text{[math]}$ for two-sample data as

$\text{[math]}$

The asymptotic probability that $\text{[math]}$ is greater than the observed value $\text{[math]}$ , under the null hypothesis of no difference between the two class levels, is computed as

$\text{[math]}$

Similarly, the asymptotic probability that $\text{[math]}$ is greater than the observed value $\text{[math]}$ is computed as

$\text{[math]}$

To request exact p-values for the Kolmogorov-Smirnov statistics, you can specify the KS option in the EXACT statement. See the section Exact Tests for more information.

Cramer-von Mises Test

The Cramer-von Mises statistic is defined as

$\text{[math]}$

where $\text{[math]}$ is the number of ties at the $\text{[math]}$ th distinct value and $\text{[math]}$ is the number of distinct values. The asymptotic value is computed as

$\text{[math]}$

PROC NPAR1WAY displays the contribution of each class level to the sum $\text{[math]}$ .

Kuiper Test

For data with two class levels, PROC NPAR1WAY computes the Kuiper statistic, its scaled value for the asymptotic distribution, and the asymptotic p-value. The Kuiper statistic is computed as

$\text{[math]}$

The asymptotic value is

$\text{[math]}$

PROC NPAR1WAY displays the value of $\text{[math]}$ for each class level.

The p-value for the Kuiper test is the probability of observing a larger value of $\text{[math]}$ under the null hypothesis of no difference between the two classes. PROC NPAR1WAY computes this p-value according to Owen (1962, p. 441).

Top of Page