Explicit bounds are given for the Kolmogorov and Wasserstein distances between a mixture of normal distributions, by which we mean that the conditional distribution given some $\sigma$ -algebra is normal, and a normal distribution with properly chosen parameter values. The bounds depend only on the first two moments of the first two conditional moments given the $\sigma$ -algebra. The proof is based on Stein’s method. As an application, we consider the Yule–Ornstein–Uhlenbeck model, used in the field of phylogenetic comparative methods. We obtain bounds for both distances between the distribution of the average value of a phenotypic trait over n related species, and a normal distribution. The bounds imply and extend earlier limit theorems by Bartoszek and Sagitov.

A spatio-temporal model of particle or star growth is defined, whereby new unit masses arrive sequentially in discrete time. These unit masses are referred to as candidate stars, which tend to arrive in mass-dense regions and then either form a new star or are absorbed by some neighbouring star of high mass. We analyse the system as time increases, and derive the asymptotic growth rate of the number of stars as well as the size of a randomly chosen star. We also prove that the size-biased mass distribution converges to a Poisson–Dirichlet distribution. This is achieved by embedding our model into a continuous-time Markov process, so that new stars arrive according to a marked Poisson process, with locations as marks, whereas existing stars grow as independent Yule processes. Our approach can be interpreted as a Hoppe-type urn scheme with a spatial structure. We discuss its relevance for and connection to models of population genetics, particle aggregation, image segmentation, epidemic spread, and random graphs with preferential attachment.

In this paper we study a generalized Pólya urn with balls of two colors and a random triangular replacement matrix. We extend some results of Janson (2004), (2005) to the case where the largest eigenvalue of the mean of the replacement matrix is not in the dominant class. Using some useful martingales and the embedding method introduced in Athreya and Karlin (1968), we describe the asymptotic composition of the urn after the nth draw, for large n.

Let Y0, Y1, Y2, … be an i.i.d. sequence of random variables with continuous distribution function, and let P be a simple point process on 0≦t≦∞, independent of the Yj's. We assume that P has a point at t = 0; we associate Yj with the jth point of j≧0, and we say that the Yj's occur at the arrival times of P. Y0 is considered a ‘reference value'. The first Yj (j≧1) to exceed all previous ones is called the first ‘record value', and the time of its occurrence is the first ‘record time'. Subsequent record values and times are defined analogously. We give an infinite series representation for the joint characteristic function of the first n record times, for general P; in some cases the series can be summed. We find the intensity of the record process when P is a general birth process, and when P is a linear birth process with m immigration sources we find the distribution of the number of records in (0, t]. For m = 0 (the Yule process) we give moments of record times and a compact form for the record process intensity. We show that the records occur according to a homogeneous Poisson process when m = 1, and we display a different model with the same behavior, leading to statistical non-identifiability if only the record times are observed. For m = 2, the records occur according to a semi-Markov process; again we display a different model with the same behavior. Finally we give a new derivation of the joint distribution of the interrecord times when P is an arbitrary Poisson process. We relate this result to existing work and to the classical record model. We also obtain a new characterization of the exponential distribution.

We provide a probabilistic proof of the characterization of point processes (on the real line) with the order statistic property. The characterization is used to investigate the homogeneity of such processes and is also related to the martingale theory associated with point processes.

It has been conjectured that a certain transformation of gene frequency space due to Fisher and Bhattacharyya will map the random genetic drift process, or its diffusion approximation, into one with radial symmetry. This paper proves rigorously that the Fisher–Bhattacharyya mapping does not do this. This implies that the initial state of an evolving ensemble can only be unbiasedly estimated from the means of a sample if we weight by the proper divergence times. If the ensemble is known not to have begun at the centroid of frequency space, the estimate of the initial state vector is not simply the arithmetic average, as symmetry analysis of the Christoffel velocity field shows.

Individuals arrive at a geographical area that is initially empty, in accordance with a pure birth process with birth parameters λj = jλ + θ, j ≧ 0. Due to contamination, this geographical area is unsafe for its population and at some fixed time T in the future everyone in the area will be evacuated and no further immigration will be allowed.

Suppose now than an intermediate evacuation time τ, 0 ≦ τ ≦ T, at which time everyone present in the area would be evacuated, is to be chosen. The area would then again fill up with individuals between times τ and T, and, at T, the final evacuation would be made and the area would be permanently sealed off. The problem is to choose τ so as to minimise the total expected cost incurred by time T, where a cost g(x) is incurred for each individual that spends a time x in the area.