Hermite distribution

In probability theory and statistics, the Hermite distribution, named after Charles Hermite, is a discrete probability distribution used to model count data with more than one parameter.

This distribution is flexible in terms of its ability to allow a moderate over-dispersion in the data.

The authors Kemp and Kemp [1] have called it "Hermite distribution" from the fact its probability function and the moment generating function can be expressed in terms of the coefficients of (modified) Hermite polynomials.

The distribution first appeared in the paper Applications of Mathematics to Medical Problems,[2] by Anderson Gray McKendrick in 1926.

In this work the author explains several mathematical methods that can be applied to medical research.

As a practical application, McKendrick considered the distribution of counts of bacteria in leucocytes.

Using the method of moments he fitted the data with the Hermite distribution and found the model more satisfactory than fitting it with a Poisson distribution.

The distribution was formally introduced and published by C. D. Kemp and Adrienne W. Kemp in 1965 in their work Some Properties of ‘Hermite’ Distribution.

The work is focused on the properties of this distribution for instance a necessary condition on the parameters and their maximum likelihood estimators (MLE), the analysis of the probability generating function (PGF) and how it can be expressed in terms of the coefficients of (modified) Hermite polynomials.

An example they have used in this publication is the distribution of counts of bacteria in leucocytes that used McKendrick but Kemp and Kemp estimate the model using the maximum likelihood method.

[3][4] The same authors published in 1966 the paper An alternative Derivation of the Hermite Distribution.

In 1971, Y. C. Patel[6] did a comparative study of various estimation procedures for the Hermite distribution in his doctoral thesis.

In 1974, Gupta and Jain[7] did a research on a generalized form of Hermite distribution.

Let X1 and X2 be two independent Poisson variables with parameters a1 and a2.

The probability distribution of the random variable Y = X1 + 2X2 is the Hermite distribution with parameters a1 and a2 and probability mass function is given by [8] where The probability generating function of the probability mass is,[8] When a random variable Y = X1 + 2X2 is distributed by an Hermite distribution, where X1 and X2 are two independent Poisson variables with parameters a1 and a2, we write The moment generating function of a random variable X is defined as the expected value of et, as a function of the real parameter t. For an Hermite distribution with parameters X1 and X2, the moment generating function exists and is equal to The cumulant generating function is the logarithm of the moment generating function and is equal to [4] If we consider the coefficient of (it)rr!

in the expansion of K(t) we obtain the r-cumulant Hence the mean and the succeeding three moments about it are The skewness is the third moment centered around the mean divided by the 3/2 power of the standard deviation, and for the hermite distribution is,[4] The kurtosis is the fourth moment centered around the mean, divided by the square of the variance, and for the Hermite distribution is,[4] The excess kurtosis is just a correction to make the kurtosis of the normal distribution equal to zero, and it is the following, In a discrete distribution the characteristic function of any real-valued random variable is defined as the expected value of

Given a sample X1, ..., Xm are independent random variables each having an Hermite distribution we wish to estimate the value of the parameters

Using these two equation, We can parameterize the probability function by μ and d Hence the log-likelihood function is,[9] where From the log-likelihood function, the likelihood equations are,[9] Straightforward calculations show that,[9] where

The likelihood equation does not always have a solution like as it shows the following proposition, Proposition:[9] Let X1, ..., Xm come from a generalized Hermite distribution with fixed n. Then the MLEs of the parameters are

indicates the empirical factorial momement of order 2.

A usual choice for discrete distributions is the zero relative frequency of the data set which is equated to the probability of zero under the assumed distribution.

Following the example of Y. C. Patel (1976) the resulting system of equations, We obtain the zero frequency and the mean estimator a1 of

When Hermite distribution is used to model a data sample is important to check if the Poisson distribution is enough to fit the data.

Following the parametrized probability mass function used to calculate the maximum likelihood estimator, is important to corroborate the following hypothesis, The likelihood-ratio test statistic [9] for hermite distribution is, Where

As d = 1 belongs to the boundary of the domain of parameters, under the null hypothesis, W does not have an asymptotic

It can be established that the asymptotic distribution of W is a 50:50 mixture of the constant 0 and the

The score statistic is,[9] where m is the number of observations.

The asymptotic distribution of the score test statistic under the null hypothesis is a

It may be convenient to use a signed version of the score test, that is,

, following asymptotically a standard normal.

Example of a multi-modal data, Hermite distribution(0.1,1.5).