Source author record

Martin Schlather

Martin Schlather 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

20works
8topics
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

20 published item(s)

preprint2022arXiv

Covariance Models for Multivariate Random Fields resulting from Pseudo Cross-Variograms

So far, the pseudo cross-variogram is primarily used as a tool for the structural analysis of multivariate random fields. Mainly applying recent theoretical results on the pseudo cross-variogram, we use it as a cornerstone in the construction of valid covariance models for multivariate random fields. In particular, we extend known univariate constructions to the multivariate case, and generalize existing multivariate models. Furthermore, we provide a general construction principle for conditionally negative definite matrix-valued kernels, which we use to reinterpret previous modeling proposals.

preprint2016arXiv

The realization problem for tail correlation functions

For a stochastic process $\{X_t\}_{t \in T}$ with identical one-dimensional margins and upper endpoint $τ_{\text{up}}$ its tail correlation function (TCF) is defined through $χ^{(X)}(s,t) = \lim_{τ\to τ_{\text{up}}} P(X_s > τ\,\mid\, X_t > τ)$. It is a popular bivariate summary measure that has been frequently used in the literature in order to assess tail dependence. In this article, we study its realization problem. We show that the set of all TCFs on $T \times T$ coincides with the set of TCFs stemming from a subclass of max-stable processes and can be completely characterized by a system of affine inequalities. Basic closure properties of the set of TCFs and regularity implications of the continuity of $χ$ are derived. If $T$ is finite, the set of TCFs on $T \times T$ forms a convex polytope of $\lvert T \rvert \times \lvert T \rvert$ matrices. Several general results reveal its complex geometric structure. Up to $\lvert T \rvert = 6$ a reduced system of necessary and sufficient conditions for being a TCF is determined. None of these conditions will become obsolete as $\lvert T \rvert\geq 3$ grows.

preprint2015arXiv

An exceptional max-stable process fully parameterized by its extremal coefficients

The extremal coefficient function (ECF) of a max-stable process $X$ on some index set $T$ assigns to each finite subset $A\subset T$ the effective number of independent random variables among the collection $\{X_t\}_{t\in A}$. We introduce the class of Tawn-Molchanov processes that is in a 1:1 correspondence with the class of ECFs, thus also proving a complete characterization of the ECF in terms of negative definiteness. The corresponding Tawn-Molchanov process turns out to be exceptional among all max-stable processes sharing the same ECF in that its dependency set is maximal w.r.t. inclusion. This entails sharp lower bounds for the finite dimensional distributions of arbitrary max-stable processes in terms of its ECF. A spectral representation of the Tawn-Molchanov process and stochastic continuity are discussed. We also show how to build new valid ECFs from given ECFs by means of Bernstein functions.

preprint2015arXiv

Maxima of independent, non-identically distributed Gaussian vectors

Let $X_{i,n},n\in \mathbb{N},1\leq i\leq n$, be a triangular array of independent $\mathbb{R}^d$-valued Gaussian random vectors with correlation matrices $Σ_{i,n}$. We give necessary conditions under which the row-wise maxima converge to some max-stable distribution which generalizes the class of Hüsler-Reiss distributions. In the bivariate case, the conditions will also be sufficient. Using these results, new models for bivariate extremes are derived explicitly. Moreover, we define a new class of stationary, max-stable processes as max-mixtures of Brown-Resnick processes. As an application, we show that these processes realize a large set of extremal correlation functions, a natural dependence measure for max-stable processes. This set includes all functions $ψ(\sqrt{γ(h)}),h\in \mathbb{R}^d$, where $ψ$ is a completely monotone function and $γ$ is an arbitrary variogram.

preprint2015arXiv

Statistical Post-Processing of Forecasts for Extremes Using Bivariate Brown-Resnick Processes with an Application to Wind Gusts

To improve the forecasts of weather extremes, we propose a joint spatial model for the observations and the forecasts, based on a bivariate Brown-Resnick process. As the class of stationary bivariate Brown-Resnick processes is fully characterized by the class of pseudo cross-variograms, we contribute to the theorical understanding of pseudo cross-variograms refining the knowledge of the asymptotic behaviour of all their components and introducing a parsimonious, but flexible parametric model. Both findings are of interest in classical geostatistics on their own. The proposed model is applied to real observation and forecast data for extreme wind gusts at 119 stations in Northern Germany.

preprint2014arXiv

A parametric variogram model bridging between stationary and intrinsically stationary processes

A simple variogram model with two parameters is presented that includes the power variogram for the fractional Brownian motion, a modified De Wijsian model, the generalized Cauchy model and the multiquadrics model. One parameter controls the smoothness of the process. The other parameter allows for a smooth parametrization between stationary and intrinsically stationary second order processes in a Gaussian framework, or between mixing and non-ergodic max-stable processes when modeling spatial extremes by a Brown-Resnick process.

preprint2014arXiv

Conditional Sampling for Max-Stable Processes with a Mixed Moving Maxima Representation

This paper deals with the question of conditional sampling and prediction for the class of stationary max-stable processes which allow for a mixed moving maxima representation. We develop an exact procedure for conditional sampling using the Poisson point process structure of such processes. For explicit calculations we restrict ourselves to the one-dimensional case and use a finite number of shape functions satisfying some regularity conditions. For more general shape functions approximation techniques are presented. Our algorithm is applied to the Smith process and the Brown-Resnick process. Finally, we compare our computational results to other approaches. Here, the algorithm for Gaussian processes with transformed marginals turns out to be surprisingly competitive.

preprint2013arXiv

Joint Extremal Behavior of Hidden and Observable Time Series with an Application to GARCH Processes

We study the behavior of a real-valued and unobservable process (Y_t) under an extreme event of a related process (X_t) that is observable. Our analysis is motivated by the well-known GARCH model which represents two such sequences, i.e. the observable log returns of an asset as well as the hidden volatility process. Our results complement the findings of Segers (2007) and Smith (1992) for a single time series. We show that under suitable assumptions their concept of a tail chain as a limiting process is also applicable to our setting. Furthermore, we discuss existence and uniqueness of a limiting process under some weaker assumptions. Finally, we apply our results to the GARCH(1,1) case.

preprint2013arXiv

On the Normalized Spectral Representation of Max-Stable Processes on a Compact Set

The normalized spectral representation of a max-stable process on a compact set is the unique representation where all spectral functions share the same supremum. Among the class of equivalent spectral representations of a process, the normalized spectral representation plays a distinctive role as a solution of two optimization problems in the context of an efficient simulation of max-stable processes. Our approach has the potential of considerably reducing the simulation time of max-stable processes.

preprint2013arXiv

The Expected Linkage Disequilibrium in Finite Populations Revisited

The expected level of linkage disequilibrium (LD) in a finite ideal population at equilibrium is of relevance for many applications in population and quantitative genetics. Several recursion formulae have been proposed during the last decades, whose derivations mostly contain heuristic parts and therefore remain mathematically questionable. We propose a more justifiable approach, including an alternative recursion formula for the expected LD. Since the exact formula depends on the distribution of allele frequencies in a very complicated manner, we suggest an approximate solution and analyze its validity extensively in a simulation study. Compared to the widely used formula of Sved, the proposed formula performs better for all parameter constellations considered. We then analyze the expected LD at equilibrium using the theory on discrete-time Markov chains based on the linear recursion formula, with equilibrium being defined as the steady-state of the chain, which finally leads to a formula for the effective population size N_e. An additional analysis considers the effect of non-exactness of a recursion formula on the steady-state, demonstrating that the resulting error in expected LD can be substantial. In an application to the HapMap data of two human populations we illustrate the dependency of the N_e-estimate on the distribution of minor allele frequencies (MAFs), showing that the estimated N_e can vary by up to 30% when a uniform instead of a skewed distribution of MAFs is taken as a basis to select SNPs for the analyses. Our analyses provide new insights into the mathematical complexity of the problem studied.

preprint2012arXiv

Characterizing extremal coefficient functions and extremal correlation functions

We focus on two dependency quantities of a max-stable random field $X$ on some space $T$: the extremal coefficient function $θ$ which we define on finite sets of $T$ and the extremal correlation function $χ(s,t)=\lim_{x \uparrow \infty} \PP(X_s \geq x \mid X_t \geq x)$. We fully characterize extremal coefficient functions $θ$ by a property called complete alternation and construct a corresponding max-stable random field. Simple properties and consequences concerning the convex geometry of extremal coefficients are derived. We study how the continuity of $X$, $θ$ and $χ$ are linked to each other, and we show that extremal correlation functions $χ$ allow for convex combinations in general, and for products and pointwise limits if the resulting function is continuous. These are operations which are well-known for positive definite functions, but the latter are non-trivial for extremal correlation functions. Finally, we regard some additional implications, when the random field $X$ on $T=\mathbb{R}^d$ is stationary.

preprint2012arXiv

Estimation of Huesler-Reiss distributions and Brown-Resnick processes

Estimation of extreme-value parameters from observations in the max-domain of attraction (MDA) of a multivariate max-stable distribution commonly uses aggregated data such as block maxima. Since we expect that additional information is contained in the non-aggregated, single "large" observations, we introduce a new approach of inference based on a multivariate peaks-over-threshold method. We show that for any process in the MDA of the frequently used Hüsler-Reiss model or its spatial extension, the Brown-Resnick process, suitably defined conditional increments asymptotically follow a multivariate Gaussian distribution. This leads to computationally efficient estimates of the Hüsler-Reiss parameter matrix. Further, the results enable parametric inference for Brown-Resnick processes. A simulation study compares the performance of the new estimators to other commonly used methods. As an application, we fit a non-isotropic Brown-Resnick process to the extremes of 12 year data of daily wind speed measurements.

preprint2012arXiv

Intrinsically Weighted Means of Marked Point Processes

For a non-stationary or non-ergodic marked point process (MPP) on $\R^d$, the definition of averages becomes ambiguous as the process might have a different stochastic behavior in different realizations (non-ergodicity) or in different areas of the observation window (non-stationarity). We investigate different definitions for the moments, including a new hierarchical definition for non-ergodic MPPs, and embed them into a family of weighted mean marks. We point out examples of application in which different weighted mean marks all have a sensible meaning. Further, asymptotic properties of the corresponding estimators are investigated as well as optimal weighting procedures.

preprint2012arXiv

Random Marked Sets

We aim to link random fields and marked point processes and therefore introduce a new class of stochastic processes which are defined on a random set in R^d. Unlike for random fields, the mark covariance function of a marked random set is in general not positive definite. This implies that in many situations the use of simple geostatistical methods appears to be questionable. Surprisingly, for a special class of processes based on Gaussian random fields, we do have positive definiteness for the corresponding mark covariance function and mark correlation function.

preprint2012arXiv

Representations of max-stable processes based on single extreme events

This paper provides the basis for new methods of inference for max-stable processes ξ on general spaces that admit a certain incremental representation, which, in important cases, has a much simpler structure than the max-stable process itself. A corresponding peaks-over-threshold approach will incorporate all single events that are extreme in some sense and will therefore rely on a substantially larger amount of data in comparison to estimation procedures based on block maxima. Conditioning a process η in the max-domain of attraction of ξ on being extremal, several convergence results for the increments of η are proved. In a similar way, the shape functions of mixed moving maxima (M3) processes can be extracted from suitably conditioned single events η. Connecting the two approaches, transformation formulae for processes that admit both an incremental and an M3 representation are identified.

preprint2011arXiv

Some covariance models based on normal scale mixtures

Modelling spatio-temporal processes has become an important issue in current research. Since Gaussian processes are essentially determined by their second order structure, broad classes of covariance functions are of interest. Here, a new class is described that merges and generalizes various models presented in the literature, in particular models in Gneiting (J. Amer. Statist. Assoc. 97 (2002) 590--600) and Stein (Nonstationary spatial covariance functions (2005) Univ. Chicago). Furthermore, new models and a multivariate extension are introduced.

preprint2001arXiv

Stochastic models which separate fractal dimension and Hurst effect

Fractal behavior and long-range dependence have been observed in an astonishing number of physical systems. Either phenomenon has been modeled by self-similar random functions, thereby implying a linear relationship between fractal dimension, a measure of roughness, and Hurst coefficient, a measure of long-memory dependence. This letter introduces simple stochastic models which allow for any combination of fractal dimension and Hurst exponent. We synthesize images from these models, with arbitrary fractal properties and power-law correlations, and propose a test for self-similarity.