Misplaced Pages

Birnbaum–Saunders distribution

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.

The Birnbaum–Saunders distribution, also known as the fatigue life distribution, is a probability distribution used extensively in reliability applications to model failure times. There are several alternative formulations of this distribution in the literature. It is named after Z. W. Birnbaum and S. C. Saunders.

Theory

This distribution was developed to model failures due to cracks. A material is placed under repeated cycles of stress. The j cycle leads to an increase in the crack by Xj amount. The sum of the Xj is assumed to be normally distributed with mean and variance . The probability that the crack does not exceed a critical length ω is

P ( X ω ) = Φ ( ω n μ σ n ) {\displaystyle P(X\leq \omega )=\Phi \left({\frac {\omega -n\mu }{\sigma {\sqrt {n}}}}\right)}

where Φ() is the cdf of normal distribution.

If T is the number of cycles to failure then the cumulative distribution function (cdf) of T is

P ( T t ) = 1 Φ ( ω t μ σ t ) = Φ ( t μ ω σ t ) = Φ ( μ t σ ω σ t ) = Φ ( μ ω σ [ ( t ω / μ ) 0.5 ( ω / μ t ) 0.5 ] ) {\displaystyle P(T\leq t)=1-\Phi \left({\frac {\omega -t\mu }{\sigma {\sqrt {t}}}}\right)=\Phi \left({\frac {t\mu -\omega }{\sigma {\sqrt {t}}}}\right)=\Phi \left({\frac {\mu {\sqrt {t}}}{\sigma }}-{\frac {\omega }{\sigma {\sqrt {t}}}}\right)=\Phi \left({\frac {\sqrt {\mu \omega }}{\sigma }}\left\right)}

The more usual form of this distribution is:

F ( x ; α , β ) = Φ ( 1 α [ ( x β ) 0.5 ( β x ) 0.5 ] ) {\displaystyle F(x;\alpha ,\beta )=\Phi \left({\frac {1}{\alpha }}\left\right)}

Here α is the shape parameter and β is the scale parameter.

Properties

The Birnbaum–Saunders distribution is unimodal with a median of β.

The mean (μ), variance (σ), skewness (γ) and kurtosis (κ) are as follows:

μ = β ( 1 + α 2 2 ) {\displaystyle \mu =\beta \left(1+{\frac {\alpha ^{2}}{2}}\right)}
σ 2 = ( α β ) 2 ( 1 + 5 α 2 4 ) {\displaystyle \sigma ^{2}=(\alpha \beta )^{2}\left(1+{\frac {5\alpha ^{2}}{4}}\right)}
γ = 4 α ( 11 α 2 + 6 ) ( 5 α 2 + 4 ) 3 2 {\displaystyle \gamma ={\frac {4\alpha (11\alpha ^{2}+6)}{(5\alpha ^{2}+4)^{\frac {3}{2}}}}}
κ = 3 + 6 α 2 ( 93 α 2 + 40 ) ( 5 α 2 + 4 ) 2 {\displaystyle \kappa =3+{\frac {6\alpha ^{2}(93\alpha ^{2}+40)}{(5\alpha ^{2}+4)^{2}}}}

Given a data set that is thought to be Birnbaum–Saunders distributed the parameters' values are best estimated by maximum likelihood.

If T is Birnbaum–Saunders distributed with parameters α and β then T is also Birnbaum-Saunders distributed with parameters α and β.

Transformation

Let T be a Birnbaum-Saunders distributed variate with parameters α and β. A useful transformation of T is

X = 1 2 [ ( T β ) 0.5 ( T β ) 0.5 ] {\displaystyle X={\frac {1}{2}}\left} .

Equivalently

T = β ( 1 + 2 X 2 + 2 X ( 1 + X 2 ) 0.5 ) {\displaystyle T=\beta \left(1+2X^{2}+2X(1+X^{2})^{0.5}\right)} .

X is then distributed normally with a mean of zero and a variance of α / 4.

Probability density function

The general formula for the probability density function (pdf) is

f ( x ) = x μ β + β x μ 2 γ ( x μ ) ϕ ( x μ β β x μ γ ) x > μ ; γ , β > 0 {\displaystyle f(x)={\frac {{\sqrt {\frac {x-\mu }{\beta }}}+{\sqrt {\frac {\beta }{x-\mu }}}}{2\gamma \left(x-\mu \right)}}\phi \left({\frac {{\sqrt {\frac {x-\mu }{\beta }}}-{\sqrt {\frac {\beta }{x-\mu }}}}{\gamma }}\right)\quad x>\mu ;\gamma ,\beta >0}

where γ is the shape parameter, μ is the location parameter, β is the scale parameter, and ϕ {\displaystyle \phi } is the probability density function of the standard normal distribution.

Standard fatigue life distribution

The case where μ = 0 and β = 1 is called the standard fatigue life distribution. The pdf for the standard fatigue life distribution reduces to

f ( x ) = x + 1 x 2 γ x ϕ ( x 1 x γ ) x > 0 ; γ > 0 {\displaystyle f(x)={\frac {{\sqrt {x}}+{\sqrt {\frac {1}{x}}}}{2\gamma x}}\phi \left({\frac {{\sqrt {x}}-{\sqrt {\frac {1}{x}}}}{\gamma }}\right)\quad x>0;\gamma >0}

Since the general form of probability functions can be expressed in terms of the standard distribution, all of the subsequent formulas are given for the standard form of the function.

Cumulative distribution function

The formula for the cumulative distribution function is

F ( x ) = Φ ( x 1 x γ ) x > 0 ; γ > 0 {\displaystyle F(x)=\Phi \left({\frac {{\sqrt {x}}-{\sqrt {\frac {1}{x}}}}{\gamma }}\right)\quad x>0;\gamma >0}

where Φ is the cumulative distribution function of the standard normal distribution.

Quantile function

The formula for the quantile function is

G ( p ) = 1 4 [ γ Φ 1 ( p ) + 4 + ( γ Φ 1 ( p ) ) 2 ] 2 {\displaystyle G(p)={\frac {1}{4}}\left^{2}}

where Φ is the quantile function of the standard normal distribution.

References

External links

Public Domain This article incorporates public domain material from the National Institute of Standards and Technology

Probability distributions (list)
Discrete
univariate
with finite
support
with infinite
support
Continuous
univariate
supported on a
bounded interval
supported on a
semi-infinite
interval
supported
on the whole
real line
with support
whose type varies
Mixed
univariate
continuous-
discrete
Multivariate
(joint)
Directional
Univariate (circular) directional
Circular uniform
Univariate von Mises
Wrapped normal
Wrapped Cauchy
Wrapped exponential
Wrapped asymmetric Laplace
Wrapped Lévy
Bivariate (spherical)
Kent
Bivariate (toroidal)
Bivariate von Mises
Multivariate
von Mises–Fisher
Bingham
Degenerate
and singular
Degenerate
Dirac delta function
Singular
Cantor
Families
Category: