User:MarkSweep/work



The Pearson distribution is a family of continuous probability distributions. It was first published by Karl Pearson in 1895 and subsequently extended by him in 1901 and 1916 in a series of articles on biostatistics.

The Pearson system was originally devised in an effort to model visibly skewed observations. It was well known at the time how to adjust a theoretical model to fit the first two cumulants or moments of observed data: Any probability distribution can be extended straightforwardly to form a location-scale family. Except in pathological cases, a location-scale family can be made to fit the observed mean (first cumulant) and variance (second cumulant) arbitrarily well. However, it was not known how to construct probability distributions in which the skewness (standardized third cumulant) and kurtosis (standardized fourth cumulant) could be adjusted equally freely. This need became apparent when trying to fit known theoretical models to observed data that exhibited skewness. Pearson's examples include survival data, which are usually asymmetric.

In his original paper, Pearson (1895, p. 360) identified four types of distributions (numbered I through IV) in addition to the normal distribution (which was originally known as type V). The classification depended on whether the distributions were supported on a bounded interval, on a half-line, or on the whole real line; and whether they were potentially skewed or necessarily symmetric. A second paper (Pearson 1901) fixed two omissions: it redefined the type V distribution (originally just the normal distribution, but now the inverse-gamma distribution) and introduced the type VI distribution. Together the first two papers cover the five main types of the Pearson system (I, III, VI, V, and IV). In a third paper, Pearson (1916) introduced further special cases and subtypes (VII through XII).

Rhind (1909, pp. 430–432) devised a simple way of visualizing the parameter space of the Pearson system, which was subsequently adopted by Pearson (1916, plate 1 and pp. 430ff., 448ff.). The Pearson types are characterized by two quantities, commonly referred to as $$\beta_1$$ and $$\beta_2$$. The first is the square of the skewness: $$\beta_1 = \gamma_1^2$$ where $$\gamma_1$$ is the skewness, or third standardized moment. The second is the traditional kurtosis, or fourth standardized moment: $$\beta_2 = \gamma_2 + 3$$. (Modern treatments define kurtosis $$\gamma_2$$ in terms of cumulants instead of moments, so that for a normal distribution we have $$\gamma_2=0$$ and $$\beta_2=3$$. Here we follow the historical precedent and use $$\beta_2$$.) The diagram on the right shows which Pearson type a given concrete distribution (identified by a point $$(\beta_1, \beta_2)$$) belongs to.

XXXX

Many of the skewed and/or non-mesokurtic distributions familiar to us today were still unknown in the early 1890s. What is now known as the beta distribution had been used by Thomas Bayes as a posterior distribution of the parameter of a Bernoulli distribution in his 1763 work on inverse probability. The Beta distribution gained prominence due to its membership in Pearson's system and was known until the 1940s as the Pearson type I distribution. (Pearson's type II distribution is a special case of type I, but is usually no longer singled out.) The gamma distribution originated from Pearson's work (Pearson 1893, p. 331; Pearson 1895, pp. 357, 360, 373–376) and was known as the Pearson type III distribution, before acquiring its modern name in the 1930s and 1940s. Pearson's 1895 paper introduced the type IV distribution, which contains Student's t-distribution as a special case, predating William Gosset's subsequent use by several years. His 1901 paper introduced the inverse-gamma distribution (type V) and the beta prime distribution (type VI).

Definition
A Pearson density p is defined to be any valid solution to the differential equation (cf. Pearson 1895, p. 381)


 * $$p'(x) + \frac{x-a_0}{b_2 x^2 + b_1 x + b_0}\;p(x) = 0.

\qquad (1) \!$$

The Pearson family has four parameters (a0, b0, b1, b2), which can be used to freely adjust the first four moments of the distribution, subject to very few constraints. The parameter a0 determines a stationary point, and hence under some conditions a mode of the distribution, since


 * $$p'(a_0) = 0 \!$$

follows directly from the differential equation.

Since we are confronted with a linear differential equation with variable coefficients, its solution is straightforward:


 * $$p(x) \propto \exp\left( -\!\int\!\!\frac{x-a_0}{b_2 x^2 + b_1 x + b_0} \,\mathrm{d}x \right).

\!$$

The integral in this solution simplifies considerably when certain special cases of the integrand are considered. Pearson (1895, p. 367) distinguished two main cases, determined by the sign of the discriminant (and hence the number of real roots) of the quadratic function


 * $$f(x) = b_2\,x^2 + b_1\,x + b_0.

\qquad (2)\!$$

Case 1, negative discriminant: The Pearson type IV distribution
If the discriminant of the quadratic function (2) is negative ($$b_1^2 - 4 b_2 b_0 < 0$$), it has no real roots. Then define


 * $$y = x + \frac{b_1}{2\,b_2} \!$$ and


 * $$\alpha = \frac{\sqrt{4\,b_2\,b_0 - b_1^2\,}}{2\,b_2}. \!$$

Observe that $$\alpha$$ is a well-defined real number and $$\alpha \neq 0$$, because by assumption $$4 b_2 b_0 - b_1^2 > 0$$ and therefore $$b_2 \neq 0$$. Applying these substitutions, the quadratic function (2) is transformed into


 * $$f(x) = b_2\,(y^2 + \alpha^2). \!$$

The absence of real roots is obvious from this formulation, because $$\alpha^2$$ is necessarily positive.

We now express the solution to the differential equation (1) as a function of y:


 * $$p(y) \propto \exp\left(-

\frac{1}{b_2}\, \int\frac{y - \frac{b_1}{2\,b_2} - a_0}{y^2 + \alpha^2} \,\mathrm{d}y \right). \!$$

Pearson (1895, p. 362) called this the "trigonometrical case", because the integral


 * $$\int\frac{y - \frac{2\,b_2\,a_0 + b_1}{2\,b_2}}{y^2 + \alpha^2} \,\mathrm{d}y

= \frac{1}{2} \ln(y^2 + \alpha^2) - \frac{2\,b_2\,a_0 + b_1}{2\,b_2\,\alpha} \arctan\left(\frac{y}{\alpha}\right) + C_0 \!$$

involves the inverse trigonometic arctan function. Then


 * $$p(y) \propto \exp\left[

-\frac{1}{2\,b_2} \ln\!\left(1+\frac{y^2}{\alpha^2}\right) -\frac{\ln\alpha}{2\,b_2} +\frac{2\,b_2\,a_0 + b_1}{2\,b_2^2\,\alpha} \arctan\left(\frac{y}{\alpha}\right) + C_1 \right] \!$$

Finally, let


 * $$m = \frac{1}{2\,b_2} \!$$ and


 * $$\nu = -\frac{2\,b_2\,a_0 + b_1}{2\,b_2^2\,\alpha} \!$$

Applying these substitutions, we obtain the parametric function:


 * $$p(y) \propto

\left[1 + \frac{y^2}{\alpha^2}\right]^{-m} \exp\left[-\nu \arctan\left(\frac{y}{\alpha}\right)\right] \!$$

This unnormalized density has support on the entire real line. It depends on a scale parameter $$\alpha > 0$$ and shape parameters $$m>1/2$$ and $$\nu$$. One parameter was lost when we chose to find the solution to the differential equation (1) as a function of y rather than x. We therefore reintroduce a fourth parameter, namely the location parameter &lambda;. We have thus derived the density of the Pearson type IV distribution:


 * $$p(x) =

\frac{\left|\frac{\Gamma\!\left(m+\frac{\nu}{2}i\right)}{\Gamma(m)}\right|^2} {\alpha\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)} \left[1 + \left(\frac{x-\lambda}{\alpha}\right)^{\!2\,} \right]^{-m} \exp\left[-\nu \arctan\left(\frac{x-\lambda}{\alpha}\right)\right]. \!$$

The normalizing constant involves the complex Gamma function (&Gamma;) and the Beta function (B).

The Pearson type VII distribution


The shape parameter &nu; of the Pearson type IV distribution controls its skewness. If we fix its value at zero, we obtain a symmetric three-parameter family. This special case is known as the Pearson type VII distribution (cf. Pearson 1916, p. 450). Its density is


 * $$p(x) =

\frac{1}{\alpha\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)} \left[1 + \left(\frac{x-\lambda}{\alpha}\right)^{\!2\,} \right]^{-m}, \!$$

where B is the Beta function.

An alternative parameterization (and slight specialization) of the type VII distribution is obtained by letting


 * $$\alpha = \sigma\,\sqrt{2\,m-3}, \!$$

which requires $$m > 3/2$$. This entails a minor loss of generality but ensures that the variance of the distribution exists and is equal to $$\sigma^2$$. Now the parameter m only controls the kurtosis of the distribution. If m approaches infinity as &lambda; and &sigma; are held constant, the normal distribution arises as a special case:


 * $$\lim_{m\to\infty}

\frac{1}{\sigma\,\sqrt{2\,m-3}\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)} \left[1 + \left(\frac{x-\lambda}{\sigma\,\sqrt{2\,m-3}}\right)^{\!2\,} \right]^{-m} \!$$


 * $$= \frac{1}{\sigma\,\sqrt{2}\,\Gamma\!\left(\frac12\right)}

\times \lim_{m\to\infty} \frac{\Gamma(m)}{\Gamma\!\left(m-\frac12\right) \sqrt{m-\frac32}} \times \lim_{m\to\infty} \left[1 + \frac{\left(\frac{x-\lambda}{\sigma}\right)^2}{2\,m-3} \right]^{-m} \!$$


 * $$= \frac{1}{\sigma\sqrt{2\,\pi}}

\times 1 \times \exp\!\left[-\frac12 \left(\frac{x-\lambda}{\sigma}\right)^{\!2\,} \right] \!$$

This is the density of a normal distribution with mean &lambda; and standard deviation &sigma;.

It is convenient to require that $$m > 5/2$$ and to let


 * $$m = \frac52 + \frac{3}{\gamma_2}. \!$$

This is another specialization, and it guarantees that the first four moments of the distribution exist. More specifically, the Pearson type VII distribution parameterized in terms of $$(\lambda,\sigma,\gamma_2)$$ has a mean of &lambda;, standard deviation of &sigma;, skewness of zero, and excess kurtosis of $$\gamma_2$$.

Student's t-distribution
The Pearson type VII distribution subsumes Student's t-distribution, and hence also the Cauchy distribution. Student's t-distribution arises as the result of applying the following substitutions to its original parameterization:


 * $$\lambda = 0, \!$$


 * $$\alpha = \sqrt{\nu}, \!$$ and


 * $$m = \frac{\nu+1}{2}, \!$$

where $$\nu > 0$$. Observe that the constraint $$m > 1/2$$ is satisfied. The density of this restricted one-parameter family is


 * $$p(x) =

\frac{1}{\sqrt{\nu}\,\mathrm{\Beta}\!\left(\frac{\nu}{2}, \frac12\right)} \left[1 + \frac{x^2}{\nu} \right]^{-\frac{\nu+1}{2}}, \!$$

which is easily recognized as the density of Student's t-distribution.

Case 2, non-negative discriminant
If the quadratic function (2) has a non-negative discriminant ($$b_1^2 - 4 b_2 b_0 \geq 0$$), it has real roots r1 and r2 (not necessarily distinct):


 * $$r_{1,2} = \frac{-b_1 \pm \sqrt{b_1^2 - 4 b_2 b_0}}{2 b_2}, \!$$

XXXX We will assume without loss of generality that $$r_1 \leq r_2$$. XXXX (is that really necessary?)

In the presence of real roots the quadratic function (2) can be written as


 * $$f(x) = b_2\,(x-r_1)(x-r_2), \!$$

and the solution to the differential equation is therefore


 * $$p(x) \propto \exp\left( -\frac{1}{b_2} \int\!\!\frac{x-a_0}{(x - r_1) (x - r_2)} \,\mathrm{d}x \right). \!$$

Pearson (1895, p. 362) called this the "logarithmic case", because the integral


 * $$\int\!\!\frac{x-a_0}{(x - r_1) (x - r_2)} \,\mathrm{d}x

= \frac{(r_1-a_0)\ln(x-r_1) - (r_2-a_0)\ln(x-r_2)}{r_1-r_2} + C \!$$

involves only the logarithm function, and not the arctan function as in the previous case.

Using the substitution


 * $$\nu = \frac{1}{b_2\,(r_1-r_2)} \!$$

we obtain the following solution to the differential equation (1):


 * $$p(x) \propto

(x-r_1)^{-\nu (r_1-a_0)} (x-r_2)^{\nu (r_2-a_0)}. \!$$

Since this density is only known up to a hidden constant of proportionality, that constant can be changed and the density written as follows:


 * $$p(x) \propto

\left(1-\frac{x}{r_1}\right)^{-\nu (r_1-a_0)} \left(1-\frac{x}{r_2}\right)^{ \nu (r_2-a_0)} \!$$

XXXXXXXXXXXXXXXXXXXXXXXXX

The Pearson type I distribution arises when the roots of the quadratic equation (2) are of opposite sign, that is, $$r_1 < 0 < r_2$$. Then the solution p is supported on the interval $$(r_1, r_2)$$. Apply the substition


 * $$x = r_1 + y (r_2 - r_1) \qquad \mbox{where}\ 0-1 \land n>-1$$ is necessary and sufficient for p to be a proper probability density function.

beta distribution

The Pearson type II distribution is a special case of the Pearson type I family restricted to symmetric distributions.

Pearson type III distribution gamma distribution, chi-squared distribution

Pearson type V distribution inverse-gamma distribution

Pearson type VI distribution beta prime distribution, F-distribuion

Relation to other distributions
The Pearson family subsumes the following distributions, among others:


 * beta distribution (type I)
 * beta prime distribution (type VI)
 * Cauchy distribution (type IV)
 * chi-squared distribution (type III)
 * continuous uniform distribution (limit of type I)
 * exponential distribution (type III)
 * gamma distribution (type III)
 * F-distribuion (type VI)
 * inverse-chi-squared distribution (type V)
 * inverse-gamma distribution (type V)
 * normal distribution (limit of type I, III, IV, V, or VI)
 * Student's t-distribution (type IV)

OLD
The Pearson family is divided into several types whose shape depends on the values of its parameters. The following cases are distinguished:  $$b_2 = 0$$. This is type III.  $$b_2 \neq 0$$. Then
 * $$b_2 = b_1 = 0$$. This yields the normal distribution with mean a and variance b0.
 * $$b_2 x^2 + b_1 x + b_0 = 0 \qquad (1)\!$$

is a proper quadratic equation and can have the following kinds of roots:  
 * No real roots. This is type IV.
 * No real roots and furthermore $$a = b_1 = 0$$. This is type VII, a special case of type IV.
 * Exactly one real root. This is type V.
 * Two real roots.
 * Two real roots of the opposite sign. This is type I.
 * Real roots r and &minus;r. This is type II, a special case of type I.
 * Two real roots of the same sign. This is type VI.

More old stuff
Pearson described five types of continuous distributions:


 * I. Limited range in both directions, with skewness.
 * II. Limited range in both directions, symmetric.
 * III. Limited range in one direction only (must be skew).
 * IV. Unlimited range in both directions, with skewness.
 * V. Unlimited range in both directions, symmetric.

Generalizing the hypergeometric distribution, Pearson proposed a probability density proportional to:

$$(1+x/a_1)^{\nu a_1} (1-x/a_2)^{\nu a_2}$$

for $$-a_1 < x < a_2 $$, and by taking various limits of this derived forms for his Types I, II, III, and V. For Type IV he derived a related form:

$$exp(-\nu \arctan(x/a))/(1+x^2/a^2)^m $$

These can be shifted to the desired location. Some of his forms correspond to other named distributions.

Applications
These models are used in financial markets, given their ability to be parametrised in a way that has intuitive meaning for market traders. A number of models are in current use that capture the stochastic nature of the volatility of rates, stocks etc. and this family of distributions may prove to be one of the more important.

Pearson type III distribution
The Pearson Type III distribution is given by the probability density function


 * $$ f(x) = \frac{1}{\beta\,\Gamma(p)} \left(\frac{x-\alpha}{\beta}\right)^{p-1} e^{-(x-\alpha)/\beta}, \!$$

where x &isin; [&alpha;,&infin;) and &alpha;, &beta; and p are parameters of the distribution with &beta; > 0 and p > 0 (Abramowitz and Stegun 1954, p. 930). Here, &Gamma; denotes the Gamma function.


 * Characteristic function:
 * $$e^{i\alpha t}(1 - i\beta t)^{-p}$$


 * Mean:
 * $$\alpha + p\beta$$


 * Variance:
 * $$p\beta^2$$


 * Skewness:
 * $$\frac{2}{\sqrt{p}}$$


 * Kurtosis:
 * $$\frac{6}{p}$$

When α=0, β=2, and p is half-integer, the Pearson Type III distribution becomes the chi squared distribution of 2p degrees of freedom.