Hanski's incidence function model is one of the most widely used metapopulation models in ecology. It models the presence/absence of a species at spatially distinct habitat patches as a discrete-time Markov chain whose transition probabilities are determined by the physical landscape. In this analysis, the limiting behaviour of the model is studied as the number of patches increases and the size of the patches decreases. Two different limiting cases are identified depending on whether or not the metapopulation is initially near extinction. Basic properties of the limiting models are derived.

We characterise the class of distributions of random stochastic matrices X with the property that the products X(n)X(n − 1) · · · X(1) of independent and identically distributed copies X(k) of X converge almost surely as n → ∞ and the limit is Dirichlet distributed. This extends a result by Chamayou and Letac (1994) and is illustrated by several examples that are of interest in applications.

The generalized Pólya urn has been extensively studied and is widely applied in many disciplines. An important application of urn models is in the development of randomized treatment allocation schemes in clinical studies. The randomly reinforced urn was recently proposed, but, although the model has some intuitively desirable properties, it lacks theoretical justification. In this paper we obtain important asymptotic properties for multicolor reinforced urn models. We derive results for the rate of convergence of the number of patients assigned to each treatment under a set of minimum required conditions and provide the distributions of the limits. Furthermore, we study the asymptotic behavior for the nonhomogeneous case.

As the name suggests, the family of general error distributions has been used to model nonnormal errors in a variety of situations. In this article we show that the asymptotic distribution of linearly normalized partial maxima of random observations from the general error distributions is Gumbel when the parameter of these distributions lies in the interval (0, 1). Our result fills a gap in the literature. We also establish the corresponding density convergence, obtain an asymptotic distribution of the partial maxima under power normalization, and state and prove a strong law. We also study the asymptotic behaviour of observations near the partial maxima and the sum of such observations.

Let {Zn}n≥0 be a random walk with a negative drift and independent and identically distributed increments with heavy-tailed distribution, and let M = supn≥0Zn be its supremum. Asmussen and Klüppelberg (1996) considered the behavior of the random walk given that M > x for large x, and obtained a limit theorem, as x → ∞, for the distribution of the quadruple that includes the time τ = τ(x) to exceed level x, position Zτ at this time, position Zτ-1 at the prior time, and the trajectory up to it (similar results were obtained for the Cramér-Lundberg insurance risk process). We obtain here several extensions of this result to various regenerative-type models and, in particular, to the case of a random walk with dependent increments. Particular attention is given to describing the limiting conditional behavior of τ. The class of models includes Markov-modulated models as particular cases. We also study fluid models, the Björk-Grandell risk process, give examples where the order of τ is genuinely different from the random walk case, and discuss which growth rates are possible. Our proofs are purely probabilistic and are based on results and ideas from Asmussen, Schmidli and Schmidt (1999), Foss and Zachary (2002), and Foss, Konstantopoulos and Zachary (2007).

We consider a broad class of fair leader election algorithms, and study the duration of contestants (the number of rounds a randomly selected contestant stays in the competition) and the overall cost of the algorithm. We give sufficient conditions for the duration to have a geometric limit distribution (a perpetuity built from Bernoulli random variables), and for the limiting distribution of the total cost (after suitable normalization) to be a perpetuity. For the duration, the proof is established via convergence (to 0) of the first-order Wasserstein distance from the geometric limit. For the normalized overall cost, the method of proof is also convergence of the first-order Wasserstein distance, augmented with an argument based on a contraction mapping in the first-order Wasserstein metric space to show that the limit approaches a unique fixed-point solution of a perpetuity distributional equation. The use of these two steps is commonly referred to as the contraction method.

In the present paper, we study the evolution of an overloaded cyclic polling model that starts empty. Exploiting a connection with multitype branching processes, we derive fluid asymptotics for the joint queue length process. Under passage to the fluid dynamics, the server switches between the queues infinitely many times in any finite time interval causing frequent oscillatory behavior of the fluid limit in the neighborhood of zero. Moreover, the fluid limit is random. In addition, we suggest a method that establishes finiteness of moments of the busy period in an M/G/1 queue.

We study an urn model introduced in the paper of Chen and Wei (2005), where at each discrete time step m balls are drawn at random from the urn containing colors white and black. Balls are added to the urn according to the inspected colors, generalizing the well known Pólya-Eggenberger urn model, case m = 1. We provide exact expressions for the expectation and the variance of the number of white balls after n draws, and determine the structure of higher moments. Moreover, we discuss extensions to more than two colors. Furthermore, we introduce and discuss a new urn model where the sampling of the m balls is carried out in a step-by-step fashion, and also introduce a generalized Friedman's urn model.

Using a characterizing equation for the beta distribution, Stein's method is applied to obtain bounds of the optimal order for the Wasserstein distance between the distribution of the scaled number of white balls drawn from a Pólya-Eggenberger urn and its limiting beta distribution. The bound is computed by making a direct comparison between characterizing operators of the target and the beta distribution, the former derived by extending Stein's density approach to discrete distributions. In addition, refinements are given to Döbler's (2012) result for the arcsine approximation for the fraction of time a simple random walk of even length spends positive, and so also to the distributions of its last return time to 0 and its first visit to its terminal point, by supplying explicit constants to the present Wasserstein bound and also demonstrating that its rate is of the optimal order.

We consider possibly nonlinear distributional fixed-point equations on weighted branching trees, which include the well-known linear branching recursion. In Jelenković and Olvera-Cravioto (2012), an implicit renewal theorem was developed that enables the characterization of the power-tail asymptotics of the solutions to many equations that fall into this category. In this paper we complement the analysis in our 2012 paper to provide the corresponding rate of convergence.

Self-exciting point processes (SEPPs), or Hawkes processes, have found applications in a wide range of fields, such as epidemiology, seismology, neuroscience, engineering, and more recently financial econometrics and social interactions. In the traditional SEPP models, the baseline intensity is assumed to be a constant. This has restricted the application of SEPPs to situations where there is clearly a self-exciting phenomenon, but a constant baseline intensity is inappropriate. In this paper, to model point processes with varying baseline intensity, we introduce SEPP models with time-varying background intensities (SEPPVB, for short). We show that SEPPVB models are competitive with autoregressive conditional SEPP models (Engle and Russell 1998) for modeling ultra-high frequency data. We also develop asymptotic theory for maximum likelihood estimation based inference of parametric SEPP models, including SEPPVB. We illustrate applications to ultra-high frequency financial data analysis, and we compare performance with the autoregressive conditional duration models.

Many regenerative arguments in stochastic processes use random times which are akin to stopping times, but which are determined by the future as well as the past behaviour of the process of interest. Such arguments based on ‘conditioning on the future’ are usually developed in an ad-hoc way in the context of the application under consideration, thereby obscuring the underlying structure. In this paper we give a simple, unified, and more general treatment of such conditioning theory. We further give a number of novel applications to various particle system models, in particular to various flavours of contact processes and to infinite-bin models. We give a number of new results for existing and new models. We further make connections with the theory of Harris ergodicity.

The goal of this paper is to prove a result conjectured in Föllmer and Schachermayer (2007) in a slightly more general form. Suppose that S is a continuous semimartingale and satisfies a large deviations estimate; this is a particular growth condition on the mean-variance tradeoff process of S. We show that S then allows asymptotic exponential arbitrage with exponentially decaying failure probability, which is a strong and quantitative form of long-term arbitrage. In contrast to Föllmer and Schachermayer (2007), our result does not assume that S is a diffusion, nor does it need any ergodicity assumption.

For a family of linear preferential attachment graphs, we provide rates of convergence for the total variation distance between the degree of a randomly chosen vertex and an appropriate power law distribution as the number of vertices tends to ∞. Our proof uses a new formulation of Stein's method for the negative binomial distribution, which stems from a distributional transformation that has the negative binomial distributions as the only fixed points.

We consider the problem of efficient simulation estimation of the density function at the tails, and the probability of large deviations for a sum of independent, identically distributed (i.i.d.), light-tailed, and nonlattice random vectors. The latter problem besides being of independent interest, also forms a building block for more complex rare event problems that arise, for instance, in queueing and financial credit risk modeling. It has been extensively studied in the literature where state-independent, exponential-twisting-based importance sampling has been shown to be asymptotically efficient and a more nuanced state-dependent exponential twisting has been shown to have a stronger bounded relative error property. We exploit the saddle-point-based representations that exist for these rare quantities, which rely on inverting the characteristic functions of the underlying random vectors. These representations reduce the rare event estimation problem to evaluating certain integrals, which may via importance sampling be represented as expectations. Furthermore, it is easy to identify and approximate the zero-variance importance sampling distribution to estimate these integrals. We identify such importance sampling measures and show that they possess the asymptotically vanishing relative error property that is stronger than the bounded relative error property. To illustrate the broader applicability of the proposed methodology, we extend it to develop an asymptotically vanishing relative error estimator for the practically important expected overshoot of sums of i.i.d. random variables.

We discuss tail behaviors, subexponentiality, and the extreme value distribution of logarithmic skew-normal random variables. With optimal normalized constants, the asymptotic expansion of the distribution of the normalized maximum of logarithmic skew-normal random variables is derived. We show that the convergence rate of the distribution of the normalized maximum to the Gumbel extreme value distribution is proportional to 1/(log n)1/2.

In this paper we introduce a simple risk model with delayed claims, an extension of the classical Poisson model. The claims are assumed to arrive according to a Poisson process and claims follow a light-tailed distribution, and each loss payment of the claims will be settled with a random period of delay. We obtain asymptotic expressions for the ruin probability by exploiting a connection to Poisson models that are not time homogeneous. A finer asymptotic formula is obtained for the special case of exponentially delayed claims and an exact formula is obtained when the claims are also exponentially distributed.

Consider a circle with perimeter N > 1 on which k < N segments of length 1 are sampled in an independent and identically distributed manner. In this paper we study the probability π (k,N) that these k segments do not overlap; the density φ(·) of the position of the disks on the circle is arbitrary (that is, it is not necessarily assumed uniform). Two scaling regimes are considered. In the first we set k≡ a√N, and it turns out that the probability of interest converges (N→ ∞) to an explicitly given positive constant that reflects the impact of the density φ(·). In the other regime k scales as aN, and the nonoverlap probability decays essentially exponentially; we give the associated decay rate as the solution to a variational problem. Several additional ramifications are presented.

The Hawkes process is a self-exciting point process with clustering effect whose intensity depends on its entire past history. It has wide applications in neuroscience, finance, and many other fields. In this paper we obtain a functional central limit theorem for the nonlinear Hawkes process. Under the same assumptions, we also obtain a Strassen's invariance principle, i.e. a functional law of the iterated logarithm.

The asymptotic variance of the sample mean of a homogeneous Poisson marked point process has been studied in the literature, but confusion has arisen as to the correct expression due to some technical intricacies. This note sets the record straight with regards to the variance of the sample mean. In addition, a central limit theorem in the general d-dimensional case is also established.