A bookmaker takes bets on a two-horse race, attempting to minimise expected loss over all possible outcomes of the race. Profits are controlled by manipulation of customers' betting behaviour, which is assumed to be determined uniquely by the price quoted for each horse. We consider different strategies for choosing these prices as bets accumulate, and examine the penalty incurred by the use of a strategy other than the optimal one, in the general case where the distribution of the customers' betting probabilities is unknown.

We consider an epidemic model where the spread of the epidemic can be described by a discrete-time Galton-Watson branching process. Between times n and n + 1, any infected individual is detected with unknown probability π and the numbers of these detected individuals are the only observations we have. Detected individuals produce a reduced number of offspring in the time interval of detection, and no offspring at all thereafter. If only the generation sizes of a Galton-Watson process are observed, it is known that one can only estimate the first two moments of the offspring distribution consistently on the explosion set of the process (and, apart from some lattice parameters, no parameters that are not determined by those moments). Somewhat surprisingly, in our context, where we observe a binomially distributed subset of each generation, we are able to estimate three functions of the parameters consistently. In concrete situations, this often enables us to estimate π consistently, as well as the mean number of offspring. We apply the estimators to data for a real epidemic of classical swine fever.

In this paper we develop a constructive structure theory for a class of exponential functionals of Brownian motion which includes Asian option values. This is done in two stages of differing natures. As a first step, the functionals are represented as Laguerre reduction series obtained from main results of Schröder (2006), this paper's companion paper. These reduction series are new and given in terms of the negative moments of the integral of geometric Brownian motion, whose structure theory is developed in a second step. Providing a new angle on these processes, this is done by establishing connections with theta functions. Integral representations and computable formulae for the negative moments are thus derived and then shown to furnish highly efficient ways for computing the negative moments. Application of this paper's Laguerre reduction series in numerical examples suggests that one of the most efficient methods for the explicit valuation of Asian options is obtained. The paper also provides mathematical background results referred to in Schröder (2005c).

In this paper we develop methods for reducing the study, the computation, and the construction of stochastic functionals to those of standard concepts such as the moments of the pertinent random variables. Principally, our methods are based on the notion of ladder height densities and their Laguerre expansions, and our results provide a unifying framework for the distinct approaches of Dufresne (2000) and Schröder (2005).

Let F(x) denote a distribution function in Rd and let F*n(x) denote the nth convolution power of F(x). In this paper we discuss the asymptotic behaviour of 1 - F*n(x) as x tends to ∞ in a certain prescribed way. It turns out that in many cases 1 - F*n(x) ∼ n(1 - F(x)). To obtain results of this type, we introduce and use a form of subexponential behaviour, thereby extending the notion of multivariate regular variation. We also discuss subordination, in which situation the index n is replaced by a random index N.

We extend the boson process first to a large class of Cox processes and second to an even larger class of infinitely divisible point processes. Density and moment results are studied in detail. These results are obtained in closed form as weighted permanents, so the extension is called a permanental process. Temporal extensions and a particularly tractable case of the permanental process are also studied. Extensions of the fermion process along similar lines, leading to so-called determinantal processes, are discussed.

We analyse SIS epidemics among populations partitioned into households. The analysis considers both the stochastic and deterministic models and, unlike in previous analyses, we consider general infectious period distributions. For the deterministic model, we prove the existence of an endemic equilibrium for the epidemic if and only if the threshold parameter, R*, is greater than 1. Furthermore, by utilising Markov chains we show that the total number of infectives converges to the endemic equilibrium as t → ∞. For the stochastic model, we prove a law of large numbers result for the convergence, to the deterministic limit, of the mean number of infectives per household. This is followed by the derivation of a Gaussian limit process for the fluctuations of the stochastic model.

We study a random field obtained by counting the number of balls containing a given point when overlapping balls are thrown at random according to a Poisson random measure. We describe a microscopic process which exhibits multifractional behavior. We are particularly interested in the local asymptotic self-similarity (LASS) properties of the field, as well as in its X-ray transform. We obtain two different LASS properties when considering the asymptotics either in law or in the sense of second-order moments, and prove a relationship between the LASS behavior of the field and the LASS behavior of its X-ray transform. These results can be used to model and analyze porous media, images, or connection networks.

The maximization of the long-term growth rate of expected utility is considered under drawdown constraints. In a general situation, the value and the optimal strategy of the problem are related to those of another ‘standard’ risk-sensitive-type portfolio optimization problem. Furthermore, an upside-chance maximization problem of a large deviation probability is stated as a ‘dual’ optimization problem. As an example, a ‘linear-quadratic’ model is studied in detail: the conditions to ensure the solvabilities of the problems are discussed and explicit expressions for the solutions are presented.

We present a new iterative procedure for solving the multiple stopping problem in discrete time and discuss the stability of the algorithm. The algorithm produces monotonically increasing approximations of the Snell envelope which coincide with the Snell envelope after finitely many steps. Unlike backward dynamic programming, the algorithm allows us to calculate approximative solutions with only a few nestings of conditional expectations and is, therefore, tailor-made for a plain Monte Carlo implementation.

Consider the random variable Ln defined as the length of a longest common subsequence of two random strings of length n and whose random characters are independent and identically distributed over a finite alphabet. Chvátal and Sankoff showed that the limit γ=limn→∞E[Ln]/n is well defined. The exact value of this constant is not known, but various methods for the computation of upper and lower bounds have been discussed in the literature. Even so, high-precision bounds are hard to come by. In this paper we discuss how large deviation theory can be used to derive a consistent sequence of upper bounds, (qm)m∈ℕ, on γ, and how Monte Carlo simulation can be used in theory to compute estimates, q̂m, of the qm such that, for given Ξ > 0 and Λ ∈ (0,1), we have P[γ < q̂ < γ + Ξ] ≥ Λ. In other words, with high probability the result is an upper bound that approximates γ to high precision. We establish O((1 − Λ)−1Ξ−(4+ε)) as a theoretical upper bound on the complexity of computing q̂m to the given level of accuracy and confidence. Finally, we discuss a practical heuristic based on our theoretical approach and discuss its empirical behavior.

We study the problem of dynamically allocating flexible workers to stations in tandem or serial manufacturing systems. Workers are trained to do a subset of consecutive tasks. We show that the optimal policy is often LBFS (last buffer first-served) or FBFS (first buffer first-served). These results generalize earlier results on the optimality of the pick-and-run, expedite, and bucket brigade-type policies. We also show that, for exponential processing times and general manufacturing networks, the optimal policy will tend to have several workers assigned to the same station.

We give three applications of the Pecherskii-Rogozin-Spitzer identity for Lévy processes. First, we find the joint distribution of the supremum and the epoch at which it is ‘attained’ if a Lévy process has phase-type upward jumps. We also find the characteristics of the ladder process. Second, we establish general properties of perturbed risk models, and obtain explicit fluctuation identities in the case that the Lévy process is spectrally positive. Third, we study the tail asymptotics for the supremum of a Lévy process under different assumptions on the tail of the Lévy measure.

In this paper we discuss the complementary theorem applied to the typical n-tuple of a Poisson point process. The theorem was first presented by Miles in 1970 and discussed by Santaló in 1976 and, within a Palm measure framework, by Møller and Zuyev in 1996. The theorems put forward by these authors are not correct for all the examples that they present, suggesting that further consideration of their work is needed if one wishes to bring all those examples within the ambit of the complementary theorem. We give alternative analyses of the errant examples and, with a modification of the technicalities in the work of the above authors, move toward a more comprehensive complementary theorem. Some open issues still remain.

In this paper we study the busy period in GI/G/1 work-conserving queues. We give the exact asymptotics of the tail distribution of the busy period under the light tail assumptions. We also study the workload process in the M/G/1 system conditioned to stay positive.