We consider a host-parasite model for a population of cells that can be of two types, A or B, and exhibits unilateral reproduction: while a B-cell always splits into two cells of the same type, the two daughter cells of an A-cell can be of any type. The random mechanism that describes how parasites within a cell multiply and are then shared into the daughter cells is allowed to depend on the hosting mother cell as well as its daughter cells. Focusing on the subpopulation of A-cells and its parasites, our model differs from the single-type model recently studied by Bansaye (2008) in that the sharing mechanism may be biased towards one of the two types. Our main results are concerned with the nonextinctive case and provide information on the behavior, as n → ∞, of the number of A-parasites in generation n and the relative proportion of A- and B-cells in this generation which host a given number of parasites. As in Bansaye (2008), proofs will make use of a so-called random cell line which, when conditioned to be of type A, behaves like a branching process in a random environment.

Both small-world models of random networks with occasional long-range connections and gossip processes with occasional long-range transmission of information have similar characteristic behaviour. The long-range elements appreciably reduce the effective distances, measured in space or in time, between pairs of typical points. In this paper we show that their common behaviour can be interpreted as a product of the locally branching nature of the models. In particular, it is shown that both typical distances between points and the proportion of space that can be reached within a given distance or time can be approximated by formulae involving the limit random variable of the branching process.

Consider a generic data unit of random size L that needs to be transmitted over a channel of unit capacity. The channel availability dynamic is modeled as an independent and identically distributed sequence {A, Ai}i≥1 that is independent of L. During each period of time that the channel becomes available, say Ai, we attempt to transmit the data unit. If L <Ai, the transmission is considered successful; otherwise, we wait for the next available period Ai+1 and attempt to retransmit the data from the beginning. We investigate the asymptotic properties of the number of retransmissions N and the total transmission time T until the data is successfully transmitted. In the context of studying the completion times in systems with failures where jobs restart from the beginning, it was first recognized by Fiorini, Sheahan and Lipsky (2005) and Sheahan, Lipsky, Fiorini and Asmussen (2006) that this model results in power-law and, in general, heavy-tailed delays. The main objective of this paper is to uncover the detailed structure of this class of heavy-tailed distributions induced by retransmissions. More precisely, we study how the functional relationship ℙ[L>x]-1 ≈ Φ (ℙ[A>x]-1) impacts the distributions of N and T; the approximation ‘≈’ will be appropriately defined in the paper based on the context. Depending on the growth rate of Φ(·), we discover several criticality points that separate classes of different functional behaviors of the distribution of N. For example, we show that if log(Φ(n)) is slowly varying then log(1/ℙ[N>n]) is essentially slowly varying as well. Interestingly, if log(Φ(n)) grows slower than e√(logn) then we have the asymptotic equivalence log(ℙ[N>n]) ≈ - log(Φ(n)). However, if log(Φ(n)) grows faster than e√(logn), this asymptotic equivalence does not hold and admits a different functional form. Similarly, different types of distributional behavior are shown for moderately heavy tails (Weibull distributions) where log(ℙ[N>n]) ≈ -(log Φ(n))1/(β+1), assuming that log Φ(n) ≈ nβ, as well as the nearly exponential ones of the form log(ℙ[N>n]) ≈ -n/(log n)1/γ, γ>0, when Φ(·) grows faster than two exponential scales log log (Φ(n)) ≈ nγ.

In the paper we consider the density functions of random variables that can be written as integrals of exponential functions of Gaussian random fields. In particular, we provide closed-form asymptotic bounds for the density functions and, under smoothness conditions, we derive exact tail approximations of the density functions.

Multivariate regular variation plays a role in assessing tail risk in diverse applications such as finance, telecommunications, insurance, and environmental science. The classical theory, being based on an asymptotic model, sometimes leads to inaccurate and useless estimates of probabilities of joint tail regions. This problem can be partly ameliorated by using hidden regular variation (see Resnick (2002) and Mitra and Resnick (2011)). We offer a more flexible definition of hidden regular variation that provides improved risk estimates for a larger class of tail risk regions.

When the search algorithm QuickSelect compares keys during its execution in order to find a key of target rank, it must operate on the keys' representations or internal structures, which were ignored by the previous studies that quantified the execution cost for the algorithm in terms of the number of required key comparisons. In this paper we analyze running costs for the algorithm that take into account not only the number of key comparisons, but also the cost of each key comparison. We suppose that keys are represented as sequences of symbols generated by various probabilistic sources and that QuickSelect operates on individual symbols in order to find the target key. We identify limiting distributions for the costs, and derive integral and series expressions for the expectations of the limiting distributions. These expressions are used to recapture previously obtained results on the number of key comparisons required by the algorithm.

The signature of a system is defined as the vector whose ith element is the probability that the system fails concurrently with the ith component failure. The signature vector is known to be a distribution-free measure and a representation of the system's survival function has been developed in terms of the system's signature. The present work is devoted to the study of the joint distribution of lifetimes of pairs of systems with shared components. Here, a new distribution-free measure, the ‘joint bivariate signature’, of a pair of systems with shared components is defined, and a new representation theorem for the joint survival function of the system lifetimes is established. The theorem is shown to facilitate the study of the dependence between systems and the comparative performance of two pairs of such systems.

We propose and study a number of layer methods for Navier‒Stokes equations (NSEs) with spatial periodic boundary conditions. The methods are constructed using probabilistic representations of solutions to NSEs and exploiting ideas of the weak sense numerical integration of stochastic differential equations. Despite their probabilistic nature, the layer methods are nevertheless deterministic.

We study the Bayesian problems of detecting a change in the drift rate of an observable diffusion process with linear and exponential penalty costs for a detection delay. The optimal times of alarms are found as the first times at which the weighted likelihood ratios hit stochastic boundaries depending on the current observations. The proof is based on the reduction of the initial problems into appropriate three-dimensional optimal stopping problems and the analysis of the associated parabolic-type free-boundary problems. We provide closed-form estimates for the value functions and the boundaries, under certain nontrivial relations between the coefficients of the observable diffusion.

We study a risk process with dividend barrier b where the claims arrive according to a Markovian additive process (MAP). For spectrally negative MAPs, we present linear equations for the expected discounted dividends and the expected discounted penalty function. We apply results for the first exit times of spectrally negative Lévy processes and change-of-measure techniques. Explicit expressions are given when there are positive and negative claims, with phase-type distribution.

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.

As a version of the secretary problem, Ferguson, Hardwick and Tamaki (1992) considered an optimal stopping problem called the duration problem. The basic duration problem is the classical duration problem, in which the objective is to maximize the time of possession of a relatively best object when a known number of rankable objects appear in random order. In this paper we generalize this classical problem in two directions by allowing the number N (of available objects) to be a random variable with a known upper bound n and also allowing the objects to appear in accordance with Bernoulli trials. Two models can be considered for our random horizon duration problem according to whether the planning horizon is N or n. Since the form of the optimal rule is in general complicated, our main concern is to give to each model a sufficient condition for the optimal rule to be simple. For N having a uniform, generalized uniform, or curtailed geometric distribution, the optimal rule is shown to be simple in the so-called secretary case. The asymptotic results, as n → ∞, will also be given for these priors.

In this paper we adopt the perturbation approach of Landriault, Renaud and Zhou (2011) to find expressions for the joint Laplace transforms of occupation times for time-homogeneous diffusion processes. The expressions are in terms of solutions to the associated differential equations. These Laplace transforms are applied to study ruin-related problems for several classes of diffusion risk processes.

In this paper we study absorbing continuous-time Markov decision processes in Polish state spaces with unbounded transition and cost rates, and history-dependent policies. The performance measure is the expected total undiscounted costs. For the unconstrained problem, we show the existence of a deterministic stationary optimal policy, whereas, for the constrained problems with N constraints, we show the existence of a mixed stationary optimal policy, where the mixture is over no more than N+1 deterministic stationary policies. Furthermore, the strong duality result is obtained for the associated linear programs.

An asymptotic model for the extreme behavior of certain Markov chains is the ‘tail chain’. Generally taking the form of a multiplicative random walk, it is useful in deriving extremal characteristics, such as point process limits. We place this model in a more general context, formulated in terms of extreme value theory for transition kernels, and extend it by formalizing the distinction between extreme and nonextreme states. We make the link between the update function and transition kernel forms considered in previous work, and we show that the tail chain model leads to a multivariate regular variation property of the finite-dimensional distributions under assumptions on the marginal tails alone.

We present a new algorithm to discretize a decoupled forward‒backward stochastic differential equation driven by a pure jump Lévy process (FBSDEL for short). The method consists of two steps. In the first step we approximate the FBSDEL by a forward‒backward stochastic differential equation driven by a Brownian motion and Poisson process (FBSDEBP for short), in which we replace the small jumps by a Brownian motion. Then, we prove the convergence of the approximation when the size of small jumps ε goes to 0. In the second step we obtain the Lp-Hölder continuity of the solution of the FBSDEBP and we construct two numerical schemes for this FBSDEBP. Based on the Lp-Hölder estimate, we prove the convergence of the scheme when the number of time steps n goes to ∞. Combining these two steps leads to the proof of the convergence of numerical schemes to the solution of FBSDEs driven by pure jump Lévy processes.

In this paper we develop an asymptotic theory of aggregated linear processes, and determine in particular the limit distribution of a large class of linear and nonlinear functionals of such processes. Given a sample {Y1(N),…,Yn(N)} of the normalized N-fold aggregated process, we describe the limiting behavior of statistics TN,n= TN,n(Y1(N),…, Yn(N)) in both of the cases n/N(n) → 0 and N(n)/n → 0, assuming either a ‘limiting long- or short-memory’ condition on the underlying linear process.

We introduce two stochastic chemostat models consisting of a coupled population-nutrient process reflecting the interaction between the nutrient and the bacteria in the chemostat with finite volume. The nutrient concentration evolves continuously but depends on the population size, while the population size is a birth-and-death process with coefficients depending on time through the nutrient concentration. The nutrient is shared by the bacteria and creates a regulation of the bacterial population size. The latter and the fluctuations due to the random births and deaths of individuals make the population go almost surely to extinction. Therefore, we are interested in the long-time behavior of the bacterial population conditioned to nonextinction. We prove the global existence of the process and its almost-sure extinction. The existence of quasistationary distributions is obtained based on a general fixed-point argument. Moreover, we prove the absolute continuity of the nutrient distribution when conditioned to a fixed number of individuals and the smoothness of the corresponding densities.

In this paper we introduce discrete-time semi-Markov random evolutions (DTSMREs) and study asymptotic properties, namely, averaging, diffusion approximation, and diffusion approximation with equilibrium by the martingale weak convergence method. The controlled DTSMREs are introduced and Hamilton–Jacobi–Bellman equations are derived for them. The applications here concern the additive functionals (AFs), geometric Markov renewal chains (GMRCs), and dynamical systems (DSs) in discrete time. The rates of convergence in the limit theorems for DTSMREs and AFs, GMRCs, and DSs are also presented.

We present two iterative methods for computing the global and partial extinction probability vectors for Galton-Watson processes with countably infinitely many types. The probabilistic interpretation of these methods involves truncated Galton-Watson processes with finite sets of types and modified progeny generating functions. In addition, we discuss the connection of the convergence norm of the mean progeny matrix with extinction criteria. Finally, we give a sufficient condition for a population to become extinct almost surely even though its population size explodes on the average, which is impossible in a branching process with finitely many types. We conclude with some numerical illustrations for our algorithmic methods.