Source author record

Clément Dombry

Clément Dombry appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2021arXiv

Hidden regular variation for point processes and the single/multiple large point heuristic

We consider regular variation for marked point processes with independent heavy-tailed marks and prove a single large point heuristic: the limit measure is concentrated on the cone of point measures with one single point. We then investigate successive hidden regular variation removing the cone of point measures with at most $k$ points, $k\geq 1$, and prove a multiple large point phenomenon: the limit measure is concentrated on the cone of point measures with $k+1$ points. We show how these results imply hidden regular variation in Skorokhod space of the associated risk process. Finally, we provide an application to risk theory in a reinsurance model where the $k$ largest claims are covered and we study the asymptotic behavior of the residual risk.

preprint2020arXiv

Behavior of linear L2-boosting algorithms in the vanishing learning rate asymptotic

We investigate the asymptotic behaviour of gradient boosting algorithms when the learning rate converges to zero and the number of iterations is rescaled accordingly. We mostly consider L2-boosting for regression with linear base learner as studied in B{ü}hlmann and Yu (2003) and analyze also a stochastic version of the model where subsampling is used at each step (Friedman 2002). We prove a deterministic limit in the vanishing learning rate asymptotic and characterize the limit as the unique solution of a linear differential equation in an infinite dimensional function space. Besides, the training and test error of the limiting procedure are thoroughly analyzed. We finally illustrate and discuss our result on a simple numerical experiment where the linear L2-boosting operator is interpreted as a smoothed projection and time is related to its number of degrees of freedom.

preprint2020arXiv

Extreme quantile regression in a proportional tail framework

We revisit the model of heteroscedastic extremes initially introduced by Einmahl et al. (JRSSB, 2016) to describe the evolution of a non stationary sequence whose extremes evolve over time and adapt it into a general extreme quantile regression framework. We provide estimates for the extreme value index and the integrated skedasis function and prove their asymptotic normality. Our results are quite similar to those developed for heteroscedastic extremes but with a different proof approach emphasizing coupling arguments. We also propose a pointwise estimator of the skedasis function and a Weissman estimator of the conditional extreme quantile and prove the asymptotic normality of both estimators.

preprint2016arXiv

A Semi-Definite Programming approach to low dimensional embedding for unsupervised clustering

This paper proposes a variant of the method of Guédon and Verhynin for estimating the cluster matrix in the Mixture of Gaussians framework via Semi-Definite Programming. A clustering oriented embedding is deduced from this estimate. The procedure is suitable for very high dimensional data because it is based on pairwise distances only. Theoretical garantees are provided and an eigenvalue optimisation approach is proposed for computing the embedding. The performance of the method is illustrated via Monte Carlo experiements and comparisons with other embeddings from the literature.

preprint2016arXiv

Ergodic decompositions of stationary max-stable processes in terms of their spectral functions

We revisit conservative/dissipative and positive/null decompositions of stationary max-stable processes. Originally, both decompositions were defined in an abstract way based on the underlying non-singular flow representation. We provide simple criteria which allow to tell whether a given spectral function belongs to the conservative/dissipative or positive/null part of the de Haan spectral representation. Specifically, we prove that a spectral function is null-recurrent iff it converges to $0$ in the Cesàro sense. For processes with locally bounded sample paths we show that a spectral function is dissipative iff it converges to $0$. Surprisingly, for such processes a spectral function is integrable a.s. iff it converges to $0$ a.s. Based on these results, we provide new criteria for ergodicity, mixing, and existence of a mixed moving maximum representation of a stationary max-stable process in terms of its spectral functions. In particular, we study a decomposition of max-stable processes which characterizes the mixing property.

preprint2016arXiv

High order expansions for renewal functions and applications to ruin theory

A high order expansion of the renewal function is provided under the assumption that the inter-renewal time distribution is light tailed with finite moment generating function g on a neighborhood of 0. This expansion relies on complex analysis and is expressed in terms of the residues of the function 1/(1 -- g). Under the assumption that g can be extended into a meromorphic function on the complex plane and some technical conditions, we obtain even an exact expansion of the renewal function. An application to risk theory is given where we consider high order expansion of the ruin probability for the standard compound Poisson risk model. This precises the well known Crámer-Lundberg approximation of the ruin probability when the initial reserve is large.

preprint2016arXiv

Random tessellations associated with max-stable random fields

With any max-stable random process $η$ on $\mathcal{X}=\mathbb{Z}^d$ or $\mathbb{R}^d$, we associate a random tessellation of the parameter space $\mathcal{X}$. The construction relies on the Poisson point process representation of the max-stable process $η$ which is seen as the pointwise maximum of a random collection of functions $Φ=\{ϕ\_i, i\geq 1\}$. The tessellation is constructed as follows: two points $x,y\in \mathcal{X}$ are in the same cell if and only if there exists a function $ϕ\inΦ$ that realizes the maximum $η$ at both points $x$ and $y$, i.e. $ϕ(x)=η(x)$ and $ϕ(y)=η(y)$. We characterize the distribution of cells in terms of coverage and inclusion probabilities. Most interesting is the stationary case where the asymptotic properties of the cells are strongly related to the ergodic properties of the non-singular flow generating the max-stable process. For example, we show that: i) the cells are bounded almost surely if and only if $η$ is generated by a dissipative flow, ii) the cells have positive asymptotic density almost surely if and only if $η$ is generated by a positive flow.

preprint2015arXiv

Exact simulation of max-stable processes

Max-stable processes play an important role as models for spatial extreme events. Their complex structure as the pointwise maximum over an infinite number of random functions makes simulation highly nontrivial. Algorithms based on finite approximations that are used in practice are often not exact and computationally inefficient. We will present two algorithms for exact simulation of a max-stable process at a finite number of locations. The first algorithm generalizes the approach by \citet{DM-2014} for Brown--Resnick processes and it is based on simulation from the spectral measure. The second algorithm relies on the idea to simulate only the extremal functions, that is, those functions in the construction of a max-stable process that effectively contribute to the pointwise maximum. We study the complexity of both algorithms and prove that the second procedure is always more efficient. Moreover, we provide closed expressions for their implementation that cover the most popular models for max-stable processes and extreme value copulas. For simulation on dense grids, an adaptive design of the second algorithm is proposed.

preprint2015arXiv

On functional records and champions

Records among a sequence of iid random variables $X_1,X_2,\dotsc$ on the real line have been investigated extensively over the past decades. A record is defined as a random variable $X_n$ such that $X_n>\max(X_1,\dotsc,X_{n-1})$. Trying to generalize this concept to the case of random vectors, or even stochastic processes with continuous sample paths, the question arises how to define records in higher dimensions. We introduce two different concepts: A simple record is meant to be a stochastic process (or a random vector) $X_n$ that is larger than $ X_1,\dotsc, X_{n-1}$ in at least one component, whereas a complete record has to be larger than its predecessors in all components. The behavior of records is investigated. In particular, the probability that a stochastic process $ X_n$ is a record as $n$ tends to infinity is studied, assuming that the processes are in the max-domain of attraction of a max-stable process. Furthermore, the distribution of $ X_n$, given that $ X_n$ is a record is derived.

preprint2015arXiv

Probabilities of concurrent extremes

The statistical modelling of spatial extremes has recently made major advances. Much of its focus so far has been on the modelling of the magnitudes of extreme events but little attention has been paid on the timing of extremes. To address this gap, this paper introduces the notion of extremal concurrence. Suppose that one measures precipitation at several synoptic stations over multiple days. We say that extremes are concurrent if the maximum precipitation over time at each station is achieved simultaneously, e.g., on a single day. Under general conditions, we show that the finite sample concurrence probability converges to an asymptotic quantity, deemed extremal concurrence probability. Using Palm calculus, we establish general expressions for the extremal concurrence probability through the max-stable process emerging in the limit of the componentwise maxima of the sample. Explicit forms of the extremal concurrence probabilities are obtained for various max-stable models and several estimators are introduced. In particular, we prove that the pairwise extremal concurrence probability for max-stable vectors is precisely equal to the Kendall's $τ$. The estimators are evaluated by using simulations and applied to study the concurrence patterns of temperature extremes in the United States. The results demonstrate that concurrence probability can provide a powerful new perspective and tools for the analysis of the spatial structure and impact of extremes.

preprint2013arXiv

Approximation of stable random measures and applications to linear fractional stable integrals

Using lattice approximations of Euclidean space, we develop a way to approximate stable processes that are represented by stochastic integrals over Euclidean space. Via a stable version of the Lindeberg-Feller Theorem we show that the approximations weakly converge as the mesh-size goes to zero. As an application, we improve upon previous approximation schemes for integrals with respect to linear fractional stable motions.

preprint2013arXiv

Maximum likelihood estimators for the extreme value index based on the block maxima method

The maximum likelihood method offers a standard way to estimate the three parameters of a generalized extreme value (GEV) distribution. Combined with the block maxima method, it is often used in practice to assess the extreme value index and normalization constants of a distribution satisfying a first order extreme value condition, assuming implicitely that the block maxima are exactly GEV distributed. This is unsatisfactory since the GEV distribution is a good approximation of the block maxima distribution only for blocks of large size. The purpose of this paper is to provide a theoretical basis for this methodology. Under a first order extreme value condition only, we prove the existence and consistency of the maximum likelihood estimators for the extreme value index and normalization constants within the framework of the block maxima method.

preprint2012arXiv

Conditional simulation of extremal Gaussian processes

Recently the regular conditional distributions of max-infinitely divisible processes were derived by \citet{Dombry2011} and although these conditional distributions have complicated closed forms, \citet{Dombry2011b} introduce an algorithm to get conditional realizations of Brown-Resnick processes. In this paper we derive the regular conditional distributions of the max-stable process introduced by \citet{Schlather2002} and adapt the framework of \citet{Dombry2011b} to this specific process. We test the methods on simulated data and give an application to extreme temperatures in Switzerland. Results show that the proposed sampling scheme provide accurate conditional simulations and can handle real-sized problems.

preprint2012arXiv

Conditional simulation of max-stable processes

Since many environmental processes such as heat waves or precipitation are spatial in extent, it is likely that a single extreme event affects several locations and the areal modelling of extremes is therefore essential if the spatial dependence of extremes has to be appropriately taken into account. This paper proposes a framework for conditional simulations of max-stable processes and give closed forms for Brown-Resnick and Schlather processes. We test the method on simulated data and give an application to extreme rainfall around Zurich and extreme temperature in Switzerland. Results show that the proposed framework provides accurate conditional simulations and can handle real-sized problems.

preprint2012arXiv

Conditional simulations of Brown-Resnick processes

Since many environmental processes such as heat waves or precipitation are spatial in extent, it is likely that a single extreme event affects several locations and the areal modeling of extremes is therefore essential if the spatial dependence of extremes has to be appropriately taken into account. Although some progress has been made to develop a geostatistic of extremes, conditional simulation of max-stable processes is still in its early stage. This paper proposes a framework to get conditional simulations of Brown-Resnick processes. Although closed forms for the regular conditional distribution of Brown-Resnick processes were recently found, sampling from this conditional distribution is a considerable challenge as it leads quickly to a combinatorial explosion. To bypass this computational burden, a Markov chain Monte-Carlo algorithm is presented. We test the method on simulated data and give an application to extreme rainfall around Zurich. Results show that the proposed framework provides accurate conditional simulations of Brown-Resnick processes and can handle real-sized problems.

preprint2012arXiv

Regular conditional distributions of max infinitely divisible processes

This paper is devoted to the prediction problem in extreme value theory. Our main result is an explicit expression of the regular conditional distribution of a max-stable (or max-infinitely divisible) process $\{η(t)\}_{t\in T}$ given observations $\{η(t_i)=y_i,\ 1\leq i\leq k\}$. Our starting point is the point process representation of max-infinitely divisible processes by Giné, Hahn and Vatan (1990). We carefully analyze the structure of the underlying point process, introduce the notions of extremal function, sub-extremal function and hitting scenario associated to the constraints and derive the associated distributions. This allows us to explicit the conditional distribution as a mixture over all hitting scenarios compatible with the conditioning constraints. This formula extends a recent related result by Wang and Stoev (2011) dealing with the case of spectrally discrete max-stable random fields. We believe this work offers new tools and perspective for prediction in extreme value theory together with numerous potential applications.

preprint2012arXiv

Strong mixing properties of max-infinitely divisible random fields

Let $η=(η(t))_{t\in T}$ be a sample continuous max-infinitely random field on a locally compact metric space $T$. For a closed subset $S\in T$, we note $η_{S}$ the restriction of $η$ to $S$. We consider $β(S_1,S_2)$ the absolute regularity coefficient between $η_{S_1}$ and $η_{S_2}$, where $S_1,S_2$ are two disjoint closed subsets of $T$. Our main result is a simple upper bound for $β(S_1,S_2)$ involving the exponent measure $μ$ of $η$: we prove that $β(S_1,S_2)\leq 2\int \bbP[η\not<_{S_1} f,\ η\not <_{S_2} f]\,μ(df)$, where $f\not<_{S} g$ means that there exists $s\in S$ such that $f(s)\geq g(s)$. If $η$ is a simple max-stable random field, the upper bound is related to the so-called extremal coefficients: for countable disjoint sets $S_1$ and $S_2$, we obtain $β(S_1,S_2)\leq 4\sum_{(s_1,s_2)\in S_1\times S_2}(2-θ(s_1,s_2))$, where $θ(s_1,s_2)$ is the pair extremal coefficient. As an application, we show that these new estimates entail a central limit theorem for stationary max-infinitely divisible random fields on $\bbZ^d$. In the stationary max-stable case, we derive the asymptotic normality of three simple estimators of the pair extremal coefficient.

preprint2011arXiv

Extremes of independent stochastic processes: a point process approach

For each $n\geq 1$, let $ {X_{in}, \quad i \geq 1} $ be independent copies of a nonnegative continuous stochastic process $X_{n}=(X_n(t))_{t\in T}$ indexed by a compact metric space $T$. We are interested in the process of partial maxima [\tilde M_n(u,t) =\max \{X_{in}(t), 1 \leq i\leq [nu]},\quad u\geq 0,\ t\in T.] where the brackets $[\,\cdot\,]$ denote the integer part. Under a regular variation condition on the sequence of processes $X_n$, we prove that the partial maxima process $\tilde M_n$ weakly converges to a superextremal process $\tilde M$ as $n\to\infty$. We use a point process approach based on the convergence of empirical measures. Properties of the limit process are investigated: we characterize its finite-dimensional distributions, prove that it satisfies an homogeneous Markov property, and show in some cases that it is max-stable and self-similar. Convergence of further order statistics is also considered. We illustrate our results on the class of log-normal processes in connection with some recent results on the extremes of Gaussian processes established by Kabluchko.

preprint2011arXiv

Stochastic algorithms for computing means of probability measures

Consider a probability measure supported by a regular geodesic ball in a manifold. For any p larger than or equal to 1 we define a stochastic algorithm which converges almost surely to the p-mean of the measure. Assuming furthermore that the functional to minimize is regular around the p-mean, we prove that a natural renormalization of the inhomogeneous Markov chain converges in law into an inhomogeneous diffusion process. We give an explicit expression of this process, as well as its local characteristic.

preprint2010arXiv

The on-off network traffic model under intermediate scaling

The result provided in this paper helps complete a unified picture of the scaling behavior in heavy-tailed stochastic models for transmission of packet traffic on high-speed communication links. Popular models include infinite source Poisson models, models based on aggregated renewal sequences, and models built from aggregated on-off sources. The versions of these models with finite variance transmission rate share the following pattern: if the sources connect at a fast rate over time the cumulative statistical fluctuations are fractional Brownian motion, if the connection rate is slow the traffic fluctuations are described by a stable Lévy process, while the limiting fluctuations for the intermediate scaling regime are given by fractional Poisson motion.