Phase-type distribution

From Wikipedia, the free encyclopedia

A phase-type distribution is a probability distribution that results from a system 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 Chain with one absorbing state. Each of the states of the Markov Chain represent one of the phases.

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

Contents

Assume there exists a continuous-time Markov process with m + 1 states, where m\geq1, 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 (\vec\boldsymbol{\alpha},\alpha_{m+1}).

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

\mathbf{Q}=\left[\begin{matrix}\mathbf{S}&\vec\mathbf{S}^0\\\vec\mathbf{0}&0\end{matrix}\right],

where \mathbf{S} is a m\times m matrix and \vec\mathbf{S}^0=-\mathbf{S}\vec\mathbf{1}. Here \vec\mathbf{1} represents an m\times 1 vector with every element being 1.

The distribution of time X until the process reaches the absorbing state is said to be phase-type distributed and is denoted PH(\vec\boldsymbol{\alpha},\mathbf{S}).

The distribution function of X is given by,

F(x)=1-\vec\boldsymbol{\alpha}\exp(\mathbf{S}x)\vec\mathbf{1},

and the density function,

f(x)=\vec\boldsymbol{\alpha}\exp(\mathbf{S}x)\vec\mathbf{S^{0}},

for all x > 0. If it is 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!\vec\boldsymbol{\alpha}\mathbf{S}^{-n}\vec\mathbf{1}.

The following probability distributions are all considered special cases of a continuous phase-type distribution:

  • 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 - 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.)

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 \mathbf{S} a k\times k matrix with diagonal elements − λ and super-diagonal elements λ, with the probability of starting in state 1 equal to 1. For example E(5,λ),

\vec\boldsymbol{\alpha}=(1,0,0,0,0),

and

\mathbf{S}=\left[\begin{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}\right].

The discrete phase-type distribution is defined similar to the continuous time.

Assume there exists a discrete-time Markov chain with m + 1 states, where m\geq1, 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 (\vec\boldsymbol{\tau},\tau_{m+1}).

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

\mathbf{P}=\left[\begin{matrix}\mathbf{T}&\vec\mathbf{T}^0\\\vec\mathbf{0}&1\end{matrix}\right],

where \mathbf{T} is a m\times m matrix and \vec\mathbf{T}^0+\mathbf{T}\vec\mathbf{1}=\vec\mathbf{1}.

The distribution of the number of steps K until the process reaches the absorbing state is said to be discretely phase-type distributed and is denoted PH_{d}(\vec\boldsymbol{\tau},\mathbf{T}).

The distribution function of K is given by,

F(k)=1-\vec\boldsymbol{\tau}\mathbf{T}^{k}\vec\mathbf{1},

for k = 0,1,2,..., and the density function,

f(k)=\vec\boldsymbol{\tau}\mathbf{T}^{k-1}\vec\mathbf{T^{0}},

for k = 1,2,.... If it is assumed the probability of process starting in the absorbing state is zero, the factorial moments of the distribution function are given by,

E[K(K-1)...(K-n+1)]=n!\vec\boldsymbol{\tau}(I-\mathbf{T})^{-n}\mathbf{T}^{n-1}\vec\mathbf{1},

where I is the appropriate dimension identity matrix.

Just as the continuous time distribution is a generalisation of exponential, the discrete time is a genearlisation of the geometric distribution, for example:

  • Sheldon M Ross. Introduction to Probability Models, 8th edition. Chapter 6: Continuous Time Markov Chains, Academic Press; December, 2002
  • G. Latouche, V. Ramaswami. Introduction to Matrix Analytic Methods in Stochastic Modelling, 1st edition. Chapter 2: PH Distributions; ASA SIAM, 1999.
Image:Bvn-small.png Probability distributionsview  talk  edit ]
Univariate Multivariate
Discrete: BenfordBernoullibinomialBoltzmanncategoricalcompound PoissondegenerateGauss-Kuzmingeometrichypergeometriclogarithmicnegative binomialparabolic fractalPoissonRademacherSkellamuniformYule-SimonzetaZipfZipf-Mandelbrot Ewensmultinomialmultivariate Polya
Continuous: BetaBeta primeCauchychi-squareDirac delta functionErlangexponentialexponential powerFfadingFisher's zFisher-TippettGammageneralized extreme valuegeneralized hyperbolicgeneralized inverse GaussianHalf-LogisticHotelling's T-squarehyperbolic secanthyper-exponentialhypoexponentialinverse chi-square (scaled inverse chi-square)• inverse Gaussianinverse gamma (scaled inverse gamma) • KumaraswamyLandauLaplaceLévyLévy skew alpha-stablelogisticlog-normalMaxwell-BoltzmannMaxwell speednormal (Gaussian)normal inverse GaussianParetoPearsonpolarraised cosineRayleighrelativistic Breit-WignerRiceshifted GompertzStudent's ttriangulartype-1 Gumbeltype-2 GumbeluniformVariance-GammaVoigtvon MisesWeibullWigner semicircleWilks' lambda Dirichletinverse-WishartKentmatrix normalmultivariate normalmultivariate Studentvon Mises-FisherWigner quasiWishart
Miscellaneous: Cantorconditionalexponential familyinfinitely divisiblelocation-scale familymarginalmaximum entropyphase-typeposteriorpriorquasisamplingsingular
Advanced Search
Included Web Search Engines


Safe Search

close

Top Matching Results

Occasionally Search.com will highlight specialized results that are based on the context of your query. Examples of specialized results include specific links to news, images, or video.

Top Matching Results may highlight information from other Search.com pages, content from the CNET Network of sites, or third party content. The listings are based purely on relevance. Search.com does not receive payment for listings in this section but our partners that provide this data may get paid for listing these products.

Sponsored Links

This section contains paid listings which have been purchased by companies that want to have their sites appear for specific search terms and related content. These listings are administered, sorted and maintained by a third party and are not endorsed by Search.com.

Search Results

Search.com sends your search query to several search engines at one time and integrates the results into one list which has been sorted by relevance using Search.com's proprietary algorithm. You can customize the list of search engines included in your metasearch from the preferences.

The search engines that are used in your metasearch may allow companies to pay to have their Web sites included within the results. To view the Paid Inclusion policy for a specific search engine, please visit their Web site. Search.com does not accept payment or share revenue with any search engine partner for listings in this section.