A sequence of random variables η0, η1, …, ηn, … is defined by the recurrence formula ηn = max (ηn–1 + ξn, 0) where η0 is a discrete random variable taking on non-negative integers only and ξ1, ξ2, … ξn, … is a semi-Markov sequence of discrete random variables taking on integers only. Define Δ as the smallest n = 1, 2, … for which ηn = 0. The random variable ηn can be interpreted as the content of a dam at time t = n(n = 0, 1, 2, …) and Δ as the time of first emptiness. This paper deals with the determination of the distributions of ηn and Δ by using the method of matrix factorisation.

The structured coalescent is a continuous-time Markov chain which describes the genealogy of a sample of homologous genes from a subdivided population. Assuming this model, some results are proved relating to the genealogy of a pair of genes and the extent of subpopulation differentiation, which are valid under certain graph-theoretic symmetry and regularity conditions on the structure of the population. We first review and extend earlier results stating conditions under which the mean time since the most recent common ancestor of a pair of genes from any single subpopulation is independent of the migration rate and equal to that of two genes from an unstructured population of the same total size. Assuming the infinite alleles model of neutral mutation with a small mutation rate, we then prove a simple relationship between the migration rate and the value of Wright's coefficient FST for a pair of neighbouring subpopulations, which does not depend on the precise structure of the population provided that this is sufficiently symmetric.

A continuous-time threshold autoregressive process of order two (CTAR(2)) is constructed as the first component of the unique (in law) weak solution of a stochastic differential equation. The Cameron–Martin–Girsanov formula and a random time-change are used to overcome the difficulties associated with possible discontinuities and degeneracies in the coefficients of the stochastic differential equation. A sequence of approximating processes that are well-suited to numerical calculations is shown to converge in distribution to a solution of this equation, provided the initial state vector has finite second moments. The approximating sequence is used to fit a CTAR(2) model to percentage relative daily changes in the Australian All Ordinaries Index of share prices by maximization of the ‘Gaussian likelihood'. The advantages of non-linear relative to linear time series models are briefly discussed and illustrated by means of the forecasting performance of the model fitted to the All Ordinaries Index.

We consider the Poisson Boolean model of percolation where the percolating shapes are convex regions. By an enhancement argument we strengthen a result of Jonasson (2000) to show that the critical intensity of percolation in two dimensions is minimized among the class of convex shapes of unit area when the percolating shapes are triangles, and, for any other shape, the critical intensity is strictly larger than this minimum value. We also obtain a partial generalization to higher dimensions. In particular, for three dimensions, the critical intensity of percolation is minimized among the class of regular polytopes of unit volume when the percolating shapes are tetrahedrons. Moreover, for any other regular polytope, the critical intensity is strictly larger than this minimum value.

The main characteristics of the time-nonhomogeneous Johnson–Mehl tessellation with specified nucleation intensity α (t) = αtß, β > −1 and α > 0 being constants and t time, are investigated as functions of β.

The paper deals with asymptotic properties of the transition probabilities of a countable non-homogeneous Markov chain, the main concept used in the proofs being that of the tail σ-field of the chain. A state classification similar to that existing in the homogeneous case is given and a strong ratio limit property is shown to parallel the basic limit theorem for positive homogeneous chains. Some global asymptotic properties for null chains are also derived.

Galton-Watson forests consisting of N roots (or trees) and n nonroot vertices are studied. The limit distributions of the number of leaves in such a forest are obtained. By a leaf we mean a vertex from which there are no arcs emanating.

Intersection densities are introduced for a large class of nonstationary Poisson processes of hypersurfaces and inequalities for them are proved. In doing so, similar results from both Wieacker (1986) and Schneider (2003) are summarized in one theorem and the concept of an associated zonoid of a Poisson process of hypersurfaces is generalized to a nonstationary setting.

Given a pure-jump subordinator (i.e. nondecreasing Lévy process with no drift) with continuous Lévy measure ν, we derive a formula for the distribution function Fs (x; t) at time t of the associated subordinator whose Lévy measure is the restriction of ν to (0,s]. It will be expressed in terms of ν and the marginal distribution function F (⋅; t) of the original process. A generalization concerning an arbitrary truncation of ν will follow. Under certain conditions, an analogous formula will be obtained for the nth derivative, ∂nFs (x; t) ∂ xn. The requirement that ν is continuous is shown to have no intrinsic meaning. A number of interesting results involving the size ordered jumps of subordinators will be derived. An appropriate approximation for the small jumps of a gamma process will be considered, leading to a revisiting of the generalized Dickman distribution.

We consider sequential decision processes in which the posterior distribution of some unobservable must be included among the sufficient statistics. The technique of an earlier paper (Whittle (1964)) is applied to show that the residual loss function is necessarily homogeneous of degree one in this distribution, if certain conventions are adopted. This point leads to considerable simplification. Examples are given of the great variety of problems which can advantageously be formulated in this manner.

Suppose that in a ballot candidate A scores a votes and candidate B scores b votes, and that all the possible voting records are equally probable. Corresponding to the first r votes, let α r and β r be the numbers of votes registered for A and B, respectively. Let p be an arbitrary positive real number. Denote by δ (a, b, p)[δ *(a, b, ρ)] the number of values of r for which the inequality , r = 1, ···, a + b, holds. Heretofore the probability distributions of δand δ* have been derived for only a restricted set of values of a, b, and ρ, although, as pointed out here, they are obtainable for all values of (a, b, ρ) by using a result of Takács (1964). In this paper we present a derivation of the distribution of δ [δ *] whose development, for any (a, b, ρ), leads to both necessary and sufficient conditions for δ [δ *] to have a discrete uniform distribution.

Suppose you have u units of ammunition and want to destroy as many as possible of a sequence of attacking enemy aircraft. If you fire v = v(u), 0 , units of your ammunition at the first enemy, it survives with probability qv, where 0 < q < 1 is given, and then kills you. With the complementary probability, 1 – qv, you destroy the aircraft and you live to face the next enemy with only u – v units of ammunition remaining. It seems almost obvious that any strategy which maximizes the expected number of enemies destroyed before you die will fire more units at the first enemy as u increases, i.e., it seems obvious that v′(u) 0 under optimal play. We show this to be false, thereby disproving an appealing conjecture proposed by Weber. We also consider a variant of this problem and find that the counterpart of Weber's conjecture holds in some cases.

A random spatial coverage process whose generating point process is homogeneous Poisson, and whose attached random sets are independent and identically distributed, is called a Boolean model. Motivated by Błaszczyszyn et al. [1], distributional and higher moment properties of the size of clumps (connected clusters of overlapping sets) in this model are derived. This provides some complements to the result on the finiteness of the first moment presented in Hall [3]. The key idea is to construct a certain coupling process for a multitype branching process that dominates the clump size.

We deal with the problem of seating an airplane's passengers optimally, namely in the fastest way. Under several simplifying assumptions, whereby the passengers are infinitely thin and react within a constant time to boarding announcements, we are able to rewrite the asymptotic problem as a calculus of variations problem with constraints. This problem is solved in turn using elementary methods. While the optimal policy is not unique, we identify a rigid discrete structure which is common to all solutions. We also compare the (nontrivial) optimal solutions we find with some simple boarding policies, one of which is shown to be near-optimal.

The classical models in risk theory consider a single type of claim. In the insurance business, however, several business lines with separate claim arrival processes appear naturally, and the individual claim processes may not be independent. We introduce a new class of models for such situations, where the underlying counting process is a multivariate continuous-time Markov chain of pure-birth type and the dependency of the components arises from the fact that the birth rate for a specific claim type may depend on the number of claims in the other component processes. Under certain conditions, we obtain a fluid limit, i.e. a functional law of large numbers for these processes. We also investigate the consequences of such results for questions of interest in insurance applications. Several specific subclasses of the general model are discussed in detail and the Cramér asymptotics of the ruin probabilities are derived in particular cases.

We give a functional central limit theorem for spatial birth and death processes based on the representation of such processes as solutions of stochastic equations. For any bounded and integrable function in Euclidean space, we define a family of processes which is obtained by integrals of this function with respect to the centered and scaled spatial birth and death process with constant death rate. We prove that this family converges weakly to a Gaussian process as the scale parameter goes to infinity. We do not need the birth rates to have a finite range of interaction. Instead, we require that the birth rates have a range of interaction that decays polynomially. In order to show the convergence of the finite-dimensional distributions of the above processes, we extend Penrose's multivariate spatial central limit theorem. An example of the asymptotic normalities of the time-invariance estimators for the birth rates of spatial point processes is given.

The Halfin–Whitt regime, or the quality-and-efficiency-driven (QED) regime, for multiserver systems refers to a situation with many servers, a critical load, and yet favorable system performance. We apply this regime to the classical multiserver loss system with slow retrials. We derive nondegenerate limiting expressions for the main steady-state performance measures, including the retrial rate and the blocking probability. It is shown that the economies of scale associated with the QED regime persist for systems with retrials, although in situations when the load becomes extremely critical the retrials cause deteriorated performance. Most of our results are obtained by a detailed analysis of Cohen's equation that defines the retrial rate in an implicit way. The limiting expressions are established by studying prelimit behavior and exploiting the connection between Cohen's equation and Mills' ratio for the Gaussian and Poisson distributions.

Let Sn, k denote the random graph obtained by placing points in a square box of area n according to a Poisson process of intensity 1 and joining each point to its k nearest neighbours. Balister, Bollobás, Sarkar and Walters (2005) conjectured that, for every 0 < ε < 1 and all sufficiently large n, there exists C = C(ε) such that, whenever the probability that Sn, k is connected is at least ε, then the probability that Sn, k+C is connected is at least 1 - ε. In this paper we prove this conjecture. As a corollary, we prove that there exists a constant C' such that, whenever k(n) is a sequence of integers such that the probability Sn,k(n) is connected tends to 1 as n → ∞, then, for any integer sequence s(n) with s(n) = o(logn), the probability Sn,k(n)+⌊C'slog logn⌋ is s-connected (i.e. remains connected after the deletion of any s − 1 vertices) tends to 1 as n → ∞. This proves another conjecture given in Balister, Bollobás, Sarkar and Walters (2009).

Consider first-passage percolation on the square lattice. Hammersley and Welsh, who introduced the subject in 1965, conjectured that the expected minimum travel time from (0, 0) to (n, 0) along paths contained in the cylinder is always non-decreasing in n. However, when the bonds have time-coordinate 1 with probability p and 0 with probability 1 – p (0 < P < 1), then, for p sufficiently small, we get a counterexample.

We consider decay properties including the decay parameter, invariant measures, invariant vectors, and quasistationary distributions of a Markovian bulk-arriving queue that stops immediately after hitting the zero state. Investigating such behavior is crucial in realizing the busy period and some other related properties of Markovian bulk-arriving queues. The exact value of the decay parameter λC is obtained and expressed explicitly. The invariant measures, invariant vectors, and quasistationary distributions are then presented. We show that there exists a family of invariant measures indexed by λ ∈ [0, λC]. We then show that, under some conditions, there exists a family of quasistationary distributions, also indexed by λ ∈ [0, λC]. The generating functions of these invariant measures and quasistationary distributions are presented. We further show that a stopped Markovian bulk-arriving queue is always λC-transient and some deep properties are revealed. The clear geometric interpretation of the decay parameter is explained. A few examples are then provided to illustrate the results obtained in this paper.