We consider the stochastic volatility model obtained by adding a compound Hawkes process to the volatility of the well-known Heston model. A Hawkes process is a self-exciting counting process with many applications in mathematical finance, insurance, epidemiology, seismology, and other fields. We prove a general result on the existence of a family of equivalent (local) martingale measures. We apply this result to a particular example where the sizes of the jumps are exponentially distributed. Finally, a practical application to efficient computation of exposures is discussed.

We construct a simple example, surely known to Harry Kesten, of an R-transient Markov chain on a countable state space S ∪ {δ}, where δ is absorbing. The transition matrix K on S is irreducible and strictly substochastic. We determine the Yaglom limit, that is, the limiting conditional behavior given nonabsorption. Each starting state x ∈ S results in a different Yaglom limit. Each Yaglom limit is an R-1-invariant quasi-stationary distribution, where R is the convergence parameter of K. Yaglom limits that depend on the starting state are related to a nontrivial R-1-Martin boundary.

We study how to sample paths of a random walk up to the first time it crosses a fixed barrier, in the setting where the step sizes are independent and identically distributed with negative mean and have a regularly varying right tail. We introduce a desirable property for a change of measure to be suitable for exact simulation. We study whether the change of measure of Blanchet and Glynn (2008) satisfies this property and show that it does so if and only if the tail index α of the right tail lies in the interval (1, 3/2).

A critical branching process {Zk, k = 0, 1, 2, ...} in a random environment is considered. A conditional functional limit theorem for the properly scaled process {log Zpu, 0 ≤ u < ∞} is established under the assumptions that Zn > 0 and p ≪ n. It is shown that the limiting process is a Lévy process conditioned to stay nonnegative. The proof of this result is based on a limit theorem describing the distribution of the initial part of the trajectories of a driftless random walk conditioned to stay nonnegative.

The (conditional or unconditional) distribution of the continuous scan statistic in a one-dimensional Poisson process may be approximated by that of a discrete analogue via time discretization (to be referred to as the discrete approximation). Using a change of measure argument, we derive the first-order term of the discrete approximation which involves some functionals of the Poisson process. Richardson's extrapolation is then applied to yield a corrected (second-order) approximation. Numerical results are presented to compare various approximations.

In this paper we examine the extremal tail probabilities of moving sums in a marked Poisson random field. These sums are computed by adding up the weighted occurrences of events lying within a scanning set of fixed shape and size. We also provide an alternative representation of the constants of the asymptotic formulae in terms of the occupation measure of the conditional local random field at zero, and extend these representations to the constants of asymptotic tail probabilities of Gaussian random fields.

We study the expected time to ruin in a risk process in which dividends are paid when the surplus is above the barrier. We consider the case in which the dividend rate is smaller than the premium rate. We obtain results for the classical compound Poisson risk process with phase-type claim size. When the ruin probability is 1, we derive the expected time to ruin and the expected dividends paid. When the ruin probability is less than 1, these quantities are derived conditioning on the event that ruin occurs.

We consider a failure-prone system operating in continuous time. Condition monitoring is conducted at discrete time epochs. The state of the system is assumed to evolve as a continuous-time Markov process with a finite state space. The observation process with continuous-range values is stochastically related to the state process, which, except for the failure state, is unobservable. Combining the failure information and the condition monitoring information, we derive a general recursive filter, and, as special cases, we obtain recursive formulae for the state estimation and other quantities of interest. Updated parameter estimates are obtained using the expectation-maximization (EM) algorithm. Some practical prediction problems are discussed and finally an illustrative example is given using a real dataset.

We consider a failure-prone system which operates in continuous time and is subject to condition monitoring at discrete time epochs. It is assumed that the state of the system evolves as a continuous-time Markov process with a finite state space. The observation process is stochastically related to the state process which is unobservable, except for the failure state. Combining the failure information and the information obtained from condition monitoring, and using the change of measure approach, we derive a general recursive filter, and, as special cases, we obtain recursive formulae for the state estimation and other quantities of interest. Up-dated parameter estimates are obtained using the EM algorithm. Some practical prediction problems are discussed and an illustrative example is given using a real dataset.

We discuss a simple new approach to calculating expectations of a specific form used for the pricing of derivative assets in financial mathematics. We show that in the ‘vanilla case’, the expectations can be found by simply integrating the respective moment generating function with a certain weight. In situations corresponding to barrier-type options, we just need to carry out one more integration. The suggested approach appears to be the first (and, apart from Monte Carlo simulation, the only) one to allow the pricing of discretely monitored exotic options when the underlying asset is modelled by a general Lévy process. We illustrate the method numerically by calculating the price of a discretely monitored lookback call option in the cases when the underlying follows the geometric Brownian and variance-gamma processes.

Let Xt be a Feller (branching) diffusion with drift αx. We consider new processes, the probability measures of which are obtained from that of X via changes of measure involving suitably normalized exponential functions of with λ > 0. The new processes can be thought of as ‘self-reinforcing’ versions of the old.

Depending on the values of α, T and λ, the process under the new measure is shown to exhibit explosion in finite time. We also obtain a number of other results related to the new processes.

Since the Feller diffusion is also the total mass process of a superprocess, we relate the finite-time explosion property to the behaviour of superprocesses with local self-interaction, and raise some interesting questions for these.

Consider a classical compound Poisson model. The safety loading can be positive, negative or zero. Explicit expressions for the distributions of the surplus prior and at ruin are given in terms of the ruin probability. Moreover, the asymptotic behaviour of these distributions as the initial capital tends to infinity are obtained. In particular, for positive safety loading the Cramer case, the case of subexponential distributions and some intermediate cases are discussed.

Numerical evaluation of compound distributions is one of the central numerical tasks in insurance mathematics. Two widely used techniques are Panjer recursion and transform methods. Many authors have pointed out that aliasing errors imply the need to consider the whole distribution if transform methods are used, a potential drawback especially for heavy-tailed distributions. We investigate the magnitude of aliasing errors and show that this problem can be solved by a suitable change of measure.

This paper examines a problem of importance to the telecommunications industry. In the design of modern ATM switches, it is necessary to use simulation to estimate the probability that a queue within the switch exceeds a given large value. Since these are extremely small probabilities, importance sampling methods must be used. Here we obtain a change of measure for a broad class of models with direct applicability to ATM switches.

We consider a model with A independent sources of cells where each source is modeled by a Markov renewal point process with batch arrivals. We do not assume the sources are necessarily identically distributed, nor that batch sizes are independent of the state of the Markov process. These arrivals join a queue served by multiple independent servers, each with service times also modeled as a Markov renewal process. We only discuss a time-slotted system. The queue is viewed as the additive component of a Markov additive chain subject to the constraint that the additive component remains non-negative. We apply the theory in McDonald (1999) to obtain the asymptotics of the tail of the distribution of the queue size in steady state plus the asymptotics of the mean time between large deviations of the queue size.

A Cox risk process with a piecewise constant intensity is considered where the sequence (Li) of successive levels of the intensity forms a Markov chain. The duration σi of the level Li is assumed to be only dependent via Li. In the small-claim case a Lundberg inequality is obtained via a martingale approach. It is shown furthermore by a Lundberg bound from below that the resulting adjustment coefficient gives the best possible exponential bound for the ruin probability. In the case where the stationary distribution of Li contains a discrete component, a Cramér–Lundberg approximation can be obtained. By way of example we consider the independent jump intensity model (Björk and Grandell 1988) and the risk model in a Markovian environment (Asmussen 1989).

The theory of finite security markets developed by Harrison and Pliska [1] used the separating hyperplane theorem to establish the relationship between the lack of arbitrage opportunities and the existence of a certain martingale measure. In this paper we treat this theory by examining certain geometric properties of the sample paths of the price process, that is, we focus on the price increments of the stocks between one time period to the next and convert them to martingale differences through an equivalent change of measure. Thus, in contrast to Harrison and Pliska&s functional analytic derivation, our approach is based on probabilistic methods and allows a geometric interpretation which not only provides a connection to linear programming but also yields an algorithm for analyzing finite security markets. Moreover, we can make precise the connection between diverse expressions of economic equilibrium such as ‘absence of arbitrage’, ‘martingale property, and ‘complementary slackness property’.