Menu
Home Explore People Places Arts History Plants & Animals Science Life & Culture Technology
On this page
Log-normal distribution
Probability distribution

In probability theory, a log-normal distribution is a continuous probability distribution of a random variable whose logarithm is normally distributed. If Y = ln(X) has a normal distribution, then X = exp(Y) follows a log-normal distribution, which takes only positive values. This model is widely used in engineering, medicine, and economics for variables like energies, concentrations, and prices. It is also known as the Galton distribution after Francis Galton and linked to names such as McAlister and Gibrat. The log-normal distribution emerges from the multiplicative product of independent positive variables and relates to the central limit theorem in the log domain.

Related Image Collections Add Image
We don't have any YouTube videos related to Log-normal distribution yet.
We don't have any PDF documents related to Log-normal distribution yet.
We don't have any Books related to Log-normal distribution yet.

Definitions

Generation and parameters

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

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

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

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

In order to produce a distribution with desired mean μ X {\displaystyle \mu _{X}} and variance   σ X 2   , {\displaystyle \ \sigma _{X}^{2}\ ,} one uses   μ = ln ⁡ ( μ X 2   μ X 2 + σ X 2     )   {\displaystyle \ \mu =\ln \left({\frac {\mu _{X}^{2}}{\ {\sqrt {\mu _{X}^{2}+\sigma _{X}^{2}\ }}\ }}\right)\ } and   σ 2 = ln ⁡ ( 1 +   σ X 2   μ X 2 )   . {\displaystyle \ \sigma ^{2}=\ln \left(1+{\frac {\ \sigma _{X}^{2}\ }{\mu _{X}^{2}}}\right)~.}

Alternatively, the "multiplicative" or "geometric" parameters   μ ∗ = e μ   {\displaystyle \ \mu ^{*}=e^{\mu }\ } and   σ ∗ = e σ   {\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   {\displaystyle \ X\ } is log-normally distributed (i.e.,   X ∼ Lognormal ⁡ (   μ , σ 2   )   {\displaystyle \ X\sim \operatorname {Lognormal} \left(\ \mu ,\sigma ^{2}\ \right)\ } ), if the natural logarithm of   X   {\displaystyle \ X\ } is normally distributed with mean μ {\displaystyle \mu } and variance   σ 2   : {\displaystyle \ \sigma ^{2}\ :}

ln ⁡ ( X ) ∼ N ( μ , σ 2 ) {\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   )   {\displaystyle \ {\mathcal {N}}(\ 0,1\ )\ } standard normal distribution, then we have that67 the probability density function of the log-normal distribution is given by:

f X ( x ) = d d x   P X [   X ≤ x   ] = d d x   P X [   ln ⁡ X ≤ ln ⁡ x   ] = d d x Φ (   ln ⁡ x − μ   σ ) = φ ( ln ⁡ x − μ σ ) d d x (   ln ⁡ x − μ   σ ) = φ (   ln ⁡ x − μ   σ ) 1   σ   x   = 1   x   σ 2   π     exp ⁡ ( −   ( ln ⁡ x − μ ) 2   2   σ 2 )   . {\displaystyle {\begin{aligned}f_{X}(x)&={\frac {\rm {d}}{{\rm {d}}x}}\ \operatorname {\mathbb {P} _{\mathit {X}}} \,\!{\bigl [}\ X\leq x\ {\bigr ]}\\[6pt]&={\frac {\rm {d}}{{\rm {d}}x}}\ \operatorname {\mathbb {P} _{\mathit {X}}} \,\!{\bigl [}\ \ln X\leq \ln x\ {\bigr ]}\\[6pt]&={\frac {\rm {d}}{{\rm {d}}x}}\operatorname {\Phi } \!\!\left({\frac {\ \ln x-\mu \ }{\sigma }}\right)\\[6pt]&=\operatorname {\varphi } \!\left({\frac {\ln x-\mu }{\sigma }}\right){\frac {\rm {d}}{{\rm {d}}x}}\left({\frac {\ \ln x-\mu \ }{\sigma }}\right)\\[6pt]&=\operatorname {\varphi } \!\left({\frac {\ \ln x-\mu \ }{\sigma }}\right){\frac {1}{\ \sigma \ x\ }}\\[6pt]&={\frac {1}{\ x\ \sigma {\sqrt {2\ \pi \ }}\ }}\exp \left(-{\frac {\ (\ln x-\mu )^{2}\ }{2\ \sigma ^{2}}}\right)~.\end{aligned}}}

Cumulative distribution function

The cumulative distribution function is

F X ( x ) = Φ ( ( ln ⁡ x ) − μ σ ) {\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 )   {\displaystyle \ \operatorname {\mathcal {N}} (\ 0,\ 1)\ } ).

This may also be expressed as follows:8

1 2 [ 1 + erf ⁡ ( ln ⁡ x − μ σ 2 ) ] = 1 2 erfc ⁡ ( − ln ⁡ x − μ σ 2 ) {\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 X ∼ N ( μ , Σ ) {\displaystyle {\boldsymbol {X}}\sim {\mathcal {N}}({\boldsymbol {\mu }},\,{\boldsymbol {\Sigma }})} is a multivariate normal distribution, then Y i = exp ⁡ ( X i ) {\displaystyle Y_{i}=\exp(X_{i})} has a multivariate log-normal distribution.910 The exponential is applied elementwise to the random vector X {\displaystyle {\boldsymbol {X}}} . The mean of Y {\displaystyle {\boldsymbol {Y}}} is

E ⁡ [ Y ] i = e μ i + 1 2 Σ i i , {\displaystyle \operatorname {E} [{\boldsymbol {Y}}]_{i}=e^{\mu _{i}+{\frac {1}{2}}\Sigma _{ii}},}

and its covariance matrix is

Var ⁡ [ Y ] i j = e μ i + μ j + 1 2 ( Σ i i + Σ j j ) ( e Σ i j − 1 ) . {\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

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

This can be derived by letting z = ln ⁡ ( x ) − ( μ + n σ 2 ) σ {\displaystyle z={\tfrac {\ln(x)-(\mu +n\sigma ^{2})}{\sigma }}} within the integral. However, the log-normal distribution is not determined by its moments.11 This implies that it cannot have a defined moment generating function in a neighborhood of zero.12 Indeed, the expected value E ⁡ [ e t X ] {\displaystyle \operatorname {E} [e^{tX}]} is not defined for any positive value of the argument t {\displaystyle t} , since the defining integral diverges.

The characteristic function E ⁡ [ e i t X ] {\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 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.13 In particular, its Taylor formal series diverges:

∑ n = 0 ∞ ( i t ) n n ! e n μ + n 2 σ 2 / 2 {\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.14151617

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

φ ( t ) ≈ exp ⁡ ( − W 2 ( − i t σ 2 e μ ) + 2 W ( − i t σ 2 e μ ) 2 σ 2 ) 1 + W ( − i t σ 2 e μ ) {\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 W {\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

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.19 (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.20 (Matlab code)

Geometric or multiplicative moments

The geometric or multiplicative mean of the log-normal distribution is GM ⁡ [ X ] = e μ = μ ∗ {\displaystyle \operatorname {GM} [X]=e^{\mu }=\mu ^{*}} . It equals the median. The geometric or multiplicative standard deviation is GSD ⁡ [ X ] = e σ = σ ∗ {\displaystyle \operatorname {GSD} [X]=e^{\sigma }=\sigma ^{*}} .2122

By analogy with the arithmetic statistics, one can define a geometric variance, GVar ⁡ [ X ] = e σ 2 {\displaystyle \operatorname {GVar} [X]=e^{\sigma ^{2}}} , and a geometric coefficient of variation,23 GCV ⁡ [ X ] = e σ − 1 {\displaystyle \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 CV {\displaystyle \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,

E ⁡ [ X ] = e μ + 1 2 σ 2 = e μ ⋅ e σ 2 = GM ⁡ [ X ] ⋅ GVar ⁡ [ X ] . {\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]}}.} 24

In finance, the term e − 1 2 σ 2 {\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 by25

E ⁡ [ X n ] = e n μ + 1 2 n 2 σ 2 . {\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 respectively given by:26

E ⁡ [ X ] = e μ + 1 2 σ 2 , E ⁡ [ X 2 ] = e 2 μ + 2 σ 2 , Var ⁡ [ X ] = E ⁡ [ X 2 ] − E ⁡ [ X ] 2 = ( E ⁡ [ X ] ) 2 ( e σ 2 − 1 ) = e 2 μ + σ 2 ( e σ 2 − 1 ) , SD ⁡ [ X ] = Var ⁡ [ X ] = E ⁡ [ X ] e σ 2 − 1 = e μ + 1 2 σ 2 e σ 2 − 1 , {\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}}}

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

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

This estimate is sometimes referred to as the "geometric CV" (GCV),2829 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:

μ = ln ⁡ ( E ⁡ [ X ] 2 E ⁡ [ X 2 ] ) = ln ⁡ ( E ⁡ [ X ] 2 Var ⁡ [ X ] + E ⁡ [ X ] 2 ) , σ 2 = ln ⁡ ( E ⁡ [ X 2 ] E ⁡ [ X ] 2 ) = ln ⁡ ( 1 + Var ⁡ [ X ] E ⁡ [ X ] 2 ) . {\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/2⁠n2σ2 for n ≥ 1. That is, there exist other distributions with the same set of moments.30 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 ( ln ⁡ f ) ′ = 0 {\displaystyle (\ln f)'=0} , we get that:

Mode ⁡ [ X ] = e μ − σ 2 . {\displaystyle \operatorname {Mode} [X]=e^{\mu -\sigma ^{2}}.}

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

q X ( α ) = e μ + σ q Φ ( α ) = μ ∗ ( σ ∗ ) q Φ ( α ) , {\displaystyle q_{X}(\alpha )=e^{\mu +\sigma q_{\Phi }(\alpha )}=\mu ^{*}(\sigma ^{*})^{q_{\Phi }(\alpha )},}

where q Φ ( α ) {\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,31

Med ⁡ [ X ] = e μ = μ ∗   . {\displaystyle \operatorname {Med} [X]=e^{\mu }=\mu ^{*}~.}

Partial expectation

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

g ( k ) = ∫ k ∞ x f X ( x ∣ X > k ) d x . {\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x\mid X>k)\,dx.}

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

g ( k ) = ∫ k ∞ x f X ( x ∣ X > k ) d x = e μ + 1 2 σ 2 Φ ( μ + σ 2 − ln ⁡ k σ ) {\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x\mid X>k)\,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 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 {\displaystyle X} —with respect to a threshold k {\displaystyle k} —is its partial expectation divided by the cumulative probability of being in that range:

E [ X ∣ X < k ] = e μ + σ 2 2 ⋅ Φ [ ln ⁡ ( k ) − μ − σ 2 σ ] Φ [ ln ⁡ ( k ) − μ σ ] E [ X ∣ X ⩾ k ] = e μ + σ 2 2 ⋅ Φ [ μ + σ 2 − ln ⁡ ( k ) σ ] 1 − Φ [ ln ⁡ ( k ) − μ σ ] E [ X ∣ X ∈ [ k 1 , k 2 ] ] = e μ + σ 2 2 ⋅ Φ [ ln ⁡ ( k 2 ) − μ − σ 2 σ ] − Φ [ ln ⁡ ( k 1 ) − μ − σ 2 σ ] Φ [ ln ⁡ ( k 2 ) − μ σ ] − Φ [ ln ⁡ ( k 1 ) − μ σ ] {\displaystyle {\begin{aligned}E[X\mid X<k]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\ln(k)-\mu -\sigma ^{2}}{\sigma }}\right]}{\Phi \left[{\frac {\ln(k)-\mu }{\sigma }}\right]}}\\[8pt]E[X\mid X\geqslant k]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\mu +\sigma ^{2}-\ln(k)}{\sigma }}\right]}{1-\Phi \left[{\frac {\ln(k)-\mu }{\sigma }}\right]}}\\[8pt]E[X\mid X\in [k_{1},k_{2}]]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\ln(k_{2})-\mu -\sigma ^{2}}{\sigma }}\right]-\Phi \left[{\frac {\ln(k_{1})-\mu -\sigma ^{2}}{\sigma }}\right]}{\Phi \left[{\frac {\ln(k_{2})-\mu }{\sigma }}\right]-\Phi \left[{\frac {\ln(k_{1})-\mu }{\sigma }}\right]}}\end{aligned}}}

Alternative parameterizations

In addition to the characterization by μ , σ {\displaystyle \mu ,\sigma } or μ ∗ , σ ∗ {\displaystyle \mu ^{*},\sigma ^{*}} , here are multiple ways how the log-normal distribution can be parameterized. ProbOnto, the knowledge base and ontology of probability distributions3233 lists seven such forms:

  • LogNormal1(μ,σ) with mean, μ, and standard deviation, σ, both on the log-scale 34 P ( x ; μ , σ ) = 1 x σ 2 π exp ⁡ [ − ( ln ⁡ x − μ ) 2 2 σ 2 ] {\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\sigma }})={\frac {1}{x\sigma {\sqrt {2\pi }}}}\exp \left[-{\frac {(\ln x-\mu )^{2}}{2\sigma ^{2}}}\right]}
  • LogNormal2(μ,υ) with mean, μ, and variance, υ, both on the log-scale P ( x ; μ , v ) = 1 x v 2 π exp ⁡ [ − ( ln ⁡ x − μ ) 2 2 v ] {\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {v}})={\frac {1}{x{\sqrt {v}}{\sqrt {2\pi }}}}\exp \left[-{\frac {(\ln x-\mu )^{2}}{2v}}\right]}
  • LogNormal3(m,σ) with median, m, on the natural scale and standard deviation, σ, on the log-scale35 P ( x ; m , σ ) = 1 x σ 2 π exp ⁡ [ − ln 2 ⁡ ( x / m ) 2 σ 2 ] {\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma }})={\frac {1}{x\sigma {\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\sigma ^{2}}}\right]}
  • LogNormal4(m,cv) with median, m, and coefficient of variation, cv, both on the natural scale P ( x ; m , c v ) = 1 x ln ⁡ ( c v 2 + 1 ) 2 π exp ⁡ [ − ln 2 ⁡ ( x / m ) 2 ln ⁡ ( c v 2 + 1 ) ] {\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {cv}})={\frac {1}{x{\sqrt {\ln(cv^{2}+1)}}{\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\ln(cv^{2}+1)}}\right]}
  • LogNormal5(μ,τ) with mean, μ, and precision, τ, both on the log-scale36 P ( x ; μ , τ ) = τ 2 π 1 x exp ⁡ [ − τ 2 ( ln ⁡ x − μ ) 2 ] {\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\tau }})={\sqrt {\frac {\tau }{2\pi }}}{\frac {1}{x}}\exp \left[-{\frac {\tau }{2}}(\ln x-\mu )^{2}\right]}
  • LogNormal6(m,σg) with median, m, and geometric standard deviation, σg, both on the natural scale37 P ( x ; m , σ g ) = 1 x ln ⁡ ( σ g ) 2 π exp ⁡ [ − ln 2 ⁡ ( x / m ) 2 ln 2 ⁡ ( σ g ) ] {\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma _{g}}})={\frac {1}{x\ln(\sigma _{g}){\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\ln ^{2}(\sigma _{g})}}\right]}
  • LogNormal7(μN,σN) with mean, μN, and standard deviation, σN, both on the natural scale38 P ( x ; μ N , σ N ) = 1 x 2 π ln ⁡ ( 1 + σ N 2 / μ N 2 ) exp ⁡ ( − [ ln ⁡ x − ln ⁡ μ N 1 + σ N 2 / μ N 2 ] 2 2 ln ⁡ ( 1 + σ N 2 / μ N 2 ) ) {\displaystyle P(x;{\boldsymbol {\mu _{N}}},{\boldsymbol {\sigma _{N}}})={\frac {1}{x{\sqrt {2\pi \ln \left(1+\sigma _{N}^{2}/\mu _{N}^{2}\right)}}}}\exp \left(-{\frac {{\Big [}\ln x-\ln {\frac {\mu _{N}}{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}}{\Big ]}^{2}}{2\ln(1+\sigma _{N}^{2}/\mu _{N}^{2})}}\right)}

Examples for re-parameterization

Consider the situation when one would like to run a model using two different optimal design tools, for example PFIM39 and PopED.40 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 LN2 ⁡ ( μ , v ) → LN7 ⁡ ( μ N , σ N ) {\displaystyle \operatorname {LN2} (\mu ,v)\to \operatorname {LN7} (\mu _{N},\sigma _{N})} following formulas hold μ N = exp ⁡ ( μ + v / 2 ) {\textstyle \mu _{N}=\exp(\mu +v/2)} and σ N = exp ⁡ ( μ + v / 2 ) exp ⁡ ( v ) − 1 {\textstyle \sigma _{N}=\exp(\mu +v/2){\sqrt {\exp(v)-1}}} .

For the transition LN7 ⁡ ( μ N , σ N ) → LN2 ⁡ ( μ , v ) {\displaystyle \operatorname {LN7} (\mu _{N},\sigma _{N})\to \operatorname {LN2} (\mu ,v)} following formulas hold μ = ln ⁡ ( μ N / 1 + σ N 2 / μ N 2 ) {\textstyle \mu =\ln \left(\mu _{N}/{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}\right)} and v = ln ⁡ ( 1 + σ N 2 / μ N 2 ) {\textstyle 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.41

Multiple, reciprocal, power

  • Multiplication by a constant: If X ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} then a X ∼ Lognormal ⁡ ( μ + ln ⁡ a ,   σ 2 ) {\displaystyle aX\sim \operatorname {Lognormal} (\mu +\ln a,\ \sigma ^{2})} for a > 0. {\displaystyle a>0.}
  • Reciprocal: If X ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} then 1 X ∼ Lognormal ⁡ ( − μ ,   σ 2 ) . {\displaystyle {\tfrac {1}{X}}\sim \operatorname {Lognormal} (-\mu ,\ \sigma ^{2}).}
  • Power: If X ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} then X a ∼ Lognormal ⁡ ( a μ ,   a 2 σ 2 ) {\displaystyle X^{a}\sim \operatorname {Lognormal} (a\mu ,\ a^{2}\sigma ^{2})} for a ≠ 0. {\displaystyle a\neq 0.}

Multiplication and division of independent, log-normal random variables

If two independent, log-normal variables X 1 {\displaystyle X_{1}} and X 2 {\displaystyle X_{2}} are multiplied [divided], the product [ratio] is again log-normal, with parameters μ = μ 1 + μ 2 {\displaystyle \mu =\mu _{1}+\mu _{2}} [ μ = μ 1 − μ 2 {\displaystyle \mu =\mu _{1}-\mu _{2}} ] and σ {\displaystyle \sigma } , where σ 2 = σ 1 2 + σ 2 2 {\displaystyle \sigma ^{2}=\sigma _{1}^{2}+\sigma _{2}^{2}} .

More generally, if X j ∼ Lognormal ⁡ ( μ j , σ j 2 ) {\displaystyle X_{j}\sim \operatorname {Lognormal} (\mu _{j},\sigma _{j}^{2})} are n {\displaystyle n} independent, log-normally distributed variables, then Y = ∏ j = 1 n X j ∼ Lognormal ⁡ ( ∑ j = 1 n μ j ,   ∑ j = 1 n σ j 2 ) . {\displaystyle Y=\textstyle \prod _{j=1}^{n}X_{j}\sim \operatorname {Lognormal} {\Big (}\textstyle \sum _{j=1}^{n}\mu _{j},\ \sum _{j=1}^{n}\sigma _{j}^{2}{\Big )}.}

Multiplicative central limit theorem

See also: Gibrat's law

The geometric or multiplicative mean of n {\displaystyle n} independent, identically distributed, positive random variables X i {\displaystyle X_{i}} shows, for n → ∞ {\displaystyle n\to \infty } , approximately a log-normal distribution with parameters μ = E [ ln ⁡ ( X i ) ] {\displaystyle \mu =E[\ln(X_{i})]} and σ 2 = var [ ln ⁡ ( X i ) ] / n {\displaystyle \sigma ^{2}={\mbox{var}}[\ln(X_{i})]/n} , assuming σ 2 {\displaystyle \sigma ^{2}} is finite.

In fact, the random variables do not have to be identically distributed. It is enough for the distributions of ln ⁡ ( X i ) {\displaystyle \ln(X_{i})} 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.

Heavy-tailness of the Log-Normal

Whether a Log-Normal can be considered or not a true heavy-tail distribution is still debated. The main reason is that its variance is always finite, differently from what happen with certain Pareto distributions, for instance. However a recent study has shown how it is possible to create a Log-Normal distribution with infinite variance using Robinson Non-Standard Analysis.42

Other

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

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

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

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

  • If X ∼ N ( μ , σ 2 ) {\displaystyle X\sim {\mathcal {N}}(\mu ,\sigma ^{2})} is a normal distribution, then exp ⁡ ( X ) ∼ Lognormal ⁡ ( μ , σ 2 ) . {\displaystyle \exp(X)\sim \operatorname {Lognormal} (\mu ,\sigma ^{2}).}
  • If X ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} is distributed log-normally, then ln ⁡ ( X ) ∼ N ( μ , σ 2 ) {\displaystyle \ln(X)\sim {\mathcal {N}}(\mu ,\sigma ^{2})} is a normal random variable.
  • Let X j ∼ Lognormal ⁡ ( μ j , σ j 2 ) {\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 Y = ∑ j = 1 n X j {\textstyle Y=\sum _{j=1}^{n}X_{j}} . The distribution of Y {\displaystyle Y} has no closed-form expression, but can be reasonably approximated by another log-normal distribution Z {\displaystyle Z} at the right tail.47 Its probability density function at the neighborhood of 0 has been characterized48 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 Marlow49) is obtained by matching the mean and variance of another log-normal distribution: σ Z 2 = ln [ ∑ e 2 μ j + σ j 2 ( e σ j 2 − 1 ) ( ∑ e μ j + σ j 2 / 2 ) 2 + 1 ] , μ Z = ln [ ∑ e μ j + σ j 2 / 2 ] − σ Z 2 2 . {\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 X j {\displaystyle X_{j}} have the same variance parameter σ j = σ {\displaystyle \sigma _{j}=\sigma } , these formulas simplify to σ Z 2 = ln [ ( e σ 2 − 1 ) ∑ e 2 μ j ( ∑ e μ j ) 2 + 1 ] , μ Z = ln [ ∑ e μ j ] + σ 2 2 − σ Z 2 2 . {\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 the right tail.5051

The sum of correlated log-normally distributed random variables can also be approximated by a log-normal distribution S + = E ⁡ [ ∑ i X i ] = ∑ i E ⁡ [ X i ] = ∑ i e μ i + σ i 2 / 2 σ Z 2 = 1 / S + 2 ∑ i , j cor i j ⁡ σ i σ j E ⁡ [ X i ] E ⁡ [ X j ] = 1 / S + 2 ∑ i , j cor i j ⁡ σ i σ j e μ i + σ i 2 / 2 e μ j + σ j 2 / 2 μ Z = ln ⁡ ( S + ) − σ Z 2 / 2 {\displaystyle {\begin{aligned}S_{+}&=\operatorname {E} \left[\sum _{i}X_{i}\right]=\sum _{i}\operatorname {E} [X_{i}]=\sum _{i}e^{\mu _{i}+\sigma _{i}^{2}/2}\\\sigma _{Z}^{2}&=1/S_{+}^{2}\,\sum _{i,j}\operatorname {cor} _{ij}\sigma _{i}\sigma _{j}\operatorname {E} [X_{i}]\operatorname {E} [X_{j}]=1/S_{+}^{2}\,\sum _{i,j}\operatorname {cor} _{ij}\sigma _{i}\sigma _{j}e^{\mu _{i}+\sigma _{i}^{2}/2}e^{\mu _{j}+\sigma _{j}^{2}/2}\\\mu _{Z}&=\ln \left(S_{+}\right)-\sigma _{Z}^{2}/2\end{aligned}}}

  • If X ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} then X + c {\displaystyle X+c} is said to have a Three-parameter log-normal distribution with support x ∈ ( c , + ∞ ) {\displaystyle x\in (c,+\infty )} .52 E ⁡ [ X + c ] = E ⁡ [ X ] + c {\displaystyle \operatorname {E} [X+c]=\operatorname {E} [X]+c} , Var ⁡ [ X + c ] = Var ⁡ [ X ] {\displaystyle \operatorname {Var} [X+c]=\operatorname {Var} [X]} .
  • The log-normal distribution is a special case of the semi-bounded Johnson's SU-distribution.53
  • If X ∣ Y ∼ Rayleigh ⁡ ( Y ) {\displaystyle X\mid Y\sim \operatorname {Rayleigh} (Y)} with Y ∼ Lognormal ⁡ ( μ , σ 2 ) {\displaystyle Y\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})} , then X ∼ Suzuki ⁡ ( μ , σ ) {\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 functions54 can be obtained based on the logistic distribution to get an approximation for the CDF F ( x ; μ , σ ) = [ ( e μ x ) π / ( σ 3 ) + 1 ] − 1 . {\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

Maximum likelihood estimator

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 L ( μ , σ ) = ∏ i = 1 n 1 x i φ μ , σ ( ln ⁡ x i ) , {\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 N ( μ , σ 2 ) {\displaystyle {\mathcal {N}}(\mu ,\sigma ^{2})} . Therefore, the log-likelihood function is ℓ ( μ , σ ∣ x 1 , x 2 , … , x n ) = − ∑ i ln ⁡ x i + ℓ N ( μ , σ ∣ ln ⁡ x 1 , ln ⁡ x 2 , … , ln ⁡ x n ) . {\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 ℓ N {\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 ln ⁡ x 1 , ln ⁡ x 2 , … , ln ⁡ x n ) {\displaystyle \ln x_{1},\ln x_{2},\dots ,\ln x_{n})} , μ ^ = ∑ i ln ⁡ x i n , σ ^ 2 = ∑ i ( ln ⁡ x i − μ ^ ) 2 n . {\displaystyle {\widehat {\mu }}={\frac {\sum _{i}\ln x_{i}}{n}},\qquad {\widehat {\sigma }}^{2}={\frac {\sum _{i}\left(\ln x_{i}-{\widehat {\mu }}\right)^{2}}{n}}.}

For finite n, the estimator for μ {\displaystyle \mu } is unbiased, but the one for σ {\displaystyle \sigma } is biased. As for the normal distribution, an unbiased estimator for σ {\displaystyle \sigma } can be obtained by replacing the denominator n by n−1 in the equation for σ ^ 2 {\displaystyle {\widehat {\sigma }}^{2}} .

From this, the MLE for the expectancy of x is:55 θ ^ M L E = E [ X ] ^ M L E = e μ ^ + σ ^ 2 2 {\displaystyle {\widehat {\theta }}_{MLE}={\widehat {E[X]}}_{MLE}=e^{{\hat {\mu }}+{\frac {{\hat {\sigma }}^{2}}{2}}}}

Method of moments

When the individual values x 1 , x 2 , … , x n {\displaystyle x_{1},x_{2},\ldots ,x_{n}} are not available, but the sample's mean x ¯ {\displaystyle {\bar {x}}} 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 E ⁡ [ X ] {\displaystyle \operatorname {E} [X]} and variance Var ⁡ [ X ] {\displaystyle \operatorname {Var} [X]} for μ {\displaystyle \mu } and σ {\displaystyle \sigma } :56 μ = ln ⁡ ( x ¯ 1 + σ ^ 2 / x ¯ 2 ) , σ 2 = ln ⁡ ( 1 + σ ^ 2 / x ¯ 2 ) . {\displaystyle \mu =\ln \left({\frac {\bar {x}}{\sqrt {1+{\widehat {\sigma }}^{2}/{\bar {x}}^{2}}}}\right),\qquad \sigma ^{2}=\ln \left(1+{{\widehat {\sigma }}^{2}}/{\bar {x}}^{2}\right).}

Other estimators

Other estimators also exist, such as Finney's UMVUE estimator,57 the "Approximately Minimum Mean Squared Error Estimator", the "Approximately Unbiased Estimator" and "Minimax Estimator",58 also "A Conditional Mean Squared Error Estimator",59 and other variations as well.6061

Interval estimates

Further information: Reference range § Log-normal distribution

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 [ μ − σ , μ + σ ] {\displaystyle [\mu -\sigma ,\mu +\sigma ]} contains approximately two thirds (68%) of the probability (or of a large sample), and [ μ − 2 σ , μ + 2 σ ] {\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 [ μ ∗ / ( σ ∗ ) 2 , μ ∗ ⋅ ( σ ∗ ) 2 ] = [ μ ∗ × / ( σ ∗ ) 2 ] {\displaystyle [\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 μ {\displaystyle \mu } is [ μ ^ ± q ⋅ s e ^ ] {\displaystyle [{\widehat {\mu }}\pm q\cdot {\widehat {\mathop {se} }}]} , where s e = σ ^ / n {\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 μ ∗ = e μ {\displaystyle \mu ^{*}=e^{\mu }} (the median), is: [ μ ^ ∗ × / ( sem ∗ ) q ] {\displaystyle [{\widehat {\mu }}^{*}{}^{\times }\!\!/(\operatorname {sem} ^{*})^{q}]} with sem ∗ = ( σ ^ ∗ ) 1 / n {\displaystyle \operatorname {sem} ^{*}=({\widehat {\sigma }}^{*})^{1/{\sqrt {n}}}}

Confidence interval for E(X)

The literature discusses several options for calculating the confidence interval for μ {\displaystyle \mu } (the mean of the log-normal distribution). These include bootstrap as well as various other methods.6263

The Cox Method64 proposes to plug-in the estimators μ ^ = ∑ i ln ⁡ x i n , S 2 = ∑ i ( ln ⁡ x i − μ ^ ) 2 n − 1 {\displaystyle {\widehat {\mu }}={\frac {\sum _{i}\ln x_{i}}{n}},\qquad S^{2}={\frac {\sum _{i}\left(\ln x_{i}-{\widehat {\mu }}\right)^{2}}{n-1}}}

and use them to construct approximate confidence intervals in the following way: C I ( E ( X ) ) : e ( μ ^ + S 2 2 ± z 1 − α 2 S 2 n + S 4 2 ( n − 1 ) ) {\displaystyle CI(E(X)):e^{\left({\hat {\mu }}+{\frac {S^{2}}{2}}\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S^{2}}{n}}+{\frac {S^{4}}{2(n-1)}}}}\right)}}

Olsson 2005, proposed a "modified Cox method" by replacing z 1 − α 2 {\displaystyle z_{1-{\frac {\alpha }{2}}}} with t n − 1 , 1 − α 2 {\displaystyle t_{n-1,1-{\frac {\alpha }{2}}}} , which seemed to provide better coverage results for small sample sizes.65: Section 3.4 

Confidence interval for comparing two log normals

Comparing two log-normal distributions can often be of interest, for example, from a treatment and control group (e.g., in an A/B test). We have samples from two independent log-normal distributions with parameters ( μ 1 , σ 1 2 ) {\displaystyle (\mu _{1},\sigma _{1}^{2})} and ( μ 2 , σ 2 2 ) {\displaystyle (\mu _{2},\sigma _{2}^{2})} , with sample sizes n 1 {\displaystyle n_{1}} and n 2 {\displaystyle n_{2}} respectively.

Comparing the medians of the two can easily be done by taking the log from each and then constructing straightforward confidence intervals and transforming it back to the exponential scale.

C I ( e μ 1 − μ 2 ) : e ( μ ^ 1 − μ ^ 2 ± z 1 − α 2 S 1 2 n + S 2 2 n ) {\displaystyle CI(e^{\mu _{1}-\mu _{2}}):e^{\left({\hat {\mu }}_{1}-{\hat {\mu }}_{2}\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n}}+{\frac {S_{2}^{2}}{n}}}}\right)}}

These CI are what's often used in epidemiology for calculation the CI for relative-risk and odds-ratio.66 The way it is done there is that we have two approximately Normal distributions (e.g., p1 and p2, for RR), and we wish to calculate their ratio.67

However, the ratio of the expectations (means) of the two samples might also be of interest, while requiring more work to develop. The ratio of their means is:

E ( X 1 ) E ( X 2 ) = e μ 1 + σ 1 2 2 e μ 2 + σ 2 2 2 = e ( μ 1 − μ 2 ) + 1 2 ( σ 1 2 − σ 2 2 ) {\displaystyle {\frac {E(X_{1})}{E(X_{2})}}={\frac {e^{\mu _{1}+{\frac {\sigma _{1}^{2}}{2}}}}{e^{\mu _{2}+{\frac {\sigma _{2}^{2}}{2}}}}}=e^{(\mu _{1}-\mu _{2})+{\frac {1}{2}}(\sigma _{1}^{2}-\sigma _{2}^{2})}}

Plugin in the estimators to each of these parameters yields also a log normal distribution, which means that the Cox Method, discussed above, could similarly be used for this use-case:

C I ( E ( X 1 ) E ( X 2 ) = e μ 1 + σ 1 2 2 e μ 2 + σ 2 2 2 ) : e ( ( μ ^ 1 − μ ^ 2 + 1 2 S 1 2 − 1 2 S 2 2 ) ± z 1 − α 2 S 1 2 n 1 + S 2 2 n 2 + S 1 4 2 ( n 1 − 1 ) + S 2 4 2 ( n 2 − 1 ) ) {\displaystyle CI\left({\frac {E(X_{1})}{E(X_{2})}}={\frac {e^{\mu _{1}+{\frac {\sigma _{1}^{2}}{2}}}}{e^{\mu _{2}+{\frac {\sigma _{2}^{2}}{2}}}}}\right):e^{\left(({\hat {\mu }}_{1}-{\hat {\mu }}_{2}+{\frac {1}{2}}S_{1}^{2}-{\frac {1}{2}}S_{2}^{2})\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}\right)}}

It's worth noting that naively using the MLE in the ratio of the two expectations to create a ratio estimator will lead to a consistent, yet biased, point-estimation (we use the fact that the estimator of the ratio is a log normal distribution)68:

E [ E ^ ( X 1 ) E ^ ( X 2 ) ] = E [ e ( μ ^ 1 − μ ^ 2 ) + 1 2 ( S 1 2 − S 2 2 ) ] ≈ e ( μ 1 − μ 2 ) + 1 2 ( σ 1 2 − σ 2 2 ) + 1 2 ( σ 1 2 n 1 + σ 2 2 n 2 + σ 1 4 2 ( n 1 − 1 ) + σ 2 4 2 ( n 2 − 1 ) ) {\displaystyle E\left[{\frac {{\widehat {E}}(X_{1})}{{\widehat {E}}(X_{2})}}\right]=E\left[e^{({\widehat {\mu }}_{1}-{\widehat {\mu }}_{2})+{\frac {1}{2}}(S_{1}^{2}-S_{2}^{2})}\right]\approx e^{(\mu _{1}-\mu _{2})+{\frac {1}{2}}(\sigma _{1}^{2}-\sigma _{2}^{2})+{\frac {1}{2}}\left({\frac {\sigma _{1}^{2}}{n_{1}}}+{\frac {\sigma _{2}^{2}}{n_{2}}}+{\frac {\sigma _{1}^{4}}{2(n_{1}-1)}}+{\frac {\sigma _{2}^{4}}{2(n_{2}-1)}}\right)}}

Extremal principle of entropy to fix the free parameter σ

In applications, σ {\displaystyle \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 that69 σ = 1 6 {\displaystyle \sigma ={\frac {1}{\sqrt {6}}}}

This value can then be used to give some scaling relation between the inflexion point and maximum point of the log-normal distribution.70 This relationship is determined by the base of natural logarithm, e = 2.718 … {\displaystyle 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 σ {\displaystyle \sigma } fits well with the size of secondarily produced droplets during droplet impact 71 and the spreading of an epidemic disease.72

The value σ = 1 / 6 {\textstyle \sigma =1{\big /}{\sqrt {6}}} is used to provide a probabilistic solution for the Drake equation.73

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.74 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.75 contains a review and table of log-normal distributions from geology, biology, medicine, food, ecology, and other areas.76 is a review article on log-normal distributions in neuroscience, with annotated bibliography.

Human behavior

  • The length of comments posted in Internet discussion forums follows a log-normal distribution.77
  • Users' dwell time on online articles (jokes, news etc.) follows a log-normal distribution.78
  • The length of chess games tends to follow a log-normal distribution.79
  • Onset durations of acoustic comparison stimuli that are matched to a standard stimulus follow a log-normal distribution.80

Biology and medicine

  • Measures of size of living tissue (length, skin area, weight).81
  • Incubation period of diseases.82
  • Diameters of banana leaf spots, powdery mildew on barley.83
  • For highly communicable epidemics, such as SARS in 2003, if public intervention control policies are involved, the number of hospitalized cases is shown to satisfy the log-normal 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.84
  • 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.
  • The PacBio sequencing read length follows a log-normal distribution.85
  • Certain physiological measurements, such as blood pressure of adult humans (after separation on male/female subpopulations).86
  • Several pharmacokinetic variables, such as Cmax, elimination half-life and the elimination rate constant.87
  • In neuroscience, the distribution of firing rates across a population of neurons is often approximately log-normal. This has been first observed in the cortex and striatum 88 and later in hippocampus and entorhinal cortex,89 and elsewhere in the brain.9091 Also, intrinsic gain distributions and synaptic weight distributions appear to be log-normal92 as well.
  • Neuron densities in the cerebral cortex, due to the noisy cell division process during neurodevelopment.93
  • In operating-rooms management, the distribution of surgery duration.
  • In the size of avalanches of fractures in the cytoskeleton of living cells, showing log-normal distributions, with significantly higher size in cancer cells than healthy ones.94

Chemistry

Physical sciences

  • 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.97
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.98 The rainfall data are represented by plotting positions as part of a cumulative frequency analysis.
  • In physical oceanography, the sizes of icebergs in the midwinter Southern Atlantic Ocean were found to follow a log-normal size distribution. The iceberg sizes, measured visually and by radar from the F.S. Polarstern in 1986, were thought to be controlled by wave action in heavy seas causing them to flex and break.99
  • In atmospheric science, log-normal distributions (or distributions made by combining multiple log-normal functions) have been used to characterize both measurements and models of the sizes and concentrations of many different types of particles, from volcanic ash, to clouds and rain, to airborne microbes.100101102103 The log-normal distribution is strictly empirical, so more physically-based distributions have been adopted to better understand processes controlling size distributions of particles such as volcanic ash.104

Social sciences and demographics

  • In economics, there is evidence that the income of 97–99% of the population is distributed log-normally.105 (The distribution of higher-income individuals follows a Pareto distribution).106
  • If an income distribution follows a log-normal distribution with standard deviation σ {\displaystyle \sigma } , then the Gini coefficient, commonly use to evaluate income inequality, can be computed as G = erf ⁡ ( σ 2 ) {\displaystyle G=\operatorname {erf} \left({\frac {\sigma }{2}}\right)} where erf {\displaystyle \operatorname {erf} } is the error function, since G = 2 Φ ( σ 2 ) − 1 {\displaystyle G=2\Phi \left({\frac {\sigma }{\sqrt {2}}}\right)-1} , where Φ ( x ) {\displaystyle \Phi (x)} is the cumulative distribution function of a standard normal distribution.
  • In finance, in particular the Black–Scholes model, changes in the logarithm of exchange rates, price indices, and stock market indices are assumed normal107 (these variables behave like compound interest, not like simple interest, and so are multiplicative). However, some mathematicians such as Benoit Mandelbrot have argued 108 that log-Lévy distributions, which possesses heavy tails would be a more appropriate model, in particular for the analysis for stock market crashes. Indeed, stock price distributions typically exhibit a fat tail.109 The fat tailed distribution of changes during stock market crashes invalidate the assumptions of the central limit theorem.
  • In scientometrics, the number of citations to journal articles and patents follows a discrete log-normal distribution.110111
  • City sizes (population) satisfy Gibrat's Law.112 The growth process of city sizes is proportionate and invariant with respect to size. From the central limit theorem therefore, the log of city size is normally distributed.
  • The number of sexual partners appears to be best described by a log-normal distribution.113

Technology

  • In reliability analysis, the log-normal distribution is often used to model times to repair a maintainable system.114
  • In wireless communication, "the local-mean power expressed in logarithmic values, such as dB or neper, has a normal (i.e., Gaussian) distribution."115 Also, the random obstruction of radio signals due to large buildings and hills, called shadowing, is often modeled as a log-normal distribution.
  • Particle size distributions produced by comminution with random impacts, such as in ball milling.116
  • The file size distribution of publicly available audio and video data files (MIME types) follows a log-normal distribution over five orders of magnitude.117
  • File sizes of 140 million files on personal computers running the Windows OS, collected in 1999.118119
  • Sizes of text-based emails (1990s) and multimedia-based emails (2000s).120
  • In computer networks and Internet traffic analysis, log-normal 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 Internet 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.121
  • in physical testing when the test produces a time-to-failure of an item under specified conditions, the data is often best analyzed using a lognormal distribution.122123

See also

Notes

Further reading

Wikimedia Commons has media related to Log-normal distribution.

References

  1. Weisstein, Eric W. "Log Normal Distribution". mathworld.wolfram.com. Retrieved 2020-09-13. https://mathworld.wolfram.com/LogNormalDistribution.html

  2. "1.3.6.6.9. Lognormal Distribution". www.itl.nist.gov. U.S. National Institute of Standards and Technology (NIST). Retrieved 2020-09-13. https://www.itl.nist.gov/div898/handbook/eda/section3/eda3669.htm

  3. 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 978-0-471-58495-7

  4. 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 978-0-471-58495-7

  5. Park, Sung Y.; Bera, Anil K. (2009). "Maximum entropy autoregressive conditional heteroskedasticity model" (PDF). Journal of Econometrics. 150 (2): 219–230, esp. Table 1, p. 221. CiteSeerX 10.1.1.511.9750. doi:10.1016/j.jeconom.2008.12.014. Archived from the original (PDF) on 2016-03-07. Retrieved 2011-06-02. https://web.archive.org/web/20160307144515/http://wise.xmu.edu.cn/uploadfiles/paper-masterdownload/2009519932327055475115776.pdf

  6. Weisstein, Eric W. "Log Normal Distribution". mathworld.wolfram.com. Retrieved 2020-09-13. https://mathworld.wolfram.com/LogNormalDistribution.html

  7. 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 978-0-471-58495-7

  8. Weisstein, Eric W. "Log Normal Distribution". mathworld.wolfram.com. Retrieved 2020-09-13. https://mathworld.wolfram.com/LogNormalDistribution.html

  9. Tarmast, Ghasem (2001). Multivariate Log–Normal Distribution (PDF). ISI Proceedings: 53rd Session. Seoul. Archived (PDF) from the original on 2013-07-19. https://isi.cbs.nl/iamamember/CD2/pdf/329.PDF

  10. Halliwell, Leigh (2015). The Lognormal Random Multivariate (PDF). Casualty Actuarial Society E-Forum, Spring 2015. Arlington, VA. Archived (PDF) from the original on 2015-09-30. http://www.casact.org/pubs/forum/15spforum/Halliwell.pdf

  11. Heyde, CC. (2010), "On a Property of the Lognormal Distribution", Journal of the Royal Statistical Society, Series B, vol. 25, no. 2, pp. 392–393, doi:10.1007/978-1-4419-5823-5_6, ISBN 978-1-4419-5822-8 978-1-4419-5822-8

  12. Billingsley, Patrick (2012). Probability and Measure (Anniversary ed.). Hoboken, N.J.: Wiley. p. 415. ISBN 978-1-118-12237-2. OCLC 780289503. 978-1-118-12237-2

  13. 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. /wiki/Doi_(identifier)

  14. 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. /wiki/Doi_(identifier)

  15. Barakat, R. (1976). "Sums of independent lognormally distributed random variables". Journal of the Optical Society of America. 66 (3): 211–216. Bibcode:1976JOSA...66..211B. doi:10.1364/JOSA.66.000211. /wiki/Bibcode_(identifier)

  16. 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. http://dspace.mit.edu/bitstream/handle/1721.1/48703/onsumsoflognorma00baro.pdf

  17. 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. https://www.cambridge.org/core/services/aop-cambridge-core/content/view/F1563B5AD8918EF2CD51092F82EB0B73/S0334270000006901a.pdf/div-class-title-on-lognormal-random-variables-i-the-characteristic-function-div.pdf

  18. 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). https://link.springer.com/article/10.1007/s11009-014-9430-7

  19. Das, Abhranil (2021). "A method to integrate and classify normal distributions". Journal of Vision. 21 (10): 1. arXiv:2012.14331. doi:10.1167/jov.21.10.1. PMC 8419883. PMID 34468706. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8419883

  20. Das, Abhranil (2021). "A method to integrate and classify normal distributions". Journal of Vision. 21 (10): 1. arXiv:2012.14331. doi:10.1167/jov.21.10.1. PMC 8419883. PMID 34468706. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8419883

  21. Kirkwood, Thomas BL (Dec 1979). "Geometric means and measures of dispersion". Biometrics. 35 (4): 908–909. JSTOR 2530139. /wiki/JSTOR_(identifier)

  22. 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. https://doi.org/10.1641%2F0006-3568%282001%29051%5B0341%3ALNDATS%5D2.0.CO%3B2

  23. Kirkwood, Thomas BL (Dec 1979). "Geometric means and measures of dispersion". Biometrics. 35 (4): 908–909. JSTOR 2530139. /wiki/JSTOR_(identifier)

  24. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5216879

  25. 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 978-0-471-58495-7

  26. Weisstein, Eric W. "Log Normal Distribution". mathworld.wolfram.com. Retrieved 2020-09-13. https://mathworld.wolfram.com/LogNormalDistribution.html

  27. "1.3.6.6.9. Lognormal Distribution". www.itl.nist.gov. U.S. National Institute of Standards and Technology (NIST). Retrieved 2020-09-13. https://www.itl.nist.gov/div898/handbook/eda/section3/eda3669.htm

  28. Sawant, S.; Mohan, N. (2011) "FAQ: Issues with Efficacy Analysis of Clinical Trial Data Using SAS" Archived 24 August 2011 at the Wayback Machine, PharmaSUG2011, Paper PO08 http://pharmasug.org/proceedings/2011/PO/PharmaSUG-2011-PO08.pdf

  29. Schiff, MH; et al. (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. doi:10.1136/annrheumdis-2014-205228. PMC 4112421. PMID 24728329. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4112421

  30. 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 978-0-471-58495-7

  31. Daly, Leslie E.; Bourke, Geoffrey Joseph (2000). Interpretation and Uses of Medical Statistics. Vol. 46 (5th ed.). Oxford, UK: Wiley-Blackwell. p. 89. doi:10.1002/9780470696750. ISBN 978-0-632-04763-5. PMC 1059583; {{cite book}}: |journal= ignored (help) print edition. Online eBook ISBN 9780470696750 978-0-632-04763-5

  32. "ProbOnto". Retrieved 1 July 2017. http://www.probonto.org

  33. Swat, MJ; Grenon, P; Wimalaratne, S (2016). "ProbOnto: ontology and knowledge base of probability distributions". Bioinformatics. 32 (17): 2719–2721. doi:10.1093/bioinformatics/btw170. PMC 5013898. PMID 27153608. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5013898

  34. Forbes et al. Probability Distributions (2011), John Wiley & Sons, Inc.

  35. Forbes et al. Probability Distributions (2011), John Wiley & Sons, Inc.

  36. Lunn, D. (2012). The BUGS book: a practical introduction to Bayesian analysis. Texts in statistical science. CRC Press.

  37. Limpert, E.; Stahel, W. A.; Abbt, M. (2001). "Log-normal 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. https://doi.org/10.1641%2F0006-3568%282001%29051%5B0341%3ALNDATS%5D2.0.CO%3B2

  38. 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. PMID 22640817. /wiki/Doi_(identifier)

  39. 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. doi:10.1016/S0169-2607(00)00117-6. PMID 11275334. /wiki/Doi_(identifier)

  40. The PopED Development Team (2016). PopED Manual, Release version 2.13. Technical report, Uppsala University.

  41. ProbOnto website, URL: http://probonto.org http://probonto.org

  42. Cococcioni, Marco; Fiorini, Francesco; Pagano, Michele (2023-04-06). "Modelling Heavy Tailed Phenomena Using a LogNormal Distribution Having a Numerically Verifiable Infinite Variance". Mathematics. 11 (7): 1758. doi:10.3390/math11071758. hdl:11568/1216554. ISSN 2227-7390. https://doi.org/10.3390%2Fmath11071758

  43. 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. /wiki/Doi_(identifier)

  44. 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). /wiki/Doi_(identifier)

  45. 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. /wiki/Doi_(identifier)

  46. 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. https://doi.org/10.1155%2F2009%2F630857

  47. 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. https://hal.archives-ouvertes.fr/hal-00595951/file/PEER_stage2_10.1016%252Fj.spl.2008.03.035.pdf

  48. 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. https://doi.org/10.1155%2F2009%2F630857

  49. 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. /wiki/Doi_(identifier)

  50. 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: IEEE. pp. 1880–1890. arXiv:1705.03196. doi:10.1109/WSC.2017.8247924. ISBN 978-1-5386-3428-8. 978-1-5386-3428-8

  51. Asmussen, A.; Goffard, P.-O.; Laub, P. J. (2016). "Orthonormal polynomial expansions and lognormal sum densities". arXiv:1601.01763v1 [math.PR]. /wiki/ArXiv_(identifier)

  52. Sangal, B.; Biswas, A. (1970). "The 3-Parameter Lognormal Distribution Applications in Hydrology". Water Resources Research. 6 (2): 505–515. doi:10.1029/WR006i002p00505. /wiki/Doi_(identifier)

  53. Johnson, N. L. (1949). "Systems of Frequency Curves Generated by Methods of Translation". Biometrika. 36 (1/2): 149–176. doi:10.2307/2332539. JSTOR 2332539. PMID 18132090. /wiki/Norman_Lloyd_Johnson

  54. 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). /wiki/Doi_(identifier)

  55. Shen, Wei-Hsiung. "Estimation of parameters of a lognormal distribution." Taiwanese Journal of Mathematics 2.2 (1998): 243–250. pdf https://projecteuclid.org/journalArticle/Download?urlid=10.11650%2Ftwjm%2F1500406934

  56. Henry (https://math.stackexchange.com/users/6460/henry), Method of moments estimator for lognormal distribution, URL (version: 2022-01-12): https://math.stackexchange.com/q/4355343 https://math.stackexchange.com/users/6460/henry

  57. Finney, D. J. "On the distribution of a variate whose logarithm is normally distributed." Supplement to the Journal of the Royal Statistical Society 7.2 (1941): 155–161.

  58. Longford, Nicholas T. "Inference with the lognormal distribution." Journal of Statistical Planning and Inference 139.7 (2009): 2329–2340.

  59. Zellner, Arnold. "Bayesian and non-Bayesian analysis of the log-normal distribution and log-normal regression." Journal of the American Statistical Association 66.334 (1971): 327–330.

  60. Tang, Qi. "Comparison of different methods for estimating log-normal means". MS thesis. East Tennessee State University, 2014. linkhttps://dc.etsu.edu/cgi/viewcontent.cgi?article=3728&context=etd#page=12.13 pdf] https://www.proquest.com/docview/1547379248?pq-origsite=gscholar&fromopenview=true

  61. Kwon, Yeil. "An alternative method for estimating lognormal means." Communications for Statistical Applications and Methods 28.4 (2021): 351–368. link http://www.csam.or.kr/journal/view.html?doi=10.29220/CSAM.2021.28.4.351

  62. Olsson, Ulf. "Confidence intervals for the mean of a log-normal distribution." Journal of Statistics Education 13.1 (2005).pdf html https://www.tandfonline.com/doi/pdf/10.1080/10691898.2005.11910638

  63. 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 https://stats.stackexchange.com/q/33395

  64. The Cox Method was quoted as "personal communication" in Land, 1971,[52] and was also given in CitationZhou and Gao (1997)[53] and Olsson 2005[50]: Section 3.3

  65. Olsson, Ulf. "Confidence intervals for the mean of a log-normal distribution." Journal of Statistics Education 13.1 (2005).pdf html https://www.tandfonline.com/doi/pdf/10.1080/10691898.2005.11910638

  66. Confidence Intervals for Risk Ratios and Odds Ratios https://sphweb.bumc.bu.edu/otlt/MPH-Modules/PH717-QuantCore/PH717-Module8-CategoricalData/PH717-Module8-CategoricalData5.html?fbclid=IwY2xjawFeH3JleHRuA2FlbQIxMAABHbmxa15uyyzJuzEwh9PIUr_m2Jsc9NGiPuS6IwfA36Ca5r1wV1EoPEz3MQ_aem_03PRd_jlRfbsnr6xCPkZmw

  67. The issue is that we don't know how to do it directly, so we take their logs, and then use the delta method to say that their logs is itself (approximately) normal. This trick allows us to pretend that their exp was log normal, and use that approximation to build the CI. Notice that in the RR case, the median and the mean in the base distribution (i.e., before taking the log), is actually identical (since they are originally normal, and not log normal). E.g., p ^ 1 ∼ ˙ N ( p 1 , p 1 ( 1 − p 1 ) / n ) {\displaystyle {\hat {p}}_{1}{\dot {\sim }}N(p_{1},p_{1}(1-p1)/n)} and l o g ( p ^ 1 ) ∼ ˙ N ( l o g ( p 1 ) , ( 1 − p 1 ) / ( p 1 ∗ n ) ) {\displaystyle log({\hat {p}}_{1}){\dot {\sim }}N(log(p_{1}),(1-p1)/(p_{1}*n))} Hence, building a CI based on the log and than back-transform will give us C I ( p 1 ) : e l o g ( p ^ 1 ) ± ( 1 − p ^ 1 ) / ( p ^ 1 ∗ n ) ) {\displaystyle CI(p_{1}):e^{log({\hat {p}}_{1})\pm (1-{\hat {p}}_{1})/({\hat {p}}_{1}*n))}} . So while we expect the CI to be for the median, in this case, it's actually also for the mean in the original distribution. i.e., if the original p ^ 1 {\displaystyle {\hat {p}}_{1}} was log-normal, we'd expect that E [ p ^ 1 ] = e l o g ( p 1 ) + 1 / 2 ∗ ( 1 − p 1 ) / ( p 1 ∗ n ) {\displaystyle E[{\hat {p}}_{1}]=e^{log(p_{1})+1/2*(1-p1)/(p_{1}*n)}} . But in practice, we KNOW that E [ p ^ 1 ] = e l o g ( p 1 ) = p 1 {\displaystyle E[{\hat {p}}_{1}]=e^{log(p_{1})}=p_{1}} . Hence, the approximation we have is in the second step (of the delta method), but the CI are actually for the expectation (not just the median). This is because we are starting from a base distribution that is normal, and then using another approximation after the log again to normal. This means that a big approximation part of the CI is from the delta method. /wiki/Delta_method

  68. The formula can found by just treating the estimated means and variances as approximetly normal, which indicates the terms is itself a log-normal, enabling us to quickly get the expectation. The bias can be partially minimized by using: [ E ( X 1 ) E ( X 2 ) ] ^ = [ E ^ ( X 1 ) E ^ ( X 2 ) ] 2 ( σ 1 2 n 1 + σ 2 2 n 2 + σ 1 4 2 ( n 1 − 1 ) + σ 2 4 2 ( n 2 − 1 ) ) ^ ≈ [ e ( μ ^ 1 − μ ^ 2 ) + 1 2 ( S 1 2 − S 2 2 ) ] 2 S 1 2 n 1 + S 2 2 n 2 + S 1 4 2 ( n 1 − 1 ) + S 2 4 2 ( n 2 − 1 ) {\displaystyle {\widehat {\left[{\frac {E(X_{1})}{E(X_{2})}}\right]}}=\left[{\frac {{\widehat {E}}(X_{1})}{{\widehat {E}}(X_{2})}}\right]{\frac {2}{\widehat {\left({\frac {\sigma _{1}^{2}}{n_{1}}}+{\frac {\sigma _{2}^{2}}{n_{2}}}+{\frac {\sigma _{1}^{4}}{2(n_{1}-1)}}+{\frac {\sigma _{2}^{4}}{2(n_{2}-1)}}\right)}}}\approx \left[e^{({\widehat {\mu }}_{1}-{\widehat {\mu }}_{2})+{\frac {1}{2}}(S_{1}^{2}-S_{2}^{2})}\right]{\frac {2}{{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}}

  69. 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. https://doi.org/10.3390%2Fe19020056

  70. 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. https://doi.org/10.3390%2Fe19020056

  71. 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. Bibcode:2003PEngM..18..241W. doi:10.1016/S0266-8920(03)00028-6. /wiki/Bibcode_(identifier)

  72. 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. PMC 7111546. PMID 32288765. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7111546

  73. 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. Bibcode:2019AcAau.155..118B. doi:10.1016/j.actaastro.2018.11.033. S2CID 117598888. /wiki/Bibcode_(identifier)

  74. Sutton, John (Mar 1997). "Gibrat's Legacy". Journal of Economic Literature. 32 (1): 40–59. JSTOR 2729692. /wiki/JSTOR_(identifier)

  75. Limpert, Eckhard; Stahel, Werner A.; Abbt, Markus (2001). "Log-normal Distributions across the Sciences: Keys and Clues". BioScience. 51 (5): 341. doi:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2. ISSN 0006-3568. https://doi.org/10.1641%2F0006-3568%282001%29051%5B0341%3ALNDATS%5D2.0.CO%3B2

  76. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4051294

  77. Pawel, Sobkowicz; et al. (2013). "Lognormal distributions of user post lengths in Internet discussions – a consequence of the Weber-Fechner law?". EPJ Data Science.

  78. 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. http://mldm.ict.ac.cn/platform/pweb/academicDetail.htm?id=16

  79. "What is the average length of a game of chess?". chess.stackexchange.com. Retrieved 14 April 2018. http://chess.stackexchange.com/questions/2506/what-is-the-average-length-of-a-game-of-chess/4899#4899

  80. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5216879

  81. Huxley, Julian S. (1932). Problems of relative growth. London. ISBN 978-0-486-61114-3. OCLC 476909537. 978-0-486-61114-3

  82. Sartwell, Philip E. "The distribution of incubation periods of infectious disease." American journal of hygiene 51 (1950): 310–318.

  83. Limpert, Eckhard; Stahel, Werner A.; Abbt, Markus (2001). "Log-normal Distributions across the Sciences: Keys and Clues". BioScience. 51 (5): 341. doi:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2. ISSN 0006-3568. https://doi.org/10.1641%2F0006-3568%282001%29051%5B0341%3ALNDATS%5D2.0.CO%3B2

  84. S. K. Chan, Jennifer; Yu, Philip L. H. (2006). "Modelling SARS data using threshold geometric process". Statistics in Medicine. 25 (11): 1826–1839. doi:10.1002/sim.2376. PMID 16345017. S2CID 46599163. /wiki/Doi_(identifier)

  85. Ono, Yukiteru; Asai, Kiyoshi; Hamada, Michiaki (2013-01-01). "PBSIM: PacBio reads simulator—toward accurate genome assembly". Bioinformatics. 29 (1): 119–121. doi:10.1093/bioinformatics/bts649. ISSN 1367-4803. PMID 23129296. https://academic.oup.com/bioinformatics/article/29/1/119/273243

  86. 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. /wiki/Doi_(identifier)

  87. 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. 7 (1): 171–178. doi:10.1080/10543409708835177. ISSN 1054-3406. PMID 9056596. https://www.tandfonline.com/doi/full/10.1080/10543409708835177

  88. Scheler, Gabriele; Schumann, Johann (2006-10-08). Diversity and stability in neuronal output rates. 36th Society for Neuroscience Meeting, Atlanta.

  89. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3804159

  90. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4051294

  91. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5985929

  92. 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. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5639933

  93. 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. doi:10.1093/cercor/bhad160. PMC 10438924. PMID 37409647. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10438924

  94. 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. Bibcode:2018NJPh...20e3057P. doi:10.1088/1367-2630/aac3c7. hdl:2164/10561. ISSN 1367-2630. https://iopscience.iop.org/article/10.1088/1367-2630/aac3c7

  95. Ahrens, L. H. (1954-02-01). "The lognormal distribution of the elements (A fundamental law of geochemistry and its subsidiary)". Geochimica et Cosmochimica Acta. 5 (2): 49–73. Bibcode:1954GeCoA...5...49A. doi:10.1016/0016-7037(54)90040-X. ISSN 0016-7037. https://dx.doi.org/10.1016/0016-7037%2854%2990040-X

  96. Limpert, Eckhard; Stahel, Werner A.; Abbt, Markus (2001). "Log-normal Distributions across the Sciences: Keys and Clues". BioScience. 51 (5): 341. doi:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2. ISSN 0006-3568. https://doi.org/10.1641%2F0006-3568%282001%29051%5B0341%3ALNDATS%5D2.0.CO%3B2

  97. 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. 978-90-70754-33-4

  98. CumFreq, free software for distribution fitting https://www.waterlog.info/cumfreq.htm

  99. Wadhams, Peter (1988). "Winter observations of iceberg frequencies and sizes in the South Atlantic Ocean". Journal of Geophysical Research: Oceans. 93 (C4): 3583–3590. doi:10.1029/JC093iC04p03583. /wiki/Doi_(identifier)

  100. Heintzenberg, Jost (1994). "Properties of the Log-Normal Particle Size Distribution". Aerosol Science and Technology. 21 (1): 46–48. doi:10.1080/02786829408959695. https://doi.org/10.1080%2F02786829408959695

  101. Limpert, Eckhard; Stahel, Werner A.; Abbt, Markus (2001). "Log-normal 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. /wiki/Doi_(identifier)

  102. Di Giorgio, C; Krempff, A; Guiraud, H; Binder, P; Tiret, C; Dumenil, G (1996). "Atmospheric pollution by airborne microorganisms in the city of Marseilles". Atmospheric Environment. 30: 155–160. doi:10.1016/1352-2310(95)00143-M. /wiki/Doi_(identifier)

  103. Sheridan, M.F.; Wohletz, K.H.; Dehn, J. (1987). "Discrimination of grain-size subpopulations in pyroclastic deposits". Geology. 15 (4): 367–370. doi:10.1130/0091-7613(1987)15<367:DOGSIP>2.0.CO;2. /wiki/Doi_(identifier)

  104. Wohletz, K.H.; Sheridan, M.F.; Brown, W.K. (1989). "Particle size distributions and the sequential fragmentation/transport theory applied to volcanic ash". Journal of Geophysical Research: Solid Earth. 94 (B11): 15703–15721. doi:10.1029/JB094iB11p15703. /wiki/Doi_(identifier)

  105. Clementi, Fabio; Gallegati, Mauro (2005) "Pareto's law of income distribution: Evidence for Germany, the United Kingdom, and the United States", EconWPA /wiki/Mauro_Gallegati

  106. 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. /wiki/ArXiv_(identifier)

  107. Black, F.; Scholes, M. (1973). "The Pricing of Options and Corporate Liabilities". Journal of Political Economy. 81 (3): 637. doi:10.1086/260062. S2CID 154552078. /wiki/Doi_(identifier)

  108. Mandelbrot, Benoit (2004). The (mis-)Behaviour of Markets. Basic Books. ISBN 9780465043552. 9780465043552

  109. Bunchen, P., Advanced Option Pricing, University of Sydney coursebook, 2007

  110. Thelwall, Mike; Wilson, Paul (2014). "Regression for citation data: An evaluation of different methods". Journal of Informetrics. 8 (4): 963–971. arXiv:1510.08877. doi:10.1016/j.joi.2014.09.011. S2CID 8338485. /wiki/ArXiv_(identifier)

  111. Sheridan, Paul; Onodera, Taku (2020). "A Preferential Attachment Paradox: How Preferential Attachment Combines with Growth to Produce Networks with Log-normal In-degree Distributions". Scientific Reports. 8 (1): 2811. arXiv:1703.06645. doi:10.1038/s41598-018-21133-2. PMC 5809396. PMID 29434232. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5809396

  112. Eeckhout, Jan (2004). "Gibrat's Law for (All) Cities". American Economic Review. 94 (5): 1429–1451. doi:10.1257/0002828043052303. JSTOR 3592829. https://www.jstor.org/stable/3592829

  113. Kault, David (1996). "The Shape of the Distribution of the Number of Sexual Partners". Statistics in Medicine. 15 (2): 221–230. doi:10.1002/(SICI)1097-0258(19960130)15:2<221::AID-SIM148>3.0.CO;2-Q. PMID 8614756. https://onlinelibrary.wiley.com/doi/10.1002/(SICI)1097-0258(19960130)15:2%3C221::AID-SIM148%3E3.0.CO;2-Q

  114. O'Connor, Patrick; Kleyner, Andre (2011). Practical Reliability Engineering. John Wiley & Sons. p. 35. ISBN 978-0-470-97982-2. 978-0-470-97982-2

  115. "Shadowing". www.WirelessCommunication.NL. Archived from the original on January 13, 2012. https://web.archive.org/web/20120113201345/http://wireless.per.nl/reference/chaptr03/shadow/shadow.htm

  116. Dexter, A. R.; Tanner, D. W. (July 1972). "Packing Densities of Mixtures of Spheres with Log-normal Size Distributions". Nature Physical Science. 238 (80): 31–32. Bibcode:1972NPhS..238...31D. doi:10.1038/physci238031a0. ISSN 2058-1106. https://www.nature.com/articles/physci238031a0

  117. 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. S2CID 17404692. /wiki/ArXiv_(identifier)

  118. 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. doi:10.1145/301464.301480. ISSN 0163-5999. https://doi.org/10.1145%2F301464.301480

  119. Pawel, Sobkowicz; et al. (2013). "Lognormal distributions of user post lengths in Internet discussions – a consequence of the Weber-Fechner law?". EPJ Data Science.

  120. Pawel, Sobkowicz; et al. (2013). "Lognormal distributions of user post lengths in Internet discussions – a consequence of the Weber-Fechner law?". EPJ Data Science.

  121. 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]. /wiki/ArXiv_(identifier)

  122. ASTM D3654, Standard Test Method for Shear Adhesion on Pressure-Sensitive Tapesw

  123. ASTM D4577, Standard Test Method for Compression Resistance of a container Under Constant Load>\