Bayesian Copulae Distributions, with Application to Operational Risk Management

The aim of this paper is to introduce a new methodology for operational risk management, based on Bayesian copulae. One of the main problems related to operational risk management is understanding the complex dependence structure of the associated variables. In order to model this structure in a flexible way, we construct a method based on copulae. This allows us to split the joint multivariate probability distribution of a random vector of losses into individual components characterized by univariate marginals. Thus, copula functions embody all the information about the correlation between variables and provide a useful technique for modelling the dependency of a high number of marginals. Another important problem in operational risk modelling is the lack of loss data. This suggests the use of Bayesian models, computed via simulation methods and, in particular, Markov chain Monte Carlo. We propose a new methodology for modelling operational risk and for estimating the required capital. This methodology combines the use of copulae and Bayesian models.      

Copula-Type Estimators for Flexible Multivariate Density Modeling Using Mixtures

Copulas are popular as models for multivariate dependence because they allow the marginal densities and the joint dependence to be modeled separately. However, they usually require that the transformation from uniform marginals to the marginals of the joint dependence structure is known. This can only be done for a restricted set of copulas, for example, a normal copula. Our article introduces copula-type estimators for flexible multivariate density estimation which also allow the marginal densities to be modeled separately from the joint dependence, as in copula modeling, but overcomes the lack of flexibility of most popular copula estimators. An iterative scheme is proposed for estimating copula-type estimators and its usefulness is demonstrated through simulation and real examples. The joint dependence is modeled by mixture of normals and mixture of normal factor analyzer models, and mixture of t and mixture of t-factor analyzer models. We develop efficient variational Bayes algorithms for fitting these in which model selection is performed automatically. Based on these mixture models, we construct four classes of copula-type densities which are far more flexible than current popular copula densities, and outperform them in a simulated dataset and several real datasets. Supplementary material for this article is available online.     

Modelling the joint distribution of competing risks survival times using copula functions

The problem of modelling the joint distribution of survival times in a competing risks model, using copula functions, is considered. In order to evaluate this joint distribution and the related overall survival function, a system of non-linear differential equations is solved, which relates the crude and net survival functions of the modelled competing risks, through the copula. A similar approach to modelling dependent multiple decrements was applied by Carriere [Carriere, J., 1994. Dependent decrement theory. Transactions, Society of Actuaries XLVI, 45-65] who used a Gaussian copula applied to an incomplete double-decrement model which makes it difficult to calculate any actuarial functions and draw relevant conclusions. Here, we extend this methodology by studying the effect of complete and partial elimination of up to four competing risks on the overall survival function, the life expectancy and life annuity values. We further investigate how different choices of the copula function affect the resulting joint distribution of survival times and in particular the actuarial functions which are of importance in pricing life insurance and annuity products. For illustrative purposes, we have used a real data set and used extrapolation to prepare a complete multiple-decrement model up to age 120. Extensive numerical results illustrate the sensitivity of the model with respect to the choice of copula and its parameter(s).     

Modeling multi-country mortality dependence and its application in pricing survivor index swaps—A dynamic copula approach

This paper introduces mortality dependence in multi-country mortality modeling using a dynamic copula approach. Specifically, we use time-varying copula models to capture the mortality dependence structure across countries, examining both symmetric and asymmetric dependence structures. In addition, to capture the phenomenon of a heavy tail for the multi-country mortality index, we consider not only the setting of Gaussian innovations but also non-Gaussian innovations under the Lee–Carter framework model. As tests of the goodness of fit of different dynamic copula models, the pattern of mortality dependence, and the distribution of the innovations, we used empirical mortality data from Finland, France, the Netherlands, and Sweden. To understand the effect of mortality dependence on longevity derivatives, we also built a valuation framework for pricing a survivor index swap, then investigated the fair swap rates of a survivor swap numerically. We demonstrate that failing to consider the dynamic copula mortality model and non-Gaussian innovations would lead to serious underestimations of the swap rates and loss reserves.     

A generalized beta copula with applications in modeling multivariate long-tailed data

This work proposes a new copula class that we call the MGB2 copula. The new copula originates from extracting the dependence function of the multivariate GB2 distribution (MGB2) whose marginals follow the univariate generalized beta distribution of the second kind (GB2). The MGB2 copula can capture non-elliptical and asymmetric dependencies among marginal coordinates and provides a simple formulation for multi-dimensional applications. This new class features positive tail dependence in the upper tail and tail independence in the lower tail. Furthermore, it includes some well-known copula classes, such as the Gaussian copula, as special or limiting cases.To illustrate the usefulness of the MGB2 copula, we build a trivariate MGB2 copula model of bodily injury liability closed claims. Extended GB2 distributions are chosen to accommodate the right-skewness and the long-tailedness of the outcome variables. For the regression component, location parameters with continuous predictors are introduced using a nonlinear additive function. For comparison purposes, we also consider the Gumbel and t copulas, alternatives that capture the upper tail dependence. The paper introduces a conditional plot graphical tool for assessing the validation of the MGB2 copula. Quantitative and graphical assessment of the goodness of fit demonstrate the advantages of the MGB2 copula over the other copulas.     

Joint survival probability via truncated invariant copula

Given an intensity-based credit risk model, this paper studies dependence structure between default intensities. To model this structure, we use a multivariate shot noise intensity process, where jumps occur simultaneously and their sizes are correlated. Through very lengthy algebra, we obtain explicitly the joint survival probability of the integrated intensities by using the truncated invariant Farlie–Gumbel–Morgenstern copula with exponential marginal distributions. We also apply our theoretical result to pricing basket default swap spreads. This result can provide a useful guide for credit risk management.     

A multivariate extension of Sarhan and Balakrishnan’s bivariate distribution and its ageing and dependence properties

A new class of bivariate distributions (NBD) was recently introduced by Sarhan and Balakrishnan [A.M. Sarhan, N. Balakrishnan, A new class of bivariate distributions and its mixture, J. Multivariate Anal. 98 (2007) 1508-1527]. In this note, we give the joint survival function of a multivariate extension of the NBD, which is not an absolutely continuous multivariate distribution, and its marginal and extreme order statistics distributions are also derived. The multivariate ageing and dependence properties of the proposed n-dimensional distribution are also discussed, and then we analyze the stochastic ageing of its marginals and its minimum and maximum order statistics.     

Crisis and risk dependencies

The knowledge of the multivariate stochastic dependence between the returns of asset classes is of importance for many finance applications, such as asset allocation or risk management. By means of goodness-of-fit tests, we analyze for a multitude of portfolios consisting of different asset classes whether the stochastic dependence between the portfolios’ constituents can be adequately described by multivariate versions of some standard parametric copula functions. Furthermore, we test whether the stochastic dependence between the returns of different asset classes has changed during the recent financial crisis. The main findings are: First, whether a specific copula assumption can be rejected or not, crucially depends on the asset class and the time period considered. Second, different goodness-of-fit tests for copulas can yield very different results and these differences can vary for different asset classes and for different tested copulas. Third, even when using various goodness-of-fit tests for copulas, it is not always possible to differentiate between various copula assumptions. Fourth, during the financial crisis, copula assumptions are more frequently rejected. However, the results also raise some concerns over the suitability of goodness-of-fit tests for copulas as a diagnostic tool for identifying stressed risk dependencies.     

Dynamic Hedging of Portfolio Credit Risk in a Markov Copula Model

We devise a bottom-up dynamic model of portfolio credit risk where instantaneous contagion is represented by the possibility of simultaneous defaults. Due to a Markovian copula nature of the model, calibration of marginals and dependence parameters can be performed separately using a two-step procedure, much like in a standard static copula setup. In this sense this solves the bottom-up top-down puzzle which the CDO industry had been trying to do for a long time. This model can be used for any dynamic portfolio credit risk issue, such as dynamic hedging of CDOs by CDSs, or CVA computations on credit portfolios.     

��Ϣ��ɾʧ�����±�������ģ�͵Ĺ�������

??In survival analysis, most existing approaches for analysing right-censored failure time data assume that the censoring time is independent of the failure time. However, investigators often face problems involving dependent censoring, i.e., failure time and censoring time are possibly dependent and they may be censored one another, especially in clinical trials. Without accounting for such dependence, survival distributions cannot be estimated consistently. Numerous attempts to model this dependence have been made. Among them, copula models are of particular interest because of their simple structure. Proportional hazard model analysis for informative right-censored data has been discussed in this paper. An Archimedean copula is assumed for the joint distribution function of failure time and censoring time variables. Under the conditions of identifiability of the parameter of the Archimedean copula, the maximum likelihood estimators of the parameter of Archimedean copula, the parameters and the cumulative hazard function of PH model are worked out. Extensive simulation studies show that the feasibility of the proposed method and the consistency of the estimators.     

Superefficient estimation of the marginals by exploiting knowledge on the copula

We consider the problem of estimating the marginals in the case where there is knowledge on the copula. If the copula is smooth, it is known that it is possible to improve on the empirical distribution functions: optimal estimators still have a rate of convergence n^−1/2, but a smaller asymptotic variance. In this paper we show that for non-smooth copulas it is sometimes possible to construct superefficient estimators of the marginals: we construct both a copula and, exploiting the information our copula provides, estimators of the marginals with the rate of convergence logn/n.     

Relations among univariate aging, bivariate aging and dependence for exchangeable lifetimes

For a couple of lifetimes (X₁,X₂) with an exchangeable joint survival function , attention is focused on notions of bivariate aging that can be described in terms of properties of the level curves of . We analyze the relations existing among those notions of bivariate aging, univariate aging, and dependence. A goal and, at the same time, a method to this purpose is to define axiomatically a correspondence among those objects; in fact, we characterize notions of univariate and bivariate aging in terms of properties of dependence. Dependence between two lifetimes will be described in terms of their survival copula. The language of copulæ turns out to be generally useful for our purposes; in particular, we shall introduce the more general notion of semicopula. It will be seen that this is a natural object for our analysis. Our definitions and subsequent results will be illustrated by considering a few remarkable cases; in particular, we find some necessary or sufficient conditions for Schur-concavity of , or for IFR properties of the one-dimensional marginals. The case characterized by the condition that the survival copula of (X₁,X₂) is Archimedean will be considered in some detail. For most of our arguments, the extension to the case of n>2 is straightforward.     

Heath–Jarrow–Morton modelling of longevity bonds and the risk minimization of life insurance portfolios

This paper has two parts. In the first, we apply the Heath–Jarrow–Morton (HJM) methodology to the modelling of longevity bond prices. The idea of using the HJM methodology is not new. We can cite Cairns et al. [Cairns A.J., Blake D., Dowd K, 2006. Pricing death: framework for the valuation and the securitization of mortality risk. Astin Bull., 36 (1), 79–120], Miltersen and Persson [Miltersen K.R., Persson S.A., 2005. Is mortality dead? Stochastic force of mortality determined by arbitrage? Working Paper, University of Bergen] and Bauer [Bauer D., 2006. An arbitrage-free family of longevity bonds. Working Paper, Ulm University]. Unfortunately, none of these papers properly defines the prices of the longevity bonds they are supposed to be studying. Accordingly, the main contribution of this section is to describe a coherent theoretical setting in which we can properly define these longevity bond prices. A second objective of this section is to describe a more realistic longevity bonds market model than in previous papers. In particular, we introduce an additional effect of the actual mortality on the longevity bond prices, that does not appear in the literature. We also study multiple term structures of longevity bonds instead of the usual single term structure. In this framework, we derive a no-arbitrage condition for the longevity bond financial market. We also discuss the links between such HJM based models and the intensity models for longevity bonds such as those of Dahl [Dahl M., 2004. Stochastic mortality in life insurance: Market reserves and mortality-linked insurance contracts, Insurance: Math. Econom. 35 (1) 113–136], Biffis [Biffis E., 2005. Affine processes for dynamic mortality and actuarial valuations. Insurance: Math. Econom. 37, 443–468], Luciano and Vigna [Luciano E. and Vigna E., 2005. Non mean reverting affine processes for stochastic mortality. ICER working paper], Schrager [Schrager D.F., 2006. Affine stochastic mortality. Insurance: Math. Econom. 38, 81–97] and Hainaut and Devolder [Hainaut D., Devolder P., 2007. Mortality modelling with Lévy processes. Insurance: Math. Econom. (in press)], and suggest the standard pricing formula of these intensity models could be extended to more general settings.In the second part of this paper, we study the asset allocation problem of pure endowment and annuity portfolios. In order to solve this problem, we study the “risk-minimizing” strategies of such portfolios, when some but not all longevity bonds are available for trading. In this way, we introduce different basis risks.     

A continuous-time stochastic model for the mortality surface of multiple populations

We formulate, study and calibrate a continuous-time model for the joint evolution of the mortality surface of multiple populations. We model the mortality intensity by age and population as a mixture of stochastic latent factors, that can be either population-specific or common to all populations. These factors are described by affine time-(in)homogeneous stochastic processes. Traditional, deterministic mortality laws can be extended to multi-population stochastic counterparts within our framework. We detail the calibration procedure when factors are Gaussian, using centralized data-fusion Kalman filter. We provide an application based on the joint mortality of UK and Dutch males and females. Although parsimonious, the specification we calibrate provides a good fit of the observed mortality surface (ages 0–89) of both sexes and populations between 1960 and 2013.     

Multi-population mortality models: A factor copula approach

Modeling mortality co-movements for multiple populations have significant implications for mortality/longevity risk management. A few two-population mortality models have been proposed to date. They are typically based on the assumption that the forecasted mortality experiences of two or more related populations converge in the long run. This assumption might be justified by the long-term mortality co-integration and thus be applicable to longevity risk modeling. However, it seems too strong to model the short-term mortality dependence. In this paper, we propose a two-stage procedure based on the time series analysis and a factor copula approach to model mortality dependence for multiple populations. In the first stage, we filter the mortality dynamics of each population using an ARMA–GARCH process with heavy-tailed innovations. In the second stage, we model the residual risk using a one-factor copula model that is widely applicable to high dimension data and very flexible in terms of model specification. We then illustrate how to use our mortality model and the maximum entropy approach for mortality risk pricing and hedging. Our model generates par spreads that are very close to the actual spreads of the Vita III mortality bond. We also propose a longevity trend bond and demonstrate how to use this bond to hedge residual longevity risk of an insurer with both annuity and life books of business.     

Top-down approaches for integrated risk management: How accurate are they?

Banks and other financial institutions try to compute the necessary amount of total capital that they need for absorbing stochastically dependent losses from different risk types (e.g., credit risk and market risk). Two sophisticated procedures of this so-called integrated risk management are the top-down and the bottom-up approaches. When banks apply a more sophisticated risk integration approach at all, it is usually the top-down approach where copula functions are employed for linking the marginal distributions of profit and losses resulting from different risk types. However, it is not clear at all how accurate this approach is. Assuming that the bottom-up approach corresponds to the real-word data-generating process and using a comprehensive simulation study, it is shown that the top-down approach can underestimate the necessary amount of total capital for lower credit qualities. Furthermore, the direction and strength of the stochastic dependence between the risk types, the copula function employed, and the loss definitions all have an impact on the performance of the top-down approach. In addition, a goodness-of-fit test shows that, based on time series of loss data with realistic length, it is rather difficult to decide which copula function is the right one.     

Modelling Joint Behaviour of Asset Prices Using Stochastic Correlation

Association or interdependence of two stock prices is analyzed, and selection criteria for a suitable model developed in the present paper. The association is generated by stochastic correlation, given by a stochastic differential equation (SDE), creating interdependent Wiener processes. These, in turn, drive the SDEs in the Heston model for stock prices. To choose from possible stochastic correlation models, two goodness-of-fit procedures are proposed based on the copula of Wiener increments. One uses the confidence domain for the centered Kendall function, and the other relies on strong and weak tail dependence. The constant correlation model and two different stochastic correlation models, given by Jacobi and hyperbolic tangent transformation of Ornstein-Uhlenbeck (HtanOU) processes, are compared by analyzing daily close prices for Apple and Microsoft stocks. The constant correlation, i.e., the Gaussian copula model, is unanimously rejected by the methods, but all other two are acceptable at a 95% confidence level. The analysis also reveals that even for Wiener processes, stochastic correlation can create tail dependence, unlike constant correlation, which results in multivariate normal distributions and hence zero tail dependence. Hence models with stochastic correlation are suitable to describe more dangerous situations in terms of correlation risk.

     

Supplier default dependencies: Empirical evidence from the automotive industry

Common ways to mitigate the detrimental consequences of supplier bankruptcies are to install redundancy and to pursue a multiple-sourcing strategy. This is based on the assumption that the adverse event of one supplier going out of business is largely independent from the default of other suppliers. However, this implicit assumption does not hold in all cases. This study – based on empirical data from automotive suppliers – reveals that default dependencies among suppliers do often exist and can have significant consequences. We use copula functions, a method of representing joint distribution functions with particular marginals, to capture the default dependence between automotive suppliers and to simulate various default dependence scenarios. We also conduct a comparative static analysis illustrating the significant impact of default correlation in a supplier portfolio. Our findings should spur managers to analyze their supplier portfolios with respect to default dependencies, and to take this phenomenon into consideration when making sourcing decisions.     

Representing Sparse Gaussian DAGs as Sparse R-Vines Allowing for Non-Gaussian Dependence

Modeling dependence in high-dimensional systems has become an increasingly important topic. Most approaches rely on the assumption of a multivariate Gaussian distribution such as statistical models on directed acyclic graphs (DAGs). They are based on modeling conditional independencies and are scalable to high dimensions. In contrast, vine copula models accommodate more elaborate features like tail dependence and asymmetry, as well as independent modeling of the marginals. This flexibility comes however at the cost of exponentially increasing complexity for model selection and estimation. We show a novel connection between DAGs with limited number of parents and truncated vine copulas under sufficient conditions. This motivates a more general procedure exploiting the fast model selection and estimation of sparse DAGs while allowing for non-Gaussian dependence using vine copulas. By numerical examples in hundreds of dimensions, we demonstrate that our approach outperforms the standard method for vine structure selection. Supplementary material for this article is available online.