Ex-Gaussian Distribution#

Univariate, Continuous, Asymmetric, Unbounded, Light-tailed

The Ex-Gaussian distribution (exponentially modified Gaussian distribution or EMG) is a continuous probability distribution that results from the convolution of a normal distribution and an exponential distribution. It is characterized by three parameters: \(\mu\), \(\sigma\), and \(\nu\), which are the mean and standard deviation of the normal component, and the mean of the exponential component, respectively. It has a bell-shaped curve like the normal distribution, but with a positive skew due to the exponential component.

The Ex-Gaussian distribution is commonly used to model reaction times in psychology. It is also used to model the shape of chromatographic peaks, the intermitotic times of cell division, cluster ion beams, and other phenomena.

Key properties and parameters#

Support	\(x \in \mathbb{R}\)
Mean	\(\mu + \nu\)
Variance	\(\sigma^2 + \nu^2\)

Parameters:

\(\mu\) : (float) Mean of the normal component.
\(\sigma\) : (float) Standard deviation of the normal component, \(\sigma > 0\).
\(\nu\) : (float) Mean of the exponential component, \(\nu > 0\).

Probability Density Function (PDF)#

\[ f(x \mid \mu, \sigma, \nu) = \frac{1}{\nu}\exp\left(\frac{\mu-x}{\nu}+\frac{\sigma^2}{2\nu^2}\right)\Phi\left(\frac{x-\mu}{\sigma}-\frac{\sigma}{\nu}\right) \]

where \(\Phi\) is the standard normal CDF.

/home/docs/checkouts/readthedocs.org/user_builds/preliz/envs/stable/lib/python3.11/site-packages/pytensor/link/c/cmodule.py:2986: UserWarning: PyTensor could not link to a BLAS installation. Operations that might benefit from BLAS will be severely degraded.
This usually happens when PyTensor is installed via pip. We recommend it be installed via conda/mamba/pixi instead.
Alternatively, you can use an experimental backend such as Numba or JAX that perform their own BLAS optimizations, by setting `pytensor.config.mode == 'NUMBA'` or passing `mode='NUMBA'` when compiling a PyTensor function.
For more options and details see https://pytensor.readthedocs.io/en/latest/troubleshooting.html#how-do-i-configure-test-my-blas-library
  warnings.warn(

Cumulative Distribution Function (CDF)#

\[ F(x \mid \mu, \sigma, \nu) = \Phi(x, \mu, \sigma) - \frac{1}{2} \exp\left( \frac{1}{2\nu} \left(2\mu + \frac{\sigma^2}{\nu} - 2x\right) \right) \left(1 + \operatorname{erf}\left( \frac{x - \left(\mu + \frac{\sigma^2}{\nu}\right)}{\sigma \sqrt{2}} \right)\right) \]

where \(\Phi\) is the CDF of a Gaussian distribution, \(\text{erf}\) is the error function

References#

Wikipedia - Ex-Gaussian distribution