Phase-type distribution

Phase-type distribution

Probability distribution
name =Phase-type
type =density
pdf_

cdf_

parameters = $S,; m imes m$ subgenerator matrix
$oldsymbol{alpha}$ , probability row vector
support = $x in [0; infty)!$
pdf = $oldsymbol{alpha}e^{xS}oldsymbol{S}^{0}$
See article for details
cdf = $1-oldsymbol{alpha}e^{xS}oldsymbol{1}$

mean = $-1oldsymbol{alpha}{S}^{-1}mathbf{1}$
mode =no simple closed form
variance = $2oldsymbol{alpha}{S}^{-2}mathbf{1}$
median =no simple closed form
skewness = $-6oldsymbol{alpha}{S}^{-3}mathbf{1}/sigma^{3}$
kurtosis = $24oldsymbol{alpha}{S}^{-4}mathbf{1}/sigma^{4}-3$
entropy =
mgf = $oldsymbol{alpha}(-tI-S)^{-1}oldsymbol{S}^{0}+alpha_{m+1}$
char = $oldsymbol{alpha}(itI-S)^{-1}oldsymbol{S}^{0}+alpha_{m+1}$

A phase-type distribution is a probability distribution that results from a system of one or more inter-related Poisson processes occurring in sequence, or phases. The sequence in which each of the phases occur may itself be a stochastic process. The distribution can be represented by a random variable describing the time until absorption of a Markov process with one absorbing state. Each of the states of the Markov process represents one of the phases.

It has a discrete time equivalent the discrete phase-type distribution.

The phase-type distribution is dense in the field of all positive-valued distributions, that is, it can be used to approximate any positive valued distribution.

Definition

There exists a continuous-time Markov process with "m"+1 states, where "m" ≥ 1. The states 1,...,"m" are transient states and state "m"+1 is an absorbing state. The process has an initial probability of starting in any of the "m"+1 phases given by the probability vector (α,α_"m"+1).

The continuous phase-type distibution is the distribution of time from the processes starting until absorption in the absorbing state.

This process can be written in the form of a transition rate matrix,

: ${Q}=left [egin{matrix}{S}&mathbf{S}^0\mathbf{0}&0end{matrix}
ight] ,$

where "S" is an "m"×"m" matrix and "S"⁰ = -S1. Here 1 represents an "m"×1 vector with every element being 1.

Characterization

The distribution of time "X" until the process reaches the absorbing state is said to be phase-type distributed and is denoted PH(α,"S").

The distribution function of "X" is given by,

: $F(x)=1-oldsymbol{alpha}exp({S}x)mathbf{1},$

and the density function,

: $f(x)=oldsymbol{alpha}exp({S}x)mathbf{S^{0,$

for all "x" > 0, where exp( · ) is the matrix exponential. It is usually assumed the probability of process starting in the absorbing state is zero. The moments of the distribution function are given by

: $E [X^{n}] =(-1)^{n}n!oldsymbol{alpha}{S}^{-n}mathbf{1}.$

pecial cases

The following probability distributions are all considered special cases of a continuous phase-type distribution:
* Degenerate distribution, point mass at zero or the empty phase-type distribution - 0 phases.
* Exponential distribution - 1 phase.
* Erlang distribution - 2 or more identical phases in sequence.
* Deterministic distribution (or constant) - The limiting case of an Erlang distribution, as the number of phases become infinite, while the time in each state becomes zero.
* Coxian distribution - 2 or more (not necessarily identical) phases in sequence, with a probability of transitioning to the terminating/absorbing state after each phase.
* Hyper-exponential distribution (also called a mixture of exponential) - 2 or more non-identical phases, that each have a probability of occurring in a mutually exclusive, or parallel, manner. (Note: The exponential distribution is the degenerate situation when all the parallel phases are identical.)
* Hypoexponential distribution - 2 or more phases in sequence, can be non-identical or a mixture of identical and non-identical phases, generalises the Erlang.As the phase-type distribution is dense in the field of all positive-valued distributions, we can represent any positive valued distribution. However, the phase-type is a light-tailed or platikurtic distribution. So the representation of heavy-tailed or leptokurtic distribution by phase type is an approximation, even if the precision of the approximation can be as good as we want.

Examples

In all the following examples it is assumed that there is no probability mass at zero, that is α_"m"+1 = 0.

Exponential distribution

The simplest non-trivial example of a phase-type distribution is the exponential distribution of parameter λ. The parameter of the phase-type distribution are : "S" = -λ and α = 1.

Hyper-exponential or mixture of exponential distribution

The mixture of exponential or hyper-exponential distribution with parameter (α₁,α₂,α₃,α₄,α₅) (such that $sum alpha_i =1$ and α_"i" > 0 for all "i") and (λ₁,λ₂,λ₃,λ₄,λ₅) can be represented as a phase type distribution with

: $oldsymbol{alpha}=(alpha_1,alpha_2,alpha_3,alpha_4,alpha_5),$

and

: ${S}=left [egin{matrix}-lambda_1&0&0&0&0\0&-lambda_2&0&0&0\0&0&-lambda_3&0&0\0&0&0&-lambda_4&0\0&0&0&0&-lambda_5\end{matrix}
ight] .$

The mixture of exponential can be characterized through its density

: $f(x)=sum_{i=1}^5 alpha_i lambda_i e^{-lambda_i x}$

or its distribution function

: $F(x)=1-sum_{i=1}^5 alpha_i e^{-lambda_i x}.$

This can be generalized to a mixture of "n" exponential distributions.

Erlang distribution

The Erlang distribution has two parameters, the shape an integer "k" > 0 and the rate λ > 0. This is sometimes denoted "E"("k",λ). The Erlang distribution can be written in the form of a phase-type distribution by making "S" a "k"×"k" matrix with diagonal elements -λ and super-diagonal elements λ, with the probability of starting in state 1 equal to 1. For example "E"(5,λ),

: $oldsymbol{alpha}=(1,0,0,0,0),$ and: ${S}=left [egin{matrix}-lambda&lambda&0&0&0\0&-lambda&lambda&0&0\0&0&-lambda&lambda&0\0&0&0&-lambda&lambda\0&0&0&0&-lambda\end{matrix}
ight] .$

The hypoexponential distribution is a generalisation of the Erlang distribution by having different rates for each transition (the non-homogeneous case).

Mixture of Erlang distribution

The mixture of two Erlang distribution with parameter "E"(3,β₁), "E"(3,β₂) and (α₁,α₂) (such that α₁ + α₂ = 1 and for each "i", α_"i" ≥ 0) can be represented as a phase type distribution with

: $oldsymbol{alpha}=(alpha_1,0,0,alpha_2,0,0),$

and

: ${S}=left [egin{matrix}-eta_1&eta_1&0&0&0&0\ 0&-eta_1&eta_1&0&0&0\0&0&-eta_1&0&0&0\ 0&0&0&-eta_2&eta_2&0\0&0&0&0&-eta_2&eta_2\0&0&0&0&0&-eta_2\end{matrix}
ight] .$

Coxian distribution

The Coxian distribution is a generalisation of the hypoexponential. Instead of only being able to enter the absorbing state from state "k" it can be reached from any phase. The phase-type representation is given by,

: $S=left [egin{matrix}-lambda_{1}&p_{1}lambda_{1}&0&dots&0&0\ 0&-lambda_{2}&p_{2}lambda_{2}&ddots&0&0\ vdots&ddots&ddots&ddots&ddots&vdots\ 0&0&ddots&-lambda_{k-2}&p_{k-2}lambda_{k-2}&0\ 0&0&dots&0&-lambda_{k-1}&p_{k-1}lambda_{k-1}\ 0&0&dots&0&0&-lambda_{k}end{matrix}
ight]$

and

: $oldsymbol{alpha}=(1,0,dots,0),$

where 0 < "p"₁,...,"p"_"k"-1 ≤ 1, in the case where all "p"_"i" = 1 we have the hypoexponential distribution. The Coxian distribution is extremely important as any acyclic phase-type distribution has an equivalent Coxian representation.

The generalised Coxian distribution relaxes the condition that requires starting in the first phase.

ee also

* Discrete phase-type distribution
* Continuous-time Markov process
* Exponential distribution
* Hyper-exponential distribution
* Queueing model
* Queuing theory

References

* M. F. Neuts. Matrix-Geometric Solutions in Stochastic Models: an Algorthmic Approach, Chapter 2: Probability Distributions of Phase Type; Dover Publications Inc., 1981.
* G. Latouche, V. Ramaswami. Introduction to Matrix Analytic Methods in Stochastic Modelling, 1st edition. Chapter 2: PH Distributions; ASA SIAM, 1999.
* C. A. O'Cinneide (1990). "Characterization of phase-type distributions". Communications in Statistics: Stocahstic Models, 6(1), 1-57.
* C. A. O'Cinneide (1999). "Phase-type distribution: open problems and a few properties", Communication in Statistic: Stochastic Models, 15(4), 731-757.

Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать реферат

Look at other dictionaries:

Discrete phase-type distribution — The discrete phase type distribution is a probability distribution that results from a system of one or more inter related geometric distributions occurring in sequence, or phases. The sequence in which each of the phases occur may itself be a… … Wikipedia
Distribution transformer — Single phase distribution transformer in Canada A distribution transformer is a transformer that provides the final voltage transformation in the electric power distribution system, stepping down the voltage used in the distribution lines to the… … Wikipedia
Hypoexponential distribution — Probability distribution name =Hypoexponential type =density pdf cdf parameters =lambda {1},dots,lambda {k} > 0, rates (real) support =x in [0; infty)! pdf =Expressed as a phase type distribution oldsymbol{alpha}e^{xTheta}Thetaoldsymbol{1} Has… … Wikipedia
Negative binomial distribution — Probability mass function The orange line represents the mean, which is equal to 10 in each of these plots; the green line shows the standard deviation. notation: parameters: r > 0 number of failures until the experiment is stopped (integer,… … Wikipedia
DISTRIBUTION (économie) — Pour la Chambre de commerce internationale, la distribution est «le stade qui suit celui de la production des biens à partir du moment où ils sont commercialisés jusqu’à leur prise en possession par le consommateur ou l’utilisateur final. Elle… … Encyclopédie Universelle
Erlang distribution — Probability distribution name =Erlang type =density pdf cdf parameters =k > 0 in mathbb{Z} shape lambda > 0, rate (real) alt.: heta = 1/lambda > 0, scale (real) support =x in [0; infty)! pdf =frac{lambda^k x^{k 1} e^{ lambda x{(k 1)!,} cdf… … Wikipedia
Phase terminale (telefilm) — Phase terminale (téléfilm) Pour les articles homonymes, voir Phase terminale. Téléfilms 0 9 A B C D E F G H I J K L M N O P Q R S T U V W X Y Z par date de sortie nationalité période genre type Phase terminale … Wikipédia en Français
Phase contrast microscopy — Phase contrast image of a cheek epithelial cell Phase contrast microscopy is an optical microscopy illumination technique of great importance to biologists in which small (invisible to the human eye) phase shifts in the light passing through a… … Wikipedia
Phase dispersion minimization — (PDM) is a data analysis technique that searches for periodic components of a time series data set. It is useful for data sets with gaps, non sinusoidal variations, poor time coverage or other problems that would make Fourier techniques unusable … Wikipedia
Hyper-exponential distribution — In probability theory, a hyper exponential distribution is a continuous distribution such that the probability density function of the random variable X is given by: f X(x) = sum {i=1}^n f {Y i}(y) p i,where Y i is an exponentially distributed… … Wikipedia

Mark and share
Search through all dictionaries
Translate…
Search Internet

Share the article and excerpts