We investigate the ray-length distributions for two different rectangular versions of Gilbert's tessellation (see Gilbert (1967)). In the full rectangular version, lines extend either horizontally (east- and west-growing rays) or vertically (north- and south-growing rays) from seed points which form a Poisson point process, each ray stopping when another ray is met. In the half rectangular version, east- and south-growing rays do not interact with west and north rays. For the half rectangular tessellation, we compute analytically, via recursion, a series expansion for the ray-length distribution, whilst, for the full rectangular version, we develop an accurate simulation technique, based in part on the stopping-set theory for Poisson processes (see Zuyev (1999)), to accomplish the same. We demonstrate the remarkable fact that plots of the two distributions appear to be identical when the intensity of seeds in the half model is twice that in the full model. In this paper we explore this coincidence, mindful of the fact that, for one model, our results are from a simulation (with inherent sampling error). We go on to develop further analytic theory for the half-Gilbert model using stopping-set ideas once again, with some novel features. Using our theory, we obtain exact expressions for the first and second moments of the ray length in the half-Gilbert model. For all practical purposes, these results can be applied to the full-Gilbert model—as much better approximations than those provided by Mackisack and Miles (1996).

The almost-sure connectivity of the Euclidean minimal spanning forest MSF(X) on a homogeneous Poisson point process X ⊂ ℝd is an open problem for dimension d>2. We introduce a descending family of graphs (Gn)n≥2 that can be seen as approximations to the MSF in the sense that MSF(X)=∩n=2∞Gn(X). For n=2, one recovers the relative neighborhood graph or, in other words, the β-skeleton with β=2. We show that almost-sure connectivity of Gn(X) holds for all n≥2, all dimensions d≥2, and also point processes X more general than the homogeneous Poisson point process. In particular, we show that almost-sure connectivity holds if certain continuum percolation thresholds are strictly positive or, more generally, if almost surely X does not admit generalized descending chains.

Given two correlated Brownian motions (Xt)t≥ 0 and (Yt)t≥ 0 with constant correlation coefficient, we give the upper and lower estimations of the probability ℙ(max0 ≤s≤tXs≥ a, max0 ≤s≤tYs≥ b) for any a,b,t > 0 through explicit formulae. Our strategy is to establish a new reflection principle for two correlated Brownian motions, which can be viewed as an extension of the reflection principle for one-dimensional Brownian motion. Moreover, we also consider the nonexit probability for linear boundaries, i.e. ℙ (Xt ≤ at+c,Yt ≤ bt+d, 0≤ t≤T) for any constants a, b≥0 and c,d, T > 0.

Motivated by a wide range of applications, we consider a development of Whittle's restless bandit model in which project activation requires a state-dependent amount of a key resource, which is assumed to be available at a constant rate. As many projects may be activated at each decision epoch as resource availability allows. We seek a policy for project activation within resource constraints which minimises an aggregate cost rate for the system. Project indices derived from a Lagrangian relaxation of the original problem exist provided the structural requirement of indexability is met. Verification of this property and derivation of the related indices is greatly simplified when the solution of the Lagrangian relaxation has a state monotone structure for each constituent project. We demonstrate that this is indeed the case for a wide range of bidirectional projects in which the project state tends to move in a different direction when it is activated from that in which it moves when passive. This is natural in many application domains in which activation of a project ameliorates its condition, which otherwise tends to deteriorate or deplete. In some cases the state monotonicity required is related to the structure of state transitions, while in others it is also related to the nature of costs. Two numerical studies demonstrate the value of the ideas for the construction of policies for dynamic resource allocation, most especially in contexts which involve a large number of projects.

The pricing of options in exponential Lévy models amounts to the computation of expectations of functionals of Lévy processes. In many situations, Monte Carlo methods are used. However, the simulation of a Lévy process with infinite Lévy measure generally requires either truncating or replacing the small jumps by a Brownian motion with the same variance. We will derive bounds for the errors generated by these two types of approximation.

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γ.

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.

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.

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.

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.

In this paper we investigate the asymptotic behaviors of the finite- and infinite-time ruin probabilities for a Poisson risk model with stochastic investment returns which constitute a general adapted càdlàg process and heavy-tailed claim sizes which are bivariate upper tail independent. The results of this paper show that the asymptotic ruin probabilities are dominated by the extreme of insurance risk but not by that of investment risk. As applications of the results, we discuss four special cases when the investment returns are determined by a fractional Brownian motion, an integrated Vasicek model, an integrated Cox–Ingersoll–Ross model, and the Heston model.

The Halfin–Whitt regime, or the quality-and-efficiency-driven (QED) regime, for multiserver systems refers to a situation with many servers, a critical load, and yet favorable system performance. We apply this regime to the classical multiserver loss system with slow retrials. We derive nondegenerate limiting expressions for the main steady-state performance measures, including the retrial rate and the blocking probability. It is shown that the economies of scale associated with the QED regime persist for systems with retrials, although in situations when the load becomes extremely critical the retrials cause deteriorated performance. Most of our results are obtained by a detailed analysis of Cohen's equation that defines the retrial rate in an implicit way. The limiting expressions are established by studying prelimit behavior and exploiting the connection between Cohen's equation and Mills' ratio for the Gaussian and Poisson distributions.