We consider a stochastic evolutionary model for a phenotype developing amongst n related species with unknown phylogeny. The unknown tree is modelled by a Yule process conditioned on n contemporary nodes. The trait value is assumed to evolve along lineages as an Ornstein-Uhlenbeck process. As a result, the trait values of the n species form a sample with dependent observations. We establish three limit theorems for the sample mean corresponding to three domains for the adaptation rate. In the case of fast adaptation, we show that for large n the normalized sample mean is approximately normally distributed. Using these limit theorems, we develop novel confidence interval formulae for the optimal trait value.

In this paper we use the Siegert formula to derive alternative expressions for the moments of the first passage time of the Ornstein-Uhlenbeck process through a constant threshold. The expression for the nth moment is recursively linked to the lower-order moments and consists of only n terms. These compact expressions can substantially facilitate (numerical) applications also for higher-order moments.

We consider a multivariate continuous-time ARMA (MCARMA) process sampled at a high-frequency time grid {hn, 2hn,…, nhn}, where hn ↓ 0 and nhn → ∞ as n → ∞, or at a constant time grid where hn = h. For this model, we present the asymptotic behavior of the properly normalized partial sum to a multivariate stable or a multivariate normal random vector depending on the domain of attraction of the driving Lévy process. Furthermore, we derive the asymptotic behavior of the sample variance. In the case of finite second moments of the driving Lévy process the sample variance is a consistent estimator. Moreover, we embed the MCARMA process in a cointegrated model. For this model, we propose a parameter estimator and derive its asymptotic behavior. The results are given for more general processes than MCARMA processes and contain some asymptotic properties of stochastic integrals.

Consider a one-dimensional diffusion process on the diffusion interval I originated in x0 ∈ I. Let a(t) and b(t) be two continuous functions of t, t > t0, with bounded derivatives, a(t) < b(t), and a(t), b(t) ∈ I, for all t > t0. We study the joint distribution of the two random variables Ta and Tb, the first hitting times of the diffusion process through the two boundaries a(t) and b(t), respectively. We express the joint distribution of Ta and Tb in terms of ℙ(Ta < t, Ta < Tb) and ℙ(Tb < t, Ta > Tb), and we determine a system of integral equations verified by these last probabilities. We propose a numerical algorithm to solve this system and we prove its convergence properties. Examples and modeling motivation for this study are also discussed.

We show that simple diffusion processes are weak limits of piecewise continuous processesconstructed within a totally deterministic framework. The proofs are based on thecontinuous mapping theorem and the functional central limit theorem.

Spot prices in energy markets exhibit special features, such as price spikes, mean reversion, stochastic volatility, inverse leverage effect, and dependencies between the commodities. In this paper a multivariate stochastic volatility model is introduced which captures these features. The second-order structure and stationarity of the model are analyzed in detail. A simulation method for Monte Carlo generation of price paths is introduced and a numerical example is presented.

In this paper we study first passage times of (reflected) Ornstein-Uhlenbeck processes over compound Poisson-type boundaries. In fact, we extend the results of first rendezvous times of (reflected) Brownian motion and compound Poisson-type processes in Perry, Stadje and Zacks (2004) to the (reflected) Ornstein-Uhlenbeck case.

We prove that the first passage time density ρ(t) for an Ornstein-Uhlenbeck process X(t) obeying dX = -β Xdt + σdW to reach a fixed threshold θ from a suprathreshold initial condition x0 > θ > 0 has a lower bound of the form ρ(t) > kexp[-pe6βt] for positive constants k and p for times t exceeding some positive value u. We obtain explicit expressions for k, p, and u in terms of β, σ, x0, and θ, and discuss the application of the results to the synchronization of periodically forced stochastic leaky integrate-and-fire model neurons.

Continuous-time autoregressive moving average (CARMA) processes with a nonnegative kernel and driven by a nondecreasing Lévy process constitute a very general class of stationary, nonnegative continuous-time processes. In financial econometrics a stationary Ornstein-Uhlenbeck (or CAR(1)) process, driven by a nondecreasing Lévy process, was introduced by Barndorff-Nielsen and Shephard (2001) as a model for stochastic volatility to allow for a wide variety of possible marginal distributions and the possibility of jumps. For such processes, we take advantage of the nonnegativity of the increments of the driving Lévy process to study the properties of a highly efficient estimation procedure for the parameters when observations are available of the CAR(1) process at uniformly spaced times 0,h,…,Nh. We also show how to reconstruct the background driving Lévy process from a continuously observed realization of the process and use this result to estimate the increments of the Lévy process itself when h is small. Asymptotic properties of the coefficient estimator are derived and the results illustrated using a simulated gamma-driven Ornstein-Uhlenbeck process.

We exhibit a scaling law for the critical SIS stochastic epidemic. If at time 0 the population consists of infected and susceptible individuals, then when the time and the number currently infected are both scaled by , the resulting process converges, as N → ∞, to a diffusion process related to the Feller diffusion by a change of drift. As a consequence, the rescaled size of the epidemic has a limit law that coincides with that of a first passage time for the standard Ornstein-Uhlenbeck process. These results are the analogs for the SIS epidemic of results of Martin-Löf (1998) and Aldous (1997) for the simple SIR epidemic.

A cut-off phenomenon is shown to occur in a sample of n independent, identically distributed Ornstein-Uhlenbeck processes and its average. Their distributions stay far from equilibrium before a certain O(log(n)) time, and converge exponentially fast after. Precise estimates show that the total variation distance drops from almost 1 to almost 0 over an interval of time of length O(1) around log(n)/(2α), where α is the viscosity coefficient of the sampled process. The distribution of the hitting time of 0 by the average of the sample is computed. As n tends to infinity, the hitting time becomes concentrated around the cut-off instant, and its tails match the estimates given for the total variation distance.

We generalize a theorem due to Keilson on the approximate exponentiality of waiting times for rare events in regenerative processes. We use the result to investigate the limit distribution for a family of first entrance times in a sequence of Ehrenfest urn models. As a second application, we consider approximate pattern matching, a problem arising in molecular biology and other areas.

Let X be a one-dimensional strong Markov process with continuous sample paths. Using Volterra-Stieltjes integral equation techniques we investigate Hölder continuity and differentiability of first passage time distributions of X with respect to continuous lower and upper moving boundaries. Under mild assumptions on the transition function of X we prove the existence of a continuous first passage time density to one-sided differentiable moving boundaries and derive a new integral equation for this density. We apply our results to Brownian motion and its nonrandom Markovian transforms, in particular to the Ornstein-Uhlenbeck process.

Two methods for approximating the limiting distribution of the present value of the benefits of a portfolio of identical endowment insurance contracts are suggested. The model used assumes that both future lifetimes and interest rates are random. The first method is similar to the one presented in Parker (1994b). The second method is based on the relationship between temporary and endowment insurance contracts.

A highly symmetric loss network is considered, the symmetric star network previously considered by Whitt and Ziedins and Kelly. As the number of links, K, becomes large, the state space for this process also grows, so we consider a functional of the network, one which contains all information relevant to blocking probabilities within the network but which is easier to analyse. We show that this reduced process obeys a functional law of large numbers and a functional central limit theorem, the limit in this latter case being an Ornstein-Uhlenbeck diffusion process. Finally, by considering the network in equilibrium, we are able to prove that the well-known Erlang fixed point approximation for blocking probabilities is correct to within o(K-1/2) as K → ∞.

Two approaches used to model interest randomness are presented. They are the modeling of the force of interest accumulation function and the modeling of the force of interest. The expected value, standard deviation and coefficient of skewness of the present value of annuities-immediate are presented as illustrations. The implicit behavior of the force of interest under the two approaches is investigated by looking at a particular conditional expectation of the force of interest accumulation function.

An approximation of the distribution of the present value of the benefits of a portfolio of temporary insurance contracts is suggested for the case where the size of the portfolio tends to infinity. The model used is the one presented in Parker (1922b) and involves random interest rates and future lifetimes. Some justifications of the approximation are given. Illustrations for limiting portfolios of temporary insurance contracts are presented for an assumed Ornstein-Uhlenbeck process for the force of interest.

Let v(x, t) denote the displacement of an infinitely long, idealized string performing damped vibrations caused by white noise.

Upper and lower bounds for the distribution of maxsv(x, s) and maxxv(x, t) are presented. The results are obtained by adapting Lévy-type inequalities and exploiting a connection of v(x, t) with the Ornstein-Uhlenbeck process through Slepian's theorem.

The case of forced-damped vibrations is also analysed. Finally, a section is devoted to the case of a semi-infinite string performing damped vibrations.

A revealing alternate proof is provided for the Iglehart (1965), (1973)–Borovkov (1967) heavy-traffic limit theorem for GI/G/s queues. This kind of heavy traffic is obtained by considering a sequence of GI/G/s systems with the numbers of servers and the arrival rates going to ∞ while the service-time distributions are held fixed. The theorem establishes convergence to a Gaussian process, which in general is not Markov, for an appropriate normalization of the sequence of stochastic processes representing the number of customers in service at arbitrary times. The key idea in the new proof is to consider service-time distributions that are randomly stopped sums of exponential phases, and then work with the discrete-time vector-valued Markov chain representing the number of customers in each phase of service at arrival epochs. It is then easy to show that this sequence of Markov chains converges to a multivariate O–U (Ornstein–Uhlenbeck) diffusion process by applying simple criteria in Stroock and Varadhan (1979). The Iglehart–Borovkov limit for these special service-time distributions is the sum of the components of this multivariate O–U process. Heavy-traffic convergence is also established for the steady-state distributions of GI/M/s queues under the same conditions by exploiting stochastic-order properties.

We propose a model of a passive nerve cylinder undergoing random stimulus along its length. It is shown that this model is approximated by the solution of a stochastic partial differential equation. Numerous properties of the sample paths are derived, such as their modulus of continuity, quadratic and quartic variation, and it is shown that the solution exhibits the phenomenon of flicker noise. The first-passage problem is studied, and it is shown to be connected with a first-hitting time for an infinite-dimensional diffusion.