class

Weibull

extendsDistribution

Weibull(scale: Tensor | float, concentration: Tensor | float, validate_args: bool | None = None)

source edit

Weibull distribution — a flexible family for lifetime and survival analysis.

Weibull(scale=λ, concentration=k) generalises both the Exponential ( $k = 1$ ) and Rayleigh ( $k = 2$ ) distributions. It is widely used for modelling time-to-failure data because the hazard rate $h(t) = (k/\lambda)(t/\lambda)^{k-1}$ can increase, be constant, or decrease depending on $k$ .

Parameters

scaleTensor | float

Scale parameter

\lambda > 0

— a characteristic life value at which the CDF is

1 - e^{-1} \approx 63.2\%

concentrationTensor | float

Shape parameter

k > 0

. Values

k < 1

give a decreasing hazard (infant mortality);

k = 1

is constant (memoryless Exponential);

k > 1

gives an increasing hazard (wear-out).

validate_argsbool | None= None

If True, validate parameter constraints at construction time.

Attributes

scaleTensor

Scale parameter

\lambda

concentrationTensor

Shape parameter

k

Notes

PDF:

p(x; \lambda, k) = \frac{k}{\lambda} \left(\frac{x}{\lambda}\right)^{k-1} \exp\!\left(-\left(\frac{x}{\lambda}\right)^k\right), \quad x \geq 0

Log-PDF:

\log p(x) = \log k - k\log\lambda + (k-1)\log x - (x/\lambda)^k

Moments:

Mean: $E[X] = \lambda \,\Gamma(1 + 1/k)$
Variance: $\operatorname{Var}[X] = \lambda^2 \left[\Gamma(1+2/k) - \Gamma(1+1/k)^2\right]$

Entropy:

H[X] = \gamma\left(1 - \tfrac{1}{k}\right) + \log(\lambda/k) + 1

where $\gamma \approx 0.5772$ is the Euler–Mascheroni constant.

Reparameterised sampling uses the inverse-CDF: $X = \lambda(-\log(1-U))^{1/k}$ for $U \sim \operatorname{Uniform}(0,1)$ .

Examples

>>> import lucid
>>> from lucid.distributions import Weibull
>>> # Exponential(rate=1) as a special case
>>> dist_exp = Weibull(scale=1.0, concentration=1.0)
>>> dist = Weibull(scale=2.0, concentration=1.5)
>>> samples = dist.rsample((500,))

Used by 1

lucid.distributions

Constructors

dunder

init

→None

__init__(scale: Tensor | float, concentration: Tensor | float, validate_args: bool | None = None)

source edit

Initialise a Weibull distribution.

Parameters

scaleTensor | float

Scale parameter

\lambda > 0

concentrationTensor | float

Shape parameter

k > 0

. Values

k < 1

give a decreasing hazard;

k = 1

gives the Exponential;

k > 1

gives an increasing hazard.

validate_argsbool | None= None

If True, validate parameter constraints at construction time.

Properties

prop

mean

→Tensor

mean: Tensor

source edit

Mean of the Weibull distribution: $E[X] = \lambda \Gamma(1 + 1/k)$ .

Returns

Tensor

Mean values of shape batch_shape.

prop

support

→Constraint

support: Constraint

source edit

Support of the Weibull distribution: $[0, \infty)$ .

Returns

Constraint

The nonnegative constraint.

prop

variance

→Tensor

variance: Tensor

source edit

Variance of the Weibull distribution.

$\operatorname{Var}[X] = \lambda^2 [\Gamma(1+2/k) - \Gamma(1+1/k)^2]$

Returns

Tensor

Variance values of shape batch_shape.

Instance methods

entropy

→Tensor

entropy()

source edit

Entropy of the Weibull distribution.

H[X] = \gamma(1 - 1/k) + \log(\lambda/k) + 1

where $\gamma \approx 0.5772$ is the Euler-Mascheroni constant.

Returns

Tensor

Entropy values of shape batch_shape (nats).

log_prob

→Tensor

log_prob(value: Tensor)

source edit

Log-probability density of the Weibull distribution.

\log p(x) = \log k - k\log\lambda + (k-1)\log x - (x/\lambda)^k

Parameters

valueTensor

Non-negative points

x \geq 0

at which to evaluate.

Returns

Tensor

Log-density values of the same shape as value.

rsample

→Tensor

rsample(sample_shape: tuple[int, ...] = ())

source edit

Draw reparameterised samples via the inverse-CDF.

Uses $X = \lambda (-\log(1-U))^{1/k}$ for $U \sim \operatorname{Uniform}(0, 1)$ . Gradients propagate through both scale and concentration.

Parameters

sample_shapetuple[int, ...]= ()

Leading shape of the output sample batch.

Returns

Tensor

Non-negative samples of shape (*sample_shape, *batch_shape).

>>> import lucid >>> from lucid.distributions import Weibull >>> # Exponential(rate=1) as a special case >>> dist_exp = Weibull(scale=1.0, concentration=1.0) >>> dist = Weibull(scale=2.0, concentration=1.5) >>> samples = dist.rsample((500,))