Efficient sampling from phase-type distributions

A phase-type distributed random variable represents the time to absorption of a Markov chain with an absorbing state. In this letter, we show that the alias method can be modified to efficiently generate phase-type distributed random variables. Both initialisation and generation are fast, at the cost of larger memory requirements. Numerical experiments show that the proposed method significantly reduces the computation time compared to direct simulation.