Researcher profile

Marius Hofert

Marius Hofert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
10topics
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

8 published item(s)

preprint2023arXiv

Matrix compatibility and correlation mixture representation of generalized Gini's gamma

Representations of measures of concordance in terms of Pearson' s correlation coefficient are studied. All transforms of random variables are characterized such that the correlation coefficient of the transformed random variables is a measure of concordance. Next, Gini' s gamma is generalized and it is shown that the resulting generalized Gini' s gamma can be represented as a mixture of measures of concordance that are Pearson' s correlation coefficients of transformed random variables. As an application of this correlation mixture representation of generalized Gini' s gamma, lower and upper bounds of the compatible set of generalized Gini' s gamma, which is the collection of all possible square matrices whose entries are pairwise bivariate generalized Gini' s gammas, are derived.

preprint2022arXiv

Smooth bootstrapping of copula functionals

The smooth bootstrap for estimating copula functionals in small samples is investigated. It can be used both to gauge the distribution of the estimator in question and to augment the data. Issues arising from kernel density and distribution estimation in the copula domain are addressed, such as how to avoid the bounded domain, which bandwidth matrix to choose, and how the smoothing can be carried out. Furthermore, we investigate how the smooth bootstrap impacts the underlying dependence structure or the functionals in question and under which conditions it does not. We provide specific examples and simulations that highlight advantages and caveats of the approach.

preprint2020arXiv

Markov Chain Monte Carlo Methods for Estimating Systemic Risk Allocations

We propose a novel framework of estimating systemic risk measures and risk allocations based on Markov chain Monte Carlo (MCMC) methods. We consider a class of allocations whose jth component can be written as some risk measure of the jth conditional marginal loss distribution given the so-called crisis event. By considering a crisis event as an intersection of linear constraints, this class of allocations covers, for example, conditional Value-at-Risk (CoVaR), conditional expected shortfall (CoES), VaR contributions, and range VaR (RVaR) contributions as special cases. For this class of allocations, analytical calculations are rarely available, and numerical computations based on Monte Carlo (MC) methods often provide inefficient estimates due to the rare-event character of the crisis events. We propose an MCMC estimator constructed from a sample path of a Markov chain whose stationary distribution is the conditional distribution given the crisis event. Efficient constructions of Markov chains, such as Hamiltonian Monte Carlo and Gibbs sampler, are suggested and studied depending on the crisis event and the underlying loss distribution. The efficiency of the MCMC estimators is demonstrated in a series of numerical experiments.

preprint2020arXiv

Normal variance mixtures: Distribution, density and parameter estimation

Normal variance mixtures are a class of multivariate distributions that generalize the multivariate normal by randomizing (or mixing) the covariance matrix via multiplication by a non-negative random variable W. The multivariate t distribution is an example of such mixture, where W has an inverse-gamma distribution. Algorithms to compute the joint distribution function and perform parameter estimation for the multivariate normal and t (with integer degrees of freedom) can be found in the literature and are implemented in, e.g., the R package mvtnorm. In this paper, efficient algorithms to perform these tasks in the general case of a normal variance mixture are proposed. In addition to the above two tasks, the evaluation of the joint (logarithmic) density function of a general normal variance mixture is tackled as well, as it is needed for parameter estimation and does not always exist in closed form in this more general setup. For the evaluation of the joint distribution function, the proposed algorithms apply randomized quasi-Monte Carlo (RQMC) methods in a way that improves upon existing methods proposed for the multivariate normal and t distributions. An adaptive RQMC algorithm that similarly exploits the superior convergence properties of RQMC methods is presented for the task of evaluating the joint log-density function. This allows the parameter estimation task to be accomplished via an EM-like algorithm where all weights and log-densities are numerically estimated. It is demonstrated through numerical examples that the suggested algorithms are quite fast; even for high dimensions around 1000 the distribution function can be estimated with moderate accuracy using only a few seconds of run time. Even log-densities around -100 can be estimated accurately and quickly. An implementation of all algorithms presented in this work is available in the R package nvmix (version >= 0.0.4).

preprint2020arXiv

Outer power transformations of hierarchical Archimedean copulas: Construction, sampling and estimation

A large number of commonly used parametric Archimedean copula (AC) families are restricted to a single parameter, connected to a concordance measure such as Kendall's tau. This often leads to poor statistical fits, particularly in the joint tails, and can sometimes even limit the ability to model concordance or tail dependence mathematically. This work suggests outer power (OP) transformations of Archimedean generators to overcome these limitations. The copulas generated by OP-transformed generators can, for example, allow one to capture both a given concordance measure and a tail dependence coefficient simultaneously. For exchangeable OP-transformed ACs, a formula for computing tail dependence coefficients is obtained, as well as two feasible OP AC estimators are proposed and their properties studied by simulation. For hierarchical extensions of OP-transformed ACs, a new construction principle, efficient sampling and parameter estimation are addressed. By simulation, convergence rate and standard errors of the proposed estimator are studied. Excellent tail fitting capabilities of OP-transformed hierarchical AC models are demonstrated in a risk management application. The results show that the OP transformation is able to improve the statistical fit of exchangeable ACs, particularly of those that cannot capture upper tail dependence or strong concordance, as well as the statistical fit of hierarchical ACs, especially in terms of tail dependence and higher dimensions. Given how comparably simple it is to include OP transformations into existing exchangeable and hierarchical AC models, this transformation provides an attractive trade-off between computational effort and statistical improvement.

preprint2020arXiv

Quasi-random sampling for multivariate distributions via generative neural networks

Generative moment matching networks (GMMNs) are introduced for generating quasi-random samples from multivariate models with any underlying copula in order to compute estimates under variance reduction. So far, quasi-random sampling for multivariate distributions required a careful design, exploiting specific properties (such as conditional distributions) of the implied parametric copula or the underlying quasi-Monte Carlo (QMC) point set, and was only tractable for a small number of models. Utilizing GMMNs allows one to construct quasi-random samples for a much larger variety of multivariate distributions without such restrictions, including empirical ones from real data with dependence structures not well captured by parametric copulas. Once trained on pseudo-random samples from a parametric model or on real data, these neural networks only require a multivariate standard uniform randomized QMC point set as input and are thus fast in estimating expectations of interest under dependence with variance reduction. Numerical examples are considered to demonstrate the approach, including applications inspired by risk management practice. All results are reproducible with the demos GMMN_QMC_paper, GMMN_QMC_data and GMMN_QMC_timings as part of the R package gnn.

preprint2020arXiv

Random number generators produce collisions: Why, how many and more

It seems surprising that when applying widely used random number generators to generate one million random numbers on modern architectures, one obtains, on average, about 116 collisions. This article explains why, how to mathematically compute such a number, why they often cannot be obtained in a straightforward way, how to numerically compute them in a robust way and, among other things, what would need to be changed to bring this number below 1. The probability of at least one collision is also briefly addressed, which, as it turns out, again needs a careful numerical treatment. Overall, the article provides an introduction to the representation of floating-point numbers on a computer and corresponding implications in statistics and simulation. All computations are carried out in R and are reproducible with the code included in this article.

preprint2020arXiv

Right-truncated Archimedean and related copulas

The copulas of random vectors with standard uniform univariate margins truncated from the right are considered and a general formula for such right-truncated conditional copulas is derived. This formula is analytical for copulas that can be inverted analytically as functions of each single argument. This is the case, for example, for Archimedean and related copulas. The resulting right-truncated Archimedean copulas are not only analytically tractable but can also be characterized as tilted Archimedean copulas. This finding allows one, for example, to more easily derive analytical properties such as the coefficients of tail dependence or sampling procedures of right-truncated Archimedean copulas. As another result, one can easily obtain a limiting Clayton copula for a general vector of truncation points converging to zero; this is an important property for (re)insurance and a fact already known in the special case of equal truncation points, but harder to prove without aforementioned characterization. Furthermore, right-truncated Archimax copulas with logistic stable tail dependence functions are characterized as tilted outer power Archimedean copulas and an analytical form of right-truncated nested Archimedean copulas is also derived.