Flexible Multivariate Density Estimation With Marginal Adaptation

This article is concerned with multivariate density estimation. We discuss deficiencies in two popular multivariate density estimators—mixture and copula estimators, and propose a new class of estimators that combines the advantages of both mixture and copula modeling, while being more robust to their weaknesses. Our method adapts any multivariate density estimator using information obtained by separately estimating the marginals. We propose two marginally adapted estimators based on a multivariate mixture of normals and a mixture of factor analyzers estimators. These estimators are implemented using computationally efficient split-and-elimination variational Bayes algorithms. It is shown through simulation and real-data examples that the marginally adapted estimators are capable of improving on their original estimators and compare favorably with other existing methods. Supplementary materials for this article are available online.     

Copula models for insurance claim numbers with excess zeros and time-dependence

This paper develops two copula models for fitting the insurance claim numbers with excess zeros and time-dependence. The joint distribution of the claims in two successive periods is modeled by a copula with discrete or continuous marginal distributions. The first model fits two successive claims by a bivariate copula with discrete marginal distributions. In the second model, a copula is used to model the random effects of the conjoint numbers of successive claims with continuous marginal distributions. Zero-inflated phenomenon is taken into account in the above copula models. The maximum likelihood is applied to estimate the parameters of the discrete copula model. A two-step procedure is proposed to estimate the parameters in the second model, with the first step to estimate the marginals, followed by the second step to estimate the unobserved random effect variables and the copula parameter. Simulations are performed to assess the proposed models and methodologies.     

Modelling stochastic mortality for dependent lives

Stochastic mortality, i.e. modelling death arrival via a jump process with stochastic intensity, is gaining an increasing reputation as a way to represent mortality risk. This paper is a first attempt to model the mortality risk of couples of individuals, according to the stochastic intensity approach. Dependence between the survival times of the members of a couple is captured by an Archimedean copula.We also provide a methodology for fitting the joint survival function by working separately on the (analytical) marginals and on the (analytical) copula. First, we provide a sample-based calibration for the intensity, using a time-homogeneous, non mean-reverting, affine process: this gives the marginal survival functions. Then we calibrate and select the best fit copula according to the Wang and Wells [Wang, W., Wells, M.T., 2000b. Model selection and semiparametric inference for bivariate failure-time data. J. Amer. Statis. Assoc. 95, 62-72] methodology for censored data. By coupling the calibrated marginals with the best fit copula, we obtain a joint survival function, which incorporates the stochastic nature of mortality improvements.We apply the methodology to a well known insurance data set, using a sample generation. The best fit copula turns out to be one listed in [Nelsen, R.B., 2006. An Introduction to Copulas, Second ed. In: Springer Series], which implies not only positive dependence, but dependence increasing with age.     

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.      

Portfolio value-at-risk estimation in energy futures markets with time-varying copula-GARCH model

This paper combines copula functions with GARCH-type models to construct the conditional joint distribution, which is used to estimate Value-at-Risk (VaR) of an equally weighted portfolio comprising crude oil futures and natural gas futures in energy market. Both constant and time-varying copulas are applied to fit the dependence structure of the two assets returns. The findings show that the constant Student t copula is a good compromise for effectively fitting the dependence structure between crude oil futures and natural gas futures. Moreover, the skewed Student t distribution has a better fit than Normal and Student t distribution to the marginal distribution of each asset. Asymmetries and excess kurtosis are found in marginal distributions as well as in dependence. We estimate VaR of the underlying portfolio to be 95% and 99%, by using the Monte Carlo simulation. Then using backtesting, we compare the out-of-sample forecasting performances of VaR estimated by different models.     

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.     

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.     

A Flexible Model for Generalized Linear Regression with Measurement Error

This paper focuses on the question of specification of measurement error distribution and the distribution of true predictors in generalized linear models when the predictors are subject to measurement errors. The standard measurement error model typically assumes that the measurement error distribution and the distribution of covariates unobservable in the main study are normal. To make the model flexible enough we, instead, assume that the measurement error distribution is multivariate t and the distribution of true covariates is a finite mixture of normal densities. Likelihood–based method is developed to estimate the regression parameters. However, direct maximization of the marginal likelihood is numerically difficult. Thus as an alternative to it we apply the EM algorithm. This makes the computation of likelihood estimates feasible. The performance of the proposed model is investigated by simulation study.     

A Copula-Based Method to Build Diffusion Models with Prescribed Marginal and Serial Dependence

This paper investigates the probabilistic properties that determine the existence of space-time transformations between diffusion processes. We prove that two diffusions are related by a monotone space-time transformation if and only if they share the same serial dependence. The serial dependence of a diffusion process is studied by means of its copula density and the effect of monotone and non-monotone space-time transformations on the copula density is discussed. This approach provides a methodology to build diffusion models by freely combining prescribed marginal behaviors and temporal dependence structures. Explicit expressions of copula densities are provided for tractable models.     

Bivariate nonparametric estimation of the Pickands dependence function using Bernstein copula with kernel regression approach

In this study, a new nonparametric approach using Bernstein copula approximation is proposed to estimate Pickands dependence function. New data points obtained with Bernstein copula approximation serve to estimate the unknown Pickands dependence function. Kernel regression method is then used to derive an intrinsic estimator satisfying the convexity. Some extreme-value copula models are used to measure the performance of the estimator by a comprehensive simulation study. Also, a real-data example is illustrated. The proposed Pickands estimator provides a flexible way to have a better fit and has a better performance than the conventional estimators.     

Efficient risk simulations for linear asset portfolios in the t-copula model

We consider the problem of calculating tail probabilities of the returns of linear asset portfolios. As a flexible and accurate model for the logarithmic returns we use the t-copula dependence structure and marginals following the generalized hyperbolic distribution. Exact calculation of the tail-loss probabilities is not possible and even simulation leads to challenging numerical problems. Applying a new numerical inversion method for the generation of the marginals and importance sampling with carefully selected mean shift we develop an efficient simulation algorithm. Numerical results for a variety of realistic portfolio examples show an impressive performance gain.     

On the strong approximation of bootstrapped empirical copula processes with applications

The purpose of the present paper is to provide a strong invariance principle for the generalized bootstrapped empirical copula processwith the rate of the approximation for multivariate empirical processes. As a by-product, we obtain a uniform-in-bandwidth consistency result for kernel-type estimators of copula derivatives, which is of its own interest. We introduce also the delta-sequence estimators of the copula derivatives. The applications discussed here are change-point detection in multivariate copula models, nonparametric tests of stochastic vectorial independence and the law of iterated logarithm for the generalized bootstrapped empirical copula process. Finally, a general notion of bootstrapped empirical copula process constructed by exchangeably weighting the sample is presented.     

Modeling volatility and dependency of agricultural price and production indices of Thailand: Static versus time-varying copulas

Volatility and dependence structure are two main sources of uncertainty in many economic issues, such as exchange rates, future prices and agricultural product prices etc. who fully embody uncertainty among relationship and variation. This paper aims at estimating the dependency between the percentage changes of the agricultural price and agricultural production indices of Thailand and also their conditional volatilities using copula-based GARCH models. The motivation of this paper is twofold. First, the strategic department of agriculture of Thailand would like to have reliable empirical models for the dependency and volatilities for use in policy strategy. Second, this paper provides less restrictive models for dependency and the conditional volatility GARCH. The copula-based multivariate analysis used in this paper nested the traditional multivariate as a special case (Tae-Hwy and Xiangdong, 2009) [13]. Appropriate marginal distributions for both, the percentage changes of the agricultural price and agricultural production indices were selected for their estimation. Static as well as time varying copulas were estimated. The empirical results were found that the suitable margins were skew t distribution and the time varying copula i.e., the time varying rotate Joe copula (270°) was the choice for the policy makers to follow. The one-period ahead forecasted-growth rate of agricultural price index conditional on growth rate of agricultural production index was also provided as an example of forecasting it using the resulted margins and time-varying copula based GARCH model.     

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

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

Penalized marginal likelihood estimation of finite mixtures of Archimedean copulas

This paper proposes finite mixtures of different Archimedean copula families as a flexible tool for modelling the dependence structure in multivariate data. A novel approach to estimating the parameters in this mixture model is presented by maximizing the penalized marginal likelihood via iterative quadratic programming. The motivation for the penalized marginal likelihood stems from an underlying Bayesian model that imposes a prior distribution on the parameter of each Archimedean copula family. An approximative marginal likelihood is obtained by a classical quadrature discretization of the integral w.r.t. each family-specific prior distribution, thus yielding a finite mixture model. Family-specific smoothness penalties are added and the penalized marginal likelihood is maximized using an iterative quadratic programming routine. For comparison purposes, we also present a fully Bayesian approach via simulation-based posterior computation. The performance of the novel estimation approach is evaluated by simulations and two examples involving the modelling of the interdependence of exchange rates and of wind speed measurements, respectively. For these examples, penalized marginal likelihood estimates are compared to the corresponding Bayesian estimates.     

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.     

Rank-based inference tools for copula regression,with property and casualty insurance applications

Rank-based procedures are commonly used for inference in copula models for continuous responses whose behavior does not depend on covariates. This paper describes how these procedures can be adapted to the broader framework in which (possibly non-linear) regression models for the marginal responses are linked by a copula that does not depend on covariates. The validity of many of these techniques can be derived from the asymptotic equivalence between the classical empirical copula process and its analog based on suitable residuals from the marginal models. Moment-based parameter estimation and copula goodness-of-fit tests are shown to remain valid under weak conditions on the marginal error term distributions, even when the residual-based empirical copula process fails to converge weakly. The performance of these procedures is evaluated through simulation in the context of two general insurance applications: micro-level multivariate insurance claims, and dependent loss triangles.     

Improving efficient marginal estimators in bivariate models with parametric marginals

Suppose we have data from a bivariate model with parametric marginals. Efficient estimators of the parameters in the marginal models are generally not efficient in the bivariate model. In this article, we propose a method of improving these marginal estimators and demonstrate that the magnitude of this improvement can be as large as 100 percent in some cases.     

Computationally efficient learning of multivariate <Emphasis Type="Italic">t</Emphasis> mixture models with missing information

A finite mixture model using the multivariate t distribution has been well recognized as a robust extension of Gaussian mixtures. This paper presents an efficient PX-EM algorithm for supervised learning of multivariate t mixture models in the presence of missing values. To simplify the development of new theoretic results and facilitate the implementation of the PX-EM algorithm, two auxiliary indicator matrices are incorporated into the model and shown to be effective. The proposed methodology is a flexible mixture analyzer that allows practitioners to handle real-world multivariate data sets with complex missing patterns in a more efficient manner. The performance of computational aspects is investigated through a simulation study and the procedure is also applied to the analysis of real data with varying proportions of synthetic missing values.