This paper considers a germ-grain model for a random system of non-overlapping spheres in ℝd for d = 1, 2 and 3. The centres of the spheres (i.e. the ‘germs’ for the ‘grains’) form a stationary Poisson process; the spheres result from a uniform growth process which starts at the same instant in all points in the radial direction and stops for any sphere when it touches any other sphere. Upper and lower bounds are derived for the volume fraction of space occupied by the spheres; simulation yields the values 0.632, 0.349 and 0.186 for d = 1, 2 and 3. The simulations also provide an estimate of the tail of the distribution function of the volume of a randomly chosen sphere; these tails are compared with those of two exponential distributions, of which one is a lower bound and is an asymptote at the origin, and the other has the same mean as the simulated distribution. An upper bound on the tail of the distribution is also an asymptote at the origin but has a heavier tail than either of these exponential distributions. More detailed information for the one-dimensional case has been found by Daley, Mallows and Shepp; relevant information is summarized, including the volume fraction 1 - e-1 = 0.63212 and the tail of the grain volume distribution e-y exp(e-y - 1), which is closer to the simulated tails for d = 2 and 3 than the exponential bounds.

We establish the uniform almost sure convergence of the partitioning estimate, which is a histogram-like mean regression function estimate, under ergodic conditions for a stationary and unbounded process. The main application of our results concerns time series analysis and prediction in the Markov processes case.

Recently in the paper by Møller and Zuyev (1996), the following Gamma-type result was established. Given n points of a homogeneous Poisson process defining a random figure, its volume is Γ(n,λ) distributed, where λ is the intensity of the process. In this paper we give an alternative description of the class of random sets for which the Gamma-type results hold. We show that it corresponds to the class of stopping sets with respect to the natural filtration of the point process with certain scaling properties. The proof uses the martingale technique for directed processes, in particular, an analogue of Doob's optional sampling theorem proved in Kurtz (1980). As well as being compact, this approach provides a new insight into the nature of geometrical objects constructed with respect to a Poisson point process. We show, in particular, that in this framework the probability that a point is covered by a stopping set does not depend on whether it is a point of the process or not.

We analyse a class of diffusion models that (i) allows an explicit expression for the likelihood function of discrete time observation, (ii) allows the possibility of heavy-tailed observations, and (iii) allows an analysis of the tails of the increments. The class simply consists of transformed Ornstein–Uhlenbeck processes and is of relevance for heavy-tailed time series. We also treat the question of the existence of an equivalent martingale measure for the class of models considered.

We consider a class of random point and germ-grain processes, obtained using a rather natural weighting procedure. Given a Poisson point process, on each point one places a grain, a (possibly random) compact convex set. Let Ξ be the union of all grains. One can now construct new processes whose density is derived from an exponential of a linear combination of quermass functionals of Ξ. If only the area functional is used, then the area-interaction point process is recovered. New point processes arise if we include the perimeter length functional, or the Euler functional (number of components minus number of holes). The main question addressed by the paper is that of when the resulting point process is well-defined: geometric arguments are used to establish conditions for the point process to be stable in the sense of Ruelle.

The mathematical model we consider here is the classical Bienaymé–Galton–Watson branching process modified with immigration in the state zero.

We study properties of the waiting time to explosion of the supercritical modified process, i.e. that time until all beginning cycles which die out have disappeared. We then derive the expected total progeny of a cycle and show how higher moments can be computed. With a view to applications the main goal is to show that any statistical inference from observed cycle lengths or estimates of total progeny on the fertility rate of the process must be treated with care. As an example we discuss population experiments with trout.

We consider a random measure for which distribution is invariant under the action of a standard transformation group. The reduced moments are defined by applying classical theorems on invariant measure decomposition. We present a general method for constructing unbiased estimators of reduced moments. Several asymptotic results are established under an extension of the Brillinger mixing condition. Examples related to stochastic geometry are given.

We consider the problem of predicting integrals of second order processes whose covariances satisfy some Hölder regularity condition of order α > 0. When α is an odd integer, linear estimators based on regular sampling designs were constructed and asymptotic results for the approximation error were derived. We extend this result to any α > 0. When 2K < α ≤ 2K + 2, K a non-negative integer, we use an appropriate predictor based on the Euler-MacLaurin formula of order K with regular sampling designs. We give the corresponding result for the mean square error.

The filtering problem concerns the estimation of a stochastic process X from its noisy partial information Y. With the notable exception of the linear-Gaussian situation, general optimal filters have no finitely recursive solution. The aim of this work is the design of a Monte Carlo particle system approach to solve discrete time and nonlinear filtering problems. The main result is a uniform convergence theorem. We introduce a concept of regularity and we give a simple ergodic condition on the signal semigroup for the Monte Carlo particle filter to converge in law and uniformly with respect to time to the optimal filter, yielding what seems to be the first uniform convergence result for a particle approximation of the nonlinear filtering equation.

An unbiased stereological estimator for surface area density is derived for gradient surface processes which form a particular class of non-stationary spatial surface processes. Vertical planar sections are used for the estimation. The variance of the estimator is studied and found to be infinite for certain types of surface processes. A modification of the estimator is presented which exhibits finite variance.

We prove a central limit theorem for conditionally centred random fields, under a moment condition and strict positivity of the empirical variance per observation. We use a random normalization, which fits non-stationary situations. The theorem applies directly to Markov random fields, including the cases of phase transition and lack of stationarity. One consequence is the asymptotic normality of the maximum pseudo-likelihood estimator for Markov fields in complete generality.

We derive formulas for the first- and higher-order derivatives of the steady state performance measures for changes in transition matrices of irreducible and aperiodic Markov chains. Using these formulas, we obtain a Maclaurin series for the performance measures of such Markov chains. The convergence range of the Maclaurin series can be determined. We show that the derivatives and the coefficients of the Maclaurin series can be easily estimated by analysing a single sample path of the Markov chain. Algorithms for estimating these quantities are provided. Markov chains consisting of transient states and multiple chains are also studied. The results can be easily extended to Markov processes. The derivation of the results is closely related to some fundamental concepts, such as group inverse, potentials, and realization factors in perturbation analysis. Simulation results are provided to illustrate the accuracy of the single sample path based estimation. Possible applications to engineering problems are discussed.

The spectral factorization problem, i.e. the problem of obtaining all possible MA representations of a process with given autocovariance function, is considered for univariate, d-periodic MA(1) (equivalently, 1-dependent in the second-order sense) processes. The solutions are provided explicitly, and their invertibility properties are investigated. A characterization, in terms of their autocovariance functions, of non-invertible d-periodic 1-dependent processes, extending to the periodic case the traditional unit root condition, is provided.

We consider a class of stationary infinite-order moving average processes with margins in the class of infinitely divisible exponential dispersion models. The processes are constructed by means of the thinning operation of Joe (1996), generalizing the binomial thinning used by McKenzie (1986, 1988) and Al-Osh and Alzaid (1987) for integer-valued time series. As a special case we obtain a class of autoregressive moving average processes that are different from the ARMA models proposed by Joe (1996). The range of possible marginal distributions for the new models is extensive and includes all infinitely divisible distributions with finite moment generating functions, hereunder many known discrete, continuous and mixed distributions.

A common problem in Bayesian object recognition using marked point process models is to produce a point estimate of the true underlying object configuration: the number of objects and the size, location and shape of each object. We use decision theory and the concept of loss functions to design a more reasonable estimator for this purpose, rather than using the common zero-one loss corresponding to the maximum a posteriori estimator. We propose to use the squared Δ-metric of Baddeley (1992) as our loss function and demonstrate that the corresponding optimal Bayesian estimator can be well approximated by combining Markov chain Monte Carlo methods with simulated annealing into a two-step algorithm. The proposed loss function is tested using a marked point process model developed for locating cells in confocal microscopy images.

In this paper, we consider the question of which convergence properties of Markov chains are preserved under small perturbations. Properties considered include geometric ergodicity and rates of convergence. Perturbations considered include roundoff error from computer simulation. We are motivated primarily by interest in Markov chain Monte Carlo algorithms.

In time series analysis, it is well-known that the differencing operator ∇d may transform a non-stationary series, {Z(t)} say, to a stationary one, {W(t)} = ∇d Z(t)}; and there are many procedures for analysing and modelling {Z(t)} which exploit this transformation. Rather differently, Matheron (1973) introduced a set of measures on R n that transform an appropriate non-stationary spatial process to stationarity, and Cressie (1988) then suggested that specialized low-order analogues of these measures, called increment-vectors, be used in time series analysis. This paper develops a general theory of increment-vectors which provides a more powerful transformation tool than mere simple differencing. The methodology gives a handle on the second-moment structure and divergence behaviour of homogeneously non-stationary series which leads to many important applications such as determining the correct degree of differencing, forecasting and interpolation.

Long-range dependence has usually been defined in terms of covariance properties relevant only to second-order stationary processes. Here we provide new definitions, almost equivalent to the original ones in that domain of applicability, which are useful for processes which may not be second-order stationary, or indeed have infinite variances. The ready applicability of this formulation for categorizing the behaviour for various infinite variance models is shown.

It is becoming increasingly recognized that some long series of data can be adequately and parsimoniously modelled by stationary processes with long-range dependence. Some new discrete-time models for long-range dependence or slow decay, defined by their correlation structures, are discussed. The exact power-law correlation structure is examined in detail.

This paper deals with a Bienaymé-Galton-Watson process having a random number of ancestors. Its asymptotic properties are studied when both the number of ancestors and the number of generations tend to infinity. This yields consistent and asymptotically normal estimators of the mean and the offspring distribution of the process. By exhibiting a connection with the BGW process with immigration, all results can be transported to the immigration case, under an appropriate sampling scheme. A key feature of independent interest is a new limit theorem for sums of a random number of random variables, which extends the Gnedenko and Fahim (1969) transfer theorem.