Biostatistics Study NotesFrom concepts and derivations to interpretation

CHAPTER 02 · TOPIC 07

Exponential Distribution

Study waiting times under a constant event rate, memorylessness, survival functions, and the Poisson process.

On this page
  1. In what cases should you use it?
  2. Definition of exponential distribution
  3. How does λ change the curve?
  4. Find probability from density
  5. Average waiting time and variation
  6. Memoryless
  7. Relationship with Poisson process
  8. Relationship with survival time analysis

The exponential distribution is a continuous probability distribution that describes waiting times between independent events occurring at a constant average rate.

In what cases should you use it?#

It is used to measure the time until the next call or customer arrives, or the failure time of equipment assuming a certain failure rate. Check that your data is not only latency, but also that the average incidence is roughly constant and that events in different time intervals can be considered independent.

Definition of exponential distribution#

If the random variable X is the waiting time until the next event and the average occurrence rate is λ, then X follows an exponential distribution with a parameter λ.

Distribution notationXExp(λ),λ>0X\sim\operatorname{Exp}(\lambda),\qquad \lambda>0
probability density functionf(x)={λeλx,x0,0,x<0.f(x)=\begin{cases}\lambda e^{-\lambda x},&x\ge 0,\\0,&x<0.\end{cases}
symbolmeaning
XWait time until next event
xPossible values ​​of X
λAverage incidence rate per unit time
f(x)Probability density of waiting time around x
ebase of natural logarithm

How does λ change the curve?#

The density is f(0)=λ at x=0 and then decreases with latency. The larger λ, the more frequent the event, the faster the curve drops from a higher point, and the shorter the typical waiting time.

Exponential density curves with rates 0.5, 1.0, and 1.5; a larger rate gives a higher starting point and a faster decline
A larger rate means events occur more frequently on average, so waiting times concentrate closer to 0. The total area under each curve is 1.Source: Created for this site from the exponential density function

Find probability from density#

Since it is a continuous distribution, the probability of a single point is 0. The probability that an event will occur within a certain amount of time, or the probability that it will not occur even after that period of time, is calculated using the area under the curve.

Cumulative Distribution FunctionF(x)=P(Xx)=1eλx,x0F(x)=P(X\le x)=1-e^{-\lambda x},\qquad x\ge0
survival functionS(x)=P(X>x)=eλx,x0S(x)=P(X>x)=e^{-\lambda x},\qquad x\ge0
Probability of falling in the intervalP(a<Xb)=F(b)F(a)=eλaeλbP(a<X\le b)=F(b)-F(a)=e^{-\lambda a}-e^{-\lambda b}
symbolmeaning
f(x)probability density at position x
F(x)Cumulative probability that waiting time is less than or equal to x
S(x)Probability that the event will not occur even if x is exceeded
a、bEndpoints of the interval where 0≤a<b

Average waiting time and variation#

Occurrence rate and average waiting time are inverses of each other. The larger λ, the smaller the average and variation of the waiting time.

Expected valueE(X)=1λE(X)=\frac{1}{\lambda}
variance and standard deviationVar(X)=1λ2,SD(X)=1λ\operatorname{Var}(X)=\frac{1}{\lambda^2},\qquad \operatorname{SD}(X)=\frac{1}{\lambda}

Memoryless#

Exponential distributions are memoryless. Given that no event has occurred after waiting s hours, the probability of waiting for more than t is the same as the probability of waiting more than t from the beginning.

P(X>s+tX>s)=P(X>t)=eλtP(X>s+t\mid X>s)=P(X>t)=e^{-\lambda t}
Supplement: Deriving exponential distribution from short time intervals

Divide a time of length x into n short intervals Δt=x/n. For a constant occurrence rate λ, the probability that an event does not occur in each interval is approximately 1−λx/n.

No events in each sectionP(no event in Δt)1λxnP(\text{no event in }\Delta t)\approx 1-\lambda\frac{x}{n}
No event in all n intervalsP(X>x)(1λxn)nP(X>x)\approx\left(1-\lambda\frac{x}{n}\right)^n
Make the interval as small as possibleS(x)=P(X>x)=limn(1λxn)n=eλxS(x)=P(X>x)=\lim_{n\to\infty}\left(1-\lambda\frac{x}{n}\right)^n=e^{-\lambda x}

This limit first gives the survival function. Subtracting it from 1 gives you the cumulative distribution function, and differentiating it with x gives you the probability density function.

From survival function to cumulative probabilityF(x)=1S(x)=1eλxF(x)=1-S(x)=1-e^{-\lambda x}
From cumulative probability to densityf(x)=F(x)=S(x)=λeλxf(x)=F'(x)=-S'(x)=\lambda e^{-\lambda x}
Supplement: Why the average waiting time is 1/λ

The expected value of a non-negative random variable can be expressed as the area under the survival function. Integrating S(x)=e^(−λx) yields 1/λ.

E(X)=0S(x)dx=0eλxdx=1λE(X)=\int_0^{\infty}S(x)\,dx=\int_0^{\infty}e^{-\lambda x}\,dx=\frac{1}{\lambda}

Relationship with Poisson process#

If the number of occurrences N(t) in time t follows a Poisson distribution with mean λt, then the time until the first event X follows an exponential distribution with parameter λ. Waiting more than t is the same as not having the event occur at all within time t.

Number of occurrencesN(t)Poisson(λt)N(t)\sim\operatorname{Poisson}(\lambda t)
waiting timeP(X>t)=P(N(t)=0)=eλtP(X>t)=P(N(t)=0)=e^{-\lambda t}

Relationship with survival time analysis#

In survival time analysis, the exponential distribution is a model with a constant λ hazard rate. The hazard rate represents the instantaneous probability that an event will occur in a short period of time, given that you have survived until time x.

h(x)=f(x)S(x)=λh(x)=\frac{f(x)}{S(x)}=\lambda