Source author record

Claudia Klüppelberg

Claudia Klüppelberg 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

28works
17topics
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

28 published item(s)

preprint2022arXiv

Max-linear models in random environment

We extend previous work of max-linear models on finite directed acyclic graphs to infinite graphs as well as random graphs, and investigate their relations to classical percolation theory, more particularly the impact of Bernoulli bond percolation on such models. We show that the critical probability of percolation on the oriented square lattice graph $\mathbb{Z}^2$ describes a phase transition in the obtained model. Focus is on the dependence introduced by this graph into the max-linear model. We discuss natural applications in communication networks, in particular, concerning the propagation of influences.

preprint2021arXiv

Indirect Inference for Time Series Using the Empirical Characteristic Function and Control Variates

We estimate the parameter of a stationary time series process by minimizing the integrated weighted mean squared error between the empirical and simulated characteristic function, when the true characteristic functions cannot be explicitly computed. Motivated by Indirect Inference, we use a Monte Carlo approximation of the characteristic function based on iid simulated blocks. As a classical variance reduction technique, we propose the use of control variates for reducing the variance of this Monte Carlo approximation. These two approximations yield two new estimators that are applicable to a large class of time series processes. We show consistency and asymptotic normality of the parameter estimators under strong mixing, moment conditions, and smoothness of the simulated blocks with respect to its parameter. In a simulation study we show the good performance of these new simulation based estimators, and the superiority of the control variates based estimator for Poisson driven time series of counts.

preprint2020arXiv

Conditional Independence in Max-linear Bayesian Networks

Motivated by extreme value theory, max-linear Bayesian networks have been recently introduced and studied as an alternative to linear structural equation models. However, for max-linear systems the classical independence results for Bayesian networks are far from exhausting valid conditional independence statements. We use tropical linear algebra to derive a compact representation of the conditional distribution given a partial observation, and exploit this to obtain a complete description of all conditional independence relations. In the context-specific case, where conditional independence is queried relative to a specific value of the conditioning variables, we introduce the notion of a source DAG to disclose the valid conditional independence relations. In the context-free case we characterize conditional independence through a modified separation concept, $\ast$-separation, combined with a tropical eigenvalue condition. We also introduce the notion of an impact graph which describes how extreme events spread deterministically through the network and we give a complete characterization of such impact graphs. Our analysis opens up several interesting questions concerning conditional independence and tropical geometry.

preprint2020arXiv

Recursive max-linear models with propagating noise

Recursive max-linear vectors model causal dependence between node variables by a structural equation model, expressing each node variable as a max-linear function of its parental nodes in a directed acyclic graph (DAG) and some exogenous innovation. For such a model, there exists a unique minimum DAG, represented by the Kleene star matrix of its edge weight matrix, which identifies the model and can be estimated. For a more realistic statistical modeling we introduce some random observational noise. A probabilistic analysis of this new noisy model reveals that the unique minimum DAG representing the distribution of the non-noisy model remains unchanged and identifiable. Moreover, the distribution of the minimum ratio estimators of the model parameters at their left limits are completely determined by the distribution of the noise variables up to a positive constant. Under a regular variation condition on the noise variables we prove that the estimated Kleene star matrix converges to a matrix of independent Weibull entries after proper centering and scaling.

preprint2020arXiv

Ruin probabilities for risk processes in a bipartite network

This paper studies risk balancing features in an insurance market by evaluating ruin probabilities for single and multiple components of a multivariate compound Poisson risk process. The dependence of the components of the process is induced by a random bipartite network. In analogy with the non-network scenario, a network ruin parameter is introduced. This random parameter, which depends on the bipartite network, is crucial for the ruin probabilities. Under certain conditions on the network and for light-tailed claim size distributions we obtain Lundberg bounds and, for exponential claim size distributions, exact results for the ruin probabilities. For large sparse networks, the network ruin parameter is approximated by a function of independent Poisson variables. T

preprint2020arXiv

Tail probabilities of random linear functions of regularly varying random vectors

We provide a new extension of Breiman's Theorem on computing tail probabilities of a product of random variables to a multivariate setting. In particular, we give a complete characterization of regular variation on cones in $[0,\infty)^d$ under random linear transformations. This allows us to compute probabilities of a variety of tail events, which classical multivariate regularly varying models would report to be asymptotically negligible. We illustrate our findings with applications to risk assessment in financial systems and reinsurance markets under a bipartite network structure.

preprint2017arXiv

Partial mean field limits in heterogeneous networks

We investigate systems of interacting stochastic differential equations with two kinds of heterogeneity: one originating from different weights of the linkages, and one concerning their asymptotic relevance when the system becomes large. To capture these effects we define a partial mean field system, and prove a law of large numbers with explicit bounds on the mean squared error. Furthermore, a large deviation result is established under reasonable assumptions. The theory will be illustrated by several examples: on the one hand, we recover the classical results of chaos propagation for homogeneous systems, and on the other hand, we demonstrate the validity of our assumptions for quite general heterogeneous networks including those arising from preferential attachment random graph models.

preprint2016arXiv

Anisotropic Brown-Resnick space-time processes: estimation and model assessment

Spatially isotropic max-stable processes have been used to model extreme spatial or space-time observations. One prominent model is the Brown-Resnick process, which has been successfully fitted to time series, spatial data and space-time data. This paper extends the process to possibly anisotropic spatial structures. For regular grid observations we prove strong consistency and asymptotic normality of pairwise maximum likelihood estimates for fixed and increasing spatial domain, when the number of observations in time tends to infinity. We also present a statistical test for isotropy versus anisotropy. We apply our test to precipitation data in Florida, and present some diagnostic tools for model assessment. Finally, we present a method to predict conditional probability fields and apply it to the data.

preprint2016arXiv

Conditional loss probabilities for systems of economic agents sharing light-tailed claims with analysis of portfolio diversification benefits

We analyze systems of agents sharing light-tailed risky claims issued by different financial objects. Assuming exponentially distributed claims, we obtain that both agents' and system's losses follow generalized exponential mixture distributions. We show that this leads to qualitatively different results on individual and system risks compared to heavy-tailed claims previously studied in the literature. By deducing conditional loss distributions we investigate the impact of stress situations on agents' and system's losses. Moreover, we present a criterion for agents to decide whether holding few objects or portfolio diversification minimizes their risks in system crisis situations.

preprint2016arXiv

Passage time and fluctuation calculations for subexponential Lévy processes

We consider the passage time problem for Lévy processes, emphasising heavy tailed cases. Results are obtained under quite mild assumptions, namely, drift to $-\infty$ a.s. of the process, possibly at a linear rate (the finite mean case), but possibly much faster (the infinite mean case), together with subexponential growth on the positive side. Local and functional versions of limit distributions are derived for the passage time itself, as well as for the position of the process just prior to passage, and the overshoot of a high level. A significant connection is made with extreme value theory via regular variation or maximum domain of attraction conditions imposed on the positive tail of the canonical measure, which are shown to be necessary for the kind of convergence behaviour we are interested in.

preprint2016arXiv

Time-consistency of risk measures with GARCH volatilities and their estimation

In this paper we study time-consistent risk measures for returns that are given by a GARCH(1,1) model. We present a construction of risk measures based on their static counterparts that overcomes the lack of time-consistency. We then study in detail our construction for the risk measures Value-at-Risk (VaR) and Average Value-at-Risk (AVaR). While in the VaR case we can derive an analytical formula for its time-consistent counterpart, in the AVaR case we derive lower and upper bounds to its time-consistent version. Furthermore, we incorporate techniques from Extreme Value Theory (EVT) to allow for a more tail-geared statistical analysis of the corresponding risk measures. We conclude with an application of our results to a data set of stock prices.

preprint2015arXiv

Conditional risk measures in a bipartite market structure

In this paper we study the effect of network structure between agents and objects on measures for systemic risk. We model the influence of sharing large exogeneous losses to the financial or (re)insuance market by a bipartite graph. Using Pareto-tailed losses and multivariate regular variation we obtain asymptotic results for systemic conditional risk measures based on the Value-at-Risk and the Conditional Tail Expectation. These results allow us to assess the influence of an individual institution on the systemic or market risk and vice versa through a collection of conditional systemic risk measures. For large markets Poisson approximations of the relevant constants are provided in the example of an insurance market. The example of an underlying homogeneous random graph is analysed in detail, and the results are illustrated through simulations.

preprint2015arXiv

Integrability conditions for space-time stochastic integrals: Theory and applications

We derive explicit integrability conditions for stochastic integrals taken over time and space driven by a random measure. Our main tool is a canonical decomposition of a random measure which extends the results from the purely temporal case. We show that the characteristics of this decomposition can be chosen as predictable strict random measures, and we compute the characteristics of the stochastic integral process. We apply our conditions to a variety of examples, in particular to ambit processes, which represent a rich model class.

preprint2015arXiv

Simulation of stochastic Volterra equations driven by space--time Lévy noise

In this paper we investigate two numerical schemes for the simulation of stochastic Volterra equations driven by space--time Lévy noise of pure-jump type. The first one is based on truncating the small jumps of the noise, while the second one relies on series representation techniques for infinitely divisible random variables. Under reasonable assumptions, we prove for both methods $L^p$- and almost sure convergence of the approximations to the true solution of the Volterra equation. We give explicit convergence rates in terms of the Volterra kernel and the characteristics of the noise. A simulation study visualizes the most important path properties of the investigated processes.

preprint2014arXiv

Asymmetric COGARCH processes

Financial data are as a rule asymmetric, although most econometric models are symmetric. This applies also to continuous-time models for high-frequency and irregularly spaced data. We discuss some asymmetric versions of the continuous-time GARCH model, concentrating then on the GJR-COGARCH. We calculate higher order moments and extend the first jump approximation. These results are prerequisites for moment estimation and pseudo maximum likelihood estimation of the GJR-COGARCH parameters, respectively, which we derive in detail.

preprint2014arXiv

Superposition of COGARCH processes

We suggest three superpositions of COGARCH (supCOGARCH) volatility processes driven by Lévy processes or Lévy bases. We investigate second-order properties, jump behaviour, and prove that they exhibit Pareto-like tails. Corresponding price processes are defined and studied. We find that the supCOGARCH models allow for more flexible autocovariance structures than the COGARCH. Moreover, other than most financial volatility models, the supCOGARCH processes do not exhibit a deterministic relationship between price and volatility jumps. Furthermore, in one supCOGARCH model not all volatility jumps entail a price jump, while in another supCOGARCH model not all price jumps necessarily lead to volatility jumps.

preprint2014arXiv

Systemic risk through contagion in a core-periphery structured banking network

We contribute to the understanding of how systemic risk arises in a network of credit-interlinked agents. Motivated by empirical studies we formulate a network model which, despite its simplicity, depicts the nature of interbank markets better than a homogeneous model. The components of a vector Ornstein-Uhlenbeck process living on the vertices of the network describe the financial robustnesses of the agents. For this system, we prove a LLN for growing network size leading to a propagation of chaos result. We state properties, which arise from such a structure, and examine the effect of inhomogeneity on several risk management issues and the possibility of contagion.

preprint2013arXiv

High-frequency sampling and kernel estimation for continuous-time moving average processes

Interest in continuous-time processes has increased rapidly in recent years, largely because of high-frequency data available in many applications. We develop a method for estimating the kernel function $g$ of a second-order stationary Lévy-driven continuous-time moving average (CMA) process $Y$ based on observations of the discrete-time process $Y^Δ$ obtained by sampling $Y$ at $Δ, 2Δ,...,nΔ$ for small $Δ$. We approximate $g$ by $g^Δ$ based on the Wold representation and prove its pointwise convergence to $g$ as $Δ\rightarrow 0$ for $\CARMA(p,q)$ processes. Two non-parametric estimators of $g^Δ$, based on the innovations algorithm and the Durbin-Levinson algorithm, are proposed to estimate $g$. For a Gaussian CARMA process we give conditions on the sample size $n$ and the grid-spacing $Δ(n)$ under which the innovations estimator is consistent and asymptotically normal as $n\rightarrow\infty$. The estimators can be calculated from sampled observations of {\it any} CMA process and simulations suggest that they perform well even outside the class of CARMA processes. We illustrate their performance for simulated data and apply them to the Brookhaven turbulent wind speed data. Finally we extend results of \citet{bfk:2011:1} for sampled CARMA processes to a much wider class of CMA processes.

preprint2013arXiv

N-density representability and the optimal transport limit of the Hohenberg-Kohn functional

We derive and analyze a hierarchy of approximations to the strongly correlated limit of the Hohenberg-Kohn functional. These "density representability approximations" are obtained by first noting that in the strongly correlated limit, N-representability of the pair density reduces to the requirement that the pair density must come from a symmetric N-point density. One then relaxes this requirement to the existence of a representing symmetric k-point density with k < N. The approximate energy can be computed by simulating a fictitious k-electron system. We investigate the approximations by deriving analytically exact results for a 2-site model problem, and by incorporating them into a self-consistent Kohn-Sham calculation for small atoms. We find that the low order representability conditions already capture the main part of the correlations.

preprint2013arXiv

Two-step estimation of a multivariate Lévy process

Based on the concept of a Lévy copula to describe the dependence structure of a multivariate Lévy process we present a new estimation procedure. We consider a parametric model for the marginal Lévy processes as well as for the Lévy copula and estimate the parameters by a two-step procedure. We first estimate the parameters of the marginal processes, and then estimate in a second step only the dependence structure parameter. For infinite Lévy measures we truncate the small jumps and base our statistical analysis on the large jumps of the model. Prominent example will be a bivariate stable \lp, which allows for analytic calculations and, hence, for a comparison of different methods. We prove asymptotic normality of the parameter estimates from the two-step procedure and, in particular, we derive the Godambe information matrix, whose inverse is the covariance matrix of the normal limit law. A simulation study investigates the loss of efficiency because of the two-step procedure and the truncation.

preprint2012arXiv

Futures pricing in electricity markets based on stable CARMA spot models

We present a new model for the electricity spot price dynamics, which is able to capture seasonality, low-frequency dynamics and the extreme spikes in the market. Instead of the usual purely deterministic trend we introduce a non-stationary independent increments process for the low-frequency dynamics, and model the large fluctuations by a non-Gaussian stable CARMA process. The model allows for analytic futures prices, and we apply these to model and estimate the whole market consistently. Besides standard parameter estimation, an estimation procedure is suggested, where we fit the non-stationary trend using futures data with long time until delivery, and a robust $L^1$-filter to find the states of the CARMA process. The procedure also involves the empirical and theoretical risk premiums which -- as a by-product -- are also estimated. We apply this procedure to data from the German electricity exchange EEX, where we split the empirical analysis into base load and peak load prices. We find an overall negative risk premium for the base load futures contracts, except for contracts close to delivery, where a small positive risk premium is detected. The peak load contracts, on the other hand, show a clear positive risk premium, when they are close to delivery, while the contracts in the longer end also have a negative premium.

preprint2012arXiv

Statistical inference for max-stable processes in space and time

Max-stable processes have proved to be useful for the statistical modelling of spatial extremes. Several representations of max-stable random fields have been proposed in the literature. One such representation is based on a limit of normalized and scaled pointwise maxima of stationary Gaussian processes that was first introduced by Kabluchko, Schlather and de Haan (2009). This paper deals with statistical inference for max-stable space-time processes that are defined in an analogous fashion. We describe pairwise likelihood estimation, where the pairwise density of the process is used to estimate the model parameters and prove strong consistency and asymptotic normality of the parameter estimates for an increasing space-time dimension, i.e., as the joint number of spatial locations and time points tends to infinity. A simulation study shows that the proposed method works well for these models.

preprint2012arXiv

Turbulence modeling by time-series methods

A general model for stationary, time-wise turbulent velocity is presented and discussed. This approach, inspired by modeling ideas of Barndorff-Nielsen and Schimgel, is coherent with the K41 hypothesis of local isotropy, and it allows us to separate second-order statistics from higher order ones. The model can be motivated by Taylor's hypothesis and a relation between time and spatial spectra. Second order statistics are used to separate the deterministic kernel function and the weakly stationary driving noise. A non-parametric estimation method for the turbulence intermittency is suggested.

preprint2011arXiv

Density functional theory and optimal transportation with Coulomb cost

We present here novel insight into exchange-correlation functionals in density functional theory, based on the viewpoint of optimal transport. We show that in the case of two electrons and in the semiclassical limit, the exact exchange-correlation functional reduces to a very interesting functional of novel form, which depends on an optimal transport map $T$ associated with a given density $ρ$. Since the above limit is strongly correlated, the limit functional yields insight into electron correlations. We prove the existence and uniqueness of such an optimal map for any number of electrons and each $ρ$, and determine the map explicitly in the case when $ρ$ is radially symmetric.

preprint2011arXiv

Fractional Lévy-driven Ornstein--Uhlenbeck processes and stochastic differential equations

Using Riemann-Stieltjes methods for integrators of bounded $p$-variation we define a pathwise integral driven by a fractional Lévy process (FLP). To explicitly solve general fractional stochastic differential equations (SDEs) we introduce an Ornstein-Uhlenbeck model by a stochastic integral representation, where the driving stochastic process is an FLP. To achieve the convergence of improper integrals, the long-time behavior of FLPs is derived. This is sufficient to define the fractional Lévy-Ornstein-Uhlenbeck process (FLOUP) pathwise as an improper Riemann-Stieltjes integral. We show further that the FLOUP is the unique stationary solution of the corresponding Langevin equation. Furthermore, we calculate the autocovariance function and prove that its increments exhibit long-range dependence. Exploiting the Langevin equation, we consider SDEs driven by FLPs of bounded $p$-variation for $p<2$ and construct solutions using the corresponding FLOUP. Finally, we consider examples of such SDEs, including various state space transforms of the FLOUP and also fractional Lévy-driven Cox-Ingersoll-Ross (CIR) models.

preprint2011arXiv

High frequency sampling of a continuous-time ARMA process

Continuous-time autoregressive moving average (CARMA) processes have recently been used widely in the modeling of non-uniformly spaced data and as a tool for dealing with high-frequency data of the form $Y_{nΔ}, n=0,1,2,...$, where $Δ$ is small and positive. Such data occur in many fields of application, particularly in finance and the study of turbulence. This paper is concerned with the characteristics of the process $(Y_{nΔ})_{n\in\bbz}$, when $Δ$ is small and the underlying continuous-time process $(Y_t)_{t\in\bbr}$ is a specified CARMA process.

preprint2011arXiv

Max-stable processes for modelling extremes observed in space and time

Max-stable processes have proved to be useful for the statistical modelling of spatial extremes. Several representations of max-stable random fields have been proposed in the literature. For statistical inference it is often assumed that there is no temporal dependence, i.e., the observations at spatial locations are independent in time. We use two representations of stationary max-stable spatial random fields and extend the concepts to the space-time domain. In a first approach, we extend the idea of constructing max-stable random fields as limits of normalized and rescaled pointwise maxima of independent Gaussian random fields, which was introduced by Kabluchko, Schlather and de Haan [2009], who construct max-stable random fields associated to a class of variograms. We use a similar approach based on a well-known result by Hüsler and Reiss and apply specific spatio-temporal covariance models for the underlying Gaussian random field, which satisfy weak regularity assumptions. Furthermore, we extend Smith's storm profile model to a space-time setting and provide explicit expressions for the bivariate distribution functions. The tail dependence coefficient is an important measure of extremal dependence. We show how the spatio-temporal covariance function underlying the Gaussian random field can be interpreted in terms of the tail dependence coefficient. Within this context, we examine different concepts for constructing spatio-temporal covariance models and analyse several specific examples, including Gneiting's class of nonseparable stationary covariance functions.

preprint2011arXiv

On the Ruin Probability of the Generalised Ornstein-Uhlenbeck Process in the Cramér Case

For a bivariate \Levy process $(ξ_t,η_t)_{t\ge 0}$ and initial value $V_0$ define the Generalised Ornstein-Uhlenbeck (GOU) process \[ V_t:=e^{ξ_t}\Big(V_0+\int_0^t e^{-ξ_{s-}}\ud η_s\Big),\quad t\ge0,\] and the associated stochastic integral process \[Z_t:=\int_0^t e^{-ξ_{s-}}\ud η_s,\quad t\ge0.\] Let $T_z:=\inf\{t>0:V_t<0\mid V_0=z\}$ and $ψ(z):=P(T_z<\infty)$ for $z\ge 0$ be the ruin time and infinite horizon ruin probability of the GOU. Our results extend previous work of Nyrhinen (2001) and others to give asymptotic estimates for $ψ(z)$ and the distribution of $T_z$ as $z\to\infty$, under very general, easily checkable, assumptions, when $ξ$ satisfies a Cramér condition.