Latouche and Nguyen (2015b) constructed a sequence of stochastic fluid processes and showed that it converges weakly to a Markov-modulated Brownian motion (MMBM). Here, we construct a different sequence of stochastic fluid processes and show that it converges strongly to an MMBM. To the best of our knowledge, this is the first result on strong convergence to a Markov-modulated Brownian motion. Besides implying weak convergence, such a strong approximation constitutes a powerful tool for developing deep results for sophisticated models. Additionally, we prove that the rate of this almost sure convergence is $o(n^{-1/2} \log n)$ . When reduced to the special case of standard Brownian motion, our convergence rate is an improvement over that obtained by a different approximation in Gorostiza and Griego (1980), which is $o(n^{-1/2}(\log n)^{5/2})$ .

We analyze the discounted probability of exponential Parisian ruin for the so-called scaled classical Cramér–Lundberg risk model. As in Cohen and Young (2020), we use the comparison method from differential equations to prove that the discounted probability of exponential Parisian ruin for the scaled classical risk model converges to the corresponding discounted probability for its diffusion approximation, and we derive the rate of convergence.

Network dynamics with point-process-based interactions are of paramount modeling interest. Unfortunately, most relevant dynamics involve complex graphs of interactions for which an exact computational treatment is impossible. To circumvent this difficulty, the replica-mean-field approach focuses on randomly interacting replicas of the networks of interest. In the limit of an infinite number of replicas, these networks become analytically tractable under the so-called ‘Poisson hypothesis’. However, in most applications this hypothesis is only conjectured. In this paper we establish the Poisson hypothesis for a general class of discrete-time, point-process-based dynamics that we propose to call fragmentation-interaction-aggregation processes, and which are introduced here. These processes feature a network of nodes, each endowed with a state governing their random activation. Each activation triggers the fragmentation of the activated node state and the transmission of interaction signals to downstream nodes. In turn, the signals received by nodes are aggregated to their state. Our main contribution is a proof of the Poisson hypothesis for the replica-mean-field version of any network in this class. The proof is obtained by establishing the propagation of asymptotic independence for state variables in the limit of an infinite number of replicas. Discrete-time Galves–Löcherbach neural networks are used as a basic instance and illustration of our analysis.

We introduce a notion of finite sampling consistency for phylogenetic trees and show that the set of finitely sampling-consistent and exchangeable distributions on n-leaf phylogenetic trees is a polytope. We use this polytope to show that the set of all exchangeable and sampling-consistent distributions on four-leaf phylogenetic trees is exactly Aldous’ beta-splitting model, and give a description of some of the vertices for the polytope of distributions on five leaves. We also introduce a new semialgebraic set of exchangeable and sampling consistent models we call the multinomial model and use it to characterize the set of exchangeable and sampling-consistent distributions. Using this new model, we obtain a finite de Finetti-type theorem for rooted binary trees in the style of Diaconis’ theorem on finite exchangeable sequences.

This paper studies the joint tail asymptotics of extrema of the multi-dimensional Gaussian process over random intervals defined as $P(u)\;:\!=\; \mathbb{P}\{\cap_{i=1}^n (\sup_{t\in[0,\mathcal{T}_i]} ( X_{i}(t) +c_i t )>a_i u )\}$ , $u\rightarrow\infty$ , where $X_i(t)$ , $t\ge0$ , $i=1,2,\ldots,n$ , are independent centered Gaussian processes with stationary increments, $\boldsymbol{\mathcal{T}}=(\mathcal{T}_1, \ldots, \mathcal{T}_n)$ is a regularly varying random vector with positive components, which is independent of the Gaussian processes, and $c_i\in \mathbb{R}$ , $a_i>0$ , $i=1,2,\ldots,n$ . Our result shows that the structure of the asymptotics of P(u) is determined by the signs of the drifts $c_i$ . We also discuss a relevant multi-dimensional regenerative model and derive the corresponding ruin probability.

We present a new and straightforward algorithm that simulates exact sample paths for a generalized stress-release process. The computation of the exact law of the joint inter-arrival times is detailed and used to derive this algorithm. Furthermore, the martingale generator of the process is derived, and induces theoretical moments which generalize some results of [3] and are used to demonstrate the validity of our simulation algorithm.

We investigate properties of random mappings whose core is composed of derangements as opposed to permutations. Such mappings arise as the natural framework for studying the Screaming Toes game described, for example, by Peter Cameron. This mapping differs from the classical case primarily in the behaviour of the small components, and a number of explicit results are provided to illustrate these differences.

We consider a variant of a classical coverage process, the Boolean model in $\mathbb{R}^d$ . Previous efforts have focused on convergence of the unoccupied region containing the origin to a well-studied limit C. We study the intersection of sets centered at points of a Poisson point process confined to the unit ball. Using a coupling between the intersection model and the original Boolean model, we show that the scaled intersection converges weakly to the same limit C. Along the way, we present some tools for studying statistics of a class of intersection models.

We study weighted ensemble, an interacting particle method for sampling distributions of Markov chains that has been used in computational chemistry since the 1990s. Many important applications of weighted ensemble require the computation of long time averages. We establish the consistency of weighted ensemble in this setting by proving an ergodic theorem for time averages. As part of the proof, we derive explicit variance formulas that could be useful for optimizing the method.

This paper is devoted to the study of regime-switching jump diffusion processes with countable regimes. It aims to establish Foster–Lyapunov-type criteria for exponential ergodicity of such processes. After recalling results concerning the petiteness of compact sets, this paper presents sufficient conditions for the existence of a Foster–Lyapunov function; this, in turn, helps to establish sufficient conditions for the desired exponential ergodicity for regime-switching jump diffusion processes. Finally, an application to feedback control problems is presented.

For a non-negative separable random field Z(t), $t\in \mathbb{R}^d$ , satisfying some mild assumptions, we show that $ H_Z^\delta =\lim_{{T} \to \infty} ({1}/{T^d}) \mathbb{E}\{{\sup_{ t\in [0,T]^d \cap \delta \mathbb{Z}^d } Z(t) }\} <\infty$ for $\delta \ge 0$ , where $0 \mathbb{Z}^d\,:\!=\,\mathbb{R}^d$ , and prove that $H_Z^0$ can be approximated by $H_Z^\delta$ if $\delta$ tends to 0. These results extend the classical findings for Pickands constants $H_{Z}^\delta$ , defined for $Z(t)= \exp( \sqrt{ 2} B_\alpha (t)- \lvert {t} \rvert^{2\alpha })$ , $t\in \mathbb{R}$ , with $B_\alpha$ a standard fractional Brownian motion with Hurst parameter $\alpha \in (0,1]$ . The continuity of $H_{Z}^\delta$ at $\delta=0$ is additionally shown for two particular extensions of Pickands constants.

A continuum of stochastic dominance rules, also referred to as fractional stochastic dominance (SD), was introduced by Müller, Scarsini, Tsetlin, and Winkler (2017) to cover preferences from first- to second-order SD. Fractional SD can be used to explain many individual behaviors in economics. In this paper we introduce the concept of fractional pure SD, a special case of fractional SD. We investigate further properties of fractional SD, for example the generating processes of fractional pure SD via $\gamma$ -transfers of probability, Yaari’s dual characterization by utilizing the special class of distortion functions, the separation theorem in terms of first-order SD and fractional pure SD, Strassen’s representation, and bivariate characterization. We also establish several closure properties of fractional SD under quantile truncation, under comonotonic sums, and under distortion, as well as its equivalence characterization. Examples of distributions ordered in the sense of fractional SD are provided.

Consider two-type linear-fractional branching processes in varying environments with asymptotically constant mean matrices. Let $\nu$ be the extinction time. Under certain conditions, we show that both $\mathbb{P}(\nu=n)$ and $\mathbb{P}(\nu>n)$ are asymptotically the same as some functions of the products of spectral radii of the mean matrices. We also give an example for which $\mathbb{P}(\nu=n)$ decays with various speeds such as ${c}/({n^{1/2}\log n)^2}$ , ${c}/{n^\beta}$ , $\beta >1$ , which are very different from those of homogeneous multitype Galton–Watson processes.

Let T be the regular tree in which every vertex has exactly $d\ge 3$ neighbours. Run a branching random walk on T, in which at each time step every particle gives birth to a random number of children with mean d and finite variance, and each of these children moves independently to a uniformly chosen neighbour of its parent. We show that, starting with one particle at some vertex 0 and conditionally on survival of the process, the time it takes for every vertex within distance r of 0 to be hit by a particle of the branching random walk is $r + ({2}/{\log(3/2)})\log\log r + {\mathrm{o}}(\log\log r)$ .

We investigate vertex levels of containment in a random hypergraph grown in the spirit of a recursive tree. We consider a local profile tracking the evolution of the containment of a particular vertex over time, and a global profile concerned with counts of the number of vertices of a particular containment level.

For the local containment profile, we obtain the exact mean, variance, and probability distribution in terms of standard combinatorial quantities such as generalized harmonic numbers and Stirling numbers of the first kind. Asymptotically, we observe phases: the early vertices have an asymptotically normal distribution, intermediate vertices have a Poisson distribution, and late vertices have a degenerate distribution.

As for the global containment profile, we establish an asymptotically normal distribution for the number of vertices at the smallest containment level as well as their covariances with the number of vertices at the second smallest containment level and the variances of these numbers. The orders in the variance–covariance matrix establish concentration laws.