Gaussian random polytopes have received a lot of attention, especially in the case where the dimension is fixed and the number of points goes to infinity. Our focus is on the less-studied case where the dimension goes to infinity and the number of points is proportional to the dimension d. We study several natural quantities associated with Gaussian random polytopes in this setting. First, we show that the expected number of facets is equal to $C(\alpha)^{d+o(d)}$, where $C(\alpha)$ is some constant which depends on the constant of proportionality $\alpha$. We also extend this result to the expected number of k-facets. We then consider the more difficult problem of the asymptotics of the expected number of pairs of estranged facets of a Gaussian random polytope. When the number of points is 2d, we determine the constant C such that the expected number of pairs of estranged facets is equal to $C^{d+o(d)}$.

We study the last exit time that a spectrally negative Lévy process is below zero until it reaches a positive level b, denoted by $g_{\tau_b^+}$. We generalize the results of the infinite-horizon last exit time explored by Chiu and Yin (2005) by incorporating a random horizon $\tau_b^+$, which represents the first passage time above b. We derive an explicit expression for the joint Laplace transform of $g_{\tau_b^+}$ and $\tau_b^+$ by utilizing a hybrid observation scheme approach proposed by Li, Willmot, and Wong (2018). We further study the optimal prediction of $g_{\tau_b^+}$ in the $L_1$ sense, and find that the optimal stopping time is the first passage time above a level $y_b^{\ast}$, with an explicit characterization of the stopping boundary $y_b^{\ast}$. As examples, Brownian motion with drift and the Cramér–Lundberg model with exponential jumps are considered.

We consider time-inhomogeneous ordinary differential equations (ODEs) whose parameters are governed by an underlying ergodic Markov process. When this underlying process is accelerated by a factor $\varepsilon^{-1}$, an averaging phenomenon occurs and the solution of the ODE converges to a deterministic ODE as $\varepsilon$ vanishes. We are interested in cases where this averaged flow is globally attracted to a point. In that case, the equilibrium distribution of the solution of the ODE converges to a Dirac mass at this point. We prove an asymptotic expansion in terms of $\varepsilon$ for this convergence, with a somewhat explicit formula for the first-order term. The results are applied in three contexts: linear Markov-modulated ODEs, randomized splitting schemes, and Lotka–Volterra models in a random environment. In particular, as a corollary, we prove the existence of two matrices whose convex combinations are all stable but are such that, for a suitable jump rate, the top Lyapunov exponent of a Markov-modulated linear ODE switching between these two matrices is positive.

In a recent paper, the authors studied the distribution properties of a class of exchangeable processes, called measure-valued Pólya sequences (MVPSs), which arise as the observation process in a generalized urn sampling scheme. Here we present several results in the form of ‘sufficientness’ postulates that characterize their predictive distributions. In particular, we show that exchangeable MVPSs are the unique exchangeable models whose predictive distributions are a mixture of the marginal distribution and the average of a probability kernel evaluated at past observations. When the latter coincides with the empirical measure, we recover a well-known result for the exchangeable model with a Dirichlet process prior. In addition, we provide a ‘pure’ sufficientness postulate for exchangeable MVPSs that does not assume a particular analytic form for the predictive distributions. Two other sufficientness postulates consider the case when the state space is finite.

The basic question in perturbation analysis of Markov chains is: how do small changes in the transition kernels of Markov chains translate to chains in their stationary distributions? Many papers on the subject have shown, roughly, that the change in stationary distribution is small as long as the change in the kernel is much less than some measure of the convergence rate. This result is essentially sharp for generic Markov chains. We show that much larger errors, up to size roughly the square root of the convergence rate, are permissible for many target distributions associated with graphical models. The main motivation for this work comes from computational statistics, where there is often a tradeoff between the per-step error and per-step cost of approximate MCMC algorithms. Our results show that larger perturbations (and thus less-expensive chains) still give results with small error.

We consider interacting urns on a finite directed network, where both sampling and reinforcement processes depend on the nodes of the network. This extends previous research by incorporating node-dependent sampling and reinforcement. We classify the sampling and reinforcement schemes, as well as the networks on which the proportion of balls of either colour in each urn converges almost surely to a deterministic limit. We also investigate conditions for achieving synchronisation of the colour proportions across the urns and analyse fluctuations under specific conditions on the reinforcement scheme and network structure.

Gut and Stadmüller (2021, 2022) initiated the study of the elephant random walk with limited memory. Aguech and El Machkouri (2024) published a paper in which they discuss an extension of the results by Gut and Stadtmüller (2022) for an ‘increasing memory’ version of the elephant random walk without stops. Here we present a formal definition of the process that was hinted at by Gut and Stadtmüller. This definition is based on the triangular array setting. We give a positive answer to the open problem in Gut and Stadtmüller (2022) for the elephant random walk, possibly with stops. We also obtain the central limit theorem for the supercritical case of this model.

This paper defines and studies a broad class of shock models by assuming that a Markovian arrival process models the arrival pattern of shocks. Under the defined class, we show that the system’s lifetime follows the well-known phase-type distribution. Further, we examine the age replacement policy for systems with a continuous phase-type distribution, identifying sufficient conditions for determining the optimal replacement time. Since phase-type distributions are dense in the class of lifetime distributions, our findings for the age replacement policy are widely applicable. We include numerical examples and graphical illustrations to support our results.

The Wiener–Hopf factors of a Lévy process are the maximum and the displacement from the maximum at an independent exponential time. The majority of explicit solutions assume the upward jumps to be either phase-type or to have a rational Laplace transform, in which case the traditional expressions are lengthy expansions in terms of roots located by means of Rouché’s theorem. As an alternative, compact matrix formulas are derived, with parameters computable by iteration schemes.

We show that for $\lambda\in[0,{m_1}/({1+\sqrt{1-{1}/{m_1}}})]$, the biased random walk’s speed on a Galton–Watson tree without leaves is strictly decreasing, where $m_1\geq 2$. Our result extends the monotonic interval of the speed on a Galton–Watson tree.

We establish a number of results concerning the limiting behaviour of the longest edges in the genealogical tree generated by a continuous-time Galton–Watson process. Separately, we consider the large-time behaviour of the longest pendant edges, the longest (strictly) interior edges, and the longest of all the edges. These results extend the special case of long pendant edges of birth–death processes established in Bocharov et al. (2023).

We study the Markov chain Monte Carlo estimator for numerical integration for functions that do not need to be square integrable with respect to the invariant distribution. For chains with a spectral gap we show that the absolute mean error for $L^p$ functions, with $p \in (1,2)$, decreases like $n^{({1}/{p}) -1}$, which is known to be the optimal rate. This improves currently known results where an additional parameter $\delta \gt 0$ appears and the convergence is of order $n^{(({1+\delta})/{p})-1}$.

We develop explicit bounds for the tail of the distribution of the all-time supremum of a random walk with negative drift, where the increments have a truncated heavy-tailed distribution. As an application, we consider a ruin problem in the presence of reinsurance.

We prove an ergodic theorem for Markov chains indexed by the Ulam–Harris–Neveu tree over large subsets with arbitrary shape under two assumptions: (i) with high probability, two vertices in the large subset are far from each other, and (ii) with high probability, those two vertices have their common ancestor close to the root. The assumption on the common ancestor can be replaced by some regularity assumption on the Markov transition kernel. We verify that these assumptions are satisfied for some usual trees. Finally, with Markov chain Monte Carlo considerations in mind, we prove that when the underlying Markov chain is stationary and reversible, the Markov chain, that is the line graph, yields minimal variance for the empirical average estimator among trees with a given number of nodes. In doing so, we prove that the Hosoya–Wiener polynomial is minimized over $[{-}1,1]$ by the line graph among trees of a given size.

We establish a sample path moderate deviation principle for the integrated shot noise process with Poisson arrivals and non-stationary noises. As in Pang and Taqqu (2019), we assume that the noise is conditionally independent given the arrival times, and the distribution of each noise depends on its arrival time. As applications, we derive moderate deviation principles for the workload process and the running maximum process for a stochastic fluid queue with the integrated shot noise process as the input; we also show that a steady-state distribution exists and derive the exact tail asymptotics.

This paper is concerned with the growth rate of susceptible–infectious–recovered epidemics with general infectious period distribution on random intersection graphs. This type of graph is characterised by the presence of cliques (fully connected subgraphs). We study epidemics on random intersection graphs with a mixed Poisson degree distribution and show that in the limit of large population sizes the number of infected individuals grows exponentially during the early phase of the epidemic, as is generally the case for epidemics on asymptotically unclustered networks. The Malthusian parameter is shown to satisfy a variant of the classical Euler–Lotka equation. To obtain these results we construct a coupling of the epidemic process and a continuous-time multitype branching process, where the type of an individual is (essentially) given by the length of its infectious period. Asymptotic results are then obtained via an embedded single-type Crump–Mode–Jagers branching process.

We consider a superprocess $\{X_t\colon t\geq 0\}$ in a random environment described by a Gaussian field $\{W(t,x)\colon t\geq 0,x\in \mathbb{R}^d\}$. First, we set up a representation of $\mathbb{E}[\langle g, X_t\rangle\mathrm{e}^{-\langle \,f,X_t\rangle }\mid\sigma(W)\vee\sigma(X_r,0\leq r\leq s)]$ for $0\leq s < t$ and some functions f,g, which generalizes the result in Mytnik and Xiong (2007, Theorem 2.15). Next, we give a uniform upper bound for the conditional log-Laplace equation with unbounded initial values. We then use this to establish the corresponding conditional entrance law. Finally, the excursion representation of $\{X_t\colon t\geq 0\}$ is given.

Considering a double-indexed array $(Y_{n,i:\,n\ge 1,i\ge 1})$ of non-negative regularly varying random variables, we study the random-length weighted sums and maxima from its ‘row’ sequences. These sums and maxima may have the same tail and extremal indices (Markovich and Rodionov 2020). The main constraints of the latter results are that there exists a unique series in a scheme of series with the minimum tail index and the tail of the term number is lighter than the tail of the terms. Here, a bounded random number of series are allowed to have the minimum tail index and the tail of the term number may be heavier than the tail of the terms. We derive the tail and extremal indices of the weighted non-stationary random-length sequences under a broader set of conditions than in Markovich and Rodionov (2020). We provide examples of random sequences for which the assumptions are valid. Perspectives in adopting the results in different application areas are formulated.

We focus on obtaining Block–Savits type characterizations for different ageing classes as well as some important renewal classes by using the Laplace transform. We also introduce a novel approach, based on the equilibrium distribution, to handle situations where the techniques of Block and Savits (1980) either fail or involve tedious calculations. Our approach in conjunction with the theory of total positivity yields Vinogradov’s (1973) result for the increasing failure rate class when the distribution function is continuous. We also present simple but elegant proofs for Block and Savits’ results for the decreasing mean residual life, new better than used in expectation, and harmonic new better than used in expectation classes as applications of our approach. We address several other related issues that are germane to our problem. Finally, we conclude with a short discussion on the issue of convolutions.

We introduce a comprehensive method for establishing stochastic orders among order statistics in the independent and identically distributed case. This approach relies on the assumption that the underlying distribution is linked to a reference distribution through a transform order. Notably, this method exhibits broad applicability, particularly since several well-known nonparametric distribution families can be defined using relevant transform orders, including the convex and the star transform orders. Moreover, for convex-ordered families, we show that an application of Jensen’s inequality gives bounds for the probability that a random variable exceeds the expected value of its corresponding order statistic.