Consider the process which starts with N ≥ 3 distinct points on ℝd, and fix a positive integer K < N. Of the total N points keep those N - K which minimize the energy amongst all the possible subsets of size N - K, and then replace the removed points by K independent and identically distributed points sampled according to some fixed distribution ζ. Repeat this process ad infinitum. We obtain various quite nonrestrictive conditions under which the set of points converges to a certain limit. This is a very substantial generalization of the `Keynesian beauty contest process' introduced in Grinfeld et al. (2015), where K = 1 and the distribution ζ was uniform on the unit cube.

Let V be an n-set, and let X be a random variable taking values in the power-set of V. Suppose we are given a sequence of random coupons $X_1, X_2, \ldots $ , where the $X_i$ are independent random variables with distribution given by X. The covering time T is the smallest integer $t\geq 0$ such that $\bigcup_{i=1}^t X_i=V$ . The distribution of T is important in many applications in combinatorial probability, and has been extensively studied. However the literature has focused almost exclusively on the case where X is assumed to be symmetric and/or uniform in some way.

In this paper we study the covering time for much more general random variables X; we give general criteria for T being sharply concentrated around its mean, precise tools to estimate that mean, as well as examples where T fails to be concentrated and when structural properties in the distribution of X allow for a very different behaviour of T relative to the symmetric/uniform case.

Whittle enunciated an important reduction principle in dynamic programming when he showed that under certain conditions optimal strategies for Markov decision processes (MDPs) placed in parallel to one another take actions in a way which is consistent with the optimal strategies for the individual MDPs. However, the necessary and sufficient conditions given by Whittle are by no means always satisfied. We explore the status of this computationally attractive reduction principle when these conditions fail.

The theory of geometric probabilities is concerned with randomly generated geometric objects. The aim is to compute probabilities of certain geometric events or distributions of random variables defined in a geometric way. Very often the computation even of simple expectations is too difficult, and one has to be satisfied with establishing estimates and, if possible, sharp inequalities. In geometric probabilities, convex sets play a prominent role, since often the convexity assumptions simplify the situation considerably. Extremal problems for geometric probabilities involving convex bodies can sometimes be attacked successfully by using suitable integral-geometric transformations and then applying classical inequalities from the geometry of convex bodies, or known methods for obtaining such inequalities. Examples of such results are the topic of this paper.

Let n points be randomly and independently placed in Rd according to a common probability law. It is known that the expected volume for the convex hull of these points, in the cases where n - d ≥ 2 and even, is related linearly to expected volumes of the convex hulls for j points, j < n. We show that similar identities for these volumes hold almost surely - and in contexts where independence and communality of law do not apply. New geometric and topological identities developed here provide a foundation for this result.

Dependence between coordinate extremes is a key factor in any multivariate risk assessment. Hence, it is of interest to know whether the components of a given multivariate random vector exhibit asymptotic independence or asymptotic dependence. In the latter case the structure of the asymptotic dependence has to be clarified. In the multivariate setting it is common to have an explicit form of the density rather than the distribution function. In this paper we therefore give criteria for asymptotic dependence in terms of the density. We consider distributions with light tails and restrict attention to continuous unimodal densities defined on the whole space or on an open convex cone. For simplicity, the density is assumed to be homothetic: all level sets have the same shape. Balkema and Nolde (2010) contains conditions on the shape which guarantee asymptotic independence. The situation for asymptotic dependence, treated in the present paper, is more delicate.

Dimension reduction of multivariate data was developed by Y. Guan for point processes with Gaussian random fields as covariates. The generalization to fibre and surface processes is straightforward. In inverse regression methods, we suggest slicing based on geometrical marks. An investigation of the properties of this method is presented in simulation studies of random marked sets. In a refined model for dimension reduction, the second-order central subspace is analyzed in detail. A real data pattern is tested for independence of a covariate.

Let {Xn}n∈ℕ be a Markov chain on a measurable space with transition kernel P, and let The Markov kernel P is here considered as a linear bounded operator on the weighted-supremum space associated with V. Then the combination of quasicompactness arguments with precise analysis of eigenelements of P allows us to estimate the geometric rate of convergence ρV(P) of {Xn}n∈ℕ to its invariant probability measure in operator norm on A general procedure to compute ρV(P) for discrete Markov random walks with identically distributed bounded increments is specified.

For an initially empty M/M/1 queue, it is shown that the transform of the expectation of the number of departures in the interval (0, t] is invariant under an interchange of arrival and service rates. However, in the GI/G/1 queue with an initial single customer, the corresponding transform does not have this symmetric property.

We study the growth of a population of bacteria in a dynamical hostile environment corresponding to the immune system of the colonized organism. The immune cells evolve as subcritical open clusters of oriented percolation and are perpetually reinforced by an immigration process, while the bacteria try to grow as a supercritical oriented percolation in the remaining empty space. We prove that the population of bacteria grows linearly when it survives. From this perspective, we build general tools to study dependent oriented percolation models issued from renormalization processes.

Most finite spatial point process models specified by a density are locally stable, implying that the Papangelou intensity is bounded by some integrable function β defined on the space for the points of the process. It is possible to superpose a locally stable spatial point process X with a complementary spatial point process Y to obtain a Poisson process X ⋃ Y with intensity function β. Underlying this is a bivariate spatial birth-death process (Xt, Yt) which converges towards the distribution of (X, Y). We study the joint distribution of X and Y, and their marginal and conditional distributions. In particular, we introduce a fast and easy simulation procedure for Y conditional on X. This may be used for model checking: given a model for the Papangelou intensity of the original spatial point process, this model is used to generate the complementary process, and the resulting superposition is a Poisson process with intensity function β if and only if the true Papangelou intensity is used. Whether the superposition is actually such a Poisson process can easily be examined using well-known results and fast simulation procedures for Poisson processes. We illustrate this approach to model checking in the case of a Strauss process.

From one point of view this paper adds to a previous formulation of the best-choice problem (Petruccelli (1981)) the possibility that the number of available observations, rather than being known, is a bounded random variable N with known distribution. From another perspective, it expands the formulations of Presman and Sonin (1972) and Rasmussen and Robbins (1975) to include recall and uncertainty of selection of observations. The behaviour of optimal stopping rules is examined under various assumptions on the general model. For optimal stopping rules and their probabilities of best choice relations are obtained between the bounded and unbounded N cases. Two particular classes of stopping rules which generalize the s(r) rules of Rasmussen and Robbins (1975) are considered in detail.

This paper investigates properties of the class of graphs based on exchangeable point processes. We provide asymptotic expressions for the number of edges, number of nodes, and degree distributions, identifying four regimes: (i) a dense regime, (ii) a sparse, almost dense regime, (iii) a sparse regime with power-law behaviour, and (iv) an almost extremely sparse regime. We show that, under mild assumptions, both the global and local clustering coefficients converge to constants which may or may not be the same. We also derive a central limit theorem for subgraph counts and for the number of nodes. Finally, we propose a class of models within this framework where one can separately control the latent structure and the global sparsity/power-law properties of the graph.

We introduce a formal limit, which we refer to as a fluid limit, of scaled stochastic models for a cache managed with the least-recently-used algorithm when requests are issued according to general stochastic point processes. We define our fluid limit as a superposition of dependent replications of the original system with smaller item sizes when the number of replications approaches ∞. We derive the average probability that a requested item is not in a cache (average miss probability) in the fluid limit. We show that, when requests follow inhomogeneous Poisson processes, the average miss probability in the fluid limit closely approximates that in the original system. Also, we compare the asymptotic characteristics, as the cache size approaches ∞, of the average miss probability in the fluid limit to those in the original system.

Consider a supercritical Crump‒Jagers process in which all births are at integer times (the lattice case). Let μ̂(z) be the generating function of the intensity of the offspring process, and consider the complex roots of μ̂(z)=1. The root of smallest absolute value is e-α=1∕m, where α>0 is the Malthusian parameter; let γ* be the root of second smallest absolute value. Subject to some technical conditions, the second-order fluctuations of the age distribution exhibit one of three types of behaviour: (i) when γ*>e-α∕2=m-1∕2, they are asymptotically normal; (ii) when γ*=e-α∕2, they are still asymptotically normal, but with a larger variance; and (iii) when γ*<e-α∕2, the fluctuations are in general oscillatory and (degenerate cases excluded) do not converge in distribution. This trichotomy is similar to what has been observed in related situations, such as some other branching processes and for Pólya urns. The results lead to a symbolic calculus describing the limits. The asymptotic results also apply to the total of other (random) characteristics of the population.

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.

Learning from Matheron's representation (1973), and using the increment vector (PIV) methodology introduced by Cressie (1988) and developed by Chen and Anderson (1994), this paper presents a theory for the representation and decomposition of integrated stationary time series and gives some applications.

Here a response of a non-linear oscillator of the Liénard type with a large parameter α ≥ 0 is formulated as a solution of a two-dimensional stochastic differential equation with mean-field of the McKean type. This solution is governed by a special form of the Fokker–Planck equation such as the Smoluchowski–Kramers equation, which is an equation of motion for distribution functions in position and velocity space describing the Brownian motion of particles in an external field. By a change of time and displacement we find that the velocity process converges to a one-dimensional Ornstein–Uhlenbeck process as α →∞.

The problem treated is that of controlling a process with values in [0, a]. The non-anticipative controls (µ(t), σ(t)) are selected from a set C(x) whenever X(t–) = x and the non-decreasing process A(t) is chosen by the controller subject to the condition where y is a constant representing the initial amount of fuel. The object is to maximize the probability that X(t) reaches a. The optimal process is determined when the function has a unique minimum on [0, a] and satisfies certain regularity conditions. The optimal process is a combination of ‘timid play' in which fuel is used gradually in the form of local time at 0, and ‘bold play' in which all the fuel is used at once.