Take the 2-minute tour ×
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It's 100% free, no registration required.

I want to use the Kolmogorov-Smirnov test to test if a sample is drawn from a Pareto distribution. Unfortunately, the only way to estimate the distribution's parameters is from the sample.

Does anybody know of the existence of tables with modified critical values or any other way to perform the test in this case?

share|improve this question

1 Answer 1

up vote 3 down vote accepted

A Kolmogorov-Smirnov test with estimated parameters is known as Liliiefors' test. The values of the test statistic tend to be smaller than with the KS test.

It is no longer nonparametric; you need to work out the distribution of the test statistic for each distribution type (and it differs again if you estimate a subset of the parameters rather than all of them).

Some packages offer Lilliefors test for the exponential distribution and the normal distribution (which are the cases Lilliefors discusses in his two papers).

Here's a way to use the exponential for the Pareto (I'll use the notation at that link):

Case 1: known $x_m$, unknown $\alpha$:

  • If $x_m$ isn't $1$, divide through by $x_m$, so you have a Pareto with $x_m=1$.

  • take logs, resulting in an exponential with (unknown) scale parameter $\alpha$.

  • apply Lilliefors test for an exponential r.v.

Case 2: unknown $x_m$, unknown $\alpha$. While this case can be done by simulation, we can make use of an exponential Lilliefors test as follows:

  • take logs, yielding a shifted exponential

  • subtract the minimum observation from all observations, and discard the minimum, leaving $n-1$ observations

  • test the reduced sample for exponentiality via the Lilliefors test as before.

share|improve this answer
    
Thanks for your answer! Do you have a source for Case 2? –  Tom Aug 12 at 7:39
    
In the meantime, I came across Coefficients of the asymptotic distribution of the kolmogorov-smirnov statistic when parameters are estimated. They propose the following transformation: x_i' = (n - i + 1) * ln[x_i / x_(i-1)] Would this be suitable as well? –  Tom Aug 12 at 7:47
    
Hi Tom, I can't get the article at the moment, but I believe that the proposed transformation would be for consecutive order statistics (the 'gaps') (in which case your expression is missing a set of parentheses), and also yields an exponential. It should work just fine. Indeed, it may be slightly better (in that it would possibly give better power against more interesting alternatives) –  Glen_b Aug 14 at 2:33

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.