Skip to content
Oakfield Operator Calculus Function Reference Site

Laplace Distribution

The Laplace distribution centers probability mass at a location parameter with exponential decay on both sides, controlled by a scale parameter.

For XLaplace(μ,b)X\sim\operatorname{Laplace}(\mu,b), the density is

f(x;μ,b)=12bexp ⁣(xμb)f(x;\mu,b) = \frac{1}{2b}\exp\!\left(-\frac{|x-\mu|}{b}\right)

Defined over real-valued outcomes; parameters set the mean and spread.


E[X]=μ,Var(X)=2b2\mathbb{E}[X] = \mu, \qquad \operatorname{Var}(X) = 2b^2

Tails decay exponentially (heavier than Gaussian, lighter than Cauchy). Inverse-CDF sampling can be expressed using UU(0,1)U\sim U(0,1) with a sign and logarithm transform.


(μ,b)=(0,1)(\mu,b) = (0,1) defines the standard Laplace density. As b0+b \to 0^+, the distribution collapses to a point mass at μ\mu; as bb increases, variance grows linearly in b2b^2. Absolute deviations Xμ|X-\mu| follow an exponential distribution.


  • (μ,b)=(0,1)(\mu,b)=(0,1) is the standard Laplace distribution.
  • b0+b\to 0^+ collapses to a point mass at μ\mu.

Gaussian distributions provide a lighter-tailed alternative; exponential distributions govern one-sided magnitudes; logistic and double-exponential families share similar shapes with differing normalization.


Oakfield uses Laplace noise primarily via integrator stochastic hooks:

  • Integrator noise option: setting integrator noise = "laplace" selects a unit-variance Laplace source (implemented via a uniform RNG and a log1p transform with clamping away from log(0)).
  • Useful for heavier-tailed stochastic increments compared to Gaussian noise when running stochastic integrator steps.

Laplace-type distributions arise naturally when modeling absolute-deviation penalties and exponentially decaying tails on both sides of a central location.

They are widely used in robust statistics and sparse modeling (via the 1\ell_1 penalty connection).


  • Kotz, Kozubowski, Podgórski, The Laplace Distribution and Generalizations
  • Devroye, Non-Uniform Random Variate Generation