Both sequential Monte Carlo (SMC) methods (a.k.a. ‘particle filters’) and sequential Markov chain Monte Carlo (sequential MCMC) methods constitute classes of algorithms which can be used to approximate expectations with respect to (a sequence of) probability distributions and their normalising constants. While SMC methods sample particles conditionally independently at each time step, sequential MCMC methods sample particles according to a Markov chain Monte Carlo (MCMC) kernel. Introduced over twenty years ago in [6], sequential MCMC methods have attracted renewed interest recently as they empirically outperform SMC methods in some applications. We establish an $\mathbb{L}_r$ -inequality (which implies a strong law of large numbers) and a central limit theorem for sequential MCMC methods and provide conditions under which errors can be controlled uniformly in time. In the context of state-space models, we also provide conditions under which sequential MCMC methods can indeed outperform standard SMC methods in terms of asymptotic variance of the corresponding Monte Carlo estimators.

In this paper, we solve exit problems for a one-sided Markov additive process (MAP) which is exponentially killed with a bivariate killing intensity $\omega(\cdot,\cdot)$ dependent on the present level of the process and the current state of the environment. Moreover, we analyze the respective resolvents. All identities are expressed in terms of new generalizations of classical scale matrices for MAPs. We also remark on a number of applications of the obtained identities to (controlled) insurance risk processes. In particular, we show that our results can be applied to the Omega model, where bankruptcy takes place at rate $\omega(\cdot,\cdot)$ when the surplus process becomes negative. Finally, we consider Markov-modulated Brownian motion (MMBM) as a special case and present analytical and numerical results for a particular choice of piecewise intensity function $\omega(\cdot,\cdot)$ .

Let V be an n-set, and let X be a random variable taking values in the power-set of V. Suppose we are given a sequence of random coupons $X_1, X_2, \ldots $ , where the $X_i$ are independent random variables with distribution given by X. The covering time T is the smallest integer $t\geq 0$ such that $\bigcup_{i=1}^t X_i=V$ . The distribution of T is important in many applications in combinatorial probability, and has been extensively studied. However the literature has focused almost exclusively on the case where X is assumed to be symmetric and/or uniform in some way.

In this paper we study the covering time for much more general random variables X; we give general criteria for T being sharply concentrated around its mean, precise tools to estimate that mean, as well as examples where T fails to be concentrated and when structural properties in the distribution of X allow for a very different behaviour of T relative to the symmetric/uniform case.

We study networks of interacting queues governed by utility-maximising service-rate allocations in both discrete and continuous time. For finite networks we establish stability and some steady-state moment bounds under natural conditions and rather weak assumptions on utility functions. These results are obtained using direct applications of Lyapunov–Foster-type criteria, and apply to a wide class of systems, including those for which fluid-limit-based approaches are not applicable. We then establish stability and some steady-state moment bounds for two classes of infinite networks, with single-hop and multi-hop message routes. These results are proved by considering the infinite systems as limits of their truncated finite versions. The uniform moment bounds for the finite networks play a key role in these limit transitions.

Large samples from a light-tailed distribution often have a well-defined shape. This paper examines the implications of the assumption that there is a limit shape. We show that the limit shape determines the upper quantiles for a large class of random variables. These variables may be described loosely as continuous homogeneous functionals of the underlying random vector. They play an important role in evaluating risk in a multivariate setting. The paper also looks at various coefficients of tail dependence and at the distribution of the scaled sample points for large samples. The paper assumes convergence in probability rather than almost sure convergence. This results in an elegant theory. In particular, there is a simple characterization of domains of attraction.

In this paper, a reflected stochastic differential equation (SDE) with jumps is studied for the case where the constraint acts on the law of the solution rather than on its paths. These reflected SDEs have been approximated by Briand et al. (2016) using a numerical scheme based on particles systems, when no jumps occur. The main contribution of this paper is to prove the existence and the uniqueness of the solutions to this kind of reflected SDE with jumps and to generalize the results obtained by Briand et al. (2016) to this context.

We study failure rate monotonicity and generalised convex transform stochastic ordering properties of random variables, with an emphasis on applications. We are especially interested in the effect of a tail-weight iteration procedure to define distributions, which is equivalent to the characterisation of moments of the residual lifetime at a given instant. For the monotonicity properties, we are mainly concerned with hereditary properties with respect to the iteration procedure providing counterexamples showing either that the hereditary property does not hold or that inverse implications are not true. For the stochastic ordering, we introduce a new criterion, based on the analysis of the sign variation of a suitable function. This criterion is then applied to prove ageing properties of parallel systems formed with components that have exponentially distributed lifetimes.

The Gaussian polytope $\mathcal P_{n,d}$ is the convex hull of n independent standard normally distributed points in $\mathbb{R}^d$ . We derive explicit expressions for the probability that $\mathcal P_{n,d}$ contains a fixed point $x\in\mathbb{R}^d$ as a function of the Euclidean norm of x, and the probability that $\mathcal P_{n,d}$ contains the point $\sigma X$ , where $\sigma\geq 0$ is constant and X is a standard normal vector independent of $\mathcal P_{n,d}$ . As a by-product, we also compute the expected number of k-faces and the expected volume of $\mathcal P_{n,d}$ , thus recovering the results of Affentranger and Schneider (Discr. and Comput. Geometry, 1992) and Efron (Biometrika, 1965), respectively. All formulas are in terms of the volumes of regular spherical simplices, which, in turn, can be expressed through the standard normal distribution function $\Phi(z)$ and its complex version $\Phi(iz)$ . The main tool used in the proofs is the conic version of the Crofton formula.

The equilibrium properties of allocation algorithms for networks with a large number of nodes with finite capacity are investigated. Every node receives a flow of requests. When a request arrives at a saturated node, i.e. a node whose capacity is fully utilized, an allocation algorithm may attempt to reallocate the request to a non-saturated node. For the algorithms considered, the reallocation comes at a price: either extra capacity is required in the system, or the processing time of a reallocated request is increased. The paper analyzes the properties of the equilibrium points of the associated asymptotic dynamical system when the number of nodes gets large. At this occasion the classical model of Gibbens, Hunt, and Kelly (1990) in this domain is revisited. The absence of known Lyapunov functions for the corresponding dynamical system significantly complicates the analysis. Several techniques are used. Analytic and scaling methods are used to identify the equilibrium points. We identify the subset of parameters for which the limiting stochastic model of these networks has multiple equilibrium points. Probabilistic approaches are used to prove the stability of some of them. A criterion of exponential stability with the spectral gap of the associated linear operator of equilibrium points is also obtained.

We consider the Voronoi diagram generated by n independent and identically distributed $\mathbb{R}^{d}$ -valued random variables with an arbitrary underlying probability density function f on $\mathbb{R}^{d}$ , and analyze the asymptotic behaviors of certain geometric properties, such as the measure, of the Voronoi cells as n tends to infinity. We adapt the methods used by Devroye et al. (2017) to conduct a study of the asymptotic properties of two types of Voronoi cells: (1) Voronoi cells that have a fixed nucleus; (2) Voronoi cells that contain a fixed point. We show that the geometric properties of both types of cells resemble those in the case when the Voronoi diagram is generated by a homogeneous Poisson point process. Additionally, for the second type of Voronoi cells, we determine the limiting distribution, which is universal in all choices of f, of the re-scaled measure of the cells.

We introduce a unified framework for solving first passage times of time-homogeneous diffusion processes. Using potential theory and perturbation theory, we are able to deduce closed-form truncated probability densities, as asymptotics or approximations to the original first passage time densities, for single-side level crossing problems. The framework is applicable to diffusion processes with continuous drift functions; in particular, for bounded drift functions, we show that the perturbation series converges. In the present paper, we demonstrate examples of applying our framework to the Ornstein–Uhlenbeck, Bessel, exponential-Shiryaev, and hypergeometric diffusion processes (the latter two being previously studied by Dassios and Li (2018) and Borodin (2009), respectively). The purpose of this paper is to provide a fast and accurate approach to estimating first passage time densities of various diffusion processes.