# Log-normal distribution

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 X is log-normally distributed, then Y = ln(X) has a normal distribution. Likewise, if Y has a normal distribution, then the exponential function of Y, X = exp(Y), 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 fields, e.g. for energies, concentrations, lengths, (financial) returns and other amounts.

Notation Probability density functionSome log-normal density functions with identical parameter ${\displaystyle \mu }$ but differing parameters ${\displaystyle \sigma }$ Cumulative distribution functionCumulative distribution function of the log-normal distribution (with ${\displaystyle \mu =0}$ ) ${\displaystyle \operatorname {Lognormal} (\mu ,\,\sigma ^{2})}$ ${\displaystyle \mu \in (-\infty ,+\infty )}$, ${\displaystyle \sigma >0}$ ${\displaystyle x\in (0,+\infty )}$ ${\displaystyle {\frac {1}{x\sigma {\sqrt {2\pi }}}}\ \exp \left(-{\frac {\left(\ln x-\mu \right)^{2}}{2\sigma ^{2}}}\right)}$ ${\displaystyle {\frac {1}{2}}+{\frac {1}{2}}\operatorname {erf} {\Big [}{\frac {\ln x-\mu }{{\sqrt {2}}\sigma }}{\Big ]}}$ ${\displaystyle \exp(\mu +{\sqrt {2\sigma ^{2}}}\operatorname {erf} ^{-1}(2F-1))}$ ${\displaystyle \exp \left(\mu +{\frac {\sigma ^{2}}{2}}\right)}$ ${\displaystyle \exp(\mu )}$ ${\displaystyle \exp(\mu -\sigma ^{2})}$ ${\displaystyle [\exp(\sigma ^{2})-1]\exp(2\mu +\sigma ^{2})}$ ${\displaystyle (e^{\sigma ^{2}}\!\!+2){\sqrt {e^{\sigma ^{2}}\!\!-1}}}$ ${\displaystyle \exp(4\sigma ^{2})+2\exp(3\sigma ^{2})+3\exp(2\sigma ^{2})-6}$ ${\displaystyle \log _{2}(\sigma e^{\mu +{\tfrac {1}{2}}}{\sqrt {2\pi }})}$ defined only for numbers with a non-positive real part, see text representation ${\displaystyle \sum _{n=0}^{\infty }{\frac {(it)^{n}}{n!}}e^{n\mu +n^{2}\sigma ^{2}/2}}$ is asymptotically divergent but sufficient for numerical purposes ${\displaystyle {\begin{pmatrix}1/\sigma ^{2}&0\\0&1/2\sigma ^{4}\end{pmatrix}}}$

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

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. The log-normal distribution is the maximum entropy probability distribution for a random variate X for which the mean and variance of ln(X) are specified.[2]

## Definitions

### Generation and parameters

Let ${\displaystyle Z}$ be a standard normal variable, and let ${\displaystyle \mu }$ and ${\displaystyle \sigma >0}$ be two real numbers. Then, the distribution of the random variable

${\displaystyle X=e^{\mu +\sigma Z}}$

is called the log-normal distribution with parameters ${\displaystyle \mu }$ and ${\displaystyle \sigma }$. Thus, these are, the expected value (or mean) and standard deviation of the variable's natural logarithm, not the expectation and standard deviation of ${\displaystyle X}$ itself.

This relationship is true regardless of the base of the logarithmic or exponential function. If ${\displaystyle \log _{a}(X)}$ is normally distributed, then so is ${\displaystyle \log _{b}(X)}$, for any two positive numbers ${\displaystyle a,b\neq 1}$. Likewise, if ${\displaystyle e^{Y}}$ is log-normally distributed, then so is ${\displaystyle a^{Y}}$, where ${\displaystyle 0.

Alternatively, the "multiplicative" or "geometric" parameters ${\displaystyle \mu ^{*}=e^{\mu }}$ and ${\displaystyle \sigma ^{*}=e^{\sigma }}$ can be used. They have a more direct interpretation: ${\displaystyle \mu ^{*}}$ is the median of the distribution, and ${\displaystyle \sigma ^{*}}$ is useful for determining "scatter" intervals, see below.

### Probability density function

A positive random variable X is log-normally distributed if the logarithm of X is normally distributed,

${\displaystyle \ln(X)\sim {\mathcal {N}}(\mu ,\sigma ^{2}).}$

Let ${\displaystyle \Phi }$ and ${\displaystyle \varphi }$ be respectively the cumulative probability distribution function and the probability density function of the N(0,1) distribution.

Then we have[1]

{\displaystyle {\begin{aligned}f_{X}(x)&={\frac {\rm {d}}{{\rm {d}}x}}\Pr(X\leq x)={\frac {\rm {d}}{{\rm {d}}x}}\Pr(\ln X\leq \ln x)={\frac {\rm {d}}{{\rm {d}}x}}\Phi \left({\frac {\ln x-\mu }{\sigma }}\right)\\[6pt]&=\varphi \left({\frac {\ln x-\mu }{\sigma }}\right){\frac {\rm {d}}{{\rm {d}}x}}\left({\frac {\ln x-\mu }{\sigma }}\right)=\varphi \left({\frac {\ln x-\mu }{\sigma }}\right){\frac {1}{\sigma x}}\\[6pt]&={\frac {1}{x}}\cdot {\frac {1}{\sigma {\sqrt {2\pi \,}}}}\exp \left(-{\frac {(\ln x-\mu )^{2}}{2\sigma ^{2}}}\right).\end{aligned}}}

### Cumulative distribution function

${\displaystyle F_{X}(x)=\Phi \left({\frac {(\ln x)-\mu }{\sigma }}\right)}$

where ${\displaystyle \Phi }$ is the cumulative distribution function of the standard normal distribution (i.e. N(0,1)).

This may also be expressed as follows:

${\displaystyle {\frac {1}{2}}\left[1+\operatorname {erf} \left({\frac {\ln x-\mu }{\sigma {\sqrt {2}}}}\right)\right]={\frac {1}{2}}\operatorname {erfc} \left(-{\frac {\ln x-\mu }{\sigma {\sqrt {2}}}}\right)}$

where erfc is the complementary error function.

### Multivariate log-normal

If ${\displaystyle {\boldsymbol {X}}\sim {\mathcal {N}}({\boldsymbol {\mu }},\,{\boldsymbol {\Sigma }})}$ is a multivariate normal distribution then ${\displaystyle {\boldsymbol {Y}}=\exp({\boldsymbol {X}})}$ has a multivariate log-normal distribution[3][4] with mean

${\displaystyle \operatorname {E} [{\boldsymbol {Y}}]_{i}=e^{\mu _{i}+{\frac {1}{2}}\Sigma _{ii}},}$
${\displaystyle \operatorname {Var} [{\boldsymbol {Y}}]_{ij}=e^{\mu _{i}+\mu _{j}+{\frac {1}{2}}(\Sigma _{ii}+\Sigma _{jj})}(e^{\Sigma _{ij}}-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

${\displaystyle \operatorname {E} [X^{n}]=e^{n\mu +n^{2}\sigma ^{2}/2}}$

This can be derived by letting ${\displaystyle z={\frac {\ln(x)-(\mu +n\sigma ^{2})}{\sigma }}}$ within the integral. However, the expected value ${\displaystyle \operatorname {E} [e^{tX}]}$ is not defined for any positive value of the argument ${\displaystyle t}$ as the defining integral diverges. In consequence the moment generating function is not defined.[5] The last is related to the fact that the lognormal distribution is not uniquely determined by its moments.

The characteristic function ${\displaystyle \operatorname {E} [e^{itX}]}$ is defined for real values of t but is not defined for any complex value of t that has a negative imaginary part, and therefore the characteristic function is not analytic at the origin. In consequence, the characteristic function of the log-normal distribution cannot be represented as an infinite convergent series.[6] In particular, its Taylor formal series diverges:

${\displaystyle \sum _{n=0}^{\infty }{\frac {(it)^{n}}{n!}}e^{n\mu +n^{2}\sigma ^{2}/2}}$

However, a number of alternative divergent series representations have been obtained[6][7][8][9]

A closed-form formula for the characteristic function ${\displaystyle \varphi (t)}$ with ${\displaystyle t}$ in the domain of convergence is not known. A relatively simple approximating formula is available in closed form and given by[10]

${\displaystyle \varphi (t)\approx {\frac {\exp \left(-{\frac {W^{2}(-it\sigma ^{2}e^{\mu })+2W(-it\sigma ^{2}e^{\mu })}{2\sigma ^{2}}}\right)}{\sqrt {1+W(-it\sigma ^{2}e^{\mu })}}}}$

where ${\displaystyle 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 ${\displaystyle \varphi }$.

## Properties

### Geometric or multiplicative moments

The geometric or multiplicative mean of the log-normal distribution is ${\displaystyle \operatorname {GM} [X]=e^{\mu }=\mu ^{*}}$. It equals the median. The geometric or multiplicative standard deviation is ${\displaystyle \operatorname {GSD} [X]=e^{\sigma }=\sigma ^{*}}$.[11][12] By analogy with the arithmetic statistics, one can define a geometric variance, ${\displaystyle \operatorname {GVar} [X]=e^{\sigma ^{2}}}$, and a geometric coefficient of variation,[11] ${\displaystyle \operatorname {GCV} [X]=e^{\sigma }-1}$.

Note that the geometric mean is smaller than the arithmetic mean. This is due to the AM–GM inequality, and corresponds to the logarithm being convex down. In fact,

${\displaystyle \operatorname {E} [X]=e^{\mu +{\frac {1}{2}}\sigma ^{2}}=e^{\mu }\cdot {\sqrt {e^{\sigma ^{2}}}}=\operatorname {GM} [X]\cdot {\sqrt {\operatorname {GVar} [X]}}.}$[13]

In finance the term ${\displaystyle e^{-{\frac {1}{2}}\sigma ^{2}}}$ 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 n, the n-th moment of a log-normally distributed variable X is given by[1]

${\displaystyle \operatorname {E} [X^{n}]=e^{n\mu +{\frac {1}{2}}n^{2}\sigma ^{2}}.}$

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

{\displaystyle {\begin{aligned}\operatorname {E} [X]&=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}},\\[4pt]\operatorname {E} [X^{2}]&=e^{2\mu +2\sigma ^{2}},\\[4pt]\operatorname {Var} [X]&=\operatorname {E} [X^{2}]-\operatorname {E} [X]^{2}=(\operatorname {E} [X])^{2}(e^{\sigma ^{2}}-1)=e^{2\mu +\sigma ^{2}}(e^{\sigma ^{2}}-1),\\[4pt]\operatorname {SD} [X]&={\sqrt {\operatorname {Var} [X]}}=\operatorname {E} [X]{\sqrt {e^{\sigma ^{2}}-1}}=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}}{\sqrt {e^{\sigma ^{2}}-1}},\end{aligned}}}

respectively.

The arithmetic coefficient of variation ${\displaystyle \operatorname {CV} [X]}$ is the ratio ${\displaystyle {\frac {\operatorname {SD} [X]}{\operatorname {E} [X]}}}$. For a log-normal distribution it is equal to

${\displaystyle \operatorname {CV} [X]={\sqrt {e^{\sigma ^{2}}-1}}.}$

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:

{\displaystyle {\begin{aligned}\mu &=\ln \left({\frac {\operatorname {E} [X]^{2}}{\sqrt {\operatorname {E} [X^{2}]}}}\right)=\ln \left({\frac {\operatorname {E} [X]^{2}}{\sqrt {\operatorname {Var} [X]+\operatorname {E} [X]^{2}}}}\right),\\[4pt]\sigma ^{2}&=\ln \left({\frac {\operatorname {E} [X^{2}]}{\operatorname {E} [X]^{2}}}\right)=\ln \left(1+{\frac {\operatorname {Var} [X]}{\operatorname {E} [X]^{2}}}\right).\end{aligned}}}

A probability distribution is not uniquely determined by the moments E[Xn] = e + 1/2n2σ2 for n ≥ 1. That is, there exist other distributions with the same set of moments.[1] 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, it solves the equation ${\displaystyle (\ln f)'=0}$:

${\displaystyle \operatorname {Mode} [X]=e^{\mu -\sigma ^{2}}.}$

Since the log-transformed variable ${\displaystyle Y=\ln X}$ has a normal distribution and quantiles are preserved under monotonic transformations, the quantiles of ${\displaystyle X}$ are

${\displaystyle q_{X}(\alpha )=e^{\mu +\sigma q_{\Phi }(\alpha )}=\mu ^{*}(\sigma ^{*})^{q_{\Phi }(\alpha )},}$

where ${\displaystyle 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,[14]

${\displaystyle \operatorname {Med} [X]=e^{\mu }=\mu ^{*}.}$

### Partial expectation

The partial expectation of a random variable ${\displaystyle X}$ with respect to a threshold ${\displaystyle k}$ is defined as

${\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x)\,dx.}$

Alternatively, and using the definition of conditional expectation, it can be written as ${\displaystyle g(k)=\operatorname {E} [X\mid X>k]P(X>k)}$. For a log-normal random variable the partial expectation is given by:

${\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x)\,dx=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}}\,\Phi \!\left({\frac {\mu +\sigma ^{2}-\ln k}{\sigma }}\right)}$

where ${\displaystyle \Phi }$ is the normal cumulative distribution function. The derivation of the formula is provided in the discussion of this Wikipedia entry. 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 lognormal random variable ${\displaystyle X}$ with respect to a threshold ${\displaystyle k}$ is its partial expectation divided by the cumulative probability of being in that range:

{\displaystyle {\begin{aligned}E[X\mid X

### Alternative parameterizations

In addition to the characterization by ${\displaystyle \mu ,\sigma }$ or ${\displaystyle \mu ^{*},\sigma ^{*}}$, here are multiple ways how the lognormal distribution can be parameterized. ProbOnto, the knowledge base and ontology of probability distributions[15][16] lists seven such forms:

• Normal1(μ,σ) with mean, μ, and standard deviation, σ [17]
${\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\sigma }})={\frac {1}{\sigma {\sqrt {2\pi }}}}\exp \left[-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right]}$
• LogNormal2(μ,υ) with mean, μ, and variance, υ, both on the log-scale
${\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {v}})={\frac {1}{x{\sqrt {v}}{\sqrt {2\pi }}}}\exp \left[{\frac {-(\log x-\mu )^{2}}{2v}}\right]}$
• LogNormal3(m,σ) with median, m, on the natural scale and standard deviation, σ, on the log-scale[17]
${\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma }})={\frac {1}{x\sigma {\sqrt {2\pi }}}}\exp \left[{\frac {-[\log(x/m)]^{2}}{2\sigma ^{2}}}\right]}$
• LogNormal4(m,cv) with median, m, and coefficient of variation, cv, both on the natural scale
${\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {cv}})={\frac {1}{x{\sqrt {\log(cv^{2}+1)}}{\sqrt {2\pi }}}}\exp \left[{\frac {-[\log(x/m)]^{2}}{2\log(cv^{2}+1)}}\right]}$
• LogNormal5(μ,τ) with mean, μ, and precision, τ, both on the log-scale[18]
${\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\tau }})={\sqrt {\frac {\tau }{2\pi }}}{\frac {1}{x}}\exp \left[{-{\frac {\tau }{2}}(\log x-\mu )^{2}}\right]}$
• LogNormal6(m,σg) with median, m, and geometric standard deviation, σg, both on the natural scale[19]
${\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma _{g}}})={\frac {1}{x\log(\sigma _{g}){\sqrt {2\pi }}}}\exp \left[{\frac {-[\log(x/m)]^{2}}{2\log ^{2}(\sigma _{g})}}\right]}$
• LogNormal7(μNN) with mean, μN, and standard deviation, σN, both on the natural scale[20]
${\displaystyle P(x;{\boldsymbol {\mu _{N}}},{\boldsymbol {\sigma _{N}}})={\frac {1}{x{\sqrt {2\pi \log \left(1+\sigma _{N}^{2}/\mu _{N}^{2}\right)}}}}\exp \left({\frac {-{\Big [}\log(x)-\log {\Big (}{\frac {\mu _{N}}{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}}{\Big )}{\Big ]}^{2}}{2\log {\Big (}1+\sigma _{N}^{2}/\mu _{N}^{2}{\Big )}}}\right)}$
##### Examples for re-parameterization

Consider the situation when one would like to run a model using two different optimal design tools, e.g. PFIM[21] and PopED.[22] 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 ${\displaystyle LN2(\mu ,v)\rightarrow LN7(\mu _{N},\sigma _{N})}$ following formulas hold ${\displaystyle \mu _{N}=\exp(\mu +v/2){\text{ and }}\sigma _{N}=\exp(\mu +v/2){\sqrt {\exp(v)-1}}}$.

For the transition ${\displaystyle LN7(\mu _{N},\sigma _{N})\rightarrow LN2(\mu ,v)}$ following formulas hold ${\displaystyle \mu =\log {\Big (}\mu _{N}/{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}{\Big )}{\text{ and }}v=\log(1+\sigma _{N}^{2}/\mu _{N}^{2})}$.

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

### Other

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

The harmonic ${\displaystyle H}$, geometric ${\displaystyle G}$ and arithmetic ${\displaystyle A}$ means of this distribution are related;[25] such relation is given by

${\displaystyle H={\frac {G^{2}}{A}}.}$

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

• If ${\displaystyle X\sim {\mathcal {N}}(\mu ,\sigma ^{2})}$ is a normal distribution, then ${\displaystyle \exp(X)\sim \operatorname {Lognormal} (\mu ,\sigma ^{2}).}$
• If ${\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$ is distributed log-normally, then ${\displaystyle \ln(X)\sim {\mathcal {N}}(\mu ,\sigma ^{2})}$ is a normal random variable.
• If ${\displaystyle X_{j}\sim \operatorname {Lognormal} (\mu _{j},\sigma _{j}^{2})}$ are ${\displaystyle n}$ independent log-normally distributed variables, and ${\displaystyle Y=\textstyle \prod _{j=1}^{n}X_{j}}$, then ${\displaystyle Y}$ is also distributed log-normally:
${\displaystyle Y\sim \operatorname {Lognormal} {\Big (}\textstyle \sum _{j=1}^{n}\mu _{j},\ \sum _{j=1}^{n}\sigma _{j}^{2}{\Big )}.}$
• Let ${\displaystyle X_{j}\sim \operatorname {Lognormal} (\mu _{j},\sigma _{j}^{2})\ }$ be independent log-normally distributed variables with possibly varying ${\displaystyle \sigma }$ and ${\displaystyle \mu }$ parameters, and ${\displaystyle Y=\textstyle \sum _{j=1}^{n}X_{j}}$. The distribution of ${\displaystyle Y}$ has no closed-form expression, but can be reasonably approximated by another log-normal distribution ${\displaystyle Z}$ at the right tail.[28] Its probability density function at the neighborhood of 0 has been characterized[27] 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 mathematical justified by Marlow[29]) is obtained by matching the mean and variance of another lognormal distribution:
{\displaystyle {\begin{aligned}\sigma _{Z}^{2}&=\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 {\sigma _{Z}^{2}}{2}}.\end{aligned}}}
In the case that all ${\displaystyle X_{j}}$ have the same variance parameter ${\displaystyle \sigma _{j}=\sigma }$, these formulas simplify to
{\displaystyle {\begin{aligned}\sigma _{Z}^{2}&=\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 {\sigma ^{2}}{2}}-{\frac {\sigma _{Z}^{2}}{2}}.\end{aligned}}}

For a more accurate approximation one can use the Monte Carlo method to estimate the cumulative distribution function, the pdf and right tail. [30] [31]

• If ${\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$ then ${\displaystyle X+c}$ is said to have a Three-parameter log-normal distribution with support ${\displaystyle x\in (c,+\infty )}$.[32] ${\displaystyle \operatorname {E} [X+c]=\operatorname {E} [X]+c}$, ${\displaystyle \operatorname {Var} [X+c]=\operatorname {Var} [X]}$.
• If ${\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$ then ${\displaystyle aX\sim \operatorname {Lognormal} (\mu +\ln a,\ \sigma ^{2}).}$
• If ${\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$ then ${\displaystyle {\tfrac {1}{X}}\sim \operatorname {Lognormal} (-\mu ,\ \sigma ^{2}).}$
• If ${\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$ then ${\displaystyle X^{a}\sim \operatorname {Lognormal} (a\mu ,\ a^{2}\sigma ^{2})}$ for ${\displaystyle a\neq 0.}$
• Lognormal distribution is a special case of semi-bounded Johnson distribution
• If ${\displaystyle X\mid Y\sim \operatorname {Rayleigh} (Y)\,}$ with ${\displaystyle Y\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}$, then ${\displaystyle X\sim \operatorname {Suzuki} (\mu ,\sigma )}$ (Suzuki distribution)
• A substitute for the log-normal whose integral can be expressed in terms of more elementary functions[33] can be obtained based on the logistic distribution to get an approximation for the CDF
${\displaystyle F(x;\mu ,\sigma )=\left[\left({\frac {e^{\mu }}{x}}\right)^{\pi /(\sigma {\sqrt {3}})}+1\right]^{-1}.}$
This is a log-logistic 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 that

${\displaystyle L(\mu ,\sigma )=\prod _{i=1}^{n}{\frac {1}{x_{i}}}\varphi _{\mu ,\sigma }(\ln x_{i})}$,

where ${\displaystyle \varphi _{}}$ is the density function of the normal distribution ${\displaystyle {\mathcal {N}}(\mu ,\sigma ^{2})}$. Therefore, the log-likelihood function is

${\displaystyle \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, ${\displaystyle \ell }$ and ${\displaystyle \ell _{N}}$, reach their maximum with the same ${\displaystyle \mu }$ and ${\displaystyle \sigma }$. Hence, the maximum likelihood estimators are identical to those for a normal distribution for the observations ${\displaystyle \ln x_{1},\ln x_{2},\dots ,\ln x_{n})}$,

${\displaystyle {\widehat {\mu }}={\frac {\sum _{k}\ln x_{k}}{n}},\qquad {\widehat {\sigma }}^{2}={\frac {\sum _{k}\left(\ln x_{k}-{\widehat {\mu }}\right)^{2}}{n}}.}$

For finite n, these estimators are biased. Whereas the bias for ${\displaystyle {\widehat {\mu }}}$ is negligible, a less biased estimator for ${\displaystyle \sigma }$ is obtained as for the normal distribution by replacing the denominator n by n-1 in the equation for ${\displaystyle {\widehat {\sigma }}^{2}}$.

When the individual values ${\displaystyle x_{1},x_{2},\ldots ,x_{n}}$ are not available, but the sample's mean ${\displaystyle {\bar {x}}}$ and standard deviation s is, then the corresponding parameters are determined by the following formulas, obtained from solving the equations for the expectation ${\displaystyle \operatorname {E} [X]}$ and variance ${\displaystyle \operatorname {var} [X]}$ for ${\displaystyle \mu }$ and ${\displaystyle \sigma }$:

${\displaystyle \mu =\ln \left({\bar {x}}\ {\Big /}\ {\sqrt {1+{\frac {{\widehat {\sigma }}^{2}}{{\bar {x}}^{2}}}}}\right),\qquad \sigma ^{2}=\ln \left(1+{\frac {{\widehat {\sigma }}^{2}}{{\bar {x}}^{2}}}\right)}$.

### Statistics

The most efficient way to analyze 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.

#### Scatter intervals

A basic example is given by scatter intervals: For the normal distribution, the interval ${\displaystyle [\mu -\sigma ,\mu +\sigma ]}$ contains approximately two thirds (68 %) of the probability (or of a large sample), and ${\displaystyle [\mu -2\sigma ,\mu +2\sigma ]}$ contain 95 %. Therefore, for a log-normal distribution,

${\displaystyle [\mu ^{*}/\sigma ^{*},\mu ^{*}\cdot \sigma ^{*}]=[\mu ^{*}{}^{\times }\!\!/\sigma ^{*}]}$ contains 2/3, and
${\displaystyle [\mu ^{*}/(\sigma ^{*})^{2},\mu ^{*}\cdot (\sigma ^{*})^{2}]=[\mu ^{*}{}^{\times }\!\!/(\sigma ^{*})^{2}]}$ contains 95 %

of the probability. Using estimated parameters, the approximately the same percentages of the data should be contained in these intervals.

#### Confidence interval for ${\displaystyle \mu ^{*}}$

Using the principle, note that a confidence interval for ${\displaystyle \mu }$ is ${\displaystyle [{\widehat {\mu }}\pm q\cdot {\widehat {\mathop {se} }}]}$, where ${\displaystyle \mathop {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 ${\displaystyle \mu ^{*}}$,

${\displaystyle [{\widehat {\mu }}^{*}{}^{\times }\!\!/(\operatorname {sem} ^{*})^{q}]}$ with ${\displaystyle \operatorname {sem} ^{*}=({\widehat {\sigma }}^{*})^{1/{\sqrt {n}}}}$

### Extremal principle of entropy to fix the free parameter ${\displaystyle \sigma }$

• In applications, ${\displaystyle \sigma }$ is a parameter to be determined. For growing processes balanced by production and dissipation, the use of a extremal principle of Shannon entropy shows that
{\displaystyle {\begin{aligned}\sigma =1{\big /}{\sqrt {6}}\end{aligned}}} [34]
• This value can then be used to give some scaling relation between the inflexion point and maximum point of the lognormal distribution.[34] It is shown that this relationship is determined by the base of natural logarithm, ${\displaystyle e=2.718\ldots }$, and exhibits some geometrical similarity to the minimal surface energy principle.
• These scaling relations are shown to be 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 instance, the lognormal function with such ${\displaystyle \sigma }$ fits well with the size of secondary produced droplet during droplet impact [35] and the spreading of one epidemic disease.[36]
• The value ${\displaystyle \sigma =1{\big /}{\sqrt {6}}}$ is used to provide a probabilistic solution for the Drake equation.[37]

## Occurrence and applications

The log-normal distribution is important in the description of natural phenomena. This follows, because many natural growth processes are driven by the accumulation of many small percentage changes. These become additive on a log scale. If the effect of any one change is negligible, the central limit theorem says that the distribution of their sum is more nearly normal than that of the summands. When back-transformed onto the original scale, it makes the distribution of sizes approximately log-normal (though if the standard deviation is sufficiently small, the normal distribution can be an adequate approximation).

This multiplicative version of the central limit theorem is also known as Gibrat's law, after Robert Gibrat (1904–1980) who formulated it for companies.[38] If the rate of accumulation of these small changes does not vary over time, growth becomes independent of size. Even if that's not true, the size distributions at any age of things that grow over time tends to be log-normal.

Examples include the following:

• Human behaviors
• The length of comments posted in Internet discussion forums follows a log-normal distribution.[39]
• Users' dwell time on online articles (jokes, news etc.) follows a log-normal distribution.[40]
• The length of chess games tends to follow a log normal distribution.[41]
• Onset durations of acoustic comparison stimuli that are matched to a standard stimulus follow a log-normal distribution.[13]
• Rubik's Cube solves, both general or by person, appear to be following a log-normal distribution[42]
• In biology and medicine,
• Measures of size of living tissue (length, skin area, weight);[43]
• For highly communicable epidemics, such as SARS in 2003, if publication intervention is involved, the number of hospitalized cases is shown to satisfy the lognormal distribution with no free parameters if an entropy is assumed and the standard deviation is determined by the principle of maximum rate of entropy production.
• The length of inert appendages (hair, claws, nails, teeth) of biological specimens, in the direction of growth;
• The normalised RNA-Seq readcount for any genomic region can be well approximated by log-normal distribution.
• Certain physiological measurements, such as blood pressure of adult humans (after separation on male/female subpopulations)[44]
• In neuroscience, the distribution of firing rates across a population of neurons is often approximately lognormal. This has been first observed in the cortex and striatum [45] and later in hippocampus and entorhinal cortex,[46] and elsewhere in the brain.[47][48] Also, intrinsic gain distributions and synaptic weight distributions appear to be lognormal[49] as well.
• In colloidal chemistry and polymer chemistry

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.

• In hydrology, the log-normal distribution is used to analyze extreme values of such variables as monthly and annual maximum values of daily rainfall and river discharge volumes.[50]
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.[51]
The rainfall data are represented by plotting positions as part of a cumulative frequency analysis.
• In social sciences and demographics
• In economics, there is evidence that the income of 97%–99% of the population is distributed log-normally.[52] (The distribution of higher-income individuals follows a Pareto distribution).

[53]

• Technology
• In reliability analysis, the lognormal distribution is often used to model times to repair a maintainable system.[58]
• In wireless communication, "the local-mean power expressed in logarithmic values, such as dB or neper, has a normal (i.e., Gaussian) distribution."[59] Also, the random obstruction of radio signals due to large buildings and hills, called shadowing, is often modeled as a lognormal distribution.
• Particle size distributions produced by comminution with random impacts, such as in ball milling
• The file size distribution of publicly available audio and video data files (MIME types) follows a log-normal distribution over five orders of magnitude.[60]
• In computer networks and Internet traffic analysis, lognormal is shown as a good statistical model to represent the amount of traffic per unit time. This has been shown by applying a robust statistical approach on a large groups of real Interent traces. In this context, the log-normal distribution has shown a good performance in two main use cases: (1) predicting the proportion of time traffic will exceed a given level (for service level agreement or link capacity estimation) i.e. link dimensioning based on bandwidth provisioning and (2) predicting 95th percentile pricing.[61]

## Notes

1. Johnson, Norman L.; Kotz, Samuel; Balakrishnan, N. (1994), "14: Lognormal Distributions", Continuous univariate distributions. Vol. 1, Wiley Series in Probability and Mathematical Statistics: Applied Probability and Statistics (2nd ed.), New York: John Wiley & Sons, ISBN 978-0-471-58495-7, MR 1299979
2. Park, Sung Y.; Bera, Anil K. (2009). "Maximum entropy autoregressive conditional heteroskedasticity model" (PDF). Journal of Econometrics. 150 (2): 219–230. CiteSeerX 10.1.1.511.9750. doi:10.1016/j.jeconom.2008.12.014. Retrieved 2011-06-02. Table 1, p. 221.
3. Tarmast, Ghasem (2001). Multivariate Log–Normal Distribution (PDF). ISI Proceedings: 53rd Session. Seoul.
4. Halliwell, Leigh (2015). The Lognormal Random Multivariate (PDF). Casualty Actuarial Society E-Forum, Spring 2015. Arlington, VA.
5. Heyde, CC. (1963), "On a property of the lognormal distribution", Journal of the Royal Statistical Society, Series B, 25 (2): 392–393, doi:10.1007/978-1-4419-5823-5_6, ISBN 978-1-4419-5822-8
6. Holgate, P. (1989). "The lognormal characteristic function, vol. 18, pp. 4539–4548, 1989". Communications in Statistics - Theory and Methods. 18 (12): 4539–4548. doi:10.1080/03610928908830173.
7. Barakat, R. (1976). "Sums of independent lognormally distributed random variables". Journal of the Optical Society of America. 66 (3): 211–216. doi:10.1364/JOSA.66.000211.
8. Barouch, E.; Kaufman, GM.; Glasser, ML. (1986). "On sums of lognormal random variables" (PDF). Studies in Applied Mathematics. 75 (1): 37–55. doi:10.1002/sapm198675137. hdl:1721.1/48703.
9. Leipnik, Roy B. (January 1991). "On Lognormal Random Variables: I – The Characteristic Function" (PDF). Journal of the Australian Mathematical Society Series B. 32 (3): 327–347. doi:10.1017/S0334270000006901.
10. 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).
11. Kirkwood, Thomas BL (Dec 1979). "Geometric means and measures of dispersion". Biometrics. 35 (4): 908–9. JSTOR 2530139.
12. Limpert, E; Stahel, W; Abbt, M (2001). "Lognormal distributions across the sciences: keys and clues". BioScience. 51 (5): 341–352. doi:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2.
13. Heil P, Friedrich B (2017). "Onset-Duration Matching of Acoustic Stimuli Revisited: Conventional Arithmetic vs. Proposed Geometric Measures of Accuracy and Precision". Frontiers in Psychology. 7: 2013. doi:10.3389/fpsyg.2016.02013. PMC 5216879. PMID 28111557.
14. Daly, Leslie E.; Bourke, Geoffrey Joseph (2000). Interpretation and uses of medical statistics. Journal of Epidemiology and Community Health. 46 (5th ed.). Wiley-Blackwell. p. 89. doi:10.1002/9780470696750. ISBN 978-0-632-04763-5. PMC 1059583.
15. "ProbOnto". Retrieved 1 July 2017.
16. Swat, MJ; Grenon, P; Wimalaratne, S (2016). "ProbOnto: ontology and knowledge base of probability distributions". Bioinformatics. 32 (17): 2719–21. doi:10.1093/bioinformatics/btw170. PMC 5013898. PMID 27153608.
17. Forbes et al. Probability Distributions (2011), John Wiley & Sons, Inc.
18. Lunn, D. (2012). The BUGS book: a practical introduction to Bayesian analysis. Texts in statistical science. CRC Press.
19. Limpert, E., Stahel, W. A., and Abbt, M. (2001). Log-normal distributions across the sciences: Keys and clues. BioScience, 51(5):341-352.
20. Nyberg J. et al. (2012) PopED - An extended, parallelized, population optimal design tool. Comput Methods Programs Biomed.; 108(2):789-805. doi: 10.1016/j.cmpb.2012.05.005
21. 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:141–151
22. The PopED Development Team (2016). PopED Manual, Release version 2.13. Technical report, Uppsala University.
23. ProbOnto website, URL: http://probonto.org
24. Damgaard, Christian; Weiner, Jacob (2000). "Describing inequality in plant size or fecundity". Ecology. 81 (4): 1139–1142. doi:10.1890/0012-9658(2000)081[1139:DIIPSO]2.0.CO;2.
25. Rossman, Lewis A (July 1990). "Design stream flows based on harmonic means". Journal of Hydraulic Engineering. 116 (7): 946–950. doi:10.1061/(ASCE)0733-9429(1990)116:7(946).
26. Thorin, Olof (1977). "On the infinite divisibility of the lognormal distribution". Scandinavian Actuarial Journal. 1977 (3): 121–148. doi:10.1080/03461238.1977.10405635. ISSN 0346-1238.
27. Gao, Xin (2009). "Asymptotic Behavior of Tail Density for Sum of Correlated Lognormal Variables". International Journal of Mathematics and Mathematical Sciences. 2009: 1–28. doi:10.1155/2009/630857.
28. Asmussen, S.; Rojas-Nandayapa, L. (2008). "Asymptotics of Sums of Lognormal Random Variables with Gaussian Copula" (PDF). Statistics and Probability Letters. 78 (16): 2709–2714. doi:10.1016/j.spl.2008.03.035.
29. Marlow, NA. (Nov 1967). "A normal limit theorem for power sums of independent normal random variables". Bell System Technical Journal. 46 (9): 2081–2089. doi:10.1002/j.1538-7305.1967.tb04244.x.
30. Botev, Z. I.; L'Ecuyer, P. (2017). "Accurate computation of the right tail of the sum of dependent log-normal variates". 2017 Winter Simulation Conference (WSC). 3rd–6th Dec 2017 Las Vegas, NV, USA: IEEE. pp. 1880–1890. arXiv:1705.03196. doi:10.1109/WSC.2017.8247924. ISBN 978-1-5386-3428-8.
31. Asmussen, A.; Goffard, P.-O.; Laub, P. J. (2016). "Orthonormal polynomial expansions and lognormal sum densities". arXiv:1601.01763v1 [math.PR].
32. Sangal, B.; Biswas, A. (1970). "The 3-Parameter Lognormal Distribution Applications in Hydrology". Water Resources Research. 6 (2): 505–515. doi:10.1029/WR006i002p00505.
33. Swamee, P. K. (2002). "Near Lognormal Distribution". Journal of Hydrologic Engineering. 7 (6): 441–444. doi:10.1061/(ASCE)1084-0699(2002)7:6(441).
34. Wu, Ziniu; Li, Juan; Bai, Chenyuan (2017). "Scaling Relations of Lognormal Type Growth Process with an Extremal Principle of Entropy". Entropy. 19 (56): 1–14. Bibcode:2017Entrp..19...56W. doi:10.3390/e19020056.
35. Wu, Zi-Niu (2003). "Prediction of the size distribution of secondary ejected droplets by crown splashing of droplets impinging on a solid wall". Probabilistic Engineering Mechanics. 18 (3): 241–249. doi:10.1016/S0266-8920(03)00028-6.
36. Wang, WenBin; Wu, ZiNiu; Wang, ChunFeng; Hu, RuiFeng (2013). "Modelling the spreading rate of controlled communicable epidemics through an entropy-based thermodynamic model". Science China Physics, Mechanics and Astronomy. 56 (11): 2143–2150. arXiv:1304.5603. Bibcode:2013SCPMA..56.2143W. doi:10.1007/s11433-013-5321-0. ISSN 1674-7348.
37. Bloetscher, Frederick (2019). "Using predictive Bayesian Monte Carlo- Markov Chain methods to provide a probabilistic solution for the Drake equation". Acta Astronautica. 155: 118–130. doi:10.1016/j.actaastro.2018.11.033.
38. Sutton, John (Mar 1997). "Gibrat's Legacy". Journal of Economic Literature. 32 (1): 40–59. JSTOR 2729692.
39. Pawel, Sobkowicz; et al. (2013). "Lognormal distributions of user post lengths in Internet discussions - a consequence of the Weber-Fechner law?". EPJ Data Science.
40. Yin, Peifeng; Luo, Ping; Lee, Wang-Chien; Wang, Min (2013). Silence is also evidence: interpreting dwell time for recommendation from psychological perspective. ACM International Conference on KDD.
41. "What is the average length of a game of chess?". chess.stackexchange.com. Retrieved 14 April 2018.
42. "Rubik's Cube Competitors' Mean times from 2019 competitions". reddit.com. 2019-08-21. Retrieved 2018-08-23.
43. Huxley, Julian S. (1932). Problems of relative growth. London. ISBN 978-0-486-61114-3. OCLC 476909537.
44. Makuch, Robert W.; D.H. Freeman; M.F. Johnson (1979). "Justification for the lognormal distribution as a model for blood pressure". Journal of Chronic Diseases. 32 (3): 245–250. doi:10.1016/0021-9681(79)90070-5. PMID 429469.
45. Scheler, Gabriele; Schumann, Johann (2006-10-08). Diversity and stability in neuronal output rates. 36th Society for Neuroscience Meeting, Atlanta.
46. 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. doi:10.1016/j.celrep.2013.07.039. ISSN 2211-1247. PMC 3804159. PMID 23994479.
47. 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. doi:10.1038/nrn3687. ISSN 1471-003X. PMC 4051294. PMID 24569488.
48. 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. doi:10.1016/j.pneurobio.2012.09.004. ISSN 1873-5118. PMC 5985929. PMID 23123501.
49. Scheler, Gabriele (2017-07-28). "Logarithmic distributions prove that intrinsic learning is Hebbian". F1000Research. 6: 1222. doi:10.12688/f1000research.12130.2. PMC 5639933. PMID 29071065.
50. Oosterbaan, R.J. (1994). "6: Frequency and Regression Analysis" (PDF). In Ritzema, H.P. (ed.). Drainage Principles and Applications, Publication 16. Wageningen, The Netherlands: International Institute for Land Reclamation and Improvement (ILRI). pp. 175–224. ISBN 978-90-70754-33-4.
51. CumFreq, free software for distribution fitting
52. Clementi, Fabio; Gallegati, Mauro (2005) "Pareto's law of income distribution: Evidence for Germany, the United Kingdom, and the United States", EconWPA
53. Wataru, Souma (2002-02-22). "Physics of Personal Income". In Takayasu, Hideki (ed.). Empirical Science of Financial Fluctuations: The Advent of Econophysics. Springer. arXiv:cond-mat/0202388. doi:10.1007/978-4-431-66993-7.
54. Black, F.; Scholes, M. (1973). "The Pricing of Options and Corporate Liabilities". Journal of Political Economy. 81 (3): 637. doi:10.1086/260062.
55. Mandelbrot, Benoit (2004). The (mis-)Behaviour of Markets. Basic Books. ISBN 9780465043552.
56. Bunchen, P., Advanced Option Pricing, University of Sydney coursebook, 2007
57. Thelwall, Mike; Wilson, Paul (2014). "Regression for citation data: An evaluation of different methods". Journal of Infometrics. 8 (4): 963–971. arXiv:1510.08877. doi:10.1016/j.joi.2014.09.011.
58. O'Connor, Patrick; Kleyner, Andre (2011). Practical Reliability Engineering. John Wiley & Sons. p. 35. ISBN 978-0-470-97982-2.
59. "Shadowing". www.WirelessCommunication.NL. Archived from the original on January 13, 2012.
60. Gros, C; Kaczor, G.; Markovic, D (2012). "Neuropsychological constraints to human data production on a global scale". The European Physical Journal B. 85 (28): 28. arXiv:1111.6849. Bibcode:2012EPJB...85...28G. doi:10.1140/epjb/e2011-20581-3.
61. Alamsar, Mohammed; Parisis, George; Clegg, Richard; Zakhleniuk, Nickolay (2019). "On the Distribution of Traffic Volumes in the Internet and its Implications". arXiv:1902.03853 [cs.NI].
• Crow, Edwin L.; Shimizu, Kunio, eds. (1988), Lognormal Distributions, Theory and Applications, Statistics: Textbooks and Monographs, 88, New York: Marcel Dekker, Inc., pp. xvi+387, ISBN 978-0-8247-7803-3, MR 0939191, Zbl 0644.62014
• Aitchison, J. and Brown, J.A.C. (1957) The Lognormal Distribution, Cambridge University Press.
• Limpert, E; Stahel, W; Abbt, M (2001). "Lognormal distributions across the sciences: keys and clues". BioScience. 51 (5): 341–352. doi:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2.
• Eric W. Weisstein et al. Log Normal Distribution at MathWorld. Electronic document, retrieved October 26, 2006.
• Holgate, P. (1989). "The lognormal characteristic function". Communications in Statistics - Theory and Methods. 18 (12): 4539–4548. doi:10.1080/03610928908830173.
• Brooks, Robert; Corson, Jon; Donal, Wales (1994). "The Pricing of Index Options When the Underlying Assets All Follow a Lognormal Diffusion". Advances in Futures and Options Research. 7. SSRN 5735.