Joseph Mecke, Professor of Mathematics at the Friedrich Schiller University, Jena, celebrates his 65th birthday on 18 February 2003. To commemorate this happy occasion, the signatories invited some colleagues, friends and former students to honour Joseph as a mathematician, teacher and research worker by submitting research articles to this issue of Advances in Applied Probability. We trust that the following papers will pay homage to the influence he has exerted in point process theory and stochastic geometry during his distinguished career.

We define a class of tessellation models based on perturbing or deforming standard tessellations such as the Voronoi tessellation. We show how distributions over this class of ‘deformed’ tessellations can be used to define prior distributions for models based on tessellations, and how inference for such models can be carried out using Markov chain Monte Carlo methods; stability properties of the algorithms are investigated. Our approach applies not only to fixed dimension problems, but also to variable dimension problems, in which the number of cells in the tessellation is unknown. We illustrate our methods with two real examples. The first relates to reconstructing animal territories, represented by the individual cells of a tessellation, from observation of an inhomogeneous Poisson point process. The second example involves the analysis of an image of a cross-section through a sample of metal, with the tessellation modelling the micro-crystalline structure of the metal.

Our aim is to estimate the volume-weighted mean of the volumes of three-dimensional ‘particles’ (compact, not-necessarily-convex subsets) from plane sections of the particle population. The standard stereological technique is to place test lines in the plane section, and measure cubed intercept lengths with the two-dimensional particle profiles. This paper discusses more efficient estimators obtained by integrating over all possible placements of the test line. We prove that these estimators have smaller variance than the line transect estimators, and indeed are related to them by the Rao-Blackwell process. In the improved estimators, the cubed intercept length is replaced by a moment of the distance between two points in the section profile. This can be computed as a moment of the set covariance function, which in turn is computable using the fast Fourier transform. We also derive an isoperimetric-type inequality between the improved estimator and the area-weighted 3/2th moment of the profile areas. Finally, we present two practical applications to particles of silicon carbide and to synaptic boutons in brain tissue. We estimate the variance of the technique and the gain in efficiency over line transect techniques; the efficiency improvement appears to be as much as one order of magnitude.

This paper introduces a new graph constructed from a point process. The idea is to connect a point with its nearest neighbour, then to the second nearest and continue this process until the point belongs to the interior of the convex hull of these nearest neighbours. The number of such neighbours is called the degree of a point. We derive the distribution of the degree of the typical point in a Poisson process, prove a central limit theorem for the sum of degrees, and propose an edge-corrected estimator of the distribution of the degree that is unbiased for a stationary Poisson process. Simulation studies show that this degree is a useful concept that allows the separation of clustering and repulsive behaviour of point processes.

A known gamma-type result for the Poisson process states that certain domains defined through configuration of the points (or ‘particles’) of the process have volumes which are gamma distributed. By proving the corresponding sequential gamma-type result, we show that in some cases such a domain allows for decomposition into subdomains each having independent exponentially distributed volumes. We consider other examples—based on the Voronoi and Delaunay tessellations—where a natural decomposition does not produce subdomains with exponentially distributed volumes. A simple algorithm for the construction of a typical Voronoi flower arises in this work. In our theoretical development, we generalize the classical theorem of Slivnyak, relating it to the strong Markov property of the Poisson process and to a result of Mecke and Muche (1995). This new theorem has interest beyond the specific problems being considered here.

Distance measurements are useful tools in stochastic geometry. For a Boolean model Z in ℝd, the classical contact distribution functions allow the estimation of important geometric parameters of Z. In two previous papers, several types of generalized contact distributions have been investigated and applied to stationary and nonstationary Boolean models. Here, we consider random sets Z which are generated as the union sets of Poisson processes X of k-flats, k ∈ {0, …, d-1}, and study distances from a fixed point or a fixed convex body to Z. In addition, we also consider the distances from a given flat or a flag consisting of flats to the individual members of X and investigate the associated process of nearest points in the flats of X. In particular, we discuss to which extent the directional distribution of X is determined by this point process. Some of our results are presented for more general stationary processes of flats.

We study percolation and Ising models defined on generalizations of quad-trees used in multiresolution image analysis. These can be viewed as trees for which each mother vertex has 2d daughter vertices, and for which daughter vertices are linked together in d-dimensional Euclidean configurations. Retention probabilities and interaction strengths differ according to whether the relevant bond is between mother and daughter or between neighbours. Bounds are established which locate phase transitions and show the existence of a coexistence phase for the percolation model. Results are extended to the corresponding Ising model using the Fortuin-Kasteleyn random-cluster representation.

In order to increase the variety of feasible models for random stationary tessellations (mosaics), two operations acting on tessellations are studied: superposition and iteration (the latter is also referred to as nesting). The superposition of two planar tessellations is the superposition of the edges of the cells of both tessellations. The iteration of tessellations means that one tessellation is chosen as a ‘frame’ tessellation. The single cells of this ‘frame’ are simultaneously and independently subdivided by cut-outs of tessellations of an independent and identically distributed sequence of tessellations. In the present paper, we investigate the limits for sequences of tessellations that are generated by consecutive application of superposition or iteration respectively. Sequences of (renormalised) superpositions of stationary planar tessellations converge weakly to Poisson line tessellations. For consecutive iteration the notion of stability of distributions is adapted and necessary conditions are formulated for those tessellations which may occur as limits.

Results about stationary Poisson hyperplane processes and the induced hyperplane mosaics are extended to the case where, instead of stationarity, it is only assumed that the intensity measure has a (possibly continuous) density with respect to some translation-invariant measure. Intensities and quermass densities, which are constant in the stationary case, are then replaced by functions. In a similar way, the associated zonoid (Matheron's Steiner convex set) is generalized.

This paper studies the law of any real powers of the integral of geometric Brownian motion over finite time intervals. As its main results, an apparently new integral representation is derived and its interrelations with the integral representations for these laws originating by Yor and by Dufresne are established. In fact, our representation is found to furnish what seems to be a natural bridge between these other two representations. Our results are obtained by enhancing the Hartman-Watson Ansatz of Yor, based on Bessel processes and the Laplace transform, by complex analytic techniques. Systematizing this idea in order to overcome the limits of Yor's theory seems to be the main methodological contribution of the paper.

Denote by αt(μ) the probability law of At(μ) = ∫0texp(2(Bs+μ s))ds for a Brownian motion {Bs, s ≥ 0}. It is well known that αt(μ) is of interest in a number of domains, e.g. mathematical finance, diffusion processes in random environments, stochastic analysis on hyperbolic spaces and so on, but that it has complicated expressions. Recently, Dufresne obtained some remarkably simple expressions for αt(0) and αt(1), as well as an equally remarkable relationship between αt(μ) and αt(ν) for two different drifts μ and ν. In this paper, hinging on previous results about αt(μ), we give different proofs of Dufresne's results and present extensions of them for the processes {At(μ), t ≥ 0}.

We consider a failure-prone system which operates in continuous time and is subject to condition monitoring at discrete time epochs. It is assumed that the state of the system evolves as a continuous-time Markov process with a finite state space. The observation process is stochastically related to the state process which is unobservable, except for the failure state. Combining the failure information and the information obtained from condition monitoring, and using the change of measure approach, we derive a general recursive filter, and, as special cases, we obtain recursive formulae for the state estimation and other quantities of interest. Up-dated parameter estimates are obtained using the EM algorithm. Some practical prediction problems are discussed and an illustrative example is given using a real dataset.

Sums of independent random variables concentrated on the same finite discrete, not necessarily lattice, set of points are approximated by compound Poisson distributions and signed compound Poisson measures. Such approximations can be more accurate than the normal distribution. Short asymptotic expansions are constructed.

We consider random recursive fractals and prove fine results about their local behaviour. We show that for a class of random recursive fractals the usual multifractal spectrum is trivial in that all points have the same local dimension. However, by examining the local behaviour of the measure at typical points in the set, we establish the size of fine fluctuations in the measure. The results are proved using a large deviation principle for a class of general branching processes which extends the known large deviation estimates for the supercritical Galton-Watson process.

Consider N towers each made up of a number of counters. At each step a tower is chosen at random, a counter removed which is then added to another tower also chosen at random. The probability distribution for the time needed to empty one of the towers is obtained in the case N = 3. Arguments are set forward as to why no simple formulae can be expected for N > 3. An asymptotic expression for the mean time before one of the towers becomes empty is derived in the case of four towers when they all initially contain a comparably large number of counters. We then study related problems, in particular the ruin problem for three players. Here we use simple martingale methodology as well as a solution proposed by T. S. Ferguson for a slightly modified problem. Throughout the paper it is our main objective to shed light on the reasons why the case N > 3 is so substantially different from the case N ≤ 3.