Skip to content

Weibull distribution

The Weibull distribution models positive waiting times with changing hazard. Separate shape from scale, calculate conditional survival, and distinguish duration from intensity.

The Weibull distribution is a family of continuous probability distributions for nonnegative quantities, especially waiting times and lifetimes. Its two main parameters are shape, which controls how the event rate changes as time passes, and scale, which sets the time units. Unlike an exponential distribution, it need not give a fresh clock and an already running clock the same remaining-wait distribution.

That flexibility matters when a financial operation has a clock: time to resolve a reconciliation exception, time until equipment failure, or time to default under a specified survival model. A distribution supplies probabilities for the stated event; fitting a positive curve alone does not establish the economic mechanism behind it.

Start with the survival curve

Use the unshifted, two-parameter convention

T∼Weibull⁡(k,η),k>0,η>0.T\sim\operatorname{Weibull}(k,\eta),\qquad k>0,\quad\eta>0.
Here kk is dimensionless shape and η\eta is scale, measured in the same units as TT. For t≥0t\ge0, the survival function and cumulative distribution function are
S(t)=P(T>t)=e−(t/η)k,F(t)=1−e−(t/η)k.S(t)=\mathbb{P}(T>t)=e^{-(t/\eta)^k}, \qquad F(t)=1-e^{-(t/\eta)^k}.
For t<0t<0, F(t)=0F(t)=0 and S(t)=1S(t)=1. Differentiating the CDF gives the probability density:
f(t)=kη(tη)k−1e−(t/η)k,t>0.f(t)=\frac{k}{\eta}\left(\frac{t}{\eta}\right)^{k-1}e^{-(t/\eta)^k},\qquad t>0.

At exactly one time, probability is zero: P(T=t)=0\mathbb{P}(T=t)=0. Probability over an interval is an area under this density, while S(t)S(t) directly answers whether the wait exceeds a threshold.

The hazard rate is the conditional event rate among cases still waiting:
h(t)=f(t)S(t)=kη(tη)k−1,t>0.h(t)=\frac{f(t)}{S(t)} =\frac{k}{\eta}\left(\frac{t}{\eta}\right)^{k-1},\qquad t>0.
More precisely, h(t) Δth(t)\,\Delta t approximates P(t<T≤t+Δt∣T>t)\mathbb{P}(t<T\le t+\Delta t\mid T>t) for a sufficiently short positive interval. Hazard has units of inverse time and can exceed one; it is a rate, not a probability. The cumulative hazard is H(t)=−ln⁡S(t)=(t/η)kH(t)=-\ln S(t)=(t/\eta)^k.
ShapeHazard as elapsed time increasesInterpretation of the remaining clock
0<k<10<k<1DecreasesCases that have lasted longer face a lower immediate event rate
k=1k=1Constant at 1/η1/\etaExponential; the clock is memoryless
k>1k>1IncreasesCases that have lasted longer face a higher immediate event rate

For k<1k<1, both density and hazard grow without bound as t↓0t\downarrow0. The density spike is integrable: it creates no probability atom at zero. At k=1k=1 the density's right limit is 1/η1/\eta; at k>1k>1 it is zero. An ordinary Weibull hazard is monotone, so it cannot represent a hazard that first falls and then rises.

A ten-minute scale is not a ten-minute average

Suppose an operations team uses a hypothetical model for the time until an unresolved trade reconciliation is completed, measured from when it enters a particular workflow:
k=2,η=10 minutes.k=2,\qquad \eta=10\text{ minutes}.

This is a stipulated teaching model, not a calibration to operational records. It gives

S(t)=e−(t/10)2,h(t)=210t10 min−1,S(t)=e^{-(t/10)^2},\qquad h(t)=\frac{2}{10}\frac{t}{10}\text{ min}^{-1},

where tt in the numerical expression is elapsed minutes. The model's increasing completion hazard says older unresolved cases have a greater immediate chance of completion, which could be a hypothesis about an escalation process. It does not say the hypothesis is true.

QuestionResult
Probability still unresolved after ten minutesS(10)=e−1≈0.367879S(10)=e^{-1}\approx0.367879, or 36.7879%
Time by which half are completed10ln⁡2≈8.32554610\sqrt{\ln 2}\approx8.325546 minutes
Expected completion time5π≈8.8622695\sqrt{\pi}\approx8.862269 minutes
Variance of completion time100(1−π/4)≈21.460184 min2100(1-\pi/4)\approx21.460184\text{ min}^2
Hazard at ten minutes0.2 min−10.2\text{ min}^{-1}
For every positive shape, S(η)=e−1S(\eta)=e^{-1}. Thus scale is the time by which approximately 63.2121% of cases are completed, not generally the mean or median. Nor is it an arrival intensity: ten minutes is a duration, while intensity is measured per minute.
Weibull survival and hazard curves for shapes 0.5, 1 and 2, all with scale ten minutes. All survival curves pass through 0.367879 at ten minutes; the shape-two hazard there is 0.2 per minute. The shape-0.5 hazard diverges near zero.
Weibull survival and hazard curves for shapes 0.5, 1 and 2, all with scale ten minutes. All survival curves pass through 0.367879 at ten minutes; the shape-two hazard there is 0.2 per minute. The shape-0.5 hazard diverges near zero.
Original synthetic comparison with zero location. Changing shape changes duration dependence while preserving the survival probability at the shared scale. The horizontal window ends at 25 minutes; the distributions continue beyond it.

Waiting changes the next question

Once a reconciliation is still unresolved after ten minutes, the relevant conditional probability of completing in the next five minutes is
P(10<T≤15∣T>10)=1−S(15)S(10)=1−e−1.25≈0.713495.\begin{aligned} \mathbb{P}(10<T\le15\mid T>10) &=1-\frac{S(15)}{S(10)}\\ &=1-e^{-1.25}\approx0.713495. \end{aligned}
That is 71.3495%, compared with 22.1199% for completion during the first five minutes of a fresh case, F(5)=1−e−0.25F(5)=1-e^{-0.25}. Restarting the clock at ten minutes and applying the original five-minute CDF would discard the model's increasing hazard. Only shape one has the exponential memoryless property.

This conditioning uses only the fact that the case is unresolved. If new information reveals its size, complexity or workflow stage, a single unconditional Weibull model need not describe the revised outlook. If observation ends before completion, the elapsed duration is a censored observation, not an observed completion time.

Quantiles, moments and simulation

Solving F(t)=pF(t)=p gives the quantile function:
Q(p)=η[−ln⁡(1−p)]1/k,0<p<1.Q(p)=\eta[-\ln(1-p)]^{1/k},\qquad0<p<1.

In particular, the median is η(ln⁡2)1/k\eta(\ln2)^{1/k}. The general moments are

E[T]=ηΓ(1+1k),Var(T)=η2[Γ(1+2k)−Γ(1+1k)2].\begin{aligned} \mathbb{E}[T]&=\eta\Gamma\left(1+\frac1k\right),\\ \text{Var}(T)&=\eta^2\left[\Gamma\left(1+\frac2k\right)-\Gamma\left(1+\frac1k\right)^2\right]. \end{aligned}
Here Γ(a)=∫0∞ua−1e−u du\Gamma(a)=\int_0^\infty u^{a-1}e^{-u}\,du is the Gamma function; it is distinct from the Gamma distribution. All positive-order moments exist for positive shape and scale.
A useful construction explains the family: if EE is exponential with rate one, then T=ηE1/kT=\eta E^{1/k} is Weibull. Equivalently, inverse transform sampling uses T=Q(U)T=Q(U) for a uniform draw UU in (0,1)(0,1). The Random Number Generation lesson develops inversion and the unit-scale Weibull recipe; its full content requires membership and has a public preview. Generating independent draws creates an independent-duration model because of the sampling construction, not because the family name guarantees independence.

Keep the event and parameterisation explicit

Some references use S(t)=e−λtkS(t)=e^{-\lambda t^k}. This is the same unshifted family with λ=η−k\lambda=\eta^{-k}; λ\lambda has units of time to the power −k-k and is a constant hazard only when k=1k=1. Other implementations add a location parameter, which changes the earliest possible event time. Record the actual survival formula alongside parameter estimates.

In credit applications, reduced-form default models connect hazard to default survival and distinguish physical probabilities from pricing probabilities. That full lesson requires membership and has a public preview. An operational completion model and a risk-neutral default model need different data, event definitions and interpretations even when their algebra matches.
A further financial use is duration-based backtesting of Value at Risk. Christoffersen and Pelletier use a Weibull alternative to an exponential no-hit-duration benchmark. A shape below one gives a decreasing hazard, accommodating excess very short and very long durations associated with violation clustering. This concerns the event process implied by durations; a Weibull marginal alone does not establish serial dependence between successive duration observations. Daily hit durations are discrete, so a continuous-duration approximation and its test calibration need care.
Finally, the bounded reverse-Weibull extreme-value family encountered in block-maxima theory is a different parameterisation and use of the name. The ordinary waiting-time Weibull defined here has no finite upper endpoint.

Sources

Last updated
  • continuous distributions
  • waiting times
  • survival
  • hazard rates
  • shape scale
← Back to the glossaryBack to top ↑