In this paper we calculate posterior distributions associated with a version of the Poisson–Dirichlet distribution called the GEM. The GEM has been shown (by several authors) to be the limiting stationary distribution for allele frequencies listed in age order associated with the neutral infinite alleles model. In view of this result, we use our posterior distributions to calculate Bayes estimators for the frequency of the oldest allele given a sample.

This paper presents the large-buffer asymptotics for a multiplexer which serves N types of heterogeneous sessions which have long-tailed session lengths. Specifically, the model considered is that sessions of type i ∈ {1,…,N} arrive as a Poisson process with rate λi. Each type of session (independently) remains active for a random duration, say τi, where P(τi > x) ~ αix-(1 + βi) for positive numbers αi and βi. While active, a session transmits at a rate ri. Under the assumption that the average load ρ = ∑Ni=1riλiE[τi] < C, where C denotes the server capacity, we show that both the tail distribution of the stationary buffer content and the loss asymptotics in finite buffers of size z behave approximately as z-κJ0, where κJ0 depends not only on the βi but also on the transmission rates ri; it is the ratio of βi to ri which determines κJ0. When specialized to the homogeneous case, i.e., when ri=r and βi = β for all i, the result coincides with results reported in the literature which have been shown under more restrictive hypotheses. Finally, it is a simple observation that light-tailed sessions only have the effect of reducing the available capacity for long-tailed sessions, but do not contribute otherwise to the definition of κJ0.

For a random closed set X and a compact observation window W the mean coverage fraction of W can be estimated by measuring the area of W covered by X. Jensen and Gundersen, and Baddeley and Cruz-Orive described cases where a point counting estimator is more efficient than area measurement. We give two other examples, where at first glance unnatural estimators are not only better than the area measurement but by Grenander's Theorem have minimal variance. Whittle's Theorem is used to show that the point counting estimator in the original Jensen-Gundersen paradox is optimal for large randomly translated discs.

This note is concerned with the continuing misconception that a discrete-time network of queues, with independent customer routing and the number of arrivals and services in a time interval following geometric and truncated geometric distributions respectively, has a product-form equilibrium distribution.

Detailed balance of a chemical reaction network can be defined in several different ways. Here we investigate the relationship among four types of detailed balance conditions: deterministic, stochastic, local, and zero-order local detailed balance. We show that the four types of detailed balance are equivalent when different reactions lead to different species changes and are not equivalent when some different reactions lead to the same species change. Under the condition of local detailed balance, we further show that the system has a global potential defined over the whole space, which plays a central role in the large deviation theory and the Freidlin–Wentzell-type metastability theory of chemical reaction networks. Finally, we provide a new sufficient condition for stochastic detailed balance, which is applied to construct a class of high-dimensional chemical reaction networks that both satisfies stochastic detailed balance and displays multistability.

There are two ways of speeding up Markov chain Monte Carlo algorithms: (a) construct more complex samplers that use gradient and higher-order information about the target and (b) design a control variate to reduce the asymptotic variance. While the efficiency of (a) as a function of dimension has been studied extensively, this paper provides the first results linking the efficiency of (b) with dimension. Specifically, we construct a control variate for a d-dimensional random walk Metropolis chain with an independent, identically distributed target using the solution of the Poisson equation for the scaling limit in [30]. We prove that the asymptotic variance of the corresponding estimator is bounded above by a multiple of $\log(d)/d$ over the spectral gap of the chain. The proof hinges on large deviations theory, optimal Young’s inequality and Berry–Esseen-type bounds. Extensions of the result to non-product targets are discussed.

‘Small worlds’ are large systems in which any given node has only a few connections to other points, but possessing the property that all pairs of points are connected by a short path, typically logarithmic in the number of nodes. The use of random walks for sampling a uniform element from a large state space is by now a classical technique; to prove that such a technique works for a given network, a bound on the mixing time is required. However, little detailed information is known about the behaviour of random walks on small-world networks, though many predictions can be found in the physics literature. The principal contribution of this paper is to show that for a famous small-world random graph model known as the Newman-Watts small-world model, the mixing time is of order log2n. This confirms a prediction of Richard Durrett [5, page 22], who proved a lower bound of order log2n and an upper bound of order log3n.

The Doss-Sussmann (DS) approach is used for uniform simulation of the Cox-Ingersoll-Ross (CIR) process. The DS formalism allows us to express trajectories of the CIR process through solutions of some ordinary differential equation (ODE) depending on realizations of a Wiener process involved. By simulating the first-passage times of the increments of the Wiener process to the boundary of an interval and solving the ODE, we uniformly approximate the trajectories of the CIR process. In this respect special attention is payed to simulation of trajectories near 0. From a conceptual point of view the proposed method gives a better quality of approximation (from a pathwise point of view) than standard, even exact, simulation of the stochastic differential equation at some deterministic time grid.

The number of returns to the origin for a simple random walk starting and ending at the origin is reconsidered; closed-form expressions in μ are given for μ3 and μ4.

Consider a set of n vertices, where each vertex has a location in $\mathbb{R}^d$ that is sampled uniformly from the unit cube in $\mathbb{R}^d$, and a weight associated to it. Construct a random graph by placing edges independently for each vertex pair with a probability that is a function of the distance between the locations and the vertex weights.

Under appropriate integrability assumptions on the edge probabilities that imply sparseness of the model, after appropriately blowing up the locations, we prove that the local limit of this random graph sequence is the (countably) infinite random graph on $\mathbb{R}^d$ with vertex locations given by a homogeneous Poisson point process, having weights which are independent and identically distributed copies of limiting vertex weights. Our set-up covers many sparse geometric random graph models from the literature, including geometric inhomogeneous random graphs (GIRGs), hyperbolic random graphs, continuum scale-free percolation, and weight-dependent random connection models.

We prove that the limiting degree distribution is mixed Poisson and the typical degree sequence is uniformly integrable, and we obtain convergence results on various measures of clustering in our graphs as a consequence of local convergence. Finally, as a byproduct of our argument, we prove a doubly logarithmic lower bound on typical distances in this general setting.

We present a numerical method to compute the survival function and the moments of the exit time for a piecewise-deterministic Markov process (PDMP). Our approach is based on the quantization of an underlying discrete-time Markov chain related to the PDMP. The approximation we propose is easily computable and is even flexible with respect to the exit time we consider. We prove the convergence of the algorithm and obtain bounds for the rate of convergence in the case of the moments. We give an academic example and a model from the reliability field to illustrate the results of the paper.

In this paper we develop an asymptotic theory of aggregated linear processes, and determine in particular the limit distribution of a large class of linear and nonlinear functionals of such processes. Given a sample {Y1(N),…,Yn(N)} of the normalized N-fold aggregated process, we describe the limiting behavior of statistics TN,n= TN,n(Y1(N),…, Yn(N)) in both of the cases n/N(n) → 0 and N(n)/n → 0, assuming either a ‘limiting long- or short-memory’ condition on the underlying linear process.

We consider a controlled i.i.d. process, where several i.i.d. sources are sampled sequentially. At each time instant, a controller determines from which source to obtain the next sample. Any causal sampling policy, possibly history-dependent, may be employed. The purpose is to characterize the extremal large deviations of the sample mean, namely to obtain asymptotic rate bounds (similar to and extending Cramér's theorem) which hold uniformly over all sampling policies. Lower and upper bounds are obtained, and it is shown that in many (but not all) cases stationary sampling policies are sufficient to obtain the extremal large deviations rates. These results are applied to a hypothesis testing problem, where data samples may be sequentially chosen from several i.i.d. sources (representing different types of experiments). The analysis provides asymptotic estimates for the error probabilities, corresponding both to optimal and to worst-case sampling policies.

Let (φ(Xn))n be a function of a finite-state Markov chain (Xn)n. In this article, we investigate the conditions under which the random variables φ(n) have the same distribution as Yn (for every n), where (Yn)n is a Markov chain with fixed transition probability matrix. In other words, for a deterministic function φ, we investigate the conditions under which (Xn)n is weakly lumpable for the state vector. We show that the set of all probability distributions of X0, such that (Xn)n is weakly lumpable for the state vector, can be finitely generated. The connections between our definition of lumpability and the usual one (i.e. as the proportional dynamics property) are discussed.

We study failure rate monotonicity and generalised convex transform stochastic ordering properties of random variables, with an emphasis on applications. We are especially interested in the effect of a tail-weight iteration procedure to define distributions, which is equivalent to the characterisation of moments of the residual lifetime at a given instant. For the monotonicity properties, we are mainly concerned with hereditary properties with respect to the iteration procedure providing counterexamples showing either that the hereditary property does not hold or that inverse implications are not true. For the stochastic ordering, we introduce a new criterion, based on the analysis of the sign variation of a suitable function. This criterion is then applied to prove ageing properties of parallel systems formed with components that have exponentially distributed lifetimes.

Under a fourth-order moment condition on the branching and a second-order moment condition on the immigration mechanisms, we show that an appropriately scaled projection of a supercritical and irreducible continuous-state and continuous-time branching process with immigration on certain left non-Perron eigenvectors of the branching mean matrix is asymptotically mixed normal. With an appropriate random scaling, under some conditional probability measure, we prove asymptotic normality as well. In the case of a non-trivial process, under a first-order moment condition on the immigration mechanism, we also prove the convergence of the relative frequencies of distinct types of individuals on a suitable event; for instance, if the immigration mechanism does not vanish, then this convergence holds almost surely.

The basis of this paper is the elementary observation that the n-step descendant distribution of any Galton–Watson process satisfies a discrete Smoluchowski coagulation equation with multiple coalescence. Using this we obtain simple necessary and sufficient criteria for the convergence of scaling limits of critical Galton–Watson processes in terms of scaled family-size distributions and a natural notion of convergence of Lévy triples. Our results provide a clear and natural interpretation, and an alternate proof, of the fact that the Lévy jump measure of certain continuous-state branching processes (CSBPs) satisfies a generalized Smoluchowski equation. (This result was previously proved by Bertoin and Le Gall (2006).) Our analysis shows that the nonlinear scaling dynamics of CSBPs become linear and purely dilatational when expressed in terms of the Lévy triple associated with the branching mechanism. We prove a continuity theorem for CSBPs in terms of the associated Lévy triples, and use our scaling analysis to prove the existence of universal critical Galton–Watson processes and CSBPs analogous to Doeblin's `universal laws'. Namely, these universal processes generate all possible critical and subcritical CSBPs as subsequential scaling limits. Our convergence results rely on a natural topology for Lévy triples and a continuity theorem for Bernstein transforms (Laplace exponents) which we develop in a self-contained appendix.

Suppose segments are distributed at random along a line, their locations being determined by a Poisson process. In the case where segment length is fixed, we compare efficiencies of several different estimates of Poisson intensity. The case of random segment length is also considered, and there we study estimation procedures based on empiric properties. The one-dimensional mosaic may be viewed as an M/G/∞ queue.