On the excursion area of perturbed Gaussian fields

We investigate Lipschitz-Killing curvatures for excursion sets of random fields on $\mathbb R^2$ under small spatial-invariant random perturbations. An expansion formula for mean curvatures is derived when the magnitude of the perturbation vanishes, which recovers the Gaussian Kinematic Formula at the limit by contiguity of the model. We develop an asymptotic study of the perturbed excursion area behaviour that leads to a quantitative non-Gaussian limit theorem, in Wasserstein distance, for fixed small perturbations and growing domain. When letting both the perturbation vanish and the domain grow, a standard Central Limit Theorem follows. Taking advantage of these results, we propose an estimator for the perturbation which turns out to be asymptotically normal and unbiased, allowing to make inference through sparse information on the field

Plug-in estimation of level sets in a non-compact setting with applications in multivariate risk theory

This paper deals with the problem of estimating the level sets of an unknown distribution function $F$. A plug-in approach is followed. That is, given a consistent estimator $F_n$ of $F$, we estimate the level sets of $F$ by the level sets of $F_n$. In our setting no compactness property is a priori required for the level sets to estimate. We state consistency results with respect to the Hausdorff distance and the volume of the symmetric difference. Our results are motivated by applications in multivariate risk theory. In this sense we also present simulated and real examples which illustrate our theoretical results.Level sets ; Distribution function ; Plug-in estimation ; Hausdorff distance ; Conditional Tail Expectation

Estimation of extreme L1L^1-multivariate expectiles with functional covariates

The present article is devoted to the semi-parametric estimation of multivariate expectiles for extreme levels. The considered multivariate risk measures also include the possible conditioning with respect to a functional covariate, belonging to an infinite-dimensional space. By using the first order optimality condition, we interpret these expectiles as solutions of a multidimensional nonlinear optimum problem. Then the inference is based on a minimization algorithm of gradient descent type, coupled with consistent kernel estimations of our key statistical quantities such as conditional quantiles, conditional tail index and conditional tail dependence functions. The method is valid for equivalently heavy-tailed marginals and under a multivariate regular variation condition on the underlying unknown random vector with arbitrary dependence structure. Our main result establishes the consistency in probability of the optimum approximated solution vectors with a speed rate. This allows us to estimate the global computational cost of the whole procedure according to the data sample size.Comment: 21 page

A local statistic for the spatial extent of extreme threshold exceedances

We introduce the extremal range, a local statistic for studying the spatial extent of extreme events in random fields on $\mathbb{R}^2$. Conditioned on exceedance of a high threshold at a location $s$, the extremal range at $s$ is the random variable defined as the smallest distance from $s$ to a location where there is a non-exceedance. We leverage tools from excursion-set theory to study distributional properties of the extremal range, propose parametric models and predict the median extremal range at extreme threshold levels. The extremal range captures the rate at which the spatial extent of conditional extreme events scales for increasingly high thresholds, and we relate its distributional properties with the bivariate tail dependence coefficient and the extremal index of time series in classical Extreme-Value Theory. Consistent estimation of the distribution function of the extremal range for stationary random fields is proven. For non-stationary random fields, we implement generalized additive median regression to predict extremal-range maps at very high threshold levels. An application to two large daily temperature datasets, namely reanalyses and climate-model simulations for France, highlights decreasing extremal dependence for increasing threshold levels and reveals strong differences in joint tail decay rates between reanalyses and simulations.Comment: 32 pages, 5 figure

Identifying regions of concomitant compound precipitation and wind speed extremes over Europe

The task of simplifying the complex spatio-temporal variables associated with climate modeling is of utmost importance and comes with significant challenges. In this research, our primary objective is to tailor clustering techniques to handle compound extreme events within gridded climate data across Europe. Specifically, we intend to identify subregions that display asymptotic independence concerning compound precipitation and wind speed extremes. To achieve this, we utilise daily precipitation sums and daily maximum wind speed data derived from the ERA5 reanalysis dataset spanning from 1979 to 2022. Our approach hinges on a tuning parameter and the application of a divergence measure to spotlight disparities in extremal dependence structures without relying on specific parametric assumptions. We propose a data-driven approach to determine the tuning parameter. This enables us to generate clusters that are spatially concentrated, which can provide more insightful information about the regional distribution of compound precipitation and wind speed extremes. In the process, we aim to elucidate the respective roles of extreme precipitation and wind speed in the resulting clusters. The proposed method is able to extract valuable information about extreme compound events while also significantly reducing the size of the dataset within reasonable computational timeframes.Comment: 34 pages, 15 figure

Estimation of the Multivariate Conditional-Tail-Expectation for extreme risk levels: illustrations on environmental data-sets

This paper deals with the problem of estimating the Multivariate version of the Conditional-Tail-Expectation introduced in the bivariate framework in Di Bernardino et al. [16], and generalized in Cousin and Di Bernardino [13]. We propose a new semi-parametric estimator for this risk measure, essentially based on statistical extrapolation techniques, well designed for extreme risk levels. Following Cai et al. [9], we prove a central limit theorem. We illustrate the practical properties of our estimator on simulations. The performances of our new estimator are discussed and compared to the ones of the empirical Kendall's process based estimator, previously proposed in Di Bernardino and Prieur [17]. We conclude with two applications on real data-sets: rainfall measurements recorded at three stations located in the south of Paris (France) and the analysis of strong wind gusts in the north west of France

Plug-in estimation of level sets in a non-compact setting with applications in multivariate risk theory

International audienceThis paper deals with the problem of estimating the level sets of an unknown distribution function $F$. A plug-in approach is followed. That is, given a consistent estimator $F_n$ of $F$, we estimate the level sets of $F$ by the level sets of $F_n$. In our setting no compactness property is a priori required for the level sets to estimate. We state consistency results with respect to the Hausdorff distance and the volume of the symmetric difference. Our results are motivated by applications in multivariate risk theory. In this sense we also present simulated and real examples which illustrate our theoretical results

Colorectal Cancer Stage at Diagnosis Before vs During the COVID-19 Pandemic in Italy

IMPORTANCE Delays in screening programs and the reluctance of patients to seek medical attention because of the outbreak of SARS-CoV-2 could be associated with the risk of more advanced colorectal cancers at diagnosis. OBJECTIVE To evaluate whether the SARS-CoV-2 pandemic was associated with more advanced oncologic stage and change in clinical presentation for patients with colorectal cancer. DESIGN, SETTING, AND PARTICIPANTS This retrospective, multicenter cohort study included all 17 938 adult patients who underwent surgery for colorectal cancer from March 1, 2020, to December 31, 2021 (pandemic period), and from January 1, 2018, to February 29, 2020 (prepandemic period), in 81 participating centers in Italy, including tertiary centers and community hospitals. Follow-up was 30 days from surgery. EXPOSURES Any type of surgical procedure for colorectal cancer, including explorative surgery, palliative procedures, and atypical or segmental resections. MAIN OUTCOMES AND MEASURES The primary outcome was advanced stage of colorectal cancer at diagnosis. Secondary outcomes were distant metastasis, T4 stage, aggressive biology (defined as cancer with at least 1 of the following characteristics: signet ring cells, mucinous tumor, budding, lymphovascular invasion, perineural invasion, and lymphangitis), stenotic lesion, emergency surgery, and palliative surgery. The independent association between the pandemic period and the outcomes was assessed using multivariate random-effects logistic regression, with hospital as the cluster variable. RESULTS A total of 17 938 patients (10 007 men [55.8%]; mean [SD] age, 70.6 [12.2] years) underwent surgery for colorectal cancer: 7796 (43.5%) during the pandemic period and 10 142 (56.5%) during the prepandemic period. Logistic regression indicated that the pandemic period was significantly associated with an increased rate of advanced-stage colorectal cancer (odds ratio [OR], 1.07; 95%CI, 1.01-1.13; P = .03), aggressive biology (OR, 1.32; 95%CI, 1.15-1.53; P < .001), and stenotic lesions (OR, 1.15; 95%CI, 1.01-1.31; P = .03). CONCLUSIONS AND RELEVANCE This cohort study suggests a significant association between the SARS-CoV-2 pandemic and the risk of a more advanced oncologic stage at diagnosis among patients undergoing surgery for colorectal cancer and might indicate a potential reduction of survival for these patients

On an asymmetric extension of multivariate Archimedean copulas based on quadratic form

International audienceAn important topic in Quantitative Risk Management concerns the modeling of dependence among risk sources and in this regard Archimedean copulas appear to be very useful. However, they exhibit symmetry, which is not always consistent with patterns observed in real world data. We investigate extensions of the Archimedean copula family that make it possible to deal with asymmetry. Our extension is based on the observation that when applied to the copula the inverse function of the generator of an Archimedean copula can be expressed as a linear form of generator inverses. We propose to add a distortion term to this linear part, which leads to asymmetric copulas. Parameters of this new class of copulas are grouped within a matrix, thus facilitating some usual applications as level curve determination or estimation. Some choices such as sub-model stability help associating each parameter to one bivariate projection of the copula. We also give some admissibility conditions for the considered copulas. We propose different examples as some natural multivariate extensions of Farlie-Gumbel-Morgenstern or Gumbel-Barnett

On certain transformations of Archimedean copulas: Application to the non-parametric estimation of their generators

We study the impact of certain transformations within the class of Archimedean copulas. We give some admissibility conditions for these transformations, and define some equivalence classes for both transformations and generators of Archimedean copulas. We extend the r-fold composition of the diagonal section of a copula, from r ∈ N to r ∈ R. This extension, coupled with results on equivalence classes, gives us new expressions of transformations and generators. Estimators deriving directly from these expressions are proposed and their convergence is investigated. We provide confidence bands for the estimated generators. Numerical illustrations show the empirical performance of these estimators