This paper reviews the principal properties of stationary Poisson processes of k-dimensional linear affine subvarieties in a d-dimensional vector space, and outlines basic results for the estimation of parameters from a single incomplete realization of the process.

Recent progress in microdissection and in DNA sequencing has facilitated the subsampling of multi-focal cancers in organs such as the liver in several hundred spots, helping to determine the pattern of mutations in each of these spots. This has led to the construction of genealogies of the primary, secondary, tertiary, and so forth, foci of the tumor. These studies have led to diverse conclusions concerning the Darwinian (selective) or neutral evolution in cancer. Mathematical models of the development of multi-focal tumors have been devised to support these claims. We offer a model for the development of a multi-focal tumor: it is a mathematically rigorous refinement of a model of Ling et al. (2015). Guided by numerical studies and simulations, we show that the rigorous model, in the form of an infinite-type branching process, displays distributions of tumor size which have heavy tails and moments that become infinite in finite time. To demonstrate these points, we obtain bounds on the tails of the distributions of the process and an infinite series expansion for the first moments. In addition to its inherent mathematical interest, the model is corroborated by recent literature on apparent super-exponential growth in cancer metastases.

In, Tory and Pickard show that a simple subclass of unilateral AR processes identifies with Gaussian Pickard random fields on Z2. First, we extend this result to the whole class of unilateral AR processes, by showing that they all satisfy a Pickard-type property, under which correlation matching and maximum entropy properties are assessed. Then, it is established that the Pickard property provides the ‘missing’ equations that complement the two-dimensional Yule-Walker equations, in the sense that the conjunction defines a one-to-one mapping between the set of AR parameters and a set of correlations. It also implies Markov chain conditions that allow exact evaluation of the likelihood and an exact sampling scheme on finite lattices.

Radix Sort is a sorting algorithm based on analyzing digital data. We study the number of swaps made by Radix Select (a one-sided version of Radix Sort) to find an element with a randomly selected rank. This kind of grand average provides a smoothing over all individual distributions for specific fixed-order statistics. We give an exact analysis for the grand mean and an asymptotic analysis for the grand variance, obtained by poissonization, the Mellin transform, and depoissonization. The digital data model considered is the Bernoulli(p). The distributions involved in the swaps experience a phase change between the biased cases (p ≠ ½) and the unbiased case (p = ½). In the biased cases, the grand distribution for the number of swaps (when suitably scaled) converges to that of a perpetuity built from a two-point distribution. The tool for this proof is contraction in the Wasserstein metric space, and identifying the limit as the fixed-point solution of a distributional equation. In the unbiased case the same scaling for the number of swaps gives a limiting constant in probability.

In this paper we show that the waiting and the sojourn times of a customer in a single-stage, multiple-server, G/G/c queueing system are increasing and starshaped with respect to the mean service time. Usefulness of this result in the design of the optimal service speed in the G/G/c queueing system is also demonstrated.

We study solvency of insurers in a comprehensive model where various economic factors affect the capital developments of the companies. The main interest is in the impact of real growth to ruin probabilities. The volume of the business is allowed to increase or decrease. In the latter case, the study is focused on run-off companies. Our main results give sharp asymptotic estimates for infinite-time ruin probabilities.

We consider an extension of the Poisson hail model where the service speed is either 0 or ∞ at each point of the Euclidean space. We use and develop tools pertaining to sub-additive ergodic theory in order to establish shape theorems for the growth of the ice-heap under light tail assumptions on the hailstone characteristics. The asymptotic shape depends on the statistics of the hailstones, the intensity of the underlying Poisson point process, and on the geometrical properties of the zero speed set.

The paper discussed the value of the ‘Optimum Model” in describing the inter-relationship between metric characters and natural fitness. This model represents one extreme of a whole spectrum of models. It can be used to predict changes in the mean of the character under natural selection, but the basic premises underlying such use are not well understood.

We consider a migration process whose singleton process is a time-dependent Markov replacement process. For the singleton process, which may be treated as either open or closed, we study the limiting distribution, the distribution of the time to replacement and related quantities. For a replacement process in equilibrium we obtain a version of Little's law and we provide conditions for reversibility. For the resulting linear population process we characterize exponential ergodicity for two types of environmental behaviour, i.e. either convergent or cyclic, and finally for large population sizes a diffusion approximation analysis is provided.

In Michel and Paroux (2003) the authors proposed a new proof of a well-known convergence result for the scaled elementary connected vacant component in the high intensity Boolean model towards the Crofton cell of the Poisson hyperplane process (see, e.g. Hall (1985)). In this paper we investigate the second-order term in this convergence when the two-dimensional Boolean model and the Poisson line process are coupled on the same probability space. We consider the particular case where the grains are discs with random radii. A precise coupling between the Boolean model and the Poisson line process is first established. A result of directional convergence in distribution for the difference of the two sets involved is then derived. Eventually, we show the convergence of the process, measuring the difference between the two random sets, once rescaled, as a function of the direction.

We consider homogeneous STIT tessellations Y in the ℓ-dimensional Euclidean space ℝℓ and show the triviality of the tail σ-algebra. This is a sharpening of the mixing result by Lachièze-Rey (2001).

Brown‒Resnick processes are max-stable processes that are associated to Gaussian processes. Their simulation is often based on the corresponding spectral representation which is not unique. We study to what extent simulation accuracy and efficiency can be improved by minimizing the maximal variance of the underlying Gaussian process. Such a minimization is a difficult mathematical problem that also depends on the geometry of the simulation domain. We extend Matheron's (1974) seminal contribution in two directions: (i) making his description of a minimal maximal variance explicit for convex variograms on symmetric domains, and (ii) proving that the same strategy also reduces the maximal variance for a huge class of nonconvex variograms representable through a Bernstein function. A simulation study confirms that our noncostly modification can lead to substantial improvements among Gaussian representations. We also compare it with three other established algorithms.

We consider the problem of numerical integration when the sampling nodes form a stationary point process on the real line. In previous papers it was argued that a naïve Riemann sum approach can cause a severe variance inflation when the sampling points are not equidistant. We show that this inflation can be avoided using a higher-order Newton–Cotes quadrature rule which exploits smoothness properties of the integrand. Under mild assumptions, the resulting estimator is unbiased and its variance asymptotically obeys a power law as a function of the mean point distance. If the Newton–Cotes rule is of sufficiently high order, the exponent of this law turns out to only depend on the point process through its mean point distance. We illustrate our findings with the stereological estimation of the volume of a compact object, suggesting alternatives to the well-established Cavalieri estimator.

We propose a two-parameter family of conjugate prior distributions for the number of undiscovered objects in a class of Bayesian search models. The family contains the one-parameter Euler and Heine families as special cases. The two parameters may be interpreted respectively as an overall success rate and a rate of depletion of the source of objects. The new family gives enhanced flexibility in modelling.

In this article we study the parabolic system of equations which is closely related to a multitype branching Brownian motion. Particular attention is paid to the monotone traveling wave solutions of this system. Provided with some moment conditions, we show the existence, uniqueness, and asymptotic behaviors of such waves with speed greater than or equal to a critical value c̲ and nonexistence of such waves with speed smaller than c̲.

By using large devaitions theory, we give asymptotic formulas for the transient blocking probabilities of M/M/N/N and M (with finite Poissonian sources) M/N/N queues.