The covariance C(r), r ≥ 0, of a stationary isotropic random closed set Ξ is typically complicated to evaluate. This is the reason that an exponential approximation formula for C(r) has been widely used in the literature, which matches C(0) and C(1)(0), and in many cases also limr→∞C(r). However, for 0 < r < ∞, the accuracy of this approximation is not very high in general. In the present paper, we derive representation formulae for the covariance C(r) and its derivative C(1)(r) using Palm calculus, where r ≥ 0 is arbitrary. As a consequence, an explicit expression is obtained for the second derivative C(2)(0). These results are then used to get a refined exponential approximation for C(r), which additionally matches the second derivative C(2)(0).

A continuum growth model is introduced. The state at time t, St, is a subset of ℝd and consists of a connected union of randomly sized Euclidean balls, which emerge from outbursts at their centre points. An outburst occurs somewhere in St after an exponentially distributed time with expected value |St|-1 and the location of the outburst is uniformly distributed over St. The main result is that, if the distribution of the radii of the outburst balls has bounded support, then St grows linearly and St/t has a nonrandom shape as t → ∞. Due to rotational invariance the asymptotic shape must be a Euclidean ball.

A new class of models for inhomogeneous spatial point processes is introduced. These locally scaled point processes are modifications of homogeneous template point processes, having the property that regions with different intensities differ only by a scale factor. This is achieved by replacing volume measures used in the density with locally scaled analogues defined by a location-dependent scaling function. The new approach is particularly appealing for modelling inhomogeneous Markov point processes. Distance-interaction and shot noise weighted Markov point processes are discussed in detail. It is shown that the locally scaled versions are again Markov and that locally the Papangelou conditional intensity of the new process behaves like that of a global scaling of the homogeneous process. Approximations are suggested that simplify calculation of the density, for example, in simulation. For sequential point processes, an alternative and simpler definition of local scaling is proposed.

The iteration of random tessellations in ℝd is considered, where each cell of an initial tessellation is further subdivided into smaller cells by so-called component tessellations. Sufficient conditions for stationarity and isotropy of iterated tessellations are given. Formulae are derived for the intensities of their facet processes, and for the expected intrinsic volumes of their typical facets. Particular emphasis is put on two special cases: superposition and nesting of tessellations. Bernoulli thinning of iterated tessellations is also considered.

In this article we prove local convergence for a Boolean model of shells conditioned by the noncovering of the origin towards the thick hyperplane Poisson process in the Euclidean space. The existing results of Hall as well as the convergence theorems proved by Paroux or Molchanov concerned the zero-width process and the connected component of the unfilled region of the origin. Our results deal with the convergence in any given window of the space, with the earlier results of Paroux and Molchanov as a corollary.

This paper studies path lengths in random binary search trees under the random permutation model. It is known that the total path length, when properly normalized, converges almost surely to a nondegenerate random variable Z. The limit distribution is commonly referred to as the ‘quicksort distribution’. For the class 𝒜m of finite binary trees with at most m nodes we partition the external nodes of the binary search tree according to the largest tree that each external node belongs to. Thus, the external path length is divided into parts, each part associated with a tree in 𝒜m. We show that the vector of these path lengths, after normalization, converges almost surely to a constant vector times Z.

Let T = (T1, T2,…) be a sequence of real random variables with ∑j=1∞1|Tj|>0 < ∞ almost surely. We consider the following equation for distributions μ: W ≅ ∑j=1∞TjWj, where W, W1, W2,… have distribution μ and T, W1, W2,… are independent. We show that the representation of general solutions is a mixture of certain infinitely divisible distributions. This result can be applied to investigate the existence of symmetric solutions for Tj ≥ 0: essentially under the condition that E ∑j=1∞Tj2 log+Tj2 < ∞, the existence of nontrivial symmetric solutions is exactly determined, revealing a connection with the existence of positive solutions of a related fixed-point equation. Furthermore, we derive results about a special class of canonical symmetric solutions including statements about Lebesgue density and moments.

A new treatment of second-order self-similarity and asymptotic self-similarity for stationary discrete time series is given, based on the fixed points of a renormalisation operator with normalisation factors which are not assumed to be power laws. A complete classification of fixed points is provided, consisting of the fractional noise and one other class. A convenient variance time function approach to process characterisation is used to exhibit large explicit families of processes asymptotic to particular fixed points. A natural, general definition of discrete long-range dependence is provided and contrasted with common alternatives. The closely related discrete form of regular variation is defined, its main properties given, and its connection to discrete self-similarity explained. Folkloric results on long-range dependence are proved or disproved rigorously.

We consider the asymptotic behaviour of a random walk when it exits from a symmetric region of the form {(x, n) :|x| ≤ rnb} as r → ∞. In order to be sure that this actually occurs, we treat only the case where the power b lies in the interval [0,½), and we further assume a condition that prevents the probability of exiting at either boundary tending to 0. Under these restrictions we establish necessary and sufficient conditions for the pth moment of the overshoot to be O(rq), and for the overshoot to be tight, as r → ∞.

This work is devoted to asymptotic properties of singularly perturbed Markov chains in discrete time. The motivation stems from applications in discrete-time control and optimization problems, manufacturing and production planning, stochastic networks, and communication systems, in which finite-state Markov chains are used to model large-scale and complex systems. To reduce the complexity of the underlying system, the states in each recurrent class are aggregated into a single state. Although the aggregated process may not be Markovian, its continuous-time interpolation converges to a continuous-time Markov chain whose generator is a function determined by the invariant measures of the recurrent states. Sequences of occupation measures are defined. A mean square estimate on a sequence of unscaled occupation measures is obtained. Furthermore, it is proved that a suitably scaled sequence of occupation measures converges to a switching diffusion.

Continuous-space-time branching processes (CSBP) are investigated in order to model random energy cascades. CSBPs are based on spectrally positive Lévy processes and, as such, are characterized by their corresponding Laplace exponents. Special emphasis is put on the CSBPs of Feller, Lamperti and Neveu and on their Poisson point process representations. The Neveu model (either supercritical or subcritical) is of particular interest in physics for its connection with the random energy model of Derrida, as revisited by Ruelle. Exploiting some connections between the partition functions of energy and the Poisson-Dirichlet distributions of Pitman and Yor, some information on the zero-temperature limit is extracted. Finally, for the subcritical versions of the three models, we compute the distribution of some of their interesting features: extinction time and probability, area under the profile (total energy) and width (maximal energy).

This paper studies the first passage times to flat boundaries for a double exponential jump diffusion process, which consists of a continuous part driven by a Brownian motion and a jump part with jump sizes having a double exponential distribution. Explicit solutions of the Laplace transforms, of both the distribution of the first passage times and the joint distribution of the process and its running maxima, are obtained. Because of the overshoot problems associated with general jump diffusion processes, the double exponential jump diffusion process offers a rare case in which analytical solutions for the first passage times are feasible. In addition, it leads to several interesting probabilistic results. Numerical examples are also given. The finance applications include pricing barrier and lookback options.

Generalizing the standard frailty models of survival analysis, we propose to model frailty as a weighted Lévy process. Hence, the frailty of an individual is not a fixed quantity, but develops over time. Formulae for the population hazard and survival functions are derived. The power variance function Lévy process is a prominent example. In many cases, notably for compound Poisson processes, quasi-stationary distributions of survivors may arise. Quasi-stationarity implies limiting population hazard rates that are constant, in spite of the continual increase of the individual hazards. A brief discussion is given of the biological relevance of this finding.