In this paper we consider reflected diffusions with positive and negative jumps, constrained to lie in the nonnegative orthant of ℝn. We allow for the drift and diffusion coefficients, as well as for the directions of reflection, to be random fields over time and space. We provide a boundary behavior characterization, generalizing known results in the nonrandom coefficients and constant directions of the reflection case. In particular, the regulator processes are related to semimartingale local times at the boundaries, and they are shown not to charge the times the process expends at the intersection of boundary faces. Using the boundary results, we extend the conditions for product-form distributions in the stationary regime to the case when the drift and diffusion coefficients, as well as the directions of reflection, are random fields over space.

Let S0 := 0 and Sk := ξ1 + ··· + ξk for k ∈ ℕ := {1, 2, …}, where {ξk : k ∈ ℕ} are independent copies of a random variable ξ with values in ℕ and distribution pk := P{ξ = k}, k ∈ ℕ. We interpret the random walk {Sk : k = 0, 1, 2, …} as a particle jumping to the right through integer positions. Fix n ∈ ℕ and modify the process by requiring that the particle is bumped back to its current state each time a jump would bring the particle to a state larger than or equal to n. This constraint defines an increasing Markov chain {Rk(n) : k = 0, 1, 2, …} which never reaches the state n. We call this process a random walk with barrier n. Let Mn denote the number of jumps of the random walk with barrier n. This paper focuses on the asymptotics of Mn as n tends to ∞. A key observation is that, under p1 > 0, {Mn : n ∈ ℕ} satisfies the distributional recursion M1 = 0 and for n = 2, 3, …, where In is independent of M2, …, Mn−1 with distribution P{In = k} = pk / (p1 + ··· + pn−1), k ∈ {1, …, n − 1}. Depending on the tail behavior of the distribution of ξ, several scalings for Mn and corresponding limiting distributions come into play, including stable distributions and distributions of exponential integrals of subordinators. The methods used in this paper are mainly probabilistic. The key tool is to compare (couple) the number of jumps, Mn, with the first time, Nn, when the unrestricted random walk {Sk : k = 0, 1, …} reaches a state larger than or equal to n. The results are applied to derive the asymptotics of the number of collision events (that take place until there is just a single block) for β(a, b)-coalescent processes with parameters 0 < a < 2 and b = 1.

We consider the M/M/∞ queue with m primary servers and infinitely many secondary servers. All the servers are numbered and ordered. An arriving customer takes the lowest available server. We define the wasted spaces as the difference between the highest numbered occupied server and the total number of occupied servers. Letting ρ = λ0/μ be the ratio of arrival to service rates, we study the probability distribution of the wasted spaces asymptotically for ρ → ∞. We also give some numerical results and the tail behavior for ρ = O(1).

This paper deals with censoring of infinite-state banded Markov chains. Censoring involves reducing the time spent in states outside a certain set of states to 0 without affecting the number of visits within this set. We show that, if all states are transient, there is, besides the standard censored Markov chain, a nonstandard censored Markov chain which is stochastic. Both the stochastic and the substochastic solutions are found by censoring a sequence of finite transition matrices. If all matrices in the sequence are stochastic, the stochastic solution arises in the limit, whereas the substochastic solution arises if the matrices in the sequence are substochastic. We also show that, if the Markov chain is recurrent, the only solution is the stochastic solution. Censoring is particularly fruitful when applied to quasi-birth-and-death (QBD) processes. It turns out that key matrices in such processes are not unique, a fact that has been observed by several authors. We note that the stochastic solution is important for the analysis of finite queues.

We consider a multiserver queue in the Halfin-Whitt regime: as the number of servers n grows without a bound, the utilization approaches 1 from below at the rate Assuming that the service time distribution is lattice valued with a finite support, we characterize the limiting scaled stationary queue length distribution in terms of the stationary distribution of an explicitly constructed Markov chain. Furthermore, we obtain an explicit expression for the critical exponent for the moment generating function of a limiting stationary queue length. This exponent has a compact representation in terms of three parameters: the amount of spare capacity and the coefficients of variation of interarrival and service times. Interestingly, it matches an analogous exponent corresponding to a single-server queue in the conventional heavy-traffic regime.

We consider stochastic replicator processes for games that are composed of finitely many trials. Several general results on the relation between Nash equilibria and the long-run behaviour of the stochastic processes are proved. In particular, a sufficient condition is given for almost sure convergence to a state where everyone plays in every trial a strict Nash equilibrium. The results are applied to multiple-trial conflicts based on wars of attrition and on sperm competition games with fair raffles, respectively.

Inference of evolutionary trees and rates from biological sequences is commonly performed using continuous-time Markov models of character change. The Markov process evolves along an unknown tree while observations arise only from the tips of the tree. Rate heterogeneity is present in most real data sets and is accounted for by the use of flexible mixture models where each site is allowed its own rate. Very little has been rigorously established concerning the identifiability of the models currently in common use in data analysis, although nonidentifiability was proven for a semiparametric model and an incorrect proof of identifiability was published for a general parametric model (GTR + Γ + I). Here we prove that one of the most widely used models (GTR + Γ) is identifiable for generic parameters, and for all parameter choices in the case of four-state (DNA) models. This is the first proof of identifiability of a phylogenetic model with a continuous distribution of rates.

We analyze asymptotically a differential-difference equation that arises in a Markov-modulated fluid model. We use singular perturbation methods to analyze the problem with appropriate scalings of the two state variables. In particular, the ray method and asymptotic matching are used.

In this paper a limit theorem is proved that establishes conditions under which the distribution of the difference between the size of the initial event in a random sequence, modeled as a finite point process, and the largest subsequent event approaches a limiting form independent of the size of the initial event. The underlying assumptions are that the sizes of the individual events follow an exponential distribution, that the expected total number of events increases exponentially with the size of the initial event, and that the structure of the sequence approximates that of a Poisson process. Particular cases to which the results apply include sequences of independent and identically distributed exponential variables, and the epidemic-type aftershock (ETAS) branching process model in the subcritical case. In all these cases the form of the limit distribution is shown to be that of a double exponential (type-I extreme value distribution). In sampling from a family of aftershock sequences, with possibly different underlying parameters, the limit distribution is a mixture of such double exponential distributions. The conditions for the simple limit to exist relate to the approximation of the distribution of the number of events by a Poisson distribution. One such condition requires the coefficient of variation (ratio of standard deviation to mean) of the number of events to converge to 0 as the mean increases. The results provide a statistical background to Båth's law in seismology, which asserts that in an aftershock sequence the magnitude of the main shock is commonly around 1.2 units higher than the magnitude of the largest aftershock.

The classical models in risk theory consider a single type of claim. In the insurance business, however, several business lines with separate claim arrival processes appear naturally, and the individual claim processes may not be independent. We introduce a new class of models for such situations, where the underlying counting process is a multivariate continuous-time Markov chain of pure-birth type and the dependency of the components arises from the fact that the birth rate for a specific claim type may depend on the number of claims in the other component processes. Under certain conditions, we obtain a fluid limit, i.e. a functional law of large numbers for these processes. We also investigate the consequences of such results for questions of interest in insurance applications. Several specific subclasses of the general model are discussed in detail and the Cramér asymptotics of the ruin probabilities are derived in particular cases.

We introduce a modification of the linear ARCH (LARCH) model (Giraitis, Robinson, and Surgailis (2000)) - a special case of Sentana's (1995) quadratic ARCH (QARCH) model - for which the conditional variance is a sum of a positive constant and the square of an inhomogeneous linear combination of past observations. Necessary and sufficient conditions for the existence of a stationary solution with finite variance are obtained. We give conditions under which the stationary solution with infinite fourth moment can exhibit long memory, the leverage effect, and a Lévy-stable limit behavior of partial sums of squares.

We consider a branching particle system where an individual particle gives birth to a random number of offspring at the place where it dies. The probability distribution of the number of offspring is given by pk, k = 2, 3, …. The corresponding branching process is related to the semilinear partial differential equation for x ∈ ℝd, where A is the infinitesimal generator of a multiplicative semigroup and the pks, k = 2, 3, …, are nonnegative functions such that We obtain sufficient conditions for the existence of global positive solutions to semilinear equations of this form. Our results extend previous work by Nagasawa and Sirao (1969) and others.

We describe and implement a novel methodology for Monte Carlo simulation of one-dimensional killed diffusions. The proposed estimators represent an unbiased and efficient alternative to current Monte Carlo estimators based on discretization methods for the cases when the finite-dimensional distributions of the process are unknown. For barrier option pricing in finance, we design a suitable Monte Carlo algorithm both for the single barrier case and the double barrier case. Results from numerical investigations are in excellent agreement with the theoretical predictions.

A weighted graph G is a pair (V, ℰ) containing vertex set V and edge set ℰ, where each edge e ∈ ℰ is associated with a weight We. A subgraph of G is a forest if it has no cycles. All forests on the graph G form a probability space, where the probability of each forest is proportional to the product of the weights of its edges. This paper aims to simulate forests exactly from the target distribution. Methods based on coupling from the past (CFTP) and rejection sampling are presented. Comparisons of these methods are given theoretically and via simulation.