We construct random dynamics for collections of nonintersecting planar contours, leaving invariant the distributions of length- and area-interacting polygonal Markov fields with V-shaped nodes. The first of these dynamics is based on the dynamic construction of consistent polygonal fields, as presented in the original articles by Arak (1983) and Arak and Surgailis (1989), (1991), and it provides an easy-to-implement Metropolis-type simulation algorithm. The second dynamics leads to a graphical construction in the spirit of Fernández et al. (1998), (2002) and yields a perfect simulation scheme in a finite window in the infinite-volume limit. This algorithm seems difficult to implement, yet its value lies in that it allows for theoretical analysis of the thermodynamic limit behaviour of length-interacting polygonal fields. The results thus obtained include, in the class of infinite-volume Gibbs measures without infinite contours, the uniqueness and exponential α-mixing of the thermodynamic limit of such fields in the low-temperature region. Outside this class, we conjecture the existence of an infinite number of extreme phases breaking both the translational and rotational symmetries.

Let (Xt)t≥0 be a continuous-time irreducible Markov chain on a finite state space E, let v: E→ℝ\{0}, and let (φt)t≥0 be defined by φt=∫0tv(Xs)d s. We consider the case in which the process (φt)t≥0 is oscillating and that in which (φt)t≥0 has a negative drift. In each of these cases, we condition the process (Xt,φt)t≥0 on the event that (φt)t≥0 hits level y before hitting 0 and prove weak convergence of the conditioned process as y→∞. In addition, we show the relationship between the conditioning of the process (φt)t≥0 with a negative drift to oscillate and the conditioning of it to stay nonnegative for a long time, and the relationship between the conditioning of (φt)t≥0 with a negative drift to drift to ∞ and the conditioning of it to hit large levels before hitting 0.

We consider risk processes that locally behave like Brownian motion with some drift and variance, these both depending on an underlying Markov chain that is also used to generate the claims arrival process. Thus, claims arrive according to a renewal process with waiting times of phase type. Positive claims (downward jumps) are always possible but negative claims (upward jumps) are also allowed. The claims are assumed to form an independent, identically distributed sequence, independent of everything else. As main results, the joint Laplace transform of the time to ruin and the undershoot at ruin, as well as the probability of ruin, are explicitly determined under the assumption that the Laplace transform of the positive claims is a rational function. Both the joint Laplace transform and the ruin probability are decomposed according to the type of ruin: ruin by jump or ruin by continuity. The methods used involve finding certain martingales by first finding partial eigenfunctions for the generator of the Markov process composed of the risk process and the underlying Markov chain. We also use certain results from complex function theory as important tools.

We derive necessary and sufficient conditions for the existence of bounded or summable solutions to systems of linear equations associated with Markov chains. This substantially extends a famous result of G. E. H. Reuter, which provides a convenient means of checking various uniqueness criteria for birth-death processes. Our result allows chains with much more general transition structures to be accommodated. One application is to give a new proof of an important result of M. F. Chen concerning upwardly skip-free processes. We then use our generalization of Reuter's lemma to prove new results for downwardly skip-free chains, such as the Markov branching process and several of its many generalizations. This permits us to establish uniqueness criteria for several models, including the general birth, death, and catastrophe process, extended branching processes, and asymptotic birth-death processes, the latter being neither upwardly skip-free nor downwardly skip-free.

We study the following model for a phylogenetic tree on n extant species: the origin of the clade is a random time in the past whose (improper) distribution is uniform on (0,∞); thereafter, the process of extinctions and speciations is a continuous-time critical branching process of constant rate, conditioned on there being the prescribed number n of species at the present time. We study various mathematical properties of this model as n→∞: namely the time of origin and of the most recent common ancestor, the pattern of divergence times within lineage trees, the time series of the number of species, the total number of extinct species, the total number of species ancestral to the extant ones, and the ‘local’ structure of the tree itself. We emphasize several mathematical techniques: the association of walks with trees; a point process representation of lineage trees; and Brownian limits.

Our main result is the proof of the existence of random stationary tessellations in d-dimensional Euclidean space with the following stability property: their distribution is invariant with respect to the operation of iteration (or nesting) of tessellations with an appropriate rescaling. This operation means that the cells of a given tessellation are individually and independently subdivided by independent, identically distributed tessellations, resulting in a new tessellation. It is also shown that, for any stationary tessellation, the sequence that is generated by repeated rescaled iteration converges weakly to such a stable tessellation; thus, the class of all stable stationary tessellations is fully characterized.

In this article, we review known results and present new ones concerning the power spectra of large classes of signals and random fields driven by an underlying point process, such as spatial shot noises (with random impulse response and arbitrary basic stationary point processes described by their Bartlett spectra) and signals or fields sampled at random times or points (where the sampling point process is again quite general). We also obtain the Bartlett spectrum for the general linear Hawkes spatial branching point process (with random fertility rate and general immigrant process described by its Bartlett spectrum). We then obtain the Bochner spectra of general spatial linear birth and death processes. Finally, we address the issues of random sampling and linear reconstruction of a signal from its random samples, reviewing and extending former results.

Let (Xt)t≥0 be a continuous-time irreducible Markov chain on a finite state space E, let v be a map v: E→ℝ\{0}, and let (φt)t≥0 be an additive functional defined by φt=∫0tv(Xs)d s. We consider the case in which the process (φt)t≥0 is oscillating and that in which (φt)t≥0 has a negative drift. In each of these cases, we condition the process (Xt,φt)t≥0 on the event that (φt)t≥0 is nonnegative until time T and prove weak convergence of the conditioned process as T→∞.

In this paper, we consider the asymptotic behavior of stationary probability vectors of Markov chains of GI/G/1 type. The generating function of the stationary probability vector is explicitly expressed by the R-measure. This expression of the generating function is more convenient for the asymptotic analysis than those in the literature. The RG-factorization of both the repeating row and the Wiener-Hopf equations for the boundary row are used to provide necessary spectral properties. The stationary probability vector of a Markov chain of GI/G/1 type is shown to be light tailed if the blocks of the repeating row and the blocks of the boundary row are light tailed. We derive two classes of explicit expression for the asymptotic behavior, the geometric tail, and the semigeometric tail, based on the repeating row, the boundary row, or the minimal positive solution of a crucial equation involved in the generating function, and discuss the singularity classes of the stationary probability vector.

We consider the problem of dynamic admission control in a Markovian loss system with two classes. Jobs arrive at the system in batches; each admitted job requires different service rates and brings different revenues depending on its class. We introduce the definition of a ‘preferred class’ for systems receiving mixed and single-class batches separately, and derive sufficient conditions for each system to have a preferred class. We also establish a monotonicity property of the optimal value functions, which reduces the number of possibly optimal actions.

We consider a hard-sphere model in ℝd generated by a stationary point process N and the lilypond growth protocol: at time 0, every point of N starts growing with unit speed in all directions to form a system of balls in which any particular ball ceases its growth at the instant that it collides with another ball. Some quite general conditions are given, under which it is shown that the model is well defined and exhibits no percolation. The absence of percolation is attributable to the fact that, under our assumptions, there can be no descending chains in N. The proof of this fact forms a significant part of the paper. It is also shown that, in the absence of descending chains, mutual-nearest-neighbour matching can be used to construct a bijective point map as defined by Thorisson.

Let X1,X2,… be a sequence of independent and identically distributed random variables with some continuous distribution function F. Let L(n) and X(n) denote the nth record time and the nth record value, respectively. We refer to the variables Xi as near-nth-record observations if Xi∈(X(n)-a,X(n)], with a>0, and L(n)<i<L(n+1). In this work we study asymptotic properties of the number of near-record observations. We also discuss sums of near-record observations.

Given a sample of genes taken from a large population, we consider the neutral coalescent genealogy and study the theoretical and empirical distributions of the size of the smallest clade containing a fixed gene. We show that the theoretical distribution is strongly related to a Yule distribution of parameter 2, and that the empirical count statistics are asymptotically Gaussian as the number of genes grows to infinity. Then we consider external branches of the coalescent tree, and describe their lengths. Using the infinitely many sites model of mutation, we also describe the conditional distribution of the external branch lengths, given the number of pairwise differences between a reference DNA sequence and the sequence of one closest relative in the sample.

Consider a sequence (Xn) of independent and identically distributed random variables taking nonnegative integer values, and call Xn a record if Xn> max{X1,…,Xn−1}. By means of martingale arguments it is shown that the counting process of records among the first n observations, suitably centered and scaled, is asymptotically normally distributed.

Let X1, X2,… be real-valued random variables. For u>0, define the time of ruin T = T(u) by T = inf{n: X1+⋯+Xn>u} or T=∞ if X1+⋯+Xn≤u for every n = 1,2,…. We are interested in the ruin probabilities of general processes {Xn} for large u. In the presence of heavy tails, one often finds power estimates. Our objective is to specify the associated powers and provide the crude estimate P(T≤xu)≈u−R(x) for large u, for a given x∈ℝ. The rate R(x) will be described by means of tails of partial sums and maxima of {Xn}. We also extend our results to the case of the infinite time horizon.

Our objective is to construct a perfect simulation algorithm for unmarked and marked Hawkes processes. The usual straightforward simulation algorithm suffers from edge effects, whereas our perfect simulation algorithm does not. By viewing Hawkes processes as Poisson cluster processes and using their branching and conditional independence structures, useful approximations of the distribution function for the length of a cluster are derived. This is used to construct upper and lower processes for the perfect simulation algorithm. A tail-lightness condition turns out to be of importance for the applicability of the perfect simulation algorithm. Examples of applications and empirical results are presented.

Consider a random walk S=(Sn: n≥0) that is ‘perturbed’ by a stationary sequence (ξn: n≥0) to produce the process S=(Sn+ξn: n≥0). In this paper, we are concerned with developing limit theorems and approximations for the distribution of Mn=max{Sk+ξk: 0≤k≤n} when the random walk has a drift close to 0. Such maxima are of interest in several modeling contexts, including operations management and insurance risk theory. The associated limits combine features of both conventional diffusion approximations for random walks and extreme-value limit theory.

We consider a compound renewal (Sparre Andersen) risk process with interclaim times that have a Kn distribution (i.e. the Laplace transform of their density function is a ratio of two polynomials of degree at most n ∈ N). The Laplace transform of the expected discounted penalty function at ruin is derived. This leads to a generalization of the defective renewal equations given by Willmot (1999) and Gerber and Shiu (2005). Finally, explicit results are given for rationally distributed claim severities.