We consider the inverse problem of estimating the input random field in a stochastic integral equation relating two random fields. The purpose of this paper is to present an approach to this problem using a Riesz-based or orthonormal-based series expansion of the input random field with uncorrelated random coefficients. We establish conditions under which the input series expansion induces (via the integral equation) a Riesz-based or orthonormal-based series expansion for the output random field. The estimation problem is studied considering two cases, depending on whether data are available from either the output random field alone, or from both the input and output random fields. Finally, we discuss this approach in the case of transmissivity estimation from piezometric head data, which was the original motivation of this work.

The inverse absorption distribution is shown to be a q-Pascal analogue of the Kemp and Kemp (1991) q-binomial distribution. The probabilities for the direct absorption distribution are obtained via the inverse absorption probabilities and exact expressions for its first two factorial moments are derived using q-series transformations of its probability generating function. Alternative models for the distribution are given.

We introduce the concept of the absorbing process for analysing a state process. Our aim is to show the existence of the absorbing process with probability one. This process is shown to be stationary, asymptotically stationary, periodic or a.m.s., if the input distribution has such properties. The real process is absorbed into this process so that its stability and some other properties are easily derived.

Criteria are determined for the variance to mean ratio to be greater than one (over-dispersed) or less than one (under-dispersed). This is done for random variables which are functions of a Markov chain in continuous time, and for the counts in a simple point process on the line. The criteria for the Markov chain are in terms of the infinitesimal generator and those for the point process in terms of the conditional intensity. Examples include a conjecture of Faddy (1994). The case of time-reversible point processes is particularly interesting, and here underdispersion is not possible. In particular, point processes which arise from Markov chains which are time-reversible, have finitely many states and are irreducible are always overdispersed.

The multitype discrete time indecomposable branching process with immigration is considered. Using a martingale approach a limit theorem is proved for such processes when the totality of immigrating individuals at a given time depends on evolution of the processes generating by previously immigrated individuals. Corollaries of the limit theorem are obtained for the cases of finite and infinite second moments of offspring distribution in critical processes.

The influence of bivariate extremal dependence on the limiting behaviour of the concomitant of the largest order statistic is examined. Our approach is to fix the marginal distributions and derive a general tail characterisation of the joint survivor function. From this, we identify the normalisation required to obtain the limiting distribution of the concomitant of the largest order statistic, obtain its tail form, and investigate the limiting probability that the vector of componentwise maxima occurs as an observation of the bivariate process. The results are illustrated for a range of extremal dependence forms.

A new class of Gibbsian models with potentials associated with the connected components or homogeneous parts of images is introduced. For these models the neighbourhood of a pixel is not fixed as for Markov random fields, but is given by the components which are adjacent to the pixel. The relationship to Markov random fields and marked point processes is explored and spatial Markov properties are established. Extensions to infinite lattices are also studied, and statistical inference problems including geostatistical applications and statistical image analysis are discussed. Finally, simulation studies are presented which show that the models may be appropriate for a variety of interesting patterns, including images exhibiting intermediate degrees of spatial continuity and images of objects against background.

The risk reserve process of an insurance company within a deteriorating Markov-modulated environment is considered. The company invests its capital with interest rate α; the premiums and claims are increasing with rates β and γ. The problem of stopping the process at a random time which maximizes the expected net gain in order to calculate new premiums is investigated. A semimartingale representation of the risk reserve process yields, under certain conditions, an explicit solution of the problem.

The input of a multiserver loss system is assumed to be a periodic random marked point process which has, with probability one, infinitely many construction points. It is shown that, independently of the initial distribution, there exists a unique periodic process modeling the periodic steady-state behaviour of the loss system. In addition, practical sufficient conditions for the existence of enough construction points are derived.

In this paper we consider an explicit solution of an optimal stopping problem arising in connection with a dice game. An optimal stopping rule and the maximum expected reward in this problem can easily be computed by means of the distributions involved and the specific rules of the game

Asymptotic expansions are obtained for the distribution function of a studentized estimator of the offspring mean sequence in an array branching process with immigration. The expansion result is shown to hold in a test function topology. As an application of this result, it is shown that the bootstrapping distribution of the estimator of the offspring mean in a sub-critical branching process with immigration also admits the same expansion (in probability). From these considerations, it is concluded that the bootstrapping distribution provides a better approximation asymptotically than the normal distribution.

We generalise the work of Cramér and Leadbetter, Ylvisaker and Ito on the level crossings of a stationary Gaussian process to multivariate processes. Necessary and sufficient conditions for the existence of the expected number of crossings E(C) of the boundary of a region of ℝp by a stationary vector stochastic process are obtained, along with an explicit formula for E(C) in the Gaussian case. A rigorous proof of Belyaev's integral formula for the factorial moments of the number of exits from a region of ℝp is given for a class of processes which includes Gaussian processes having a finite second order spectral moment matrix. Applications to χ2 processes are briefly considered.

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.

Let Tr be the first time at which a random walk Sn escapes from the strip [-r,r], and let |STr|-r be the overshoot of the boundary of the strip. We investigate the order of magnitude of the overshoot, as r → ∞, by providing necessary and sufficient conditions for the ‘stability’ of |STr|, by which we mean that |STr|/r converges to 1, either in probability (weakly) or almost surely (strongly), as r → ∞. These also turn out to be equivalent to requiring only the boundedness of |STr|/r, rather than its convergence to 1, either in the weak or strong sense, as r → ∞. The almost sure characterisation turns out to be extremely simple to state and to apply: we have |STr|/r → 1 a.s. if and only if EX2 < ∞ and EX = 0 or 0 < |EX| ≤ E|X| < ∞. Proving this requires establishing the equivalence of the stability of STr with certain dominance properties of the maximum partial sum Sn* = max{|Sj|: 1 ≤ j ≤ n} over its maximal increment.

The Brownian density process is a Gaussian distribution-valued process. It can be defined either as a limit of a functional over a Poisson system of independent Brownian particles or as a solution of a stochastic partial differential equation with respect to Gaussian martingale measure. We show that, with an appropriate change in the initial distribution of the infinite particle system, the limiting density process is non-Gaussian and it solves a stochastic partial differential equation where the initial measure and the driving measure are non-Gaussian, possibly having infinite second moment.

We consider a stationary germ-grain model Ξ with convex and compact grains and the distance r(x) from x ε ℝd to Ξ. For almost all points x ε ℝd there exists a unique point p(x) in the boundary of Ξ such that r(x) is the length of the vector x-p(x), which is called the spherical contact vector at x. In this paper we relate the distribution of the spherical contact vector to the times it takes a typical boundary point of Ξ to hit another grain if all grains start growing at the same time and at the same speed. The notion of a typical point is made precise by using the generalized curvature measures of Ξ. The result generalizes a well known formula for the Boolean model. Specific examples are discussed in detail.

For a stationary point process X of sets in the convex ring in ℝd, a relation is given between the mean particles of the section process X ∩ E (where E varies through the set of k-dimensional subspaces in ℝd) and a mean particle of X. In particular, it is shown that the mean bodies of all planar sections of X determine the Blaschke body of X and hence the mean normal distribution of X.

We give an alternative proof of a point-process version of the FKG–Holley–Preston inequality which provides a sufficient condition for stochastic domination of probability measures, and for positive correlations of increasing functions.

We formulate stochastic indicator parameters that characterize pollution levels in geographical regions with heterogeneous contaminant distributions. The indicator parameters are expressed in terms of the random fields representing the contaminant distributions and the critical threshold level specified by health and environmental standards. Certain theoretical results are proven regarding univariate and bivariate indicator parameters. The analytical expressions obtained are general and can be used in practice for various types of contaminant distributions. A test of ergodicity-breaking is suggested for scientific and engineering applications in terms of the indicator parameters. Fractal characteristics of the indicator parameters are discussed. The effects of modelling and observation scale on exceedance contamination analysis are examined. Indicator random field parameters are studied on both continuum and lattice domains using analytical means and numerical simulations.