In this paper we investigate the number and maximum severity of the ruin excursion of the insurance portfolio reserve process in the Cramér–Lundberg model with and without tax payments. We also provide a relation of the Cramér–Lundberg risk model with the G/G/∞ queue and use it to derive some explicit ruin probability formulae. Finally, the renewal risk model with tax is considered, and an asymptotic identity is derived that in some sense extends the tax identity of the Cramér– Lundberg risk model.

For a bivariate Lévy process (ξt,ηt)t≥ 0 and initial value V0 define the generalised Ornstein–Uhlenbeck (GOU) process Vt:=eξt (V0+∫t0 e-ξs-dηs), t≥0, and the associated stochastic integral process Zt:=∫0t e-ξs-dηs, t≥0. Let Tz:=inf{t>0: Vt<0|V0=z} and ψ(z):=P(Tz<∞) for z≥0 be the ruin time and infinite horizon ruin probability of the GOU process. Our results extend previous work of Nyrhinen (2001) and others to give asymptotic estimates for ψ(z) and the distribution of Tz as z→∞, under very general, easily checkable, assumptions, when ξ satisfies a Cramér condition.

We consider a risk model with threshold strategy, where the insurance company pays off a certain percentage of the income as dividend whenever the current surplus is larger than a given threshold. We investigate the ruin time, ruin probability, and the total dividend, using methods and results from queueing theory.

We consider an insurance model, where the underlying point process is a Cox process. Using a martingale approach applied to diffusion processes, finite-time Lundberg inequalities are obtained. By change-of-measure techniques, Cramér–Lundberg approximations are derived.

For a random walk with both downward and upward jumps (increments), the joint distribution of the exit time across a given level and the undershoot or overshoot at crossing is determined through its generating function, when assuming that the distribution of the jump in the direction making the exit possible has a Laplace transform which is a rational function. The expected exit time is also determined and the paper concludes with exact distribution results concerning exits from bounded intervals. The proofs use simple martingale techniques together with some classical expansions of polynomials and Rouché's theorem from complex function theory.

Suppose that N is a ℤ+-valued random variable and that X,X1,X2,… is a sequence of independent and identically distributed ℤ+ random variables independent of N. In this paper we are interested in properties of the conditional variable In particular, we want to know the asymptotic behavior of the conditional mean ENk or the conditional variance varNk as k→∞. We consider the cases when X is Poisson and when X is mixed Poisson. The problem is motivated by modeling loss reserves in nonlife insurance.

In this paper we extend our previous work. We consider the local-time process L of a strong Markov process X, add negative drift to L, and reflect it à la Skorokhod to obtain a process Q. The reflection of X, together with Q, is, in some sense, a macroscopic model for a service system with two priorities. We derive an expression for the joint law of the duration of an excursion, the maximum value of the process on it, and the time between successive excursions. We work with a properly constructed stationary version of the process. Examples are also given in the paper.

In this paper we consider a general Lévy process X reflected at a downward periodic barrier At and a constant upper barrier K, giving a process VKt=Xt +LAt−LKt. We find the expression for a loss rate defined by lK=ELK1 and identify its asymptotics as K→∞ when X has light-tailed jumps and EX1<0.

We consider a stochastic differential equation (SDE) with piecewise linear drift driven by a spectrally one-sided Lévy process. We show that this SDE has some connections with queueing and storage models, and we use this observation to obtain the invariant distribution.

We use the properties of the Matuszewska indices to show asymptotic inequalities for hazard rates. We discuss the relation between membership in the classes of dominatedly or extended rapidly varying tail distributions and corresponding hazard rate conditions. Convolution closure is established for the class of distributions with extended rapidly varying tails.

We prove large deviation results for Minkowski sums Sn of independent and identically distributed random compact sets where we assume that the summands have a regularly varying distribution and finite expectation. The main focus is on random convex compact sets. The results confirm the heavy-tailed large deviation heuristics: ‘large’ values of the sum are essentially due to the ‘largest’ summand. These results extend those in Mikosch, Pawlas and Samorodnitsky (2011) for generally nonconvex sets, where we assumed that the normalization of Sn grows faster than n.

We study asymptotically optimal simulation algorithms for approximating the tail probability of P(eX1+⋯+ eXd>u) as u→∞. The first algorithm proposed is based on conditional Monte Carlo and assumes that (X1,…,Xd) has an elliptical distribution with very mild assumptions on the radial component. This algorithm is applicable to a large class of models in finance, as we demonstrate with examples. In addition, we propose an importance sampling algorithm for an arbitrary dependence structure that is shown to be asymptotically optimal under mild assumptions on the marginal distributions and, basically, that we can simulate efficiently (X1,…,Xd|Xj >b) for large b. Extensions that allow us to handle portfolios of financial options are also discussed.

A stochastic perpetuity takes the form D∞=∑n=0∞ exp(Y1+⋯+Yn)Bn, where Yn:n≥0) and (Bn:n≥0) are two independent sequences of independent and identically distributed random variables (RVs). This is an expression for the stationary distribution of the Markov chain defined recursively by Dn+1=AnDn+Bn, n≥0, where An=eYn; D∞ then satisfies the stochastic fixed-point equation D∞D̳AD∞+B, where A and B are independent copies of the An and Bn (and independent of D∞ on the right-hand side). In our framework, the quantity Bn, which represents a random reward at time n, is assumed to be positive, unbounded with EBnp <∞ for some p>0, and have a suitably regular continuous positive density. The quantity Yn is assumed to be light tailed and represents a discount rate from time n to n-1. The RV D∞ then represents the net present value, in a stochastic economic environment, of an infinite stream of stochastic rewards. We provide an exact simulation algorithm for generating samples of D∞. Our method is a variation of dominated coupling from the past and it involves constructing a sequence of dominating processes.

The variance minimization (VM) and cross-entropy (CE) methods are two versatile adaptive importance sampling procedures that have been successfully applied to a wide variety of difficult rare-event estimation problems. We compare these two methods via various examples where the optimal VM and CE importance densities can be obtained analytically. We find that in the cases studied both VM and CE methods prescribe the same importance sampling parameters, suggesting that the criterion of minimizing the CE distance is very close, if not asymptotically identical, to minimizing the variance of the associated importance sampling estimator.

We apply the known formulae of the RESTART problem to Markov models of software (and many other) systems, and derive new equations. We show how checkpoints might be included, with their resultant performance under RESTART. The result is a complete procedure for finding the mean, variance, and tail behavior of the job completion time as a function of the failure rate. We also provide a detailed example.

We present an exact simulation algorithm for the stationary distribution of the customer delay D for first-in–first-out (FIFO) M/G/c queues in which ρ=λ/μ<1. We assume that the service time distribution G(x)=P(S≤x),x≥0 (with mean 0<E(S)=1/μ<∞), and its corresponding equilibrium distribution Ge(x)=μ∫0x P(S>y)dy are such that samples of them can be simulated. We further assume that G has a finite second moment. Our method involves the general method of dominated coupling from the past (DCFTP) and we use the single-server M/G/1 queue operating under the processor sharing discipline as an upper bound. Our algorithm yields the stationary distribution of the entire Kiefer–Wolfowitz workload process, the first coordinate of which is D. Extensions of the method to handle simulating generalized Jackson networks in stationarity are also remarked upon.

In this paper we generalize existing results for the steady-state distribution of growth-collapse processes with independent exponential intercollapse times to the case where they have a general distribution on the positive real line having a finite mean. In order to compute the moments of the stationary distribution, no further assumptions are needed. However, in order to compute the stationary distribution, the price that we are required to pay is the restriction of the collapse ratio distribution from a general distribution concentrated on the unit interval to minus-log-phase-type distributions. A random variable has such a distribution if the negative of its natural logarithm has a phase-type distribution. Thus, this family of distributions is dense in the family of all distributions concentrated on the unit interval. The approach is to first study a certain Markov-modulated shot noise process from which the steady-state distribution for the related growth-collapse model can be inferred via level crossing arguments.

Based on a simple model due to Dietz, it is shown that the size of a major epidemic of a vector-borne disease with basic reproduction ratio R0>1 is dominated by the size of a standard SIR (susceptible–infected–removed) epidemic with direct host-to-host transmission of disease and the same R0. Further bounds and numerical illustrations are provided, broadly spanning situations where the size of the epidemic is short of infecting almost all those susceptible. The total size is moderately sensitive to changes in the population parameters that contribute to R0, so that the fluctuating behaviour in ‘annual’ epidemics is not surprising.