Queueing networks are stochastic systems formed by interconnected resources routing and serving jobs. They induce jump processes with distinctive properties, and find widespread use in inferential tasks. Here, service rates for jobs and potential bottlenecks in the routing mechanism must be estimated from a reduced set of observations. However, this calls for the derivation of complex conditional density representations, over both the stochastic network trajectories and the rates, which is considered an intractable problem. Numerical simulation procedures designed for this purpose do not scale, because of high computational costs; furthermore, variational approaches relying on approximating measures and full independence assumptions are unsuitable. In this paper, we offer a probabilistic interpretation of variational methods applied to inference tasks with queueing networks, and show that approximating measure choices routinely used with jump processes yield ill-defined optimization problems. Yet we demonstrate that it is still possible to enable a variational inferential task, by considering a novel space expansion treatment over an analogous counting process for job transitions. We present and compare exemplary use cases with practical queueing networks, showing that our framework offers an efficient and improved alternative where existing variational or numerically intensive solutions fail.

The logistic birth and death process is perhaps the simplest stochastic population model that has both density-dependent reproduction and a phase transition, and a lot can be learned about the process by studying its extinction time, $\tau_n$, as a function of system size n. A number of existing results describe the scaling of $\tau_n$ as $n\to\infty$ for various choices of reproductive rate $r_n$ and initial population $X_n(0)$ as a function of n. We collect and complete this picture, obtaining a complete classification of all sequences $(r_n)$ and $(X_n(0))$ for which there exist rescaling parameters $(s_n)$ and $(t_n)$ such that $(\tau_n-t_n)/s_n$ converges in distribution as $n\to\infty$, and identifying the limits in each case.

We reduce the upper bound for the bond percolation threshold of the cubic lattice from 0.447 792 to 0.347 297. The bound is obtained by a growth process approach which views the open cluster of a bond percolation model as a dynamic process. A three-dimensional dynamic process on the cubic lattice is constructed and then projected onto a carefully chosen plane to obtain a two-dimensional dynamic process on a triangular lattice. We compare the bond percolation models on the cubic lattice and their projections, and demonstrate that the bond percolation threshold of the cubic lattice is no greater than that of the triangular lattice. Applying the approach to the body-centered cubic lattice yields an upper bound of 0.292 893 for its bond percolation threshold.

In this paper we analyse the limiting conditional distribution (Yaglom limit) for stochastic fluid models (SFMs), a key class of models in the theory of matrix-analytic methods. So far, only transient and stationary analyses of SFMs have been considered in the literature. The limiting conditional distribution gives useful insights into what happens when the process has been evolving for a long time, given that its busy period has not ended yet. We derive expressions for the Yaglom limit in terms of the singularity˜$s^*$ such that the key matrix of the SFM, ${\boldsymbol{\Psi}}(s)$, is finite (exists) for all $s\geq s^*$ and infinite for $s<s^*$. We show the uniqueness of the Yaglom limit and illustrate the application of the theory with simple examples.

In this note, we consider dynamic assortment optimization with incomplete information under the capacitated multinomial logit choice model. Recently, it has been shown that the regret (the cumulative expected revenue loss caused by offering suboptimal assortments) that any decision policy endures is bounded from below by a constant times $\sqrt {NT}$, where $N$ denotes the number of products and $T$ denotes the time horizon. This result is shown under the assumption that the product revenues are constant, and thus leaves the question open whether a lower regret rate can be achieved for nonconstant revenue parameters. In this note, we show that this is not the case: we show that, for any vector of product revenues there is a positive constant such that the regret of any policy is bounded from below by this constant times $\sqrt {N T}$. Our result implies that policies that achieve ${{\mathcal {O}}}(\sqrt {NT})$ regret are asymptotically optimal for all product revenue parameters.

We study propagation of avalanches in a certain excitable network. The model is a particular case of the one introduced by Larremore et al. (Phys. Rev. E, 2012) and is mathematically equivalent to an endemic variation of the Reed–Frost epidemic model introduced by Longini (Math. Biosci., 1980). Two types of heuristic approximation are frequently used for models of this type in applications: a branching process for avalanches of a small size at the beginning of the process and a deterministic dynamical system once the avalanche spreads to a significant fraction of a large network. In this paper we prove several results concerning the exact relation between the avalanche model and these limits, including rates of convergence and rigorous bounds for common characteristics of the model.

Oscillatory systems of interacting Hawkes processes with Erlang memory kernels were introduced by Ditlevsen and Löcherbach (Stoch. Process. Appl., 2017). They are piecewise deterministic Markov processes (PDMP) and can be approximated by a stochastic diffusion. In this paper, first, a strong error bound between the PDMP and the diffusion is proved. Second, moment bounds for the resulting diffusion are derived. Third, approximation schemes for the diffusion, based on the numerical splitting approach, are proposed. These schemes are proved to converge with mean-square order 1 and to preserve the properties of the diffusion, in particular the hypoellipticity, the ergodicity, and the moment bounds. Finally, the PDMP and the diffusion are compared through numerical experiments, where the PDMP is simulated with an adapted thinning procedure.

Signatures are useful in analyzing and evaluating coherent systems. However, their computation is a challenging problem, especially for complex coherent structures. In most cases the reliability of a binary coherent system can be linked to a tail probability associated with a properly defined waiting time random variable in a sequence of binary trials. In this paper we present a method for computing the minimal signature of a binary coherent system. Our method is based on matrix-geometric distributions. First, a proper matrix-geometric random variable corresponding to the system structure is found. Second, its probability generating function is obtained. Finally, the companion representation for the distribution of matrix-geometric distribution is used to obtain a matrix-based expression for the minimal signature of the coherent system. The results are also extended to a system with two types of components.

We consider a birth–death process with killing where transitions from state i may go to either state $i-1$ or state $i+1$ or an absorbing state (killing). Stochastic ordering results on the killing time are derived. In particular, if the killing rate in state i is monotone in i, then the distribution of the killing time with initial state i is stochastically monotone in i. This result is a consequence of the following one for a non-negative tri-diagonal matrix M: if the row sums of M are monotone, so are the row sums of $M^n$ for all $n\ge 2$.

We describe a process where two types of particles, marked red and blue, arrive in a domain at a constant rate. When a new particle arrives into the domain, if there are particles of the opposite color present within a distance of 1 from the new particle, then, among these particles, it matches to the one with the earliest arrival time, and both particles are removed. Otherwise, the particle is simply added to the system. Additionally, particles may lose patience and depart at a constant rate. We study the existence of a stationary regime for this process, when the domain is either a compact space or a Euclidean space. In the compact setting, we give a product-form characterization of the stationary distribution, and then prove an FKG-type inequality that establishes certain clustering properties of the particles in the steady state.

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.

In this paper we consider a new generalized finite mixture model formed by dependent and identically distributed (d.i.d.) components. We then establish results for the comparisons of lifetimes of two such generalized finite mixture models in two different cases: (i) when the two mixture models are formed from two random vectors $\textbf{X}$ and $\textbf{Y}$ but with the same weights, and (ii) when the two mixture models are formed with the same random vectors but with different weights. Because the lifetimes of k-out-of-n systems and coherent systems are special cases of the mixture model considered, we used the established results to compare the lifetimes of k-out-of-n systems and coherent systems with respect to the reversed hazard rate and hazard rate orderings.

We introduce a definition of long range dependence of random processes and fields on an (unbounded) index space $T\subseteq \mathbb{R}^d$ in terms of integrability of the covariance of indicators that a random function exceeds any given level. This definition is specifically designed to cover the case of random functions with infinite variance. We show the value of this new definition and its connection to limit theorems via some examples including subordinated Gaussian as well as random volatility fields and time series.

This paper introduces a non-linear and continuous-time opinion dynamics model with additive noise and state-dependent interaction rates between agents. The model features interaction rates which are proportional to a negative power of the opinion distances. We establish a non-local partial differential equation for the distribution of opinion distances and use Mellin transforms to provide an explicit formula for the stationary solution of the latter, when it exists. Our approach leads to new qualitative and quantitative results on this type of dynamics. To the best of our knowledge these Mellin transform results are the first quantitative results on the equilibria of opinion dynamics with distance-dependent interaction rates. The closed-form expressions for this class of dynamics are obtained for the two-agent case. However, the results can be used in mean-field models featuring several agents whose interaction rates depend on the empirical average of their opinions. The technique also applies to linear dynamics, namely with a constant interaction rate, on an interaction graph.

This paper considers a variant of the classical Cramér–Lundberg model that is particularly appropriate in the credit context, with the distinguishing feature that it corresponds to a finite number of obligors. The focus is on computing the ruin probability, i.e. the probability that the initial reserve, increased by the interest received from the obligors and decreased by the losses due to defaults, drops below zero. As well as an exact analysis (in terms of transforms) of this ruin probability, an asymptotic analysis is performed, including an efficient importance-sampling-based simulation approach.

The base model is extended in multiple dimensions: (i) we consider a model in which there may, in addition, be losses that do not correspond to defaults, (ii) then we analyze a model in which the individual obligors are coupled via a regime switching mechanism, (iii) then we extend the model so that between the losses the reserve process behaves as a Brownian motion rather than a deterministic drift, and (iv) we finally consider a set-up with multiple groups of statistically identical obligors.

We study, under mild conditions, the weak approximation constructed from a standard Poisson process for a class of Gaussian processes, and establish its sample path moderate deviations. The techniques consist of a good asymptotic exponential approximation in moderate deviations, the Besov–Lèvy modulus embedding, and an exponential martingale technique. Moreover, our results are applied to the weak approximations associated with the moving average of Brownian motion, fractional Brownian motion, and an Ornstein–Uhlenbeck process.