Let f be the density function associated to a matrix-exponential distribution of parameters $(\boldsymbol{\alpha}, T,\boldsymbol{{s}})$ . By exponentially tilting f, we find a probabilistic interpretation which generalizes the one associated to phase-type distributions. More specifically, we show that for any sufficiently large $\lambda\ge 0$ , the function $x\mapsto \left(\int_0^\infty e^{-\lambda s}f(s)\textrm{d} s\right)^{-1}e^{-\lambda x}f(x)$ can be described in terms of a finite-state Markov jump process whose generator is tied to T. Finally, we show how to revert the exponential tilting in order to assign a probabilistic interpretation to f itself.

The paper discusses the risk of ruin in insurance coverage of an epidemic in a closed population. The model studied is an extended susceptible–infective–removed (SIR) epidemic model built by Lefèvre and Simon (Methodology Comput. Appl. Prob.22, 2020) as a block-structured Markov process. A fluid component is then introduced to describe the premium amounts received and the care costs reimbursed by the insurance. Our interest is in the risk of collapse of the corresponding reserves of the company. The use of matrix-analytic methods allows us to determine the distribution of ruin time, the probability of ruin, and the final amount of reserves. The case where the reserves are subjected to a Brownian noise is also studied. Finally, some of the results obtained are illustrated for two particular standard SIR epidemic models.

In this paper we consider a class of quasi-birth-and-death processes for which explicit solutions can be obtained for the rate matrix R and the associated matrix G. The probabilistic interpretations of these matrices allow us to describe their elements in terms of paths on the two-dimensional lattice. Then determining explicit expressions for the matrices becomes equivalent to solving a lattice path counting problem, the solution of which is derived using path decomposition, Bernoulli excursions, and hypergeometric functions. A few applications are provided, including classical models for which we obtain some new results.

In this paper we introduce certain Hankel matrices that can be used to study ME (matrix exponential) distributions, in particular to compute their ME orders. The Hankel matrices for a given ME probability distribution can be constructed if one of the following five types of information about the distribution is available: (i) an ME representation, (ii) its moments, (iii) the derivatives of its distribution function, (iv) its Laplace-Stieltjes transform, or (v) its distribution function. Using the Hankel matrices, a necessary and sufficient condition for a probability distribution to be an ME distribution is found and a method of computing the ME order of the ME distribution developed. Implications for the PH (phase-type) order of PH distributions are examined. The relationship between the ME order, the PH order, and a lower bound on the PH order given by Aldous and Shepp (1987) is discussed in numerical examples.

In this paper, we study the classification of matrix GI/M/1-type Markov chains with a tree structure. We show that the Perron–Frobenius eigenvalue of a Jacobian matrix provides information for classifying these Markov chains. A fixed-point approach is utilized. A queueing application is presented to show the usefulness of the classification method developed in this paper.

The isotropic planar point processes of phase-type are natural generalizations of the Poisson process on the plane. On the one hand, those processes are isotropic and stationary for the mean count, as in the case of the Poisson process. On the other hand, they exhibit dependence of counts in disjoint sets. In a recent paper, we have proved that the number of points in a square window has a Poisson distribution asymptotically as the window is located far away from the origin of the process. We extend our work to the case of a window of arbitrary shape.

When analyzing the equilibrium behavior of M/G/1 type Markov chains by transform methods, restrictive hypotheses are often made to avoid technical problems that arise in applying results from complex analysis and linear algebra. It is shown that such restrictive assumptions are unnecessary, and an analysis of these chains using generating functions is given under only the natural hypotheses that first moments (or second moments in the null recurrent case) exist. The key to the analysis is the identification of an important subspace of the space of bounded solutions of the system of homogeneous vector-valued Wiener–Hopf equations associated with the chain. In particular, the linear equations in the boundary probabilities obtained from the transform method are shown to correspond to a spectral basis of the shift operator on this subspace. Necessary and sufficient conditions under which the chain is ergodic, null recurrent or transient are derived in terms of properties of the matrix-valued generating functions determined by transitions of the Markov chain. In the transient case, the Martin exit boundary is identified and shown to be associated with certain eigenvalues and vectors of one of these generating functions. An equilibrium analysis of the class of G/M/1 type Markov chains by similar methods is also presented.