We consider the problem of estimating a latent point process, given the realization of another point process. We establish an expression for the conditional distribution of a latent Poisson point process given the observation process when the transformation from the latent process to the observed process includes displacement, thinning, and augmentation with extra points. Our original analysis is based on an elementary and self-contained random measure theoretic approach. This simplifies and complements previous derivations given in Mahler (2003), and Singh, Vo, Baddeley and Zuyev (2009).

For stationary Poisson hyperplane tessellations in d-dimensional Euclidean space and a dimension k ∈ {1, …, d}, we investigate the typical k-face and the weighted typical k-face (weighted by k-dimensional volume), without isotropy assumptions on the tessellation. The case k = d concerns the previously studied typical cell and zero cell, respectively. For k < d, we first find the conditional distribution of the typical k-face or weighted typical k-face, given its direction. Then we investigate how the shapes of the faces are influenced by assumptions of different types: either via containment of convex bodies of given volume (including a new result for k = d), or, for weighted typical k-faces, in the spirit of D. G. Kendall's asymptotic problem, suitably generalized. In all these results on typical or weighted typical k-faces with given direction space L, the Blaschke body of the section process of the underlying hyperplane process with L plays a crucial role.

We consider the following stochastic model for a mobile service scenario. Consider a stationary Poisson process in Rd, with its points radially ordered with respect to the origin (the anchor); if d = 2, the points may correspond to locations of, e.g. restaurants. A user, with a location different from the origin, asks for the location of the first Poisson point and keeps asking for the location of the next Poisson point until the first time that he/she can be completely certain that he/she knows which Poisson point is his/her nearest neighbour. This waiting time is the communication cost, while the inferred privacy region is a random set obtained by an adversary who only knows the anchor and the points received from the server, where the adversary ‘does the best’ to infer the possible locations of the user. Probabilistic results related to the communication cost and the inferred privacy region are established for any dimension d ≥ 1. Furthermore, special results when d = 1 and particularly when d = 2 are derived.

We introduce a new 1-dependent percolation model to describe and analyze the spread of an epidemic on a general directed and locally finite graph. We assign a two-dimensional random weight vector to each vertex of the graph in such a way that the weights of different vertices are independent and identically distributed, but the two entries of the vector assigned to a vertex need not be independent. The probability for an edge to be open depends on the weights of its end vertices, but, conditionally on the weights, the states of the edges are independent of each other. In an epidemiological setting, the vertices of a graph represent the individuals in a (social) network and the edges represent the connections in the network. The weights assigned to an individual denote its (random) infectivity and susceptibility, respectively. We show that one can bound the percolation probability and the expected size of the cluster of vertices that can be reached by an open path starting at a given vertex from above by the corresponding quantities for independent bond percolation with a certain density; this generalizes a result of Kuulasmaa (1982). Many models in the literature are special cases of our general model.

Coalescents with multiple collisions (also called Λ-coalescents or simple exchangeable coalescents) are used as models of genealogies. We study a new class of Markovian coalescent processes connected to a population model with immigration. Consider an infinite population with immigration labelled at each generation by N := {1, 2, …}. Some ancestral lineages cannot be followed backwards after some time because their ancestor is outside the population. The individuals with an immigrant ancestor constitute a distinguished family and we define exchangeable distinguished coalescent processes as a model for genealogy with immigration, focusing on simple distinguished coalescents, i.e. such that when a coagulation occurs all the blocks involved merge as a single block. These processes are characterized by two finite measures on [0, 1] denoted by M = (Λ0, Λ1). We call them M-coalescents. We show by martingale arguments that the condition of coming down from infinity for the M-coalescent coincides with that obtained by Schweinsberg for the Λ-coalescent. In the same vein as Bertoin and Le Gall, M-coalescents are associated with some stochastic flows. The superprocess embedded can be viewed as a generalized Fleming-Viot process with immigration. The measures Λ0 and Λ1 respectively specify the reproduction and the immigration. The coming down from infinity of the M-coalescent will be interpreted as the initial types extinction: after a certain time all individuals are immigrant children.

Organisms adapt to fluctuating environments by regulating their dynamics, and by adjusting their phenotypes to environmental changes. We model population growth using multitype branching processes in random environments, where the offspring distribution of some organism having trait t ∈ in environment e ∈ ε is given by some (fixed) distribution ϒt,e on ℕ. Then, the phenotypes are attributed using a distribution (strategy) πt,e on the trait space . We look for the optimal strategy πt,e, t ∈ , e ∈ ε, maximizing the net growth rate or Lyapounov exponent, and characterize the set of optimal strategies. This is considered for various models of interest in biology: hereditary versus nonhereditary strategies and strategies involving or not involving a sensing mechanism. Our main results are obtained in the setting of nonhereditary strategies: thanks to a reduction to simple branching processes in a random environment, we derive an exact expression for the net growth rate and a characterization of optimal strategies. We also focus on typical genealogies, that is, we consider the problem of finding the typical lineage of a randomly chosen organism.

Mixed-level orthogonal arrays are basic structures in experimental design. We develop three algorithms that compute Rao- and Gilbert-Varshamov-type bounds for mixed-level orthogonal arrays. The computational complexity of the terms involved in the original combinatorial representations of these bounds can grow fast as the parameters of the arrays increase and this justifies the construction of these algorithms. The first is a recursive algorithm that computes the bounds exactly, the second is based on an asymptotic analysis, and the third is a simulation algorithm. They are all based on the representation of the combinatorial expressions that appear in the bounds as expectations involving a symmetric random walk. The Markov property of the underlying random walk gives the recursive formula to compute the expectations. A large deviation (LD) analysis of the expectations provides the asymptotic algorithm. The asymptotically optimal importance sampling (IS) of the same expectation provides the simulation algorithm. Both the LD analysis and the construction of the IS algorithm use a representation of these problems as a sequence of stochastic optimal control problems converging to a limit calculus of a variations problem. The construction of the IS algorithm uses a recently discovered method of using subsolutions to the Hamilton-Jacobi-Bellman equations associated with the limit problem.

We provide necessary and sufficient conditions for the asymptotic normality of Nn, the number of records among the first n observations from a sequence of independent and identically distributed random variables, with general distribution F. In the case of normality we identify the centering and scaling sequences. Also, we characterize distributions for which the limit is not normal in terms of their discrete and continuous components.

We supply some relations that establish intertwining from duality and give a probabilistic interpretation. This is carried out in the context of discrete Markov chains, fixing up the background of previous relations established for monotone chains and their Siegmund duals. We revisit the duality for birth-and-death chains and the nonneutral Moran model, and we also explore the duality relations in an ultrametric-type dual that extends the Siegmund kernel. Finally, we discuss the sharp dual, following closely the Diaconis-Fill study.

Max-stable random fields play a central role in modeling extreme value phenomena. We obtain an explicit formula for the conditional probability in general max-linear models, which include a large class of max-stable random fields. As a consequence, we develop an algorithm for efficient and exact sampling from the conditional distributions. Our method provides a computational solution to the prediction problem for spectrally discrete max-stable random fields. This work offers new tools and a new perspective to many statistical inference problems for spatial extremes, arising, for example, in meteorology, geology, and environmental applications.

In a fork-join network each incoming job is split into K tasks and the K tasks are simultaneously assigned to K parallel service stations for processing. For the distributions of response times and queue lengths of fork-join networks, no explicit formulae are available. Existing methods provide only analytic approximations for the response time and the queue length distributions. The accuracy of such approximations may be difficult to justify for some complicated fork-join networks. In this paper we propose a perfect simulation method based on coupling from the past to generate exact realisations from the equilibrium of fork-join networks. Using the simulated realisations, Monte Carlo estimates for the distributions of response times and queue lengths of fork-join networks are obtained. Comparisons of Monte Carlo estimates and theoretical approximations are also provided. The efficiency of the sampling algorithm is shown theoretically and via simulation.

The modeling of random bi-phasic, or porous, media has been, and still is, under active investigation by mathematicians, physicists, and physicians. In this paper we consider a thresholded random process X as a source of the two phases. The intervals when X is in a given phase, named chords, are the subject of interest. We focus on the study of the tails of the chord length distribution functions. In the literature concerned with real data, different types of tail behavior have been reported, among them exponential-like or power-like decay. We look for the link between the dependence structure of the underlying thresholded process X and the rate of decay of the chord length distribution. When the process X is a stationary Gaussian process, we relate the latter to the rate at which the covariance function of X decays at large lags. We show that exponential, or nearly exponential, decay of the tail of the distribution of the chord lengths is very common, perhaps surprisingly so.

Radix Sort is a sorting algorithm based on analyzing digital data. We study the number of swaps made by Radix Select (a one-sided version of Radix Sort) to find an element with a randomly selected rank. This kind of grand average provides a smoothing over all individual distributions for specific fixed-order statistics. We give an exact analysis for the grand mean and an asymptotic analysis for the grand variance, obtained by poissonization, the Mellin transform, and depoissonization. The digital data model considered is the Bernoulli(p). The distributions involved in the swaps experience a phase change between the biased cases (p ≠ ½) and the unbiased case (p = ½). In the biased cases, the grand distribution for the number of swaps (when suitably scaled) converges to that of a perpetuity built from a two-point distribution. The tool for this proof is contraction in the Wasserstein metric space, and identifying the limit as the fixed-point solution of a distributional equation. In the unbiased case the same scaling for the number of swaps gives a limiting constant in probability.

We consider a feed-forward network with a single-server station serving jobs with multiple levels of priority. The service discipline is preemptive in that the server always serves a job with the current highest level of priority. For this system with discontinuous dynamics, we establish the sample path large deviation principle using a weak convergence argument. In the special case where jobs have two different levels of priority, we also explicitly identify the exponential decay rate of the total population overflow probabilities by examining the geometry of the zero-level sets of the system Hamiltonians.

We study optimal control problems for (time-)delayed stochastic differential equations with jumps. We establish sufficient and necessary stochastic maximum principles for an optimal control of such systems. The associated adjoint processes are shown to satisfy a (time-)advanced backward stochastic differential equation (ABSDE). Several results on existence and uniqueness of such ABSDEs are shown. The results are illustrated by an application to optimal consumption from a cash flow with delay.