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.

We consider a queue where the server is the Euclidean space, and the customers are random closed sets (RACSs) of the Euclidean space. These RACSs arrive according to a Poisson rain and each of them has a random service time (in the case of hail falling on the Euclidean plane, this is the height of the hailstone, whereas the RACS is its footprint). The Euclidean space serves customers at speed 1. The service discipline is a hard exclusion rule: no two intersecting RACSs can be served simultaneously and service is in the first-in–first-out order, i.e. only the hailstones in contact with the ground melt at speed 1, whereas the others are queued. A tagged RACS waits until all RACSs that arrived before it and intersecting it have fully melted before starting its own melting. We give the evolution equations for this queue. We prove that it is stable for a sufficiently small arrival intensity, provided that the typical diameter of the RACS and the typical service time have finite exponential moments. We also discuss the percolation properties of the stationary regime of the RACS in the queue.

The limit behaviour in probability of realised quadratic variation is discussed under a relatively simple ambit process setting. The relation of this to the underlying volatility/intermittency field is in focus, especially as concerns the question of no volatility/intermittency memory.

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.

A Hawkes process is also known under the name of a self-exciting point process and has numerous applications throughout science and engineering. We derive the statistical estimation (maximum likelihood estimation) and goodness-of-fit (mainly graphical) for multivariate Hawkes processes with possibly dependent marks. As an application, we analyze two data sets from finance.

This paper provides a simple proof for the fact that the hitting time to an infrequently visited subset for a one-dependent regenerative process converges weakly to an exponential distribution. Special cases are positive recurrent Harris chains and Harris processes. The paper further extends this class of limit theorems to ‘rewards’ that are cumulated to the hitting time of such a rare set.

Let ξ be a random measure on a locally compact second countable topological group, and let X be a random element in a measurable space on which the group acts. In the compact case we give a natural definition of the concept that the origin is a typical location for X in the mass of ξ, and prove that when this holds, the same is true on sets placed uniformly at random around the origin. This new result motivates an extension of the concept of typicality to the locally compact case where it coincides with the concept of mass-stationarity. We describe recent developments in Palm theory where these ideas play a central role.

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.

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 consider the problem of estimating the probability that the maximum of a Gaussian process with negative mean and indexed by positive integers reaches a high level, say b. In great generality such a probability converges to 0 exponentially fast in a power of b. Under mild assumptions on the marginal distributions of the process and no assumption on the correlation structure, we develop an importance sampling procedure, called the target bridge sampler (TBS), which takes a polynomial (in b) number of function evaluations to achieve a small relative error. The procedure also yields samples of the underlying process conditioned on hitting b in finite time. In addition, we apply our method to the problem of estimating the tail of the maximum of a superposition of a large number, n, of independent Gaussian sources. In this situation TBS achieves a prescribed relative error with a bounded number of function evaluations as n ↗ ∞. A remarkable feature of TBS is that it is not based on exponential changes of measure. Our numerical experiments validate the performance indicated by our theoretical findings.

We investigate the conditional full support (CFS) property, introduced in Guasoni et al. (2008a), for Gaussian processes with stationary increments. We give integrability conditions on the spectral measure of such a process which ensure that the process has CFS or not. In particular, the general results imply that, for a process with spectral density f such that f(λ) ∼ c1λpe−c2λq for λ → ∞ (with necessarily p < 1 if q = 0), the CFS property holds if and only if q < 1.

We consider the following stochastic model for a mobile service scenario. Consider a stationary Poisson process in Rd, with its points radially ordered with respect to the origin (the anchor); if d = 2, the points may correspond to locations of, e.g. restaurants. A user, with a location different from the origin, asks for the location of the first Poisson point and keeps asking for the location of the next Poisson point until the first time that he/she can be completely certain that he/she knows which Poisson point is his/her nearest neighbour. This waiting time is the communication cost, while the inferred privacy region is a random set obtained by an adversary who only knows the anchor and the points received from the server, where the adversary ‘does the best’ to infer the possible locations of the user. Probabilistic results related to the communication cost and the inferred privacy region are established for any dimension d ≥ 1. Furthermore, special results when d = 1 and particularly when d = 2 are derived.

The modeling of random bi-phasic, or porous, media has been, and still is, under active investigation by mathematicians, physicists, and physicians. In this paper we consider a thresholded random process X as a source of the two phases. The intervals when X is in a given phase, named chords, are the subject of interest. We focus on the study of the tails of the chord length distribution functions. In the literature concerned with real data, different types of tail behavior have been reported, among them exponential-like or power-like decay. We look for the link between the dependence structure of the underlying thresholded process X and the rate of decay of the chord length distribution. When the process X is a stationary Gaussian process, we relate the latter to the rate at which the covariance function of X decays at large lags. We show that exponential, or nearly exponential, decay of the tail of the distribution of the chord lengths is very common, perhaps surprisingly so.

Let Xn be a sequence of integrable real random variables, adapted to a filtration (Gn). Define Cn = √{(1 / n)∑k=1nXk − E(Xn+1 | Gn)} and Dn = √n{E(Xn+1 | Gn) − Z}, where Z is the almost-sure limit of E(Xn+1 | Gn) (assumed to exist). Conditions for (Cn, Dn) → N(0, U) x N(0, V) stably are given, where U and V are certain random variables. In particular, under such conditions, we obtain √n{(1 / n)∑k=1nX_k - Z} = Cn + Dn → N(0, U + V) stably. This central limit theorem has natural applications to Bayesian statistics and urn problems. The latter are investigated, by paying special attention to multicolor randomly reinforced urns.

We provide necessary and sufficient conditions for the asymptotic normality of Nn, the number of records among the first n observations from a sequence of independent and identically distributed random variables, with general distribution F. In the case of normality we identify the centering and scaling sequences. Also, we characterize distributions for which the limit is not normal in terms of their discrete and continuous components.

In this paper we introduce conditionally independent increment point processes, that is, processes that are conditionally independent inside and outside a bounded set A given N(A), the number of points in A. We show that these point processes can be characterized by means of the avoidance function of a multinomial ‘support process’, the solution of a suitably defined linear system of equations, and, finally, the infinitesimal matrix of a continuous-time Markov chain.

This paper is motivated by relations between association and independence of random variables. It is well known that, for real random variables, independence implies association in the sense of Esary, Proschan and Walkup (1967), while, for random vectors, this simple relationship breaks. We modify the notion of association in such a way that any vector-valued process with independent increments also has associated increments in the new sense - association between blocks. The new notion is quite natural and admits nice characterization for some classes of processes. In particular, using the covariance interpolation formula due to Houdré, Pérez-Abreu and Surgailis (1998), we show that within the class of multidimensional Gaussian processes, block association of increments is equivalent to supermodularity (in time) of the covariance functions. We also define corresponding versions of weak association, positive association, and negative association. It turns out that the central limit theorem for weakly associated random vectors due to Burton, Dabrowski and Dehling (1986) remains valid, if the weak association is relaxed to the weak association between blocks.

Max-stable random fields play a central role in modeling extreme value phenomena. We obtain an explicit formula for the conditional probability in general max-linear models, which include a large class of max-stable random fields. As a consequence, we develop an algorithm for efficient and exact sampling from the conditional distributions. Our method provides a computational solution to the prediction problem for spectrally discrete max-stable random fields. This work offers new tools and a new perspective to many statistical inference problems for spatial extremes, arising, for example, in meteorology, geology, and environmental applications.

Coalescents with multiple collisions (also called Λ-coalescents or simple exchangeable coalescents) are used as models of genealogies. We study a new class of Markovian coalescent processes connected to a population model with immigration. Consider an infinite population with immigration labelled at each generation by N := {1, 2, …}. Some ancestral lineages cannot be followed backwards after some time because their ancestor is outside the population. The individuals with an immigrant ancestor constitute a distinguished family and we define exchangeable distinguished coalescent processes as a model for genealogy with immigration, focusing on simple distinguished coalescents, i.e. such that when a coagulation occurs all the blocks involved merge as a single block. These processes are characterized by two finite measures on [0, 1] denoted by M = (Λ0, Λ1). We call them M-coalescents. We show by martingale arguments that the condition of coming down from infinity for the M-coalescent coincides with that obtained by Schweinsberg for the Λ-coalescent. In the same vein as Bertoin and Le Gall, M-coalescents are associated with some stochastic flows. The superprocess embedded can be viewed as a generalized Fleming-Viot process with immigration. The measures Λ0 and Λ1 respectively specify the reproduction and the immigration. The coming down from infinity of the M-coalescent will be interpreted as the initial types extinction: after a certain time all individuals are immigrant children.