Let ${\mathrm{d}} X(t) = -Y(t) \, {\mathrm{d}} t$, where Y(t) is a one-dimensional diffusion process, and let $\tau(x,y)$ be the first time the process (X(t), Y(t)), starting from (x, y), leaves a subset of the first quadrant. The problem of computing the probability $p(x,y)\,:\!=\, \mathbb{P}[X(\tau(x,y))=0]$ is considered. The Laplace transform of the function p(x, y) is obtained in important particular cases, and it is shown that the transform can at least be inverted numerically. Explicit expressions for the Laplace transform of $\mathbb{E}[\tau(x,y)]$ and of the moment-generating function of $\tau(x,y)$ can also be derived.

A diffusion approximation to a risk process under dynamic proportional reinsurance is considered. The goal is to minimise the discounted time in drawdown; that is, the time where the distance of the present surplus to the running maximum is larger than a given level $d > 0$. We calculate the value function and determine the optimal reinsurance strategy. We conclude that the drawdown measure stabilises process paths but has a drawback as it also prevents surpassing the initial maximum. That is, the insurer is, under the optimal strategy, not interested in any more profits. We therefore suggest using optimisation criteria that do not avoid future profits.

We study the transient and limiting behavior of a queue with a Pólya arrival process. The Pólya process is interesting because it exhibits path-dependent behavior, e.g. it satisfies a non-ergodic law of large numbers: the average number of arrivals over time [0, t] converges almost surely to a nondegenerate limit as $t \rightarrow \infty$. We establish a heavy-traffic diffusion limit for the $\sum_{i=1}^{n} P_i/GI/1$ queue, with arrivals occurring exogenously according to the superposition of n independent and identically distributed Pólya point processes. That limit yields a tractable approximation for the transient queue-length distribution, because the limiting net input process is a Gaussian Markov process with stationary increments. We also provide insight into the long-run performance of queues with path-dependent arrival processes. We show how Little’s law can be stated in this context, and we provide conditions under which there is stability for a queue with a Pólya arrival process.

Our aim is to find sufficient conditions for weak convergence of stochastic integrals with respect to the state occupation measure of a Markov chain. First, we study properties of the state indicator function and the state occupation measure of a Markov chain. In particular, we establish weak convergence of the state occupation measure under a scaling of the generator matrix. Then, relying on the connection between the state occupation measure and the Dynkin martingale, we provide sufficient conditions for weak convergence of stochastic integrals with respect to the state occupation measure. We apply our results to derive diffusion limits for the Markov-modulated Erlang loss model and the regime-switching Cox–Ingersoll–Ross process.

We provide the first generic exact simulation algorithm for multivariate diffusions. Current exact sampling algorithms for diffusions require the existence of a transformation which can be used to reduce the sampling problem to the case of a constant diffusion matrix and a drift which is the gradient of some function. Such a transformation, called the Lamperti transformation, can be applied in general only in one dimension. So, completely different ideas are required for the exact sampling of generic multivariate diffusions. The development of these ideas is the main contribution of this paper. Our strategy combines techniques borrowed from the theory of rough paths, on the one hand, and multilevel Monte Carlo on the other.

This work proposes and analyzes a family of spatially inhomogeneous epidemic models. This is our first effort to use stochastic partial differential equations (SPDEs) to model epidemic dynamics with spatial variations and environmental noise. After setting up the problem, the existence and uniqueness of solutions of the underlying SPDEs are examined. Then, definitions of permanence and extinction are given, and certain sufficient conditions are provided for permanence and extinction. Our hope is that this paper will open up windows for investigation of epidemic models from a new angle.

We study the problem of stopping a Brownian bridge X in order to maximise the expected value of an exponential gain function. The problem was posed by Ernst and Shepp (2015), and was motivated by bond selling with non-negative prices.

Due to the non-linear structure of the exponential gain, we cannot rely on methods used in the literature to find closed-form solutions to other problems involving the Brownian bridge. Instead, we must deal directly with a stopping problem for a time-inhomogeneous diffusion. We develop techniques based on pathwise properties of the Brownian bridge and martingale methods of optimal stopping theory, which allow us to find the optimal stopping rule and to show the regularity of the value function.

The d-dimensional Λ-Fleming-Viot generator acting on functions g(x ), with x being a vector of d allele frequencies, can be written as a Wright-Fisher generator acting on functions g with a modified random linear argument of x induced by partitioning occurring in the Λ-Fleming-Viot process. The eigenvalues and right polynomial eigenvectors are easy to see from this representation. The two-dimensional process, which has a one-dimensional generator, is considered in detail. A nonlinear equation is found for the Green's function. In a model with genic selection a proof is given that there is a critical selection value such that if the selection coefficient is greater than or equal to the critical value then fixation, when the boundary 1 is hit, has probability 1 beginning from any nonzero frequency. This is an analytic proof different from the proofs of Der, Epstein and Plotkin (2011) and Foucart (2013). An application in the infinitely-many-alleles Λ-Fleming-Viot process is finding an interesting identity for the frequency spectrum of alleles that is based on size biasing. The moment dual process in the Fleming-Viot process is the usual Λ-coalescent tree back in time. The Wright-Fisher representation using a different set of polynomials g n (x) as test functions produces a dual death process which has a similarity to the Kingman coalescent and decreases by units of one. The eigenvalues of the process are analogous to the Jacobi polynomials when expressed in terms of g n (x), playing the role of x n . Under the stationary distribution when there is mutation, is analogous to the nth moment in a beta distribution. There is a d-dimensional version g n (X ), and even an intriguing Ewens' sampling formula analogy when d → ∞.

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.

The drawdown process of a one-dimensional regular diffusion process X is given by X reflected at its running maximum. The drawup process is given by X reflected at its running minimum. We calculate the probability that a drawdown precedes a drawup in an exponential time-horizon. We then study the law of the occupation times of the drawdown process and the drawup process. These results are applied to address problems in risk analysis and for option pricing of the drawdown process. Finally, we present examples of Brownian motion with drift and three-dimensional Bessel processes, where we prove an identity in law.