We prove strong convergence and asymptotic normality for the record and the weak record rate of observations of the form Yn = Xn + Tn, n ≥ 1, where (Xn)n ∈ Z is a stationary ergodic sequence of random variables and (Tn)n ≥ 1 is a stochastic trend process with stationary ergodic increments. The strong convergence result follows from the Dubins-Freedman law of large numbers and Birkhoff's ergodic theorem. For the asymptotic normality we rely on the approach of Ballerini and Resnick (1987), coupled with a moment bound for stationary sequences, which is used to deal with the random trend process. Examples of applications are provided. In particular, we obtain strong convergence and asymptotic normality for the number of ladder epochs in a random walk with stationary ergodic increments.

This paper considers Markov decision processes (MDPs) with unbounded rates, as a function of state. We are especially interested in studying structural properties of optimal policies and the value function. A common method to derive such properties is by value iteration applied to the uniformised MDP. However, due to the unboundedness of the rates, uniformisation is not possible, and so value iteration cannot be applied in the way we need. To circumvent this, one can perturb the MDP. Then we need two results for the perturbed sequence of MDPs: 1. there exists a unique solution to the discounted cost optimality equation for each perturbation as well as for the original MDP; 2. if the perturbed sequence of MDPs converges in a suitable manner then the associated optimal policies and the value function should converge as well. We can model both the MDP and perturbed MDPs as a collection of parametrised Markov processes. Then both of the results above are essentially implied by certain continuity properties of the process as a function of the parameter. In this paper we deduce tight verifiable conditions that imply the necessary continuity properties. The most important of these conditions are drift conditions that are strongly related to nonexplosiveness.

Nonlinear filtering is investigated in a system where both the signal system and the observation system are under non-Gaussian Lévy fluctuations. Firstly, the Zakai equation is derived, and it is further used to derive the Kushner-Stratonovich equation. Secondly, by a filtered martingale problem, uniqueness for strong solutions of the Kushner-Stratonovich equation and the Zakai equation is proved. Thirdly, under some extra regularity conditions, the Zakai equation for the unnormalized density is also derived in the case of α-stable Lévy noise.

In this paper we consider a classic problem concerning the high excursion probabilities of a Gaussian random field f living on a compact set T. We develop efficient computational methods for the tail probabilities ℙ{supTf(t) > b}. For each positive ε, we present Monte Carlo algorithms that run in constant time and compute the probabilities with relative error ε for arbitrarily large b. The efficiency results are applicable to a large class of Hölder continuous Gaussian random fields. Besides computations, the change of measure and its analysis techniques have several theoretical and practical indications in the asymptotic analysis of Gaussian random fields.

We consider a branching random walk. Biggins and Kyprianou (2004) proved that, in the boundary case, the associated derivative martingale converges almost surely to a finite nonnegative limit, whose law serves as a fixed point of a smoothing transformation (Mandelbrot's cascade). In this paper, we give a necessary and sufficient condition for the nontriviality of the limit in this boundary case.

We investigate the large deviation properties of the maximum likelihood estimators for the Ornstein-Uhlenbeck process with shift. We propose a new approach to establish large deviation principles which allows us, via a suitable transformation, to circumvent the classical nonsteepness problem. We estimate simultaneously the drift and shift parameters. On the one hand, we prove a large deviation principle for the maximum likelihood estimates of the drift and shift parameters. Surprisingly, we find that the drift estimator shares the same large deviation principle as the estimator previously established for the Ornstein-Uhlenbeck process without shift. Sharp large deviation principles are also provided. On the other hand, we show that the maximum likelihood estimator of the shift parameter satisfies a large deviation principle with a very unusual implicit rate function.

We present the first class of perfect sampling (also known as exact simulation) algorithms for the steady-state distribution of non-Markovian loss systems. We use a variation of dominated coupling from the past. We first simulate a stationary infinite server system backwards in time and analyze the running time in heavy traffic. In particular, we are able to simulate stationary renewal marked point processes in unbounded regions. We then use the infinite server system as an upper bound process to simulate the loss system. The running time analysis of our perfect sampling algorithm for loss systems is performed in the quality-driven (QD) and the quality-and-efficiency-driven regimes. In both cases, we show that our algorithm achieves subexponential complexity as both the number of servers and the arrival rate increase. Moreover, in the QD regime, our algorithm achieves a nearly optimal rate of complexity.

For a class of cell division processes in the Euclidean space ℝd, spatial consistency is investigated. This addresses the problem whether the distribution of the generated structures, restricted to a bounded set V, depends on the choice of a larger region W ⊃ V where the construction of the cell division process is performed. This can also be understood as the problem of boundary effects in the cell division procedure. It is known that the STIT tessellations are spatially consistent. In the present paper it is shown that, within a reasonable wide class of cell division processes, the STIT tessellations are the only ones that are consistent.

For any strictly positive martingale S = eX for which X has a characteristic function, we provide an expansion for the implied volatility. This expansion is explicit in the sense that it involves no integrals, but only polynomials in the log-strike. We illustrate the versatility of our expansion by computing the approximate implied volatility smile in three well-known martingale models: one finite activity exponential Lévy model, Merton (1976), one infinite activity exponential Lévy model (variance gamma), and one stochastic volatility model, Heston (1993). Finally, we illustrate how our expansion can be used to perform a model-free calibration of the empirically observed implied volatility surface.

A random k-out mapping (digraph) on [n] is generated by choosing k random images of each vertex one at a time, subject to a 'preferential attachment' rule: the current vertex selects an image i with probability proportional to a given parameter α = α(n) plus the number of times i has already been selected. Intuitively, the larger α becomes, the closer the resulting k-out mapping is to the uniformly random k-out mapping. We prove that α = Θ(n1/2) is the threshold for α growing 'fast enough' to make the random digraph approach the uniformly random digraph in terms of the total variation distance. We also determine an exact limit for this distance for the α = βn1/2 case.

This paper is concerned with the bottom-up hierarchical system and public debate model proposed by Galam (2008), as well as a spatial version of the public debate model. In all three models, there is a population of individuals who are characterized by one of two competing opinions, say opinion −1 and opinion +1. This population is further divided into groups of common size s. In the bottom-up hierarchical system, each group elects a representative candidate, whereas in the other two models, all the members of each group discuss at random times until they reach a consensus. At each election/discussion, the winning opinion is chosen according to Galam's majority rule: the opinion with the majority of representatives wins when there is a strict majority, while one opinion, say opinion −1, is chosen by default in the case of a tie. For the public debate models we also consider the following natural updating rule that we call proportional rule: the winning opinion is chosen at random with a probability equal to the fraction of its supporters in the group. The three models differ in term of their population structure: in the bottom-up hierarchical system, individuals are located on a finite regular tree, in the nonspatial public debate model, they are located on a complete graph, and in the spatial public debate model, they are located on the d-dimensional regular lattice. For the bottom-up hierarchical system and nonspatial public debate model, Galam studied the probability that a given opinion wins under the majority rule and, assuming that individuals' opinions are initially independent, making the initial number of supporters of a given opinion a binomial random variable. The first objective of this paper is to revisit Galam's result, assuming that the initial number of individuals in favor of a given opinion is a fixed deterministic number. Our analysis reveals phase transitions that are sharper under our assumption than under Galam's assumption, particularly with small population size. The second objective is to determine whether both opinions can coexist at equilibrium for the spatial public debate model under the proportional rule, which depends on the spatial dimension.

We study a multistage epidemic model which generalizes the SIR model and where infected individuals go through K ≥ 1 stages of the epidemic before being removed. An infected individual in stage k ∈ {1, …, K} may infect a susceptible individual, who directly goes to stage k of the epidemic; or it may go to the next stage k + 1 of the epidemic. For this model, we identify the critical regime in which we establish diffusion approximations. Surprisingly, the limiting diffusion exhibits an unusual form of state space collapse which we analyze in detail.

The calculation of multivariate normal probabilities is of great importance in many statistical and economic applications. In this paper we propose a spherical Monte Carlo method with both theoretical analysis and numerical simulation. We start by writing the multivariate normal probability via an inner radial integral and an outer spherical integral using the spherical transformation. For the outer spherical integral, we apply an integration rule by randomly rotating a predetermined set of well-located points. To find the desired set, we derive an upper bound for the variance of the Monte Carlo estimator and propose a set which is related to the kissing number problem in sphere packings. For the inner radial integral, we employ the idea of antithetic variates and identify certain conditions so that variance reduction is guaranteed. Extensive Monte Carlo simulations on some probabilities confirm these claims.

In this paper we consider the beta(2 − α, α)-coalescents with 1 < α < 2 and study the moments of external branches, in particular, the total external branch length of an initial sample of n individuals. For this class of coalescents, it has been proved that nα-1T(n) →DT, where T(n) is the length of an external branch chosen at random and T is a known nonnegative random variable. For beta(2 − α, α)-coalescents with 1 < α < 2, we obtain limn→+∞n3α-5 𝔼(Lext(n) − n2-α𝔼T)2 = ((α − 1)Γ(α + 1))2Γ(4 − α) / ((3 − α)Γ(4 − 2α)).

The class of restless bandits as proposed by Whittle (1988) have long been known to be intractable. This paper presents an optimality result which extends that of Weber and Weiss (1990) for restless bandits to a more general setting in which individual bandits have multiple levels of activation but are subject to an overall resource constraint. The contribution is motivated by the recent works of Glazebrook et al. (2011a), (2011b) who discussed the performance of index heuristics for resource allocation in such systems. Hitherto, index heuristics have been shown, under a condition of full indexability, to be optimal for a natural Lagrangian relaxation of such problems in which a resource is purchased rather than constrained. We find that under key assumptions about the nature of solutions to a deterministic differential equation that the index heuristics above are asymptotically optimal in a sense described by Whittle. We then demonstrate that these assumptions always hold for three-state bandits.

We provide a characterization of realisable set covariograms, bringing a rigorous yet abstract solution to the S2 problem in materials science. Our method is based on the covariogram functional for random measurable sets (RAMS) and on a result about the representation of positive operators on a noncompact space. RAMS are an alternative to the classical random closed sets in stochastic geometry and geostatistics, and they provide a weaker framework that allows the manipulation of more irregular functionals, such as the perimeter. We therefore use the illustration provided by the S2 problem to advocate the use of RAMS for solving theoretical problems of a geometric nature. Along the way, we extend the theory of random measurable sets, and in particular the local approximation of the perimeter by local covariograms.

This paper is concerned with the solution of the optimal stopping problem associated to the value of American options driven by continuous-time Markov chains. The value-function of an American option in this setting is characterised as the unique solution (in a distributional sense) of a system of variational inequalities. Furthermore, with continuous and smooth fit principles not applicable in this discrete state-space setting, a novel explicit characterisation is provided of the optimal stopping boundary in terms of the generator of the underlying Markov chain. Subsequently, an algorithm is presented for the valuation of American options under Markov chain models. By application to a suitably chosen sequence of Markov chains, the algorithm provides an approximate valuation of an American option under a class of Markov models that includes diffusion models, exponential Lévy models, and stochastic differential equations driven by Lévy processes. Numerical experiments for a range of different models suggest that the approximation algorithm is flexible and accurate. A proof of convergence is also provided.