Many prominent continuous-time stochastic volatility models exhibit certain functional relationships between price jumps and volatility jumps. We show that stochastic volatility models like the Ornstein–Uhlenbeck and other continuous-time CARMA models as well as continuous-time GARCH and EGARCH models all exhibit such functional relations. We investigate the asymptotic behaviour of certain functionals of price and volatility processes for discrete observations of the price process on a grid, which are relevant for estimation and testing problems.

The risk processes considered in this paper are generated by an underlying Markov process with a regenerative structure and an independent sequence of independent and identically distributed claims. Between the arrivals of claims the process increases at a rate which is a nonnegative function of the present value of the Markov process. The intensity for a claim to occur is another nonnegative function of the value of the Markov process. The claim arrival times are the regeneration times for the Markov process. Two-sided claims are allowed, but the distribution of the positive claims is assumed to have a Laplace transform that is a rational function. The main results describe the joint Laplace transform of the time at ruin and the deficit at ruin. The method used consists in finding partial eigenfunctions for the generator of the joint process consisting of the Markov process and the accumulated claims process, a joint process which is also Markov. These partial eigenfunctions are then used to find a martingale that directly leads to an expression for the desired Laplace transform. In the final section, three examples are given involving different types of the underlying Markov process.

Consider a general bivariate Lévy-driven risk model. The surplus process Y, starting with Y0=x > 0, evolves according to dYt= Yt- dRt -dPt for t > 0, where P and R are two independent Lévy processes respectively representing a loss process in a world without economic factors and a process describing the return on investments in real terms. Motivated by a conjecture of Paulsen, we study the finite-time and infinite-time ruin probabilities for the case in which the loss process P has a Lévy measure of extended regular variation and the stochastic exponential of R fulfills a moment condition. We obtain a simple and unified asymptotic formula as x→∞, which confirms Paulsen's conjecture.

A finite-horizon insurance model is studied where the risk/reserve process can be controlled by reinsurance and investment in the financial market. Our setting is innovative in the sense that we describe in a unified way the timing of the events, that is, the arrivals of claims and the changes of the prices in the financial market, by means of a continuous-time semi-Markov process which appears to be more realistic than, say, classical diffusion-based models. Obtaining explicit optimal solutions for the minimizing ruin probability is a difficult task. Therefore we derive a specific methodology, based on recursive relations for the ruin probability, to obtain a reinsurance and investment policy that minimizes an exponential bound (Lundberg-type bound) on the ruin probability.

We consider the question of an optimal transaction between two investors to minimize their risks. We define a dynamic entropic risk measure using backward stochastic differential equations related to a continuous-time single jump process. The inf-convolution of dynamic entropic risk measures is a key transformation in solving the optimization problem.

We study the asymptotic behaviors of estimators of the parameters in an Ornstein–Uhlenbeck process with linear drift, such as the law of the iterated logarithm (LIL) and Berry–Esseen bounds. As an application of the Berry–Esseen bounds, the precise rates in the LIL for the estimators are obtained.

Based on the explicit coupling property, the ergodicity and the exponential ergodicity of Lévy-driven Ornstein–Uhlenbeck processes are established.

In the setting of the classical Cramér–Lundberg risk insurance model, Albrecher and Hipp (2007) introduced the idea of tax payments. More precisely, if X = {Xt: t≥ 0} represents the Cramér–Lundberg process and, for all t≥ 0, St=sup_{s≤ t}Xs, then Albrecher and Hipp studied Xt - γ St,t≥ 0, where γ∈(0,1) is the rate at which tax is paid. This model has been generalised to the setting that X is a spectrally negative Lévy process by Albrecher, Renaud and Zhou (2008). Finally, Kyprianou and Zhou (2009) extended this model further by allowing the rate at which tax is paid with respect to the process S = {St: t≥ 0} to vary as a function of the current value of S. Specifically, they considered the so-called perturbed spectrally negative Lévy process, Ut:=Xt -∫(0,t]γ(S_u)dSu,t≥ 0, under the assumptions that γ:[0,∞)→ [0,1) and ∫0∞ (1-γ(s))d s =∞. In this article we show that a number of the identities in Kyprianou and Zhou (2009) are still valid for a much more general class of rate functions γ:[0,∞)→∝. Moreover, we show that, with appropriately chosen γ, the perturbed process can pass continuously (i.e. creep) into (-∞, 0) in two different ways.

We extend the class of tempered stable distributions, which were first introduced in Rosiński (2007). Our new class allows for more structure and more variety of the tail behaviors. We discuss various subclasses and the relations between them. To characterize the possible tails, we give detailed results about finiteness of various moments. We also give necessary and sufficient conditions for the tails to be regularly varying. This last part allows us to characterize the domain of attraction to which a particular tempered stable distribution belongs.

For a sequence of finite state space birth–death processes, each having a single absorbing state, we show that, under certain conditions, as the size of the state space tends to infinity, the quasistationary distributions converge to the stationary distribution of a limiting infinite state space birth–death process. This generalizes a result of Keilson and Ramaswamy by allowing birth and death rates to depend upon the size of the state space. We give sufficient conditions under which the convergence result of Keilson and Ramaswamy remains valid. The generalization allows us to apply our convergence result to examples from population biology: a Pearl–Verhulst logistic population growth model and the susceptible-infective-susceptible (SIS) model for infectious spread. The limit distributions obtained suggest new finite-population approximations to the quasistationary distributions of these models, obtained by the method of cumulant closure. The new approximations are found to be both simple in form and accurate.

We consider a single-server queue with Poisson input operating under first-come–first-served (FCFS) or last-come–first-served (LCFS) disciplines. The service times of the customers are independent and obey a general distribution. The system is subject to costs for holding a customer per unit of time, which can be customer specific or customer class specific. We give general expressions for the corresponding value functions, which have elementary compact forms, similar to the Pollaczek–Khinchine mean value formulae. The results generalize earlier work where similar expressions have been obtained for specific service time distributions. The obtained value functions can be readily applied to develop nearly optimal dispatching policies for a broad range of systems with parallel queues, including multiclass scenarios and the cases where service time estimates are available.

We are interested in continuous-time, denumerable state controlled Markov chains (CMCs), with compact Borel action sets, and possibly unbounded transition and reward rates, under the discounted reward optimality criterion. For such CMCs, we propose a definition of a sequence of control models {ℳn} converging to a given control model ℳ, which ensures that the discount optimal reward and policies of ℳn converge to those of ℳ. As an application, we propose a finite-state and finite-action truncation technique of the original control model ℳ, which is illustrated by approximating numerically the optimal reward and policy of a controlled population system with catastrophes. We study the corresponding convergence rates.

In this paper we study the functional central limit theorem (CLT) for stationary Markov chains with a self-adjoint operator and general state space. We investigate the case when the variance of the partial sum is not asymptotically linear in n, and establish that conditional convergence in distribution of partial sums implies the functional CLT. The main tools are maximal inequalities that are further exploited to derive conditions for tightness and convergence to the Brownian motion.

Let {Xn(t), t∈[0,∞)}, n∈ℕ, be standard stationary Gaussian processes. The limit distribution of t∈[0,T(n)]|Xn(t)| is established as rn(t), the correlation function of {Xn(t), t∈[0,∞)}, n∈ℕ, which satisfies the local and long-range strong dependence conditions, extending the results obtained in Seleznjev (1991).

In this paper we explore the first passage times of constant-elasticity-of-variance (CEV) processes with two-sided reflecting barriers. The explicit Laplace transforms of the first passage times are derived. Our results can include analytic formulae concerning Laplace transforms of first passage times of reflected Ornstein–Uhlenbeck processes, reflected geometric Brownian motions, and reflected square-root processes.

Upper deviation results are obtained for the split time of a supercritical continuous-time Markov branching process. More precisely, we establish the existence of logarithmic limits for the likelihood that the split times of the process are greater than an identified value and determine an expression for the limiting quantity. We also give an estimation for the lower deviation probability of the split times, which shows that the scaling is completely different from the upper deviations.

In this paper we introduce a variation on Glaser's method for determining the shape of the failure rate function of a mixture. It has often been seen that the shape of the failure rate depends on the mixing parameter q. Our method provides an explanation for this phenomenon. We then illustrate our technique with the mixture of an exponential and a gamma density for all possible cases.

Consider a random graph where the mean degree is given and fixed. In this paper we derive the maximal size of the largest connected component in the graph. We also study the related question of the largest possible outbreak size of an epidemic occurring ‘on’ the random graph (the graph describing the social structure in the community). More precisely, we look at two different classes of random graphs. First, the Poissonian random graph in which each node i is given an independent and identically distributed (i.i.d.) random weight Xi with E(Xi)=µ, and where there is an edge between i and j with probability 1-e-XiXj/(µ n), independently of other edges. The second model is the thinned configuration model in which the n vertices of the ground graph have i.i.d. ground degrees, distributed as D, with E(D) = µ. The graph of interest is obtained by deleting edges independently with probability 1-p. In both models the fraction of vertices in the largest connected component converges in probability to a constant 1-q, where q depends on X or D and p. We investigate for which distributions X and D with given µ and p, 1-q is maximized. We show that in the class of Poissonian random graphs, X should have all its mass at 0 and one other real, which can be explicitly determined. For the thinned configuration model, D should have all its mass at 0 and two subsequent positive integers.

Markovian systems with multiple interacting subsystems under the influence of a control unit are considered. The state spaces of the subsystems are countably infinite, whereas that of the control unit is finite. A recent infinite level-dependent quasi-birth-and-death model for such systems is extended by facilitating the automatic representation and generation of the nonzero blocks in its underlying infinitesimal generator matrix with sums of Kronecker products. Experiments are performed on systems of stochastic chemical kinetics having two or more countably infinite state space subsystems. Results indicate that, even though more memory is consumed, there are many cases where a matrix-analytic solution coupled with Lyapunov theory yields a faster and more accurate steady-state measure compared to that obtained with simulation.

The Asmussen–Kroese Monte Carlo estimators of P(Sn > u) and P(SN > u) are known to work well in rare event settings, where SN is the sum of independent, identically distributed heavy-tailed random variables X1,…,XN and N is a nonnegative, integer-valued random variable independent of the Xi. In this paper we show how to improve the Asmussen–Kroese estimators of both probabilities when the Xi are nonnegative. We also apply our ideas to estimate the quantity E[(SN-u)+].