Consider a sample 𝒳n={X1,…,Xn} of independent and identically distributed variables drawn with a probability distribution ℙX supported on a compact set M⊂ℝd. In this paper we mainly deal with the study of a natural estimator for the geodesic distance on M. Under rather general geometric assumptions on M, we prove a general convergence result. Assuming M to be a compact manifold of known dimension d′≤d, and under regularity assumptions on ℙX, we give an explicit convergence rate. In the case when M has no boundary, knowledge of the dimension d′ is not needed to obtain this convergence rate. The second part of the work consists in building an estimator for the Fréchet expectations on M, and proving its convergence under regularity conditions, applying the previous results.

We consider an open problem of obtaining the optimal operational sequence for the 1-out-of-n system with warm standby. Using the virtual age concept and the cumulative exposure model, we show that the components should be activated in accordance with the increasing sequence of their lifetimes. Lifetimes of the components and the system are compared with respect to the stochastic precedence order and its generalization. Only specific cases of this optimal problem were considered in the literature previously.

In their 1960 book on finite Markov chains, Kemeny and Snell established that a certain sum is invariant. The value of this sum has become known as Kemeny’s constant. Various proofs have been given over time, some more technical than others. We give here a very simple physical justification, which extends without a hitch to continuous-time Markov chains on a finite state space. For Markov chains with denumerably infinite state space, the constant may be infinite and even if it is finite, there is no guarantee that the physical argument will hold. We show that the physical interpretation does go through for the special case of a birth-and-death process with a finite value of Kemeny’s constant.

The study of finite approximations of probability measures has a long history. In Xu and Berger (2017), the authors focused on constrained finite approximations and, in particular, uniform ones in dimension d=1. In the present paper we give an elementary construction of a uniform decomposition of probability measures in dimension d≥1. We then use this decomposition to obtain upper bounds on the rate of convergence of the optimal uniform approximation error. These bounds appear to be the generalization of the ones obtained by Xu and Berger (2017) and to be sharp for generic probability measures.

Bizarrely shaped voting districts are frequently lambasted as likely instances of gerrymandering. In order to systematically identify such instances, researchers have devised several tests for so-called geographic compactness (i.e. shape niceness). We demonstrate that under certain conditions, a party can gerrymander a competitive state into geographically compact districts to win an average of over 70% of the districts. Our results suggest that geometric features alone may fail to adequately combat partisan gerrymandering.

The study of high-dimensional distributions is of interest in probability theory, statistics, and asymptotic convex geometry, where the object of interest is the uniform distribution on a convex set in high dimensions. The ℓp-spaces and norms are of particular interest in this setting. In this paper we establish a limit theorem for distributions on ℓp-spheres, conditioned on a rare event, in a high-dimensional geometric setting. As part of our proof, we establish a certain large deviation principle that is also relevant to the study of the tail behavior of random projections of ℓp-balls in a high-dimensional Euclidean space.

Introduced recently in mathematical finance by Bayer et al. (2016), the rough Bergomi model has proved particularly efficient to calibrate option markets. We investigate some of its probabilistic properties, in particular proving a pathwise large deviations principle for a small-noise version of the model. The exponential function (continuous but superlinear) as well as the drift appearing in the volatility process fall beyond the scope of existing results, and a dedicated analysis is needed.

We examine a system of interacting random walks with leftward drift on ℤ, which begins with a single active particle at the origin and some distribution of inactive particles on the positive integers. Inactive particles become activated when landed on by other particles, and all particles beginning at the same point possess equal leftward drift. Once activated, the trajectories of distinct particles are independent. This system belongs to a broader class of problems involving interacting random walks on rooted graphs, referred to collectively as the frog model. Additional conditions that we impose on our model include that the number of frogs (i.e. particles) at positive integer points is a sequence of independent random variables which is increasing in terms of the standard stochastic order, and that the sequence of leftward drifts associated with frogs originating at these points is decreasing. Our results include sharp conditions with respect to the sequence of random variables and the sequence of drifts that determine whether the model is transient (meaning the probability infinitely many frogs return to the origin is 0) or nontransient. We consider several, more specific, versions of the model described, and a cleaner, more simplified set of sharp conditions will be established for each case.

The selection model in population genetics is a dynamic system on the set of the probability distributions 𝒑=(p1,…,pn) of the alleles A1…,An, with pi(t+1) proportional to pi(t) multiplied by ∑jfi,jpj(t), and fi,j=fj,i interpreted as a fitness of the gene pair (Ai,Aj). It is known that 𝒑̂ is a locally stable equilibrium if and only if 𝒑̂ is a strict local maximum of the quadratic form 𝒑T𝒇𝒑. Usually, there are multiple local maxima and lim𝒑(t) depends on 𝒑(0). To address the question of a typical behavior of {𝒑(t)}, John Kingman considered the case when the fi,j are independent and [0,1]-uniform. He proved that with high probability (w.h.p.) no local maximum may have more than 2.49n1∕2 positive components, and reduced 2.49 to 2.14 for a nonbiological case of exponentials on [0,∞). We show that the constant 2.14 serves a broad class of smooth densities on [0,1] with the increasing hazard rate. As for a lower bound, we prove that w.h.p. for all k≤2n1∕3, there are many k-element subsets of [n] that pass a partial test to be a support of a local maximum. Still, it may well be that w.h.p. the actual supports are much smaller. In that direction, we prove that w.h.p. a support of a local maximum, which does not contain a support of a local equilibrium, is very unlikely to have size exceeding ⅔log2n and, for the uniform fitnesses, there are super-polynomially many potential supports free of local equilibriums of size close to ½log2n.

Let Xn(k) be the number of vertices at level k in a random recursive tree with n+1 vertices. We are interested in the asymptotic behavior of Xn(k) for intermediate levels k=kn satisfying kn→∞ and kn=o(logn) as n→∞. In particular, we prove weak convergence of finite-dimensional distributions for the process (Xn ([knu]))u>0, properly normalized and centered, as n→∞. The limit is a centered Gaussian process with covariance (u,v)↦(u+v)−1. One-dimensional distributional convergence of Xn(kn), properly normalized and centered, was obtained with the help of analytic tools by Fuchs et al. (2006). In contrast, our proofs, which are probabilistic in nature, exploit a connection of our model with certain Crump–Mode–Jagers branching processes.

In this paper we prove asymptotic upper bounds on the variance of the number of vertices and the missed area of inscribed random disc-polygons in smooth convex discs whose boundary is C+2. We also consider a circumscribed variant of this probability model in which the convex disc is approximated by the intersection of random circles.

In this paper we consider an infinite system of instantaneously coalescing rate 1 simple symmetric random walks on ℤ2, started from the initial condition with all sites in ℤ2 occupied. Two-dimensional coalescing random walks are a `critical' model of interacting particle systems: unlike coalescence models in dimension three or higher, the fluctuation effects are important for the description of large-time statistics in two dimensions, manifesting themselves through the logarithmic corrections to the `mean field' answers. Yet the fluctuation effects are not as strong as for the one-dimensional coalescence, in which case the fluctuation effects modify the large time statistics at the leading order. Unfortunately, unlike its one-dimensional counterpart, the two-dimensional model is not exactly solvable, which explains a relative scarcity of rigorous analytic answers for the statistics of fluctuations at large times. Our contribution is to find, for any N≥2, the leading asymptotics for the correlation functions ρN(x1,…,xN) as t→∞. This generalises the results for N=1 due to Bramson and Griffeath (1980) and confirms a prediction in the physics literature for N>1. An analogous statement holds for instantaneously annihilating random walks. The key tools are the known asymptotic ρ1(t)∼logt∕πt due to Bramson and Griffeath (1980), and the noncollision probability 𝒑NC(t), that no pair of a finite collection of N two-dimensional simple random walks meets by time t, whose asymptotic 𝒑NC(t)∼c0(logt)-(N2) was found by Cox et al. (2010). We re-derive the asymptotics, and establish new error bounds, both for ρ1(t) and 𝒑NC(t) by proving that these quantities satisfy effective rate equations; that is, approximate differential equations at large times. This approach can be regarded as a generalisation of the Smoluchowski theory of renormalised rate equations to multi-point statistics.

Different Markov chains can be used for approximate sampling of a distribution given by an unnormalized density function with respect to the Lebesgue measure. The hit-and-run, (hybrid) slice sampler, and random walk Metropolis algorithm are popular tools to simulate such Markov chains. We develop a general approach to compare the efficiency of these sampling procedures by the use of a partial ordering of their Markov operators, the covariance ordering. In particular, we show that the hit-and-run and the simple slice sampler are more efficient than a hybrid slice sampler based on hit-and-run, which, itself, is more efficient than a (lazy) random walk Metropolis algorithm.

In this paper we consider the degree-wise effect of a second step for a random walk on a graph. We prove that under the configuration model, for any fixed degree sequence the probability of exceeding a given degree threshold is smaller after two steps than after one. This builds on recent work of Kramer et al. (2016) regarding the friendship paradox under random walks.

We extend the work of Antunović et al. (2016) on competing types in preferential attachment models to include cases where the types have different fitnesses, which may be either multiplicative or additive. We show that, depending on the values of the parameters of the models, there are different possible limiting behaviours depending on the zeros of a certain function. In particular, we show the existence of choices of the parameters where one type is favoured both by having higher fitness and by the type of attachment mechanism, but the other type has a positive probability of dominating the network in the limit.

Consider a bipartite random geometric graph on the union of two independent homogeneous Poisson point processes in d-space, with distance parameter r and intensities λ,μ. For any λ>0 we consider the percolation threshold μc(λ) associated to the parameter μ. Denoting by λc the percolation threshold for the standard Poisson Boolean model with radii r, we show the lower bound μc(λ)≥clog(c∕(λ−λc)) for any λ>λc with c>0 a fixed constant. In particular, there is no phase transition in μ at the critical value of λ, that is, μc(λc) =∞.

The mean time to failure (MTTF) function in age replacement is used to evaluate the performance and effectiveness of the age replacement policy. In this paper, based on the MTTF function, we introduce two new nonparametric classes of lifetime distributions with nonmonotonic mean time to failure in age replacement; increasing then decreasing MTTF (IDMTTF) and decreasing then increasing MTTF (DIMTTF). The implications between these classes of distributions and some existing classes of nonmonotonic ageing classes are studied. The characterizations of IDMTTF and DIMTTF in terms of the scaled total time on test transform are also obtained.

In this paper we are concerned with modelling the reliability of a system subject to external shocks. In a run shock model, the system fails when a sequence of shocks above a threshold arrive in succession. Nevertheless, using a single threshold to measure the severity of a shock is too critical in real practice. To this end, we develop a generalized run shock model with two thresholds. We employ a phase-type distribution to model the damage size and the inter-arrival time of shocks, which is highly versatile and may be used to model many quantitative features of random phenomenon. Furthermore, we use the Markovian property to construct a multi-state system which degrades with the arrival of shocks. We also provide a numerical example to illustrate our results.

The Samaniego signature is a relevant tool for studying the performance of a system whose component lifetimes are exchangeable. It is well known that the stochastic ordering of the signatures of two systems implies the same for the respective system lifetimes. We prove that the reverse claim is not true when the component lifetimes are independent and identically distributed. There exist small proportions of systems with stochastically ordered lifetimes whose signatures are not ordered. We present a simple procedure in order to check whether the system lifetimes are stochastically ordered even if their signatures are not comparable.

De Finetti’s optimal dividend problem has recently been extended to the case when dividend payments can be made only at Poisson arrival times. In this paper we consider the version with bail-outs where the surplus must be nonnegative uniformly in time. For a general spectrally negative Lévy model, we show the optimality of a Parisian-classical reflection strategy that pays the excess above a given barrier at each Poisson arrival time and also reflects from below at 0 in the classical sense.