In this paper we consider two classes of reflected Ornstein–Uhlenbeck (OU) processes: the reflected OU process with jumps and the Markov-modulated reflected OU process. We prove that their stationary distributions exist. Furthermore, for the jump reflected OU process, the Laplace transform (LT) of the stationary distribution is given. As for the Markov-modulated reflected OU process, we derive an equation satisfied by the LT of the stationary distribution.

We consider a modified version of the classical optimal dividends problem of de Finetti in which the objective function is altered by adding in an extra term which takes account of the ruin time of the risk process, the latter being modeled by a spectrally negative Lévy process. We show that, with the exception of a small class, a barrier strategy forms an optimal strategy under the condition that the Lévy measure has a completely monotone density. As a prerequisite for the proof, we show that, under the aforementioned condition on the Lévy measure, the q-scale function of the spectrally negative Lévy process has a derivative which is strictly log-convex.

In this paper we study m-discount optimality (m ≥ −1) and Blackwell optimality for a general class of controlled (Markov) diffusion processes. To this end, a key step is to express the expected discounted reward function as a Laurent series, and then search certain control policies that lexicographically maximize the mth coefficient of this series for m = −1,0,1,…. This approach naturally leads to m-discount optimality and it gives Blackwell optimality in the limit as m → ∞.

The probability that two randomly selected phylogenetic trees of the same size are isomorphic is found to be asymptotic to a decreasing exponential modulated by a polynomial factor. The number of symmetrical nodes in a random phylogenetic tree of large size obeys a limiting Gaussian distribution, in the sense of both central and local limits. The probability that two random phylogenetic trees have the same number of symmetries asymptotically obeys an inverse square-root law. Precise estimates for these problems are obtained by methods of analytic combinatorics, involving bivariate generating functions, singularity analysis, and quasi-powers approximations.

We derive a moderate deviation principle for word counts (which is extended to counts of multiple patterns) in biological sequences under different models: independent and identically distributed letters, homogeneous Markov chains of order 1 and m, and, in view of the codon structure of DNA sequences, Markov chains with three different transition matrices. This enables us to approximate P-values for the number of word occurrences in DNA and protein sequences in a new manner.

We consider a discrete-time risk process driven by proportional reinsurance and an interest rate process. We assume that the interest rate process behaves as a Markov chain. To reduce the risk of ruin, we may reinsure a part or even all of the reserve. Recursive and integral equations for ruin probabilities are given. Generalized Lundberg inequalities for the ruin probabilities are derived given a stationary policy. To illustrate these results, a numerical example is included.

In this paper, a stochastic integral of Ornstein–Uhlenbeck type is represented to be the sum of two independent random variables: one has a tempered stable distribution and the other has a compound Poisson distribution. In distribution, the compound Poisson random variable is equal to the sum of a Poisson-distributed number of positive random variables, which are independent and identically distributed and have a common specified density function. Based on the representation of the stochastic integral, we prove that the transition distribution of the tempered stable Ornstein–Uhlenbeck process is self-decomposable and that the transition density is a C∞-function.

In this paper we introduce the transformed two-parameter Poisson–Dirichlet distribution on the ordered infinite simplex. Furthermore, we prove the central limit theorem related to this distribution when both the mutation rate θ and the selection rate σ become large in a specified manner. As a consequence, we find that the properly scaled homozygosities have asymptotical normal behavior. In particular, there is a certain phase transition with the limit depending on the relative strength of σ and θ.

Suppose that an unknown number of objects arrive sequentially according to a Poisson process with random intensity λ on some fixed time interval [0,T]. We assume a gamma prior density Gλ(r, 1/a) for λ. Furthermore, we suppose that all arriving objects can be ranked uniquely among all preceding arrivals. Exactly one object can be selected. Our aim is to find a stopping time (selection time) which maximizes the time during which the selected object will stay relatively best. Our main result is the following. It is optimal to select the ith object that is relatively best and arrives at some time si(r) onwards. The value of si(r) can be obtained for each r and i as the unique root of a deterministic equation.

Let X be a pure-jump subordinator (i.e. nondecreasing Lévy process with no drift) with infinite Lévy measure, let Xε be the sum of jumps not exceeding ε, and let µ(ε)=E[Xε(1)]. We study the question of weak convergence of Xε/µ(ε) as ε ↓0, in terms of the limit behavior of µ(ε)/ε. The most interesting case reduces to the weak convergence of Xε/ε to a subordinator whose marginals are generalized Dickman distributions; we give some necessary and sufficient conditions for this to hold. For a certain significant class of subordinators for which the latter convergence holds, and whose most prominent representative is the gamma process, we give some detailed analysis regarding the convergence quality (in particular, in the context of approximating X itself). This paper completes, in some respects, the study made by Asmussen and Rosiński (2001).

In this paper we consider an extension of the Sparre Andersen insurance risk model by relaxing one of its independence assumptions. The newly proposed dependence structure is introduced through the premise that the joint distribution of the interclaim time and the subsequent claim size is bivariate phase-type (see, e.g. Assaf et al. (1984) and Kulkarni (1989)). Relying on the existing connection between risk processes and fluid flows (see, e.g. Badescu et al. (2005), Badescu, Drekic and Landriault (2007), Ramaswami (2006), and Ahn, Badescu and Ramaswami (2007)), we construct an analytically tractable fluid flow that leads to the analysis of various ruin-related quantities in the aforementioned risk model. Using matrix-analytic methods, we obtain an explicit expression for the Gerber–Shiu discounted penalty function (see Gerber and Shiu (1998)) when the penalty function depends on the deficit at ruin only. Finally, we investigate how some ruin-related quantities involving the surplus immediately prior to ruin can also be analyzed via our fluid flow methodology.

Let (Xi)i∈ℕ be a sequence of independent and identically distributed random variables with values in the set ℕ0 of nonnegative integers. Motivated by applications in enumerative combinatorics and analysis of algorithms we investigate the number of gaps and the length of the longest gap in the set {X1,…,Xn} of the first n values. We obtain necessary and sufficient conditions in terms of the tail sequence (qk)k∈ℕ0, qk=P(X1≥ k), for the gaps to vanish asymptotically as n→∞: these are ∑k=0∞qk+1/qk <∞ and limk→∞qk+1/qk=0 for convergence almost surely and convergence in probability, respectively. We further show that the length of the longest gap tends to ∞ in probability if qk+1/qk→ 1. For the family of geometric distributions, which can be regarded as the borderline case between the light-tailed and the heavy-tailed situations and which is also of particular interest in applications, we study the distribution of the length of the longest gap, using a construction based on the Sukhatme–Rényi representation of exponential order statistics to resolve the asymptotic distributional periodicities.

Convergence in probability and central limit laws of bipower variation for Gaussian processes with stationary increments and for integrals with respect to such processes are derived. The main tools of the proofs are some recent powerful techniques of Wiener/Itô/Malliavin calculus for establishing limit laws, due to Nualart, Peccati, and others.

In this paper we consider the following problem. An angler buys a fishing ticket that allows him/her to fish for a fixed time. There are two locations to fish at the lake. The fish are caught according to a renewal process, which is different for each fishing location. The angler's success is defined as the difference between the utility function, which is dependent on the size of the fish caught, and the time-dependent cost function. These functions are different for each fishing location. The goal of the angler is to find two optimal stopping times that maximize his/her success: when to change fishing location and when to stop fishing. Dynamic programming methods are used to find these two optimal stopping times and to specify the expected success of the angler at these times.

Convolutions of long-tailed and subexponential distributions play a major role in the analysis of many stochastic systems. We study these convolutions, proving some important new results through a simple and coherent approach, and also showing that the standard properties of such convolutions follow as easy consequences.

A model for describing the lifetimes of coherent systems, where the failures of components may have an impact on the lifetimes of the remaining components, is proposed. The model is motivated by the definition of sequential order statistics (cf. Kamps (1995)). Sequential order statistics describe the successive failure times in a sequential k-out-of-n system, where the distribution of the remaining components' lifetimes is allowed to change after every failure of a component. In the present paper, general component lifetimes which can be influenced by failures are considered. The ordered failure times of these components can be used to extend the concept of sequential order statistics. In particular, a definition of sequential order statistics based on exchangeable components is proposed. By utilizing the system signature (cf. Samaniego (2007)), the distribution of the lifetime of a coherent system with failure-dependent exchangeable component lifetimes is shown to be given by a mixture of the distributions of sequential order statistics. Furthermore, some results on the joint distribution of sequential order statistics based on exchangeable components are given.

We approximate the distribution of the sum of independent but not necessarily identically distributed Bernoulli random variables using a shifted binomial distribution, where the three parameters (the number of trials, the probability of success, and the shift amount) are chosen to match the first three moments of the two distributions. We give a bound on the approximation error in terms of the total variation metric using Stein's method. A numerical study is discussed that shows shifted binomial approximations are typically more accurate than Poisson or standard binomial approximations. The application of the approximation to solving a problem arising in Bayesian hierarchical modeling is also discussed.

We investigate the almost sure asymptotic properties of vector martingale transforms. Assuming some appropriate regularity conditions both on the increasing process and on the moments of the martingale, we prove that normalized moments of any even order converge in the almost sure central limit theorem for martingales. A conjecture about almost sure upper bounds under wider hypotheses is formulated. The theoretical results are supported by examples borrowed from statistical applications, including linear autoregressive models and branching processes with immigration, for which new asymptotic properties are established on estimation and prediction errors.

In this paper we consider the splitting method first introduced in rare event analysis. In this technique, the sample paths are split into R multiple copies at various stages to speed up the simulation. Given the cost, the optimization of the algorithm suggests taking all the transition probabilities to be equal; nevertheless, in practice, these quantities are unknown. We address this problem by presenting an algorithm in two phases.

In this paper we study efficient simulation algorithms for estimating P(X›x), where X is the total time of a job with ideal time T that needs to be restarted after a failure. The main tool is importance sampling, where a good importance distribution is identified via an asymptotic description of the conditional distribution of T given X›x. If T≡t is constant, the problem reduces to the efficient simulation of geometric sums, and a standard algorithm involving a Cramér-type root, γ(t), is available. However, we also discuss an algorithm that avoids finding the root. If T is random, particular attention is given to T having either a gamma-like tail or a regularly varying tail, and to failures at Poisson times. Different types of conditional limit occur, in particular exponentially tilted Gumbel distributions and Pareto distributions. The algorithms based upon importance distributions for T using these asymptotic descriptions have bounded relative error as x→∞ when combined with the ideas used for a fixed t. Nevertheless, we give examples of algorithms carefully designed to enjoy bounded relative error that may provide little or no asymptotic improvement over crude Monte Carlo simulation when the computational effort is taken into account. To resolve this problem, an alternative algorithm using two-sided Lundberg bounds is suggested.