A generalization of Markov point processes is introduced in which interactions occur between connected components of the point pattern. A version of the Hammersley-Clifford characterization theorem is proved which states that a point process is a Markov interacting component process if and only if its density function is a product of interaction terms associated with cliques of connected components. Integrability and superpositional properties of the processes are shown and a pairwise interaction example is used for detailed exploration.

Based on remote sensing of a potential minefield, point locations are identified, some of which may not be mines. The mines and mine-like objects are to be distinguished based on their point patterns, although it must be emphasized that all one sees is the superposition of their locations. In this paper, we construct a hierarchical spatial point-process model that accounts for the different patterns of mines and mine-like objects and uses posterior analysis to distinguish between them. Our Bayesian approach is applied to minefield data obtained from a multispectral video remote-sensing system.

New metrics are introduced in the space of random measures and are applied, with various modifications of the contraction method, to prove existence and uniqueness results for self-similar random fractal measures. We obtain exponential convergence, both in distribution and almost surely, of an iterative sequence of random measures (defined by means of the scaling operator) to a unique self-similar random measure. The assumptions are quite weak, and correspond to similar conditions in the deterministic case.

The fixed mass case is handled in a direct way based on regularity properties of the metrics and the properties of a natural probability space. Proving convergence in the random mass case needs additional tools, such as a specially adapted choice of the space of random measures and of the space of probability distributions on measures, the introduction of reweighted sequences of random measures and a comparison technique.

Consider the Delaunay graph and the Voronoi tessellation constructed with respect to a Poisson point process. The sequence of nuclei of the Voronoi cells that are crossed by a line defines a path on the Delaunay graph. We show that the evolution of this path is governed by a Markov chain. We study the ergodic properties of the chain and find its stationary distribution. As a corollary, we obtain the ratio of the mean path length to the Euclidean distance between the end points, and hence a bound for the mean asymptotic length of the shortest path.

We apply these results to define a family of simple incremental algorithms for constructing short paths on the Delaunay graph and discuss potential applications to routeing in mobile communication networks.

Consider the following birth-growth model in ℝ. Seeds are born randomly according to an inhomogeneous space-time Poisson process. A newly formed point immediately initiates a bi-directional coverage by sending out a growing branch. Each frontier of a branch moves at a constant speed until it meets an opposing one. New seeds continue to form on the uncovered parts on the line. We are interested in the time until a bounded interval is completely covered. The exact and limiting distributions as the length of interval tends to infinity are obtained for this completion time by considering a related Markov process. Moreover, some strong limit results are also established.

The paper compares non-parametric (design-based) and parametric (model-based) approaches to the analysis of data in the form of replicated spatial point patterns in two or more experimental groups. Basic questions for data of this kind concern estimating the properties of the underlying spatial point process within each experimental group, and comparing the properties between groups. A non-parametric approach, building on work by Diggle et. al. (1991), summarizes each pattern by an estimate of the reduced second moment measure or K-function (Ripley (1977)) and compares mean K-functions between experimental groups using a bootstrap testing procedure. A parametric approach fits particular classes of parametric model to the data, uses the model parameter estimates as summaries and tests for differences between groups by comparing fits with and without the assumption of common parameter values across groups. The paper discusses how either approach can be implemented in the specific context of a single-factor replicated experiment and uses simulations to show how the parametric approach can be more efficient when the underlying model assumptions hold, but potentially misleading otherwise.

Let n random points be uniformly and independently distributed in the unit square, and count the number W of subsets of k of the points which are covered by some translate of a small square C. If n|C| is small, the number of such clusters is approximately Poisson distributed, but the quality of the approximation is poor. In this paper, we show that the distribution of W can be much more closely approximated by an appropriate compound Poisson distribution CP(λ1, λ2,…). The argument is based on Stein's method, and is far from routine, largely because the approximating distribution does not satisfy the simplifying condition that iλi be decreasing.

We study the transformation of a non-stationary point process ξ on ℝn into a weakly stationary point process ͂ξ, with ͂ξ(B) = ξ(Φ-1(B)), where B is a Borel set, via a deformation Φ of the space ℝn. When the second-order measure is regular, Φ is uniquely determined by the homogenization equations of the second-order measure. In contrast, the first-order homogenization transformation is not unique. Several examples of point processes and transformations are investigated with a particular interest to Poisson processes.

We discuss necessary and sufficient conditions for power-law and polynomial models to be correlation functions on bounded domains. These results date back to unpublished work by Matheron (1974) and generalize the findings of Gneiting (1999).

Useful approximations have been developed along the years to predict the precision of systematic sampling for measurable functions of a bounded support in ℝd. Recently, the theory of systematic sampling on ℝ has received a thrust. In geometric sampling, design based unbiased estimators exist, however, which imply systematic sampling on the circle (𝕊1) and the semicircle (ℍ1); the planimeter estimator of an area, or the Buffon-Steinhaus estimator of curve length in the plane constitute popular examples. Over the last two decades, many other estimators of geometric measures have been obtained which imply systematic sampling also on the sphere (𝕊2). In this paper we adapt the theory available for non-periodic functions of bounded support on ℝ to periodic functions, and thereby to 𝕊1 and ℍ1, and we obtain new estimators of the corresponding variance approximations. Further we consider - we believe for the first time - the problem of predicting the precision of systematic sampling in 𝕊2. The paper starts with a historical perspective, and ends with suggestions for further research.

Let G be an infinite, locally finite, connected graph with bounded degree. We show that G supports phase transition in all or none of the following five models: bond percolation, site percolation, the Ising model, the Widom-Rowlinson model and the beach model. Some, but not all, of these implications hold without the bounded degree assumption. We finally give two examples of (random) unbounded degree graphs in which phase transition in all five models can be established: supercritical Galton-Watson trees, and Poisson-Voronoi tessellations of ℝd for d ≥ 2.

Often, the statistical analysis of the shape of a random planar curve is based on a model for a polygonal approximation to the curve. In the present paper, we instead describe the curve as a continuous stochastic deformation of a template curve. The advantage of this continuous approach is that the parameters in the model do not relate to a particular polygonal approximation. A somewhat similar approach has been used in Kent et. al. (1996), who describe the limiting behaviour of a model with a first-order Markov property as the landmarks on the curve become closely spaced; see also Grenander (1993). The model studied in the present paper is an extension of this model. Our model possesses a second-order Markov property. Its geometrical characteristics are studied in some detail and an explicit expression for the covariance function is derived. The model is applied to the boundaries of profiles of cell nuclei from a benign tumour and a malignant tumour. It turns out that the model with the second-order Markov property is the most appropriate, and that it is indeed possible to distinguish between the two samples.

We study Poisson-Voronoi tessellations in three-dimensional hyperbolic spaces, and give explicit expressions for mean surface area, mean perimeter length, and mean number of vertices of their cells. Furthermore we compare these mean characteristics with those for Poisson-Voronoi tessellations in three-dimensional Euclidean spaces. It is shown that, as the absolute value of the curvature of hyperbolic spaces increases from zero to infinity, these mean characteristics increase monotonically from those for the Euclidean case to infinity.

We consider a fully connected queueing network in which customers have one direct and many alternative routes through the network, and where customer routeing is dynamic. We obtain an asymptotically optimal routeing policy, taking the limit as the number of queues of the network increases. We observe that good policies route customers directly, unless there is a danger of servers becoming idle, in which case customers should be routed alternatively so as to avoid such idleness, and this leads to policies that perform well in moderate-sized networks.

The isotropic planar point processes of phase-type are natural generalizations of the Poisson process on the plane. On the one hand, those processes are isotropic and stationary for the mean count, as in the case of the Poisson process. On the other hand, they exhibit dependence of counts in disjoint sets. In a recent paper, we have proved that the number of points in a square window has a Poisson distribution asymptotically as the window is located far away from the origin of the process. We extend our work to the case of a window of arbitrary shape.

By using an adaptation of the radial generation method, we give an integral formula for the proportion of triangles in a Poisson-Voronoi tessellation, which gives a value of 0.0112354 to 7 decimal places. We also obtain the first four moments of some characteristics of triangles.

In [8], Le showed that procrustean mean shapes of samples are consistent estimates of Fréchet means for a class of probability measures in Kendall's shape spaces. In this paper, we investigate the analogous case in Bookstein's shape space for labelled triangles and propose an estimator that is easy to compute and is a consistent estimate of the Fréchet mean, with respect to sinh(δ/√2), of any probability measure for which such a mean exists. Furthermore, for a certain class of probability measures, this estimate also tends almost surely to the Fréchet mean calculated with respect to the Riemannian distance δ.

Suppose n particles xi in a region of the plane (possibly representing biological individuals such as trees or smaller organisms) have a joint density proportional to exp{-∑i<jϕ(n(xi-xj))}, with ℝd; a specified function of compact support. We obtain a Poisson process limit for the collection of rescaled interparticle distances as n becomes large. We give corresponding results for the case of several types of particles, representing different species, and also for the area-interaction (Widom-Rowlinson) point process of interpenetrating spheres.

Haploid population models with non-overlapping generations and fixed population size N are considered. It is assumed that the family sizes ν1,…,νN within a generation are exchangeable random variables. Rates of convergence for the finite-dimensional distributions of a properly time-scaled ancestral coalescent process are established and expressed in terms of the transition probabilities of the ancestral process, i.e., in terms of the joint factorial moments of the offspring variables ν1,…,νN.

The Kingman coalescent appears in the limit as the population size N tends to infinity if and only if triple mergers are asymptotically negligible in comparison with binary mergers. In this case, a simple upper bound for the rate of convergence of the finite-dimensional distributions is derived. It depends (up to some constants) only on the three factorial moments E((ν1)2), E((ν1)3) and E((ν1)2(ν2)2), where (x)k := x(x-1)…(x-k+1). Examples are the Wright-Fisher model, where the rate of convergence is of order N-1, and the Moran model, with a convergence rate of order N-2.

We establish new multidimensional martingales for Markov additive processes and certain modifications of such processes (e.g., such processes with reflecting barriers). These results generalize corresponding one-dimensional martingale results for Lévy processes. This martingale is then applied to various storage processes, queues and Brownian motion models.