In the paper we consider the following modification of a discrete-time branching process with stationary immigration. In each generation a binomially distributed subset of the population will be observed. The number of observed individuals constitute a partially observed branching process. After inspection both observed and unobserved individuals may change their offspring distributions. In the subcritical case we investigate the possibility of using the known estimators for the offspring mean and for the mean of the stationary-limiting distribution of the process when the observation of the population sizes is restricted. We prove that, if both the population and the number of immigrants are partially observed, the estimators are still strongly consistent. We also prove that the `skipped' version of the estimator for the offspring mean is asymptotically normal and the estimator of the stationary distribution's mean is asymptotically normal under additional assumptions.

The extremal behaviour of a Markov chain is typically characterised by its tail chain. For asymptotically dependent Markov chains, existing formulations fail to capture the full evolution of the extreme event when the chain moves out of the extreme tail region, and, for asymptotically independent chains, recent results fail to cover well-known asymptotically independent processes, such as Markov processes with a Gaussian copula between consecutive values. We use more sophisticated limiting mechanisms that cover a broader class of asymptotically independent processes than current methods, including an extension of the canonical Heffernan‒Tawn normalisation scheme, and reveal features which existing methods reduce to a degenerate form associated with nonextreme states.

We are interested in a large queue in a GI/G/k queue with heterogeneous servers. For this, we consider tail asymptotics and weak limit approximations for the stationary distribution of its queue length process in continuous time under a stability condition. Here, two weak limit approximations are considered. One is when the variances of the interarrival and/or service times are bounded, and the other is when they become large. Both require a heavy-traffic condition. Tail asymptotics and heavy-traffic approximations have been separately studied in the literature. We develop a unified approach based on a martingale produced by a good test function for a Markov process to answer both problems.

We are interested in the rate of convergence of a subordinate Markov process to its invariant measure. Given a subordinator and the corresponding Bernstein function (Laplace exponent), we characterize the convergence rate of the subordinate Markov process; the key ingredients are the rate of convergence of the original process and the (inverse of the) Bernstein function. At a technical level, the crucial point is to bound three types of moment (subexponential, algebraic, and logarithmic) for subordinators as time t tends to ∞. We also discuss some concrete models and we show that subordination can dramatically change the speed of convergence to equilibrium.

We study the impact of a hard selective sweep on the genealogy of partially linked neutral loci in the vicinity of the positively selected allele. We consider a sexual population of stochastically varying size and, focusing on two neighboring loci, derive an approximate formula for the neutral genealogy of a sample of individuals taken at the end of the sweep. Individuals are characterized by ecological parameters depending on their genetic type, and governing their growth rate and interactions with other individuals (competition). As a consequence, the `fitness' of an individual depends on the population state and is not an intrinsic characteristic of individuals. We provide a deep insight into the dynamics of the mutant and wild-type populations during the different stages of a selective sweep.

In Bingham and Doney (1988) the authors presented the applied probability community with a question which is very simply stated, yet is extremely difficult to solve: what is the distribution of the quadrant occupation time of planar Brownian motion? In this paper we study an alternate formulation of this long-standing open problem: let X(t), Y(t) t≥0, be standard Brownian motions starting at x, y, respectively. Find the distribution of the total time T=Leb{t∈[0,1]: X(t)×Y(t)>0}, when x=y=0, i.e. the occupation time of the union of the first and third quadrants. If two adjacent quadrants are used, the problem becomes much easier and the distribution of T follows the arcsine law.

The (conditional or unconditional) distribution of the continuous scan statistic in a one-dimensional Poisson process may be approximated by that of a discrete analogue via time discretization (to be referred to as the discrete approximation). Using a change of measure argument, we derive the first-order term of the discrete approximation which involves some functionals of the Poisson process. Richardson's extrapolation is then applied to yield a corrected (second-order) approximation. Numerical results are presented to compare various approximations.

In this paper we aim to apply simple actuarial methods to build an insurance plan protecting against an epidemic risk in a population. The studied model is an extended SIR epidemic in which the removal and infection rates may depend on the number of registered removals. The costs due to the epidemic are measured through the expected epidemic size and infectivity time. The premiums received during the epidemic outbreak are measured through the expected susceptibility time. Using martingale arguments, a method by recursion is developed to calculate the cost components and the corresponding premium levels in this extended epidemic model. Some numerical examples illustrate the effect of removals and the premium calculation in an insurance plan.

We establish a local martingale M associate with f(X,Y) under some restrictions on f, where Y is a process of bounded variation (on compact intervals) and either X is a jump diffusion (a special case being a Lévy process) or X is some general (càdlàg metric-space valued) Markov process. In the latter case, f is restricted to the form f(x,y)=∑k=1 K ξk (x)ηk (y). This local martingale unifies both Dynkin's formula for Markov processes and the Lebesgue–Stieltjes integration (change of variable) formula for (right-continuous) functions of bounded variation. For the jump diffusion case, when further relatively easily verifiable conditions are assumed, then this local martingale becomes an L 2-martingale. Convergence of the product of this Martingale with some deterministic function ( of time ) to 0 both in L 2 and almost sure is also considered and sufficient conditions for functions for which this happens are identified.

We study a representation for the inverse transform of the generalised Fourier–Feynman transform on the function space $C_{a,b}[0,T]$ which is induced by a generalised Brownian motion process. To do this, we define a transform via the concept of the convolution product of functionals on $C_{a,b}[0,T]$ . We establish that the composition of these transforms acts like an inverse generalised Fourier–Feynman transform and that the transforms are vector space automorphisms of a vector space ${\mathcal{E}}(C_{a,b}[0,T])$ . The vector space ${\mathcal{E}}(C_{a,b}[0,T])$ consists of exponential-type functionals on $C_{a,b}[0,T]$ .

The shapes of branching trees have been linked to disease transmission patterns. In this paper we use the general Crump‒Mode‒Jagers branching process to model an outbreak of an infectious disease under mild assumptions. Introducing a new class of characteristic functions, we are able to derive a formula for the limit of the frequency of the occurrences of a given shape in a general tree. The computational challenges concerning the evaluation of this formula are in part overcome using the jumping chronological contour process. We apply the formula to derive the limit of the frequency of cherries, pitchforks, and double cherries in the constant-rate birth‒death model, and the frequency of cherries under a nonconstant death rate.

We study the distributional properties of jumps in a continuous-state branching process with immigration. In particular, a representation is given for the distribution of the first jump time of the process with jump size in a given Borel set. From this result we derive a characterization for the distribution of the local maximal jump of the process. The equivalence of this distribution and the total Lévy measure is then studied. For the continuous-state branching process without immigration, we also study similar problems for its global maximal jump.

A multiclass c-server retrial queueing system in which customers arrive according to a class-dependent Markovian arrival process (MAP) is considered. Service and retrial times follow class-dependent phase-type (PH) distributions with the further assumption that PH distributions of retrial times are acyclic. A necessary and sufficient condition for ergodicity is obtained from criteria based on drifts. The infinite state space of the model is truncated with an appropriately chosen Lyapunov function. The truncated model is described as a multidimensional Markov chain, and a Kronecker representation of its generator matrix is numerically analyzed.

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.

We describe in detail the speed of `coming down from infinity' for birth-and-death processes which eventually become extinct. Under general assumptions on the birth-and-death rates, we firstly determine the behavior of the successive hitting times of large integers. We identify two different regimes depending on whether the mean time for the process to go from n+1 to n is negligible or not compared to the mean time to reach n from ∞. In the first regime, the coming down from infinity is very fast and the convergence is weak. In the second regime, the coming down from infinity is gradual and a law of large numbers and a central limit theorem for the hitting times sequence hold. By an inversion procedure, we deduce that the process is almost surely equivalent to a nonincreasing function when the time goes to 0. Our results are illustrated by several examples including applications to population dynamics and population genetics. The particular case where the death rate varies regularly is studied in detail.

Branching processes in random environments have been widely studied and applied to population growth systems to model the spread of epidemics, infectious diseases, cancerous tumor growth, and social network traffic. However, Ebola virus, tuberculosis infections, and avian flu grow or change at rates that vary with time—at peak rates during pandemic time periods, while at low rates when near extinction. The branching processes in generalized autoregressive conditional environments we propose provide a novel approach to branching processes that allows for such time-varying random environments and instances of peak growth and near extinction-type rates. Offspring distributions we consider to illustrate the model include the generalized Poisson, binomial, and negative binomial integer-valued GARCH models. We establish conditions on the environmental process that guarantee stationarity and ergodicity of the mean offspring number and environmental processes and provide equations from which their variances, autocorrelation, and cross-correlation functions can be deduced. Furthermore, we present results on fundamental questions of importance to these processes—the survival-extinction dichotomy, growth behavior, necessary and sufficient conditions for noncertain extinction, characterization of the phase transition between the subcritical and supercritical regimes, and survival behavior in each phase and at criticality.

We give a criterion for unlimited growth with positive probability for a large class of multidimensional stochastic models. As a by-product, we recover the necessary and sufficient conditions for recurrence and transience for critical multitype Galton–Watson with immigration processes and also significantly improve some results on multitype size-dependent Galton–Watson processes.

We define and analyze a coalescent process as a recursive box-filling process whose genealogy is given by an ancestral time-reversed, time-inhomogeneous Bienyamé‒Galton‒Watson process. Special interest is on the expected size of a typical box and its probability of being empty. Special cases leading to exact asymptotic computations are investigated when the coalescing mechanisms are either linear fractional or quadratic.

Comparison results are given for time-inhomogeneous Markov processes with respect to function classes with induced stochastic orderings. The main result states the comparison of two processes, provided that the comparability of their infinitesimal generators as well as an invariance property of one process is assumed. The corresponding proof is based on a representation result for the solutions of inhomogeneous evolution problems in Banach spaces, which extends previously known results from the literature. Based on this representation, an ordering result for Markov processes induced by bounded and unbounded function classes is established. We give various applications to time-inhomogeneous diffusions, to processes with independent increments, and to Lévy-driven diffusion processes.

We propose a stochastic model describing a process of awareness, evaluation, and decision making by agents on the d-dimensional integer lattice. Each agent may be in any of the three states belonging to the set {0, 1, 2. In this model 0 stands for ignorants, 1 for aware, and 2 for adopters. Aware and adopters inform its nearest ignorant neighbors about a new product innovation at rate λ. At rate α an agent in aware state becomes an adopter due to the influence of adopters' neighbors. Finally, aware and adopters forget the information about the new product, thus becoming ignorant, at rate 1. Our purpose is to analyze the influence of the parameters on the qualitative behavior of the process. We obtain sufficient conditions under which the innovation diffusion (and adoption) either becomes extinct or propagates through the population with positive probability.