We consider residue expansions for survival and density/mass functions of first-passage distributions in finite-state semi-Markov processes (SMPs) in continuous and integer time. Conditions are given which guarantee that the residue expansions for these functions have a dominant exponential/geometric term. The key condition assumes that the relevant states for first passage contain an irreducible class, thus ensuring the same sort of dominant exponential/geometric terms as one gets for phase-type distributions in Markov processes. Essentially, the presence of an irreducible class along with some other conditions ensures that the boundary singularity b for the moment generating function (MGF) of the first-passage-time distribution is a simple pole. In the continuous-time setting we prove that b is a dominant pole, in that the MGF has no other pole on the vertical line $\{\text{Re}(s)=b\}.$ In integer time we prove that b is dominant if all holding-time mass functions for the SMP are aperiodic and non-degenerate. The expansions and pole characterisations address first passage to a single new state or a subset of new states, and first return to the starting state. Numerical examples demonstrate that the residue expansions are considerably more accurate than saddlepoint approximations and can provide a substitute for exact computation above the 75th percentile.

Transform inversions, in which density and survival functions are computed from their associated moment generating function $\mathcal{M}$, have largely been based on methods which use values of $\mathcal{M}$ in its convergence region. Prominent among such methods are saddlepoint approximations and Fourier-series inversion methods, including the fast Fourier transform. In this paper we propose inversion methods which make use of values for $\mathcal{M}$ which lie outside of its convergence region and in its analytic continuation. We focus on the simplest and perhaps richest setting for applications in which $\mathcal{M}$ is either a meromorphic function in its analytic continuation, so that all of its singularities are poles, or else the singularities are isolated essential. Asymptotic expansions of finite- and infinite-orders are developed for density and survival functions using the poles of $\mathcal{M}$ in its analytic continuation. For finite-order expansions, the expansion error is a contour integral in the analytic continuation, which we approximate using the saddlepoint method based on following the path of steepest descent. Such saddlepoint error approximations accurately determine expansion errors and, thus, provide the means for determining the order of the expansion needed to achieve some preset accuracy. They also provide an additive correction term which increases accuracy of the expansion. Further accuracy is achieved by computing the expansion errors numerically using a contour path which ultimately tracks the steepest descent direction. Important applications include Wilks’ likelihood ratio test in MANOVA, compound distributions, and the Sparre Andersen and Cramér–Lundberg ruin models.

Let (X1,...,Xn) be multivariate normal, with mean vector 𝛍 and covariance matrix 𝚺, and let Sn=eX1+⋯+eXn. The Laplace transform ℒ(θ)=𝔼e-θSn∝∫exp{-hθ(𝒙)}d𝒙 is represented as ℒ̃(θ)I(θ), where ℒ̃(θ) is given in closed form and I(θ) is the error factor (≈1). We obtain ℒ̃(θ) by replacing hθ(𝒙) with a second-order Taylor expansion around its minimiser 𝒙*. An algorithm for calculating the asymptotic expansion of 𝒙* is presented, and it is shown that I(θ)→ 1 as θ→∞. A variety of numerical methods for evaluating I(θ) is discussed, including Monte Carlo with importance sampling and quasi-Monte Carlo. Numerical examples (including Laplace-transform inversion for the density of Sn) are also given.

We consider a multi-type branching random walk on d-dimensional Euclidian space. The~uniform convergence, as n goes to infinity, of a scaled version of the Laplace transform of the point process given by the nth generation particles of each type is obtained. Similar results in the one-type case, where the transform gives a martingale, were obtained in Biggins (1992) and Barral (2001). This uniform convergence of transforms is then used to obtain limit results for numbers in the underlying point processes. Supporting results, which are of interest in their own right, are obtained on (i) ‘Perron-Frobenius theory’ for matrices that are smooth functions of a variable λ∈L and are nonnegative when λ∈L−⊂L, where L is an open set in ℂd, and (ii) saddlepoint approximations of multivariate distributions. The saddlepoint approximations developed are strong enough to give a refined large deviation theorem of Chaganty and Sethuraman (1993) as a by-product.

Prediction of future values of a shot noise process observed on a discrete lattice of points is considered. The shot magnitudes are assumed to be independent and identically distributed and to arrive via a Poisson process; the effect of each shot dissipates and/or accumulates according to a known shot function. Conditional mean and linear point predictors of future process values are developed. Distributional prediction, obtained through saddlepoint approximation of the conditional distributions of the process, is also explored.

We consider a class of stationary infinite-order moving average processes with margins in the class of infinitely divisible exponential dispersion models. The processes are constructed by means of the thinning operation of Joe (1996), generalizing the binomial thinning used by McKenzie (1986, 1988) and Al-Osh and Alzaid (1987) for integer-valued time series. As a special case we obtain a class of autoregressive moving average processes that are different from the ARMA models proposed by Joe (1996). The range of possible marginal distributions for the new models is extensive and includes all infinitely divisible distributions with finite moment generating functions, hereunder many known discrete, continuous and mixed distributions.

Shot noise processes form an important class of stochastic processes modeling phenomena which occur as shocks to a system and with effects that diminish over time. In this paper we present extreme value results for two cases — a homogeneous Poisson process of shocks and a non-homogeneous Poisson process with periodic intensity function. Shocks occur with a random amplitude having either a gamma or Weibull density and dissipate via a compactly supported impulse response function. This work continues work of Hsing and Teugels (1989) and Doney and O'Brien (1991) to the case of random amplitudes.

In each of a large number N of independent cells a breakdown mechanism is under way and proceeds until the first of the cells actually fails. At such a time, in each cell, the situation reverts to some initial state and the mechanism restarts. In this paper we consider those mechanisms for which breakdown may be modelled as the explosion of a pure birth process. Of interest is the distribution of time between failures and the possibility of estimating N and/or model parameters by observing a sequence of failure times. Saddlepoint approximation methods are used in the relevant extreme-value theory analysis for two important cases.

Saddlepoint approximations to the probabilities for a general time-homogeneous birth process are derived from the known Laplace transform with respect to time. Numerical computations for the simple (logistic) epidemic show the approximation to have a remarkably low relative error over the useful range of times for populations as small as 5. The merit or otherwise of renormalising the approximate probabilities is discussed. To study the accuracy for moderately large population sizes, a special form of birth rate is used, for which the exact probabilities can be computed directly.

It is known that, in the deterministic model of a one-dimensional spatial linear birth process, a finite initial distribution of individuals develops into an advancing wave with the minimum velocity possible for such waves. In this paper, the spatial covariance density of the process is studied in the vicinity of a point moving outwards with a given velocity. Its behaviour casts doubt on the applicability of the deterministic result to individual realisations of the process.