BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Marius Hofert

Marius Hofert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
ComputationMethodologymath.STStatistics Theorymath.PRq-fin.RMmath.NAstat.OTApplicationsMachine LearningNumerical Analysisq-fin.CPq-fin.PR

Catalog footprint

What is connected

18works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bernoulli and tail-dependence compatibility

The tail-dependence compatibility problem is introduced. It raises the question whether a given $d\times d$-matrix of entries in the unit interval is the matrix of pairwise tail-dependence coefficients of a $d$-dimensional random vector. The problem is studied together with Bernoulli-compatible matrices, that is, matrices which are expectations of outer products of random vectors with Bernoulli margins. We show that a square matrix with diagonal entries being 1 is a tail-dependence matrix if and only if it is a Bernoulli-compatible matrix multiplied by a constant. We introduce new copula models to construct tail-dependence matrices, including commonly used matrices in statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On structure, family and parameter estimation of hierarchical Archimedean copulas

Research on structure determination and parameter estimation of hierarchical Archimedean copulas (HACs) has so far mostly focused on the case in which all appearing Archimedean copulas belong to the same Archimedean family. The present work addresses this issue and proposes a new approach for estimating HACs that involve different Archimedean families. It is based on employing goodness-of-fit test statistics directly into HAC estimation. The approach is summarized in a simple algorithm, its theoretical justification is given and its applicability is illustrated by several experiments, which include estimation of HACs involving up to five different Archimedean families.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quasi-random numbers for copula models

The present work addresses the question how sampling algorithms for commonly applied copula models can be adapted to account for quasi-random numbers. Besides sampling methods such as the conditional distribution method (based on a one-to-one transformation), it is also shown that typically faster sampling methods (based on stochastic representations) can be used to improve upon classical Monte Carlo methods when pseudo-random number generators are replaced by quasi-random number generators. This opens the door to quasi-random numbers for models well beyond independent margins or the multivariate normal distribution. Detailed examples (in the context of finance and insurance), illustrations and simulations are given and software has been developed and provided in the R packages copula and qrng.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An importance sampling approach for copula models in insurance

An importance sampling approach for sampling copula models is introduced. We propose two algorithms that improve Monte Carlo estimators when the functional of interest depends mainly on the behaviour of the underlying random vector when at least one of the components is large. Such problems often arise from dependence models in finance and insurance. The importance sampling framework we propose is general and can be easily implemented for all classes of copula models from which sampling is feasible. We show how the proposal distribution of the two algorithms can be optimized to reduce the sampling error. In a case study inspired by a typical multivariate insurance application, we obtain variance reduction factors between 10 and 30 in comparison to standard Monte Carlo estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Improved Algorithms for Computing Worst Value-at-Risk: Numerical Challenges and the Adaptive Rearrangement Algorithm

Numerical challenges inherent in algorithms for computing worst Value-at-Risk in homogeneous portfolios are identified and solutions as well as words of warning concerning their implementation are provided. Furthermore, both conceptual and computational improvements to the Rearrangement Algorithm for approximating worst Value-at-Risk for portfolios with arbitrary marginal loss distributions are given. In particular, a novel Adaptive Rearrangement Algorithm is introduced and investigated. These algorithms are implemented using the R package qrmtools.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Parallel and other simulations in R made easy: An end-to-end study

It is shown how to set up, conduct, and analyze large simulation studies with the new R package simsalapar = simulations simplified and launched parallel. A simulation study typically starts with determining a collection of input variables and their values on which the study depends, such as sample sizes, dimensions, types and degrees of dependence, estimation methods, etc. Computations are desired for all com- binations of these variables. If conducting these computations sequentially is too time- consuming, parallel computing can be applied over all combinations of select variables. The final result object of a simulation study is typically an array. From this array, sum- mary statistics can be derived and presented in terms of (flat contingency or LATEX) tables or visualized in terms of (matrix-like) figures. The R package simsalapar provides several tools to achieve the above tasks. Warnings and errors are dealt with correctly, various seeding methods are available, and run time is measured. Furthermore, tools for analyzing the results via tables or graphics are pro- vided. In contrast to rather minimal examples typically found in R packages or vignettes, an end-to-end, not-so-minimal simulation problem from the realm of quantitative risk management is given. The concepts presented and solutions provided by simsalapar may be of interest to students, researchers, and practitioners as a how-to for conducting real- istic, large-scale simulation studies in R. Also, the development of the package revealed useful improvements to R itself, which are available in R 3.0.0.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Densities of nested Archimedean copulas

Nested Archimedean copulas recently gained interest since they generalize the well-known class of Archimedean copulas to allow for partial asymmetry. Sampling algorithms and strategies have been well investigated for nested Archimedean copulas. However, for likelihood based inference it is important to have the density. The present work fills this gap. A general formula for the derivatives of the nodes and inner generators appearing in nested Archimedean copulas is developed. This leads to a tractable formula for the density of nested Archimedean copulas in arbitrary dimensions if the number of nesting levels is not too large. Various examples including famous Archimedean families and transformations of such are given. Furthermore, a numerically efficient way to evaluate the log-density is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Estimators for Archimedean copulas in high dimensions

The performance of known and new parametric estimators for Archimedean copulas is investigated, with special focus on large dimensions and numerical difficulties. In particular, method-of-moments-like estimators based on pairwise Kendall's tau, a multivariate extension of Blomqvist's beta, minimum distance estimators, the maximum-likelihood estimator, a simulated maximum-likelihood estimator, and a maximum-likelihood estimator based on the copula diagonal are studied. Their performance is compared in a large-scale simulation study both under known and unknown margins (pseudo-observations), in small and high dimensions, under small and large dependencies, various different Archimedean families and sample sizes. High dimensions up to one hundred are considered for the first time and computational problems arising from such large dimensions are addressed in detail. All methods are implemented in the open source \R{} package \pkg{copula} and can thus be easily accessed and studied.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Likelihood inference for Archimedean copulas

Explicit functional forms for the generator derivatives of well-known one-parameter Archimedean copulas are derived. These derivatives are essential for likelihood inference as they appear in the copula density, conditional distribution functions, or the Kendall distribution function. They are also required for several asymmetric extensions of Archimedean copulas such as Khoudraji-transformed Archimedean copulas. Access to the generator derivatives makes maximum-likelihood estimation for Archimedean copulas feasible in terms of both precision and run time, even in large dimensions. It is shown by simulation that the root mean squared error is decreasing in the dimension. This decrease is of the same order as the decrease in sample size. Furthermore, confidence intervals for the parameter vector are derived. Moreover, extensions to multi-parameter Archimedean families are given. All presented methods are implemented in the open-source R package nacopula and can thus easily be accessed and studied.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sibuya copulas

The standard intensity-based approach for modeling defaults is generalized by making the deterministic term structure of the survival probability stochastic via a common jump process. The survival copula of the vector of default times is derived and it is shown to be explicit and of the functional form as dealt with in the work of Sibuya. Besides the parameters of the jump process, the marginal survival functions of the default times appear in the copula. Sibuya copulas therefore allow for functional parameters and asymmetries. Due to the jump process in the construction, they allow for a singular component. Depending on the parameters, they may also be extreme-value copulas or Levy-frailty copulas. Further, Sibuya copulas are easy to sample in any dimension. Properties of Sibuya copulas including positive lower orthant dependence, tail dependence, and extremal dependence are investigated. An application to pricing first-to-default contracts is outlined and further generalizations of this copula class are addressed.

0 citations0 reviewsNo code linkedOpen access