Researcher profile

Martin Schlather

Martin Schlather contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 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.

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.