Log-normal distribution explained

Log-normal distribution
Type:continuous
Pdf Image:
Identical parameter

\mu

but differing parameters

\sigma

Cdf Image:

\mu=0 

Notation:

\operatorname{Lognormal}\left(\mu,\sigma2 \right)

Parameters:

\mu\in( -infty,+infty)

(logarithm of location),

\sigma>0 

(logarithm of scale)
Support:

x\in( 0,+infty)

Pdf:
1
x\sigma\sqrt{2\pi

}\exp\left(-

\left(lnx-\mu\right)2
2\sigma2

\right)

Cdf:
 1 
2

\left[1+\operatorname{erf}\left(

 lnx-\mu
\sigma\sqrt{2 
} \right)\right] = \Phi\left(\frac \right)
Quantile:

\exp\left(\mu+\sqrt{2\sigma2}\operatorname{erf}-1(2p-1)\right)

Mean:

\exp\left(\mu+

\sigma2
2

\right)

Median:

\exp(\mu)

Mode:

\exp\left(\mu-\sigma2 \right)

Variance:

\left[\exp(\sigma2)-1 \right]\exp\left(2 \mu+\sigma2\right)

Skewness:

\left[\exp\left(\sigma2\right)+2 \right]\sqrt{\exp(\sigma2)-1}

Kurtosis:

 1 \exp\left(4 \sigma2\right)+2 \exp\left(3 \sigma2\right)+3 \exp\left(2\sigma2\right)-6 

Entropy:

 log2\left(\sqrt{2\pi}\sigmae{2}}\right)

Mgf: defined only for numbers with a
 non-positive real part, see text
Char: representation
infty
\sum
n=0
(it)n
n!
n\mu+n2\sigma2/2
e


 is asymptotically divergent, but adequate
 for most numerical purposes
Fisher:

\begin{pmatrix}

1
\sigma2 

&0\ 0&

2
\sigma2 

\end{pmatrix}

Moments:

\mu=log\left(

\operatorname{E
[X]}{\sqrt{
\operatorname{Var
[X]

~~}{\operatorname{E}[X]2 }+1 }}\right),



\sigma=\sqrt{log\left(

\operatorname{Var
[X]

~~}{\operatorname{E}[X]2 }+1 \right)}

Es:
\operatorname{erfc
\left(
s
\sqrt{2 

}-\operatorname{erf}-1(2p-1)\right)}2(1-p)

\mu+
~s2 
2
e

[1]

In probability theory, a log-normal (or lognormal) distribution is a continuous probability distribution of a random variable whose logarithm is normally distributed. Thus, if the random variable is log-normally distributed, then has a normal distribution.[2] [3] Equivalently, if has a normal distribution, then the exponential function of, has a log-normal distribution. A random variable which is log-normally distributed takes only positive real values. It is a convenient and useful model for measurements in exact and engineering sciences, as well as medicine, economics and other topics (e.g., energies, concentrations, lengths, prices of financial instruments, and other metrics).

The distribution is occasionally referred to as the Galton distribution or Galton's distribution, after Francis Galton. The log-normal distribution has also been associated with other names, such as McAlister, Gibrat and Cobb–Douglas.

A log-normal process is the statistical realization of the multiplicative product of many independent random variables, each of which is positive. This is justified by considering the central limit theorem in the log domain (sometimes called Gibrat's law). The log-normal distribution is the maximum entropy probability distribution for a random variate —for which the mean and variance of are specified.[4]

Definitions

Generation and parameters

Let

Z

be a standard normal variable, and let

\mu

and

\sigma

be two real numbers, with

\sigma>0

. Then, the distribution of the random variable

X=e\mu

is called the log-normal distribution with parameters

\mu

and

\sigma

. These are the expected value (or mean) and standard deviation of the variable's natural logarithm, not the expectation and standard deviation of

X

itself.

This relationship is true regardless of the base of the logarithmic or exponential function: If

 loga(X)

is normally distributed, then so is

 logb(X)

for any two positive numbers

a,b1~.

Likewise, if

eY

is log-normally distributed, then so is

aY,

where

0<a1

.

In order to produce a distribution with desired mean

\muX

and variance
2 ,
\sigma
X
one uses

\mu=ln\left(

2
\mu
X
\sqrt{
2
\mu
X
+
2 
\sigma
X

}\right)

and

\sigma2=ln\left(1+

2 
\sigma
X
2
\mu
X

\right)~.

Alternatively, the "multiplicative" or "geometric" parameters

\mu*=e\mu

and

\sigma*=e\sigma

can be used. They have a more direct interpretation:

\mu*

is the median of the distribution, and

\sigma*

is useful for determining "scatter" intervals, see below.

Probability density function

A positive random variable

X

is log-normally distributed (i.e.,

X\sim\operatorname{Lognormal}\left(\mu,\sigma2 \right)

), if the natural logarithm of

X

is normally distributed with mean

\mu

and variance

\sigma2 :

ln(X)\siml{N}(\mu,\sigma2)

Let

\Phi

and

\varphi

be respectively the cumulative probability distribution function and the probability density function of the

l{N}( 0,1 )

standard normal distribution, then we have that[2] the probability density function of the log-normal distribution is given by:

\beginf_X(x) & = \frac\ \operatorname\,\!\bigl[\ X \le x\ \bigr] \\[6pt]& = \frac\ \operatorname\,\!\bigl[\ \ln X \le \ln x\ \bigr] \\[6pt]& = \frac \operatorname\!\!\left(\frac \right) \\[6pt]& = \operatorname\!\left(\frac \sigma \right) \frac \left(\frac\right) \\[6pt]& = \operatorname\!\left(\frac \right) \frac \\[6pt]& = \frac \exp\left(-\frac \right) ~.\end

Cumulative distribution function

The cumulative distribution function is

FX(x)=\Phi\left(

(lnx)-\mu
\sigma

\right)

where

\Phi

is the cumulative distribution function of the standard normal distribution (i.e.,

\operatornamel{N}( 0, 1)

).

This may also be expressed as follows:

12
\left[

1+\operatorname{erf}\left(

lnx-\mu
\sigma\sqrt{2
}\right) \right] = \frac12 \operatorname \left(-\frac\right)

where is the complementary error function.

Multivariate log-normal

If

\boldsymbolX\siml{N}(\boldsymbol\mu,\boldsymbol\Sigma)

is a multivariate normal distribution, then

Yi=\exp(Xi)

has a multivariate log-normal distribution.[5] [6] The exponential is applied elementwise to the random vector

\boldsymbolX

. The mean of

\boldsymbolY

is

\operatorname{E}[\boldsymbol

\mu\Sigmaii
i+1
2
Y]
i=e

,

and its covariance matrix is

\operatorname{Var}[\boldsymbolY]ij

\muj+
1
2
(\Sigmaii+\Sigmajj)
i+\mu
=e

(

\Sigmaij
e

-1).

Since the multivariate log-normal distribution is not widely used, the rest of this entry only deals with the univariate distribution.

Characteristic function and moment generating function

All moments of the log-normal distribution exist and

\operatorname{E}[Xn]=

n\mu+n2\sigma2/2
e

This can be derived by letting

z=\tfrac{ln(x)-(\mu+n\sigma2)}{\sigma}

within the integral. However, the log-normal distribution is not determined by its moments. This implies that it cannot have a defined moment generating function in a neighborhood of zero.[7] Indeed, the expected value

\operatorname{E}[et]

is not defined for any positive value of the argument

t

, since the defining integral diverges.

\operatorname{E}[ei]

is defined for real values of, but is not defined for any complex value of that has a negative imaginary part, and hence the characteristic function is not analytic at the origin. Consequently, the characteristic function of the log-normal distribution cannot be represented as an infinite convergent series.[8] In particular, its Taylor formal series diverges:
infty
\sum
n=0
(it)n
n!
n\mu+n2\sigma2/2
e

However, a number of alternative divergent series representations have been obtained.[9] [10] [11]

A closed-form formula for the characteristic function

\varphi(t)

with

t

in the domain of convergence is not known. A relatively simple approximating formula is available in closed form, and is given by[12]
\varphi(t)
\exp\left(-W2(-it\sigma2e\mu)+2W(-it\sigma2e\mu)\right)
2\sigma2
\sqrt{1+W(-it\sigma2e\mu)
}

where

W

is the Lambert W function. This approximation is derived via an asymptotic method, but it stays sharp all over the domain of convergence of

\varphi

.

Properties

Probability in different domains

The probability content of a log-normal distribution in any arbitrary domain can be computed to desired precision by first transforming the variable to normal, then numerically integrating using the ray-trace method.[13] (Matlab code)

Probabilities of functions of a log-normal variable

Since the probability of a log-normal can be computed in any domain, this means that the cdf (and consequently pdf and inverse cdf) of any function of a log-normal variable can also be computed.[13] (Matlab code)

Geometric or multiplicative moments

The geometric or multiplicative mean of the log-normal distribution is

\operatorname{GM}[X]=e\mu=\mu*

. It equals the median. The geometric or multiplicative standard deviation is

\operatorname{GSD}[X]=e\sigma=\sigma*

.[14] [15]

By analogy with the arithmetic statistics, one can define a geometric variance,

\operatorname{GVar}[X]=

\sigma2
e
, and a geometric coefficient of variation,

\operatorname{GCV}[X]=e\sigma-1

, has been proposed. This term was intended to be analogous to the coefficient of variation, for describing multiplicative variation in log-normal data, but this definition of GCV has no theoretical basis as an estimate of

\operatorname{CV}

itself (see also Coefficient of variation).

Note that the geometric mean is smaller than the arithmetic mean. This is due to the AM–GM inequality and is a consequence of the logarithm being a concave function. In fact,

\operatorname{E}[X]=

\mu+
12
\sigma
2
e

=e\mu

\sigma2
\sqrt{e
} = \operatorname[X] \cdot \sqrt.[16]

In finance, the term

-12\sigma
2
e
is sometimes interpreted as a convexity correction. From the point of view of stochastic calculus, this is the same correction term as in Itō's lemma for geometric Brownian motion.

Arithmetic moments

For any real or complex number, the -th moment of a log-normally distributed variable is given by

\operatorname{E}[Xn]=

n\mu+
12n
2\sigma
2
e

.

Specifically, the arithmetic mean, expected square, arithmetic variance, and arithmetic standard deviation of a log-normally distributed variable are respectively given by:

\begin{align} \operatorname{E}[X]&=e\mu{2}\sigma2},\\[4pt] \operatorname{E}[X2]&=

2\mu+2\sigma2
e

,\\[4pt] \operatorname{Var}[X]&=\operatorname{E}[X2]-\operatorname{E}[X]2=(\operatorname{E}[X])2(e

\sigma2

-1)=

2\mu+\sigma2
e
\sigma2
(e

-1),\\[4pt] \operatorname{SD}[X]&=\sqrt{\operatorname{Var}[X]}=\operatorname{E}[X]

\sigma2
\sqrt{e

-1} =e\mu{2}\sigma2}\sqrt{e

\sigma2

-1}, \end{align}

\operatorname{CV}[X]

is the ratio

\tfrac{\operatorname{SD}[X]}{\operatorname{E}[X]}

. For a log-normal distribution it is equal to

\operatorname{CV}[X]=

\sigma2
\sqrt{e

-1}.

This estimate is sometimes referred to as the "geometric CV" (GCV),[17] [18] due to its use of the geometric variance. Contrary to the arithmetic standard deviation, the arithmetic coefficient of variation is independent of the arithmetic mean.

The parameters and can be obtained, if the arithmetic mean and the arithmetic variance are known:

\begin{align} \mu&=ln\left(

\operatorname{E
[X]

2}{\sqrt{\operatorname{E}[X2]}}\right)=ln\left(

\operatorname{E
[X]

2}{\sqrt{\operatorname{Var}[X]+\operatorname{E}[X]2}}\right),\\[4pt] \sigma2&=ln\left(

\operatorname{E
[X

2]}{\operatorname{E}[X]2}\right)=ln\left(1+

\operatorname{Var
[X]}{\operatorname{E}[X]

2}\right). \end{align}

A probability distribution is not uniquely determined by the moments for . That is, there exist other distributions with the same set of moments. In fact, there is a whole family of distributions with the same moments as the log-normal distribution.

Mode, median, quantiles

The mode is the point of global maximum of the probability density function. In particular, by solving the equation

(lnf)'=0

, we get that:

\operatorname{Mode}[X]=

\mu-\sigma2
e

.

Since the log-transformed variable

Y=lnX

has a normal distribution, and quantiles are preserved under monotonic transformations, the quantiles of

X

are

qX(\alpha)=

\mu+\sigmaq\Phi(\alpha)
e

=\mu*(\sigma*)

q\Phi(\alpha)

,

where

q\Phi(\alpha)

is the quantile of the standard normal distribution.

Specifically, the median of a log-normal distribution is equal to its multiplicative mean,[19]

\operatorname{Med}[X]=e\mu=\mu*~.

Partial expectation

The partial expectation of a random variable

X

with respect to a threshold

k

is defined as

g(k)=

infty
\int
k

xfX(x\midX>k)dx.

Alternatively, by using the definition of conditional expectation, it can be written as

g(k)=\operatorname{E}[X\midX>k]P(X>k)

. For a log-normal random variable, the partial expectation is given by:

g(k)=

infty
\int
k

xfX(x\midX>k)dx=e\mu+\tfrac{1{2}\sigma2}\Phi\left(

\mu+\sigma2-lnk
\sigma

\right)

where

\Phi

is the normal cumulative distribution function. The derivation of the formula is provided in the Talk page. The partial expectation formula has applications in insurance and economics, it is used in solving the partial differential equation leading to the Black–Scholes formula.

Conditional expectation

The conditional expectation of a log-normal random variable

X

—with respect to a threshold

k

—is its partial expectation divided by the cumulative probability of being in that range:

\begin{align} E[X\midX<k]&

\mu
+\sigma2
2
=e

\Phi
\left[ln(k)-\mu-\sigma2
\sigma
\right]
\Phi
\left[ln(k)-\mu
\sigma
\right]

\\[8pt] E[X\midX\geqslantk]

\mu
+\sigma2
2
&=e

\Phi
\left[\mu+\sigma2-ln(k)
\sigma
\right]
1-\Phi
\left[ln(k)-\mu
\sigma
\right]

\\ [8pt] E[X\midX\in[k1,k2]]

\mu
+\sigma2
2
&=e

\Phi
\left[ln(k2)-\mu-\sigma2
\sigma
\right]-\Phi
\left[ln(k1)-\mu-\sigma2
\sigma
\right]
\Phi
\left[ln(k2)-\mu
\sigma
\right]-\Phi
\left[ln(k1)-\mu
\sigma
\right]

\end{align}

Alternative parameterizations

In addition to the characterization by

\mu,\sigma

or

\mu*,\sigma*

, here are multiple ways how the log-normal distribution can be parameterized. ProbOnto, the knowledge base and ontology of probability distributions[20] [21] lists seven such forms:
P(x;\boldsymbol\mu,\boldsymbol\sigma)=1
x\sigma\sqrt{2\pi
} \exp\left[-\frac{(\ln x - \mu)^2}{2 \sigma^2}\right]

P(x;\boldsymbol\mu,\boldsymbol{v})=

1
x\sqrt{v

\sqrt{2\pi}}\exp\left[-

(lnx-\mu)2
2v

\right]

P(x;\boldsymbolm,\boldsymbol\sigma)=

1
x\sigma\sqrt{2\pi
} \exp\left[-\frac{\ln^2(x/m)}{2 \sigma^2}\right]

P(x;\boldsymbolm,\boldsymbol{cv})=

1
x\sqrt{ln(cv2+1)

\sqrt{2\pi}}\exp\left[-

ln2(x/m)
2ln(cv2+1)

\right]

P(x;\boldsymbol\mu,\boldsymbol\tau)=\sqrt{

\tau
2\pi
} \frac \exp\left[-\frac{\tau}{2}(\ln x-\mu)^2\right]

P(x;\boldsymbolm,\boldsymbol

{\sigma
g})=1
xln(\sigmag)\sqrt{2\pi
} \exp\left[-\frac{\ln^2(x/m)}{2 \ln^2(\sigma_g)}\right]

P(x;\boldsymbol{\muN},\boldsymbol{\sigmaN})=

1
x\sqrt{2\pi
2\right)
ln\left(1+\sigma
N
} \exp\left(-\frac\right)

Examples for re-parameterization

Consider the situation when one would like to run a model using two different optimal design tools, for example PFIM[26] and PopED.[27] The former supports the LN2, the latter LN7 parameterization, respectively. Therefore, the re-parameterization is required, otherwise the two tools would produce different results.

For the transition

\operatorname{LN2}(\mu,v)\to\operatorname{LN7}(\muN,\sigmaN)

following formulas hold \mu_N = \exp(\mu+v/2) and \sigma_N = \exp(\mu+v/2)\sqrt.

For the transition

\operatorname{LN7}(\muN,\sigmaN)\to\operatorname{LN2}(\mu,v)

following formulas hold \mu = \ln\left(\mu_N / \sqrt \right) and v = \ln(1+\sigma_N^2/\mu_N^2).

All remaining re-parameterisation formulas can be found in the specification document on the project website.[28]

Multiple, reciprocal, power

X\sim\operatorname{Lognormal}(\mu,\sigma2)

then

aX\sim\operatorname{Lognormal}(\mu+lna,\sigma2)

for

a>0.

X\sim\operatorname{Lognormal}(\mu,\sigma2)

then

\tfrac{1}{X}\sim\operatorname{Lognormal}(-\mu,\sigma2).

X\sim\operatorname{Lognormal}(\mu,\sigma2)

then

Xa\sim\operatorname{Lognormal}(a\mu,a2\sigma2)

for

a0.

Multiplication and division of independent, log-normal random variables

If two independent, log-normal variables

X1

and

X2

are multiplied [divided], the product [ratio] is again log-normal, with parameters

\mu=\mu1+\mu2

[<math>\mu=\mu_1-\mu_2</math>] and

\sigma

, where
2
\sigma
2
. This is easily generalized to the product of

n

such variables.

More generally, if

Xj\sim\operatorname{Lognormal}(\muj,

2)
\sigma
j
are

n

independent, log-normally distributed variables, then

Y=

n
style\prod
j=1

Xj\sim\operatorname{Lognormal}

n\mu
(style\sum
j,\sum
n
j=1
2
\sigma
j

).

Multiplicative central limit theorem

See also: Gibrat's law.

The geometric or multiplicative mean of

n

independent, identically distributed, positive random variables

Xi

shows, for

n\toinfty

, approximately a log-normal distribution with parameters

\mu=E[ln(Xi)]

and

\sigma2=var[ln(Xi)]/n

, assuming

\sigma2

is finite.

In fact, the random variables do not have to be identically distributed. It is enough for the distributions of

ln(Xi)

to all have finite variance and satisfy the other conditions of any of the many variants of the central limit theorem.

This is commonly known as Gibrat's law.

Other

A set of data that arises from the log-normal distribution has a symmetric Lorenz curve (see also Lorenz asymmetry coefficient).[29]

The harmonic

H

, geometric

G

and arithmetic

A

means of this distribution are related;[30] such relation is given by

H=

G2
A.

Log-normal distributions are infinitely divisible, but they are not stable distributions, which can be easily drawn from.

Related distributions

X\siml{N}(\mu,\sigma2)

is a normal distribution, then

\exp(X)\sim\operatorname{Lognormal}(\mu,\sigma2).

X\sim\operatorname{Lognormal}(\mu,\sigma2)

is distributed log-normally, then

ln(X)\siml{N}(\mu,\sigma2)

is a normal random variable.

Xj\sim\operatorname{Lognormal}(\muj,

2)
\sigma
j
be independent log-normally distributed variables with possibly varying

\sigma

and

\mu

parameters, and Y = \sum_^n X_j. The distribution of

Y

has no closed-form expression, but can be reasonably approximated by another log-normal distribution

Z

at the right tail.[31] Its probability density function at the neighborhood of 0 has been characterized and it does not resemble any log-normal distribution. A commonly used approximation due to L.F. Fenton (but previously stated by R.I. Wilkinson and mathematically justified by Marlow[32]) is obtained by matching the mean and variance of another log-normal distribution: \begin \sigma^2_Z &= \ln\!\left[\frac{\sum e^{2\mu_j+\sigma_j^2}(e^{\sigma_j^2}-1)}{(\sum e^{\mu_j+\sigma_j^2/2})^2} + 1\right], \\ \mu_Z &= \ln\!\left[\sum e^{\mu_j+\sigma_j^2/2} \right] - \frac. \end In the case that all

Xj

have the same variance parameter

\sigmaj=\sigma

, these formulas simplify to \begin \sigma^2_Z &= \ln\!\left[(e^{\sigma^2}-1)\frac{\sum e^{2\mu_j}}{(\sum e^{\mu_j})^2} + 1\right], \\ \mu_Z &= \ln\!\left[\sum e^{\mu_j} \right] + \frac - \frac. \endFor a more accurate approximation, one can use the Monte Carlo method to estimate the cumulative distribution function, the pdf and the right tail.[33] [34]

The sum of correlated log-normally distributed random variables can also be approximated by a log-normal distribution \begin S_+ &= \operatorname\left[\sum_i X_i \right] = \sum_i \operatorname[X_i] = \sum_i e^ \\ \sigma^2_ &= 1/S_+^2 \, \sum_ \operatorname_ \sigma_i \sigma_j \operatorname[X_i] \operatorname[X_j] = 1/S_+^2 \, \sum_ \operatorname_ \sigma_i \sigma_j e^ e^ \\ \mu_Z &= \ln\left(S_+ \right) - \sigma_^2/2 \end

X\sim\operatorname{Lognormal}(\mu,\sigma2)

then

X+c

is said to have a Three-parameter log-normal distribution with support

x\in(c,+infty)

.[35]

\operatorname{E}[X+c]=\operatorname{E}[X]+c

,

\operatorname{Var}[X+c]=\operatorname{Var}[X]

.

X\midY\sim\operatorname{Rayleigh}(Y)

with

Y\sim\operatorname{Lognormal}(\mu,\sigma2)

, then

X\sim\operatorname{Suzuki}(\mu,\sigma)

(Suzuki distribution).

Statistical inference

Estimation of parameters

For determining the maximum likelihood estimators of the log-normal distribution parameters μ and σ, we can use the same procedure as for the normal distribution. Note thatL(\mu, \sigma) = \prod_^n \frac 1 \varphi_ (\ln x_i),where

\varphi

is the density function of the normal distribution

lN(\mu,\sigma2)

. Therefore, the log-likelihood function is\ell (\mu,\sigma \mid x_1, x_2, \ldots, x_n) = - \sum _i \ln x_i + \ell_N (\mu, \sigma \mid \ln x_1, \ln x_2, \dots, \ln x_n).

Since the first term is constant with regard to μ and σ, both logarithmic likelihood functions,

\ell

and

\ellN

, reach their maximum with the same

\mu

and

\sigma

. Hence, the maximum likelihood estimators are identical to those for a normal distribution for the observations

lnx1,lnx2,...,lnxn)

,\widehat \mu = \frac, \qquad \widehat \sigma^2 = \frac .

For finite n, the estimator for

\mu

is unbiased, but the one for

\sigma

is biased. As for the normal distribution, an unbiased estimator for

\sigma

can be obtained by replacing the denominator n by n−1 in the equation for

\widehat\sigma2

.

When the individual values

x1,x2,\ldots,xn

are not available, but the sample's mean

\barx

and standard deviation s is, then the Method of moments can be used. The corresponding parameters are determined by the following formulas, obtained from solving the equations for the expectation

\operatorname{E}[X]

and variance

\operatorname{Var}[X]

for

\mu

and

\sigma

: \mu = \ln\left(\frac \right), \qquad \sigma^2 = \ln\left(1 + / \bar x^2 \right).

Interval estimates

The most efficient way to obtain interval estimates when analyzing log-normally distributed data consists of applying the well-known methods based on the normal distribution to logarithmically transformed data and then to back-transform results if appropriate.

Prediction intervals

A basic example is given by prediction intervals: For the normal distribution, the interval

[\mu-\sigma,\mu+\sigma]

contains approximately two thirds (68%) of the probability (or of a large sample), and

[\mu-2\sigma,\mu+2\sigma]

contain 95%. Therefore, for a log-normal distribution,[\mu^*/\sigma^*,\mu^*\cdot\sigma^*]=[\mu^* {}^\times\!\!/ \sigma^*] contains 2/3, and[\mu^*/(\sigma^*)^2,\mu^*\cdot(\sigma^*)^2] = [\mu^* {}^\times\!\!/ (\sigma^*)^2] contains 95% of the probability. Using estimated parameters, then approximately the same percentages of the data should be contained in these intervals.

Confidence interval for μ*

Using the principle, note that a confidence interval for

\mu

is

[\widehat\mu\pmq\widehatse]

, where

se=\widehat\sigma/\sqrt{n}

is the standard error and q is the 97.5% quantile of a t distribution with n-1 degrees of freedom. Back-transformation leads to a confidence interval for

\mu*

(the median), is:[\widehat\mu^* {}^\times\!\!/ (\operatorname{sem}^*)^q] with

\operatorname{sem}*=(\widehat\sigma*)1/\sqrt{n

}

Confidence interval for μ

The literature discusses several options for calculating the confidence interval for

\mu

(the mean of the log-normal distribution). These include bootstrap as well as various other methods.[38] [39]

Extremal principle of entropy to fix the free parameter σ

In applications,

\sigma

is a parameter to be determined. For growing processes balanced by production and dissipation, the use of an extremal principle of Shannon entropy shows that[40] \sigma = \frac 1 \sqrt

This value can then be used to give some scaling relation between the inflexion point and maximum point of the log-normal distribution.[40] This relationship is determined by the base of natural logarithm,

e=2.718\ldots

, and exhibits some geometrical similarity to the minimal surface energy principle.These scaling relations are useful for predicting a number of growth processes (epidemic spreading, droplet splashing, population growth, swirling rate of the bathtub vortex, distribution of language characters, velocity profile of turbulences, etc.).For example, the log-normal function with such

\sigma

fits well with the size of secondarily produced droplets during droplet impact and the spreading of an epidemic disease.[41]

The value \sigma = 1 \big/ \sqrt is used to provide a probabilistic solution for the Drake equation.[42]

Occurrence and applications

The log-normal distribution is important in the description of natural phenomena. Many natural growth processes are driven by the accumulation of many small percentage changes which become additive on a log scale. Under appropriate regularity conditions, the distribution of the resulting accumulated changes will be increasingly well approximated by a log-normal, as noted in the section above on "Multiplicative Central Limit Theorem". This is also known as Gibrat's law, after Robert Gibrat (1904–1980) who formulated it for companies.[43] If the rate of accumulation of these small changes does not vary over time, growth becomes independent of size. Even if this assumption is not true, the size distributions at any age of things that grow over time tends to be log-normal. Consequently, reference ranges for measurements in healthy individuals are more accurately estimated by assuming a log-normal distribution than by assuming a symmetric distribution about the mean.

A second justification is based on the observation that fundamental natural laws imply multiplications and divisions of positive variables. Examples are the simple gravitation law connecting masses and distance with the resulting force, or the formula for equilibrium concentrations of chemicals in a solution that connects concentrations of educts and products. Assuming log-normal distributions of the variables involved leads to consistent models in these cases.

Specific examples are given in the following subsections.[44] contains a review and table of log-normal distributions from geology, biology, medicine, food, ecology, and other areas. is a review article on log-normal distributions in neuroscience, with annotated bibliography.

Human behavior

Biology and medicine

Chemistry

Hydrology

The image on the right, made with CumFreq, illustrates an example of fitting the log-normal distribution to ranked annually maximum one-day rainfalls showing also the 90% confidence belt based on the binomial distribution.[63]

The rainfall data are represented by plotting positions as part of a cumulative frequency analysis.

Social sciences and demographics

\sigma

, then the Gini coefficient, commonly use to evaluate income inequality, can be computed as

G=\operatorname{erf}\left(

\sigma
2

\right)

where

\operatorname{erf}

is the error function, since

G=2\Phi\left(

\sigma
\sqrt{2
}\right)-1, where

\Phi(x)

is the cumulative distribution function of a standard normal distribution.

Technology

See also

Further reading

External links

Notes and References

  1. Norton . Matthew . Khokhlov . Valentyn . Uryasev . Stan . 2019 . Calculating CVaR and bPOE for common probability distributions with application to portfolio optimization and density estimation . . 299 . 1–2 . 1281–1315 . Springer . 10.1007/s10479-019-03373-1 . 1811.11301 . 254231768 . https://web.archive.org/web/20210418151110/http://uryasev.ams.stonybrook.edu/wp-content/uploads/2019/10/Norton2019_CVaR_bPOE.pdf . 2021-04-18 . live . 2023-02-27 . stonybrook.edu .
  2. Web site: Weisstein . Eric W. . Log Normal Distribution . mathworld.wolfram.com . en . 2020-09-13 .
  3. Web site: 1.3.6.6.9. Lognormal Distribution . www.itl.nist.gov . U.S. National Institute of Standards and Technology (NIST) . 2020-09-13 .
  4. Park . Sung Y. . Bera . Anil K. . 2009 . Maximum entropy autoregressive conditional heteroskedasticity model . Journal of Econometrics . 150 . 2 . 219–230, esp. Table 1, p. 221 . 10.1.1.511.9750 . 10.1016/j.jeconom.2008.12.014 . dead . https://web.archive.org/web/20160307144515/http://wise.xmu.edu.cn/uploadfiles/paper-masterdownload/2009519932327055475115776.pdf . 2016-03-07 . 2011-06-02.
  5. Tarmast . Ghasem . 2001 . Multivariate Log–Normal Distribution . https://web.archive.org/web/20130719214220/http://isi.cbs.nl/iamamember/CD2/pdf/329.PDF . 2013-07-19 . live . ISI Proceedings: 53rd Session . Seoul.
  6. Halliwell . Leigh . 2015 . The Lognormal Random Multivariate . https://web.archive.org/web/20150930111908/http://www.casact.org/pubs/forum/15spforum/Halliwell.pdf . 2015-09-30 . live . Casualty Actuarial Society E-Forum, Spring 2015 . Arlington, VA.
  7. Book: Billingsley, Patrick. Probability and Measure. 2012. Wiley. 978-1-118-12237-2. Anniversary. Hoboken, N.J.. 415. 780289503.
  8. Holgate . P. . 1989 . The lognormal characteristic function, vol. 18, pp. 4539–4548, 1989 . Communications in Statistics - Theory and Methods . 18 . 12 . 4539–4548 . 10.1080/03610928908830173.
  9. Barakat . R. . 1976 . Sums of independent lognormally distributed random variables . Journal of the Optical Society of America . 66 . 3 . 211–216 . 1976JOSA...66..211B . 10.1364/JOSA.66.000211.
  10. Barouch . E. . Kaufman . GM. . Glasser . ML. . 1986 . On sums of lognormal random variables . Studies in Applied Mathematics . 75 . 1 . 37–55 . 10.1002/sapm198675137 . 1721.1/48703. free .
  11. Leipnik . Roy B. . January 1991 . On Lognormal Random Variables: I – The Characteristic Function . Journal of the Australian Mathematical Society, Series B . 32 . 3 . 327–347 . 10.1017/S0334270000006901. free .
  12. S. Asmussen, J.L. Jensen, L. Rojas-Nandayapa (2016). "On the Laplace transform of the Lognormal distribution",Methodology and Computing in Applied Probability 18 (2), 441-458.Thiele report 6 (13).
  13. Das. Abhranil. 2012.14331. A method to integrate and classify normal distributions. Journal of Vision . 2021. 21 . 10 . 1 . 10.1167/jov.21.10.1 . 34468706 . 8419883 .
  14. Kirkwood . Thomas BL . Geometric means and measures of dispersion . Biometrics . Dec 1979 . 35 . 4 . 908–9. 2530139 .
  15. Limpert . E . Stahel . W . Abbt . M . Lognormal distributions across the sciences: keys and clues . BioScience . 2001 . 51 . 5 . 341–352 . 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2. free .
  16. Heil P. Friedrich B. Onset-Duration Matching of Acoustic Stimuli Revisited: Conventional Arithmetic vs. Proposed Geometric Measures of Accuracy and Precision. Frontiers in Psychology. 7. 2013. 10.3389/fpsyg.2016.02013. 28111557. 5216879. 2017. free.
  17. Sawant, S.; Mohan, N. (2011) "FAQ: Issues with Efficacy Analysis of Clinical Trial Data Using SAS", PharmaSUG2011, Paper PO08
  18. Schiff . MH . etal . 2014 . Head-to-head, randomised, crossover study of oral versus subcutaneous methotrexate in patients with rheumatoid arthritis: drug-exposure limitations of oral methotrexate at doses >=15 mg may be overcome with subcutaneous administration. Ann Rheum Dis . 73. 8. 1–3 . 10.1136/annrheumdis-2014-205228 . 24728329 . 4112421.
  19. Book: Leslie E. . Daly . Geoffrey Joseph . Bourke . 2000 . Interpretation and Uses of Medical Statistics . Journal of Epidemiology and Community Health,Journal of Epidemiology and Community Health . 46,46 . 3,3 --- not found in WorldCat.org --> . 5th . Oxford, UK . Wiley-Blackwell . 978-0-632-04763-5 . 89 . 10.1002/9780470696750 . 1059583 . . print edition. Online eBook
  20. Web site: ProbOnto . 1 July 2017.
  21. 27153608 . 10.1093/bioinformatics/btw170 . 5013898 . 32 . 17 . 2719–21 . ProbOnto: ontology and knowledge base of probability distributions . 2016 . Bioinformatics . Swat . MJ . Grenon . P . Wimalaratne . S.
  22. Forbes et al. Probability Distributions (2011), John Wiley & Sons, Inc.
  23. Lunn, D. (2012). The BUGS book: a practical introduction to Bayesian analysis. Texts instatistical science. CRC Press.
  24. Limpert . E. . Stahel . W. A. . Abbt . M. . 2001 . Log-normal distributions across the sciences: Keys and clues . BioScience . 51 . 5. 341–352 . 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 . free .
  25. Nyberg . J. . etal . 2012 . PopED - An extended, parallelized, population optimal design tool . Comput Methods Programs Biomed . 108 . 2. 789–805 . 10.1016/j.cmpb.2012.05.005 . 22640817 .
  26. Retout . S . Duffull . S . Mentré . F . 2001 . Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs . Comp Meth Pro Biomed . 65 . 2. 141–151 . 10.1016/S0169-2607(00)00117-6 . 11275334 .
  27. The PopED Development Team (2016). PopED Manual, Release version 2.13. Technical report, Uppsala University.
  28. ProbOnto website, URL: http://probonto.org
  29. 10.1890/0012-9658(2000)081[1139:DIIPSO]2.0.CO;2 . Damgaard . Christian . Jacob . Weiner . Describing inequality in plant size or fecundity . Ecology . 2000 . 81 . 4 . 1139–1142.
  30. Rossman . Lewis A . July 1990 . Design stream flows based on harmonic means . Journal of Hydraulic Engineering . 116 . 7 . 946–950 . 10.1061/(ASCE)0733-9429(1990)116:7(946).
  31. S. . Asmussen . L. . Rojas-Nandayapa . Asymptotics of Sums of Lognormal Random Variables with Gaussian Copula . Statistics and Probability Letters . 78 . 16 . 2709–2714 . 2008 . 10.1016/j.spl.2008.03.035 .
  32. NA. . Marlow . A normal limit theorem for power sums of independent normal random variables . Bell System Technical Journal . 46 . 9 . 2081–2089 . Nov 1967 . 10.1002/j.1538-7305.1967.tb04244.x.
  33. Accurate computation of the right tail of the sum of dependent log-normal variates . Botev . Z. I. . L'Ecuyer . P. . 2017 . IEEE. 978-1-5386-3428-8. 2017 Winter Simulation Conference (WSC), 3rd–6th Dec 2017 . 1880–1890 . Las Vegas, NV, USA. 10.1109/WSC.2017.8247924 . 1705.03196.
  34. Asmussen. A.. Goffard. P.-O. . Laub. P. J.. 2016 . Orthonormal polynomial expansions and lognormal sum densities. 1601.01763v1. math.PR.
  35. B. . Sangal . A. . Biswas . The 3-Parameter Lognormal Distribution Applications in Hydrology . Water Resources Research . 6 . 2 . 505–515 . 1970 . 10.1029/WR006i002p00505.
  36. Norman Lloyd Johnson . Johnson . N. L. . 1949 . Systems of Frequency Curves Generated by Methods of Translation. . 36 . 1/2 . 149–176 . 2332539 . 10.2307/2332539. 18132090 .
  37. Swamee . P. K. . Near Lognormal Distribution . 10.1061/(ASCE)1084-0699(2002)7:6(441) . Journal of Hydrologic Engineering . 7 . 6 . 441–444 . 2002 .
  38. Olsson, Ulf. "Confidence intervals for the mean of a log-normal distribution." Journal of Statistics Education 13.1 (2005).pdf html
  39. user10525, How do I calculate a confidence interval for the mean of a log-normal data set?, URL (version: 2022-12-18): https://stats.stackexchange.com/q/33395
  40. Wu. Ziniu. Li. Juan. Bai. Chenyuan. Scaling Relations of Lognormal Type Growth Process with an Extremal Principle of Entropy. Entropy. 19. 56. 2017. 1–14. 10.3390/e19020056. 2017Entrp..19...56W. free.
  41. Wang. WenBin. Wu. ZiNiu. Wang. ChunFeng. Hu. RuiFeng. Modelling the spreading rate of controlled communicable epidemics through an entropy-based thermodynamic model. Science China Physics, Mechanics and Astronomy. 56. 11. 2013. 2143–2150. 1674-7348. 10.1007/s11433-013-5321-0. 32288765. 7111546. 1304.5603. 2013SCPMA..56.2143W.
  42. Bloetscher. Frederick . Using predictive Bayesian Monte Carlo- Markov Chain methods to provide a probabilistic solution for the Drake equation . Acta Astronautica. 155. 2019. 118–130 . 10.1016/j.actaastro.2018.11.033 . 2019AcAau.155..118B . 117598888.
  43. 2729692. Sutton . John . Mar 1997 . Gibrat's Legacy . Journal of Economic Literature . 32 . 1 . 40–59.
  44. Limpert . Eckhard . Stahel . Werner A. . Abbt . Markus . 2001 . Log-normal Distributions across the Sciences: Keys and Clues . BioScience . en . 51 . 5 . 341 . 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 . 0006-3568. free .
  45. Pawel . Sobkowicz. Lognormal distributions of user post lengths in Internet discussions - a consequence of the Weber-Fechner law? . EPJ Data Science . 2013. etal.
  46. Yin. Peifeng . Luo. Ping . Lee. Wang-Chien . Wang. Min . Silence is also evidence: interpreting dwell time for recommendation from psychological perspective . ACM International Conference on KDD . 2013 .
  47. Web site: What is the average length of a game of chess?. chess.stackexchange.com. 14 April 2018.
  48. Book: Huxley, Julian S. . 1932 . Problems of relative growth . London . 476909537 . 978-0-486-61114-3 .
  49. Sartwell, Philip E. "The distribution of incubation periods of infectious disease." American journal of hygiene 51 (1950): 310-318.
  50. S. K. Chan . Jennifer . Yu . Philip L. H. . Modelling SARS data using threshold geometric process . Statistics in Medicine . 2006 . 25 . 11 . 1826–1839 . 10.1002/sim.2376 . 16345017 . 46599163 .
  51. Ono. Yukiteru. Asai. Kiyoshi. Hamada. Michiaki. 2013-01-01. PBSIM: PacBio reads simulator—toward accurate genome assembly. Bioinformatics. en. 29. 1. 119–121. 10.1093/bioinformatics/bts649. 23129296. 1367-4803. free.
  52. Makuch. Robert W. . D.H. Freeman . M.F. Johnson. Justification for the lognormal distribution as a model for blood pressure. Journal of Chronic Diseases. 1979. 32. 3. 245–250. 10.1016/0021-9681(79)90070-5. 429469 .
  53. Lacey. L. F.. Keene. O. N.. Pritchard. J. F.. Bye. A.. 1997-01-01. Common noncompartmental pharmacokinetic variables: are they normally or log-normally distributed?. Journal of Biopharmaceutical Statistics. en. 7. 1. 171–178. 10.1080/10543409708835177. 9056596. 1054-3406.
  54. Scheler. Gabriele. Schumann. Johann. Diversity and stability in neuronal output rates. 36th Society for Neuroscience Meeting, Atlanta. 2006-10-08.
  55. Mizuseki. Kenji. Buzsáki. György. 2013-09-12. Preconfigured, skewed distribution of firing rates in the hippocampus and entorhinal cortex. Cell Reports. 4. 5. 1010–1021. 10.1016/j.celrep.2013.07.039. 2211-1247. 3804159. 23994479.
  56. Buzsáki. György. Mizuseki. Kenji. 2017-01-06. The log-dynamic brain: how skewed distributions affect network operations. Nature Reviews. Neuroscience. 15. 4. 264–278. 10.1038/nrn3687. 1471-003X. 4051294. 24569488.
  57. Wohrer. Adrien. Humphries. Mark D.. Machens. Christian K.. 2013-04-01. Population-wide distributions of neural activity during perceptual decision-making. Progress in Neurobiology. 103. 156–193. 10.1016/j.pneurobio.2012.09.004. 1873-5118. 23123501. 5985929.
  58. Scheler. Gabriele. Logarithmic distributions prove that intrinsic learning is Hebbian. F1000Research. 10.12688/f1000research.12130.2. 2017-07-28. 29071065. 6. 5639933. 1222 . free .
  59. Morales-Gregorio . Aitor . van Meegen . Alexander . van Albada . Sacha . 2023 . Ubiquitous lognormal distribution of neuron densities in mammalian cerebral cortex . Cerebral Cortex . 33 . 16 . 9439–9449 . 10.1093/cercor/bhad160 . 37409647 . 10438924 .
  60. Polizzi . Stefano . Laperrousaz . Bastien . Perez-Reche . Francisco J . Nicolini . Franck E . Satta . Véronique Maguer . Arneodo . Alain . Argoul . Françoise . 2018-05-29 . A minimal rupture cascade model for living cell plasticity . New Journal of Physics . 20 . 5 . 053057 . 10.1088/1367-2630/aac3c7 . 2018NJPh...20e3057P . 1367-2630. 2164/10561 . free .
  61. Ahrens . L. H. . 1954-02-01 . The lognormal distribution of the elements (A fundamental law of geochemistry and its subsidiary) . Geochimica et Cosmochimica Acta . en . 5 . 2 . 49–73 . 10.1016/0016-7037(54)90040-X . 1954GeCoA...5...49A . 0016-7037.
  62. Book: Oosterbaan, R.J. . Ritzema . H.P. . 6: Frequency and Regression Analysis . 1994 . Drainage Principles and Applications, Publication 16 . International Institute for Land Reclamation and Improvement (ILRI) . Wageningen, The Netherlands . 175–224 . http://www.waterlog.info/pdf/freqtxt.pdf . 978-90-70754-33-4 .
  63. https://www.waterlog.info/cumfreq.htm CumFreq, free software for distribution fitting
  64. Clementi, Fabio; Gallegati, Mauro (2005) "Pareto's law of income distribution: Evidence for Germany, the United Kingdom, and the United States", EconWPA
  65. cond-mat/0202388 . Physics of Personal Income. Wataru. Souma. 2002-02-22. Springer . Empirical Science of Financial Fluctuations: The Advent of Econophysics . 10.1007/978-4-431-66993-7 . Takayasu . Hideki .
  66. 10.1086/260062. The Pricing of Options and Corporate Liabilities. Journal of Political Economy. 81. 3. 637. 1973. Black . F. . Scholes . M. . 154552078.
  67. Book: Mandelbrot, Benoit. The (mis-)Behaviour of Markets . 2004 . Basic Books . 9780465043552.
  68. Bunchen, P., Advanced Option Pricing, University of Sydney coursebook, 2007
  69. Thelwall. Mike. Wilson. Paul . Regression for citation data: An evaluation of different methods . Journal of Informetrics . 2014 . 8 . 4 . 963–971 . 10.1016/j.joi.2014.09.011. 1510.08877. 8338485.
  70. Sheridan. Paul. Onodera. Taku . A Preferential Attachment Paradox: How Preferential Attachment Combines with Growth to Produce Networks with Log-normal In-degree Distributions . Scientific Reports . 2020 . 8 . 1 . 2811 . 10.1038/s41598-018-21133-2. 29434232. 5809396. 1703.06645.
  71. Eeckhout. Jan. 2004. Gibrat's Law for (All) Cities. American Economic Review. 94. 5. 1429–1451. 10.1257/0002828043052303. 3592829. JSTOR.
  72. Kault . David . The Shape of the Distribution of the Number of Sexual Partners . Statistics in Medicine . 1996 . 15 . 2 . 221–230. 10.1002/(SICI)1097-0258(19960130)15:2<221::AID-SIM148>3.0.CO;2-Q . 8614756 .
  73. Book: O'Connor . Patrick . Kleyner . Andre . 2011 . Practical Reliability Engineering . John Wiley & Sons . 978-0-470-97982-2 . 35 .
  74. Web site: Shadowing . www.WirelessCommunication.NL . dead . https://web.archive.org/web/20120113201345/http://wireless.per.nl/reference/chaptr03/shadow/shadow.htm . January 13, 2012 .
  75. Dexter . A. R. . Tanner . D. W. . July 1972 . Packing Densities of Mixtures of Spheres with Log-normal Size Distributions . Nature Physical Science . en . 238 . 80 . 31–32 . 10.1038/physci238031a0 . 1972NPhS..238...31D . 2058-1106.
  76. Gros . C . Kaczor . G.. Markovic . D . Neuropsychological constraints to human data production on a global scale . The European Physical Journal B . 2012 . 85 . 28 . 28 . 10.1140/epjb/e2011-20581-3. 1111.6849 . 2012EPJB...85...28G . 17404692 .
  77. Douceur . John R. . Bolosky . William J. . 1999-05-01 . A large-scale study of file-system contents . ACM SIGMETRICS Performance Evaluation Review . 27 . 1 . 59–70 . 10.1145/301464.301480 . 0163-5999. free .
  78. Alamsar . Mohammed . Parisis . George . Clegg . Richard . Zakhleniuk . Nickolay . 2019 . On the Distribution of Traffic Volumes in the Internet and its Implications . 1902.03853 . cs.NI .
  79. ASTM D3654, Standard Test Method for Shear Adhesion on Pressure-Sensitive Tapesw
  80. ASTM D4577, Standard Test Method for Compression Resistance of a container Under Constant Load>\