We prove that, for Poisson transmission and recovery processes, the classic susceptible→infected→recovered (SIR) epidemic model of Kermack and McKendrick provides, for any given time t>0, a strict lower bound on the expected number of susceptibles and a strict upper bound on the expected number of recoveries in the general stochastic SIR epidemic. The proof is based on the recent message passing representation of SIR epidemics applied to a complete graph.

In this paper we study the asymptotic behaviour near extinction of (sub-)critical continuous-state branching processes. In particular, we establish an analogue of Khintchine's law of the iterated logarithm near extinction time for a continuous-state branching process whose branching mechanism satisfies a given condition.

In this paper we estimate quantile sensitivities for dependent sequences via infinitesimal perturbation analysis, and prove asymptotic unbiasedness, weak consistency, and a central limit theorem for the estimators under some mild conditions. Two common cases, the regenerative setting and ϕ-mixing, are analyzed further, and a new batched estimator is constructed based on regenerative cycles for regenerative processes. Two numerical examples, the G/G/1 queue and the Ornstein–Uhlenbeck process, are given to show the effectiveness of the estimator.

We analyze the dynamics of a stochastic process driven by binomial random variables, where the probability of success depends on the past realization. We study the limit behavior when the group size is fixed but the number of iterations increases. It will become apparent that the so-called policy function and its fixed point play an outstanding role. Some applications to a statistical analysis of gender bias are also briefly discussed.

The problem of inferring the distribution of a random vector given that its norm is large requires modeling a homogeneous limiting density. We suggest an approach based on graphical models which is suitable for high-dimensional vectors. We introduce the notion of one-component regular variation to describe a function that is regularly varying in its first component. We extend the representation and Karamata's theorem to one-component regularly varying functions, probability distributions and densities, and explain why these results are fundamental in multivariate extreme-value theory. We then generalize the Hammersley–Clifford theorem to relate asymptotic conditional independence to a factorization of the limiting density, and use it to model multivariate tails.

Consider a random directed graph on n vertices with independent and identically distributed outdegrees with distribution F having mean μ, and destinations of arcs selected uniformly at random. We show that if μ > 1 then for large n there is very likely to be a unique giant strong component with proportionate size given as the product of two branching process survival probabilities, one with offspring distribution F and the other with Poisson offspring distribution with mean μ. If μ ≤ 1 there is very likely to be no giant strong component. We also extend this to allow for F varying with n.

In this paper we consider a variation of the full-information secretary problem where the random variables to be observed are independent but not necessary identically distributed. The main result is a sharp lower bound for the optimal win probability. Precisely, if X 1,...,X n are independent random variables with known continuous distributions and V n (X 1,...,Xn ):=supτℙ(X τ=M n ), where M n ≔max{X 1,...,X n } and the supremum is over all stopping times adapted to X 1,...,X n then V n (X 1,...,X n )≥(1-1/n)n-1, and this bound is attained. The method of proof consists in reducing the problem to that of a sequence of random variables taking at most two possible values, and then applying Bruss' sum-the-odds theorem, Bruss (2000). In order to obtain a sharp bound for each n, we improve Bruss' lower bound, Bruss (2003), for the sum-the-odds problem.

In this paper we connect Archimedean survival processes (ASPs) with the theory of Markov copulas. ASPs were introduced by Hoyle and Mengütürk (2013) to model the realized variance of two assets. We present some new properties of ASPs related to their dependency structure. We study weak and strong Markovian consistency properties of ASPs. An ASP is weak Markovian consistent, but generally not strong Markovian consistent. Our results contain necessary and sufficient conditions for an ASP to be strong Markovian consistent. These properties are closely related to the concept of Markov copulas, which is very useful in modelling different dependence phenomena. At the end we present possible applications.

We connect known results about diffusion limits of Markov chain Monte Carlo (MCMC) algorithms to the computer science notion of algorithm complexity. Our main result states that any weak limit of a Markov process implies a corresponding complexity bound (in an appropriate metric). We then combine this result with previously-known MCMC diffusion limit results to prove that under appropriate assumptions, the random-walk Metropolis algorithm in d dimensions takes O(d) iterations to converge to stationarity, while the Metropolis-adjusted Langevin algorithm takes O(d1/3) iterations to converge to stationarity.

The longest stretch L(n) of consecutive heads in n independent and identically distributed coin tosses is seen from the prism of large deviations. We first establish precise asymptotics for the moment generating function of L(n) and then show that there are precisely two large deviation principles, one concerning the behavior of the distribution of L(n) near its nominal value log1∕p n and one away from it. We discuss applications to inference and to logarithmic asymptotics of functionals of L(n).

We consider a continuous, infinitely divisible random field in Rd given as an integral of a kernel function with respect to a Lévy basis with convolution equivalent Lévy measure. For a large class of such random fields we compute the asymptotic probability that the supremum of the field exceeds the level x as x → ∞. Our main result is that the asymptotic probability is equivalent to the right tail of the underlying Lévy measure.

The sequential stochastic assignment problem (SSAP) allocates distinct workers with deterministic values to sequentially arriving tasks with stochastic parameters to maximize the expected total reward. In this paper we study an extension of the SSAP, in which the worker values are considered to be random variables, taking on new values upon each task arrival. Several SSAP models with different assumptions on the distribution of the worker values and closed-form expressions for optimal assignment policies are presented.

We consider the stability of robust scheduling policies for multiclass queueing networks. These are open networks with arbitrary routeing matrix and several disjoint groups of queues in which at most one queue can be served at a time. The arrival and potential service processes and routeing decisions at the queues are independent, stationary, and ergodic. A scheduling policy is called robust if it does not depend on the arrival and service rates nor on the routeing probabilities. A policy is called throughput-optimal if it makes the system stable whenever the parameters are such that the system can be stable. We propose two robust policies: longest-queue scheduling and a new policy called longest-dominating-queue scheduling. We show that longest-queue scheduling is throughput-optimal for two groups of two queues. We also prove the throughput-optimality of longest-dominating-queue scheduling when the network topology is acyclic, for an arbitrary number of groups and queues.

We consider a Markov chain (M n )n≥0 on the set ℕ0 of nonnegative integers which is eventually decreasing, i.e. ℙ{M n+1<M n   |  M n ≥a}=1 for some a∈ℕ and all n≥0. We are interested in the asymptotic behavior of the law of the stopping time T=T(a)≔inf{k∈ℕ0:  M k <a} under ℙn ≔ℙ (·  |  M 0=n) as n→∞. Assuming that the decrements of (M n )n≥0 given M 0=n possess a kind of stationarity for large n, we derive sufficient conditions for the convergence in the minimal L p -distance of ℙn (T−a n)∕b n∈·) to some nondegenerate, proper law and give an explicit form of the constants a n and b n .

We introduce and study a new model that we call the matching model. Items arrive one by one in a buffer and depart from it as soon as possible but by pairs. The items of a departing pair are said to be matched. There is a finite set of classes 𝒱 for the items, and the allowed matchings depend on the classes, according to a matching graph on 𝒱. Upon arrival, an item may find several possible matches in the buffer. This indeterminacy is resolved by a matching policy. When the sequence of classes of the arriving items is independent and identically distributed, the sequence of buffer-content is a Markov chain, whose stability is investigated. In particular, we prove that the model may be stable if and only if the matching graph is nonbipartite.

We consider the drawing with and without replacement models from numerical populations. For order statistics based on the above drawing schemes, we provide sharp upper bounds for their variances, expressed in the single observation variance units. We also characterize the populations for which the bounds are attained.

The susceptible→exposed→infectious→susceptible (SEIS) model is well known in mathematical epidemiology as a model of infection in which there is a latent period between the moment of infection and the onset of infectiousness. The compartment model is well studied, but the corresponding particle system has so far received no attention. For the particle system model in one spatial dimension, we give upper and lower bounds on the critical values, prove convergence of critical values in the limit of small and large latent time, and identify a limiting process to which the SEIS model converges in the limit of large latent time.

We investigate several aspects of a self-similar evolutionary process that builds a random bipolar network from building blocks that are themselves small bipolar networks. We characterize admissible outdegrees in the history of the evolution. We obtain the limit distribution of the polar degrees (when suitably scaled) characterized by its sequence of moments. We also obtain the asymptotic joint multivariate normal distribution of the number of nodes of small admissible outdegrees. Five possible substructures arise, and each has its own parameters (mean vector and covariance matrix) in the multivariate distribution. Several results are obtained by mapping bipolar networks into Pólya urns.

By considering the extreme behavior of bivariate random vectors with a polar representation R(u(T), v(T)), it is commonly assumed that the radial component R and the angular component T are stochastically independent. We investigate how to relax this rigid independence assumption such that conditional limit theorems can still be deduced. For this purpose, we introduce a novel measure for the dependence structure and present convenient criteria for validity of limit theorems possessing a geometrical meaning. Thus, our results verify a stability of the available limit results, which is essential in applications where the independence of the polar components is not necessarily present or exactly fulfilled.

In this paper we present closed-form expressions for the distribution of the virtual (actual) queueing time for the BMAP/R/1 and BMAP/D/1 queues, where `R' represents a class of distributions having rational Laplace‒Stieltjes transforms. The closed-form analysis is based on the roots of the underlying characteristic equation. Numerical aspects have been tested for a variety of arrival and service-time distributions and results are matched with those obtained using the matrix-analytic method (MAM). Further, a comparative study of computation time of the proposed method with the MAM has been carried out. Finally, we also present closed-form expressions for the distribution of the virtual (actual) system time. The proposed method is analytically quite simple and easy to implement.