Copula credibility for aggregate loss models

This paper develops credibility predictors of aggregate losses using a longitudinal data framework. For a model of aggregate losses, the interest is in predicting both the claims number process as well as the claims amount process. In a longitudinal data framework, one encounters data from a cross-section of risk classes with a history of insurance claims available for each risk class. Further, explanatory variables for each risk class over time are available to help explain and predict both the claims number and claims amount process.For the marginal claims distributions, this paper uses generalized linear models, an extension of linear regression, to describe cross-sectional characteristics. Elliptical copulas are used to model the dependencies over time, extending prior work that used multivariate t-copulas. The claims number process is represented using a Poisson regression model that is conditioned on a sequence of latent variables. These latent variables drive the serial dependencies among claims numbers; their joint distribution is represented using an elliptical copula. In this way, the paper provides a unified treatment of both the continuous claims amount and discrete claims number processes.The paper presents an illustrative example of Massachusetts automobile claims. Estimates of the latent claims process parameters are derived and simulated predictions are provided.     

基于倾斜-t Copula的半参数马尔可夫链

在不指定时间序列结构的情况下,我们的分布模型是基于多变量离散时间的相应马尔可夫族和相关变量一维的边际分布.这样的模型可以同时处理时间序列之间的相互依赖和每个时间序列沿时间方向的依赖.具体的参数copula被指定为倾斜-t. 倾斜-t Copla能够处理不对称,偏斜和粗尾的数据分布.三个股票指数日均收益的实证研究表明,倾斜-t copula的马尔可夫模型要比以下模型更好:倾斜正态Copula马可夫, t-copula马可夫, 倾斜-t copula但无马尔可夫特性.     

Long-tail longitudinal modeling of insurance company expenses

The insurance industry is known to have high operating expenses in the financial services sector. Insurers, investors and regulators are interested in models to understand the behavior of expenses. However, the current practice ignores skewness, occasional negative values as well as their temporal dependence.Addressing these three features, this paper develops a longitudinal model of insurance company expenses that can be used for prediction, to identify unusual behavior, and to measure firm efficiency. Specifically, we use a three-parameter asymmetric Laplace density for the marginal distribution of insurers’ expenses in each year. Copula functions are employed to accommodate their temporal dependence. As a function of explanatory variables, the location parameter allows us to analyze an insurer’s expenses in light of the firm’s characteristics. Our model can be interpreted as a longitudinal quantile regression.The analysis is performed using property-casualty insurance company data from the National Association of Insurance Commissioners of years 2001-2006. Due to the long-tailed nature of insurers’ expenses, two alternative approaches are proposed to improve the performance of the longitudinal quantile regression model: rescaling and transformation. Predictive densities are derived that allow one to compare the predictions for individual insurers in a hold-out-sample. Both predictive models are shown to be reasonable with the rescaling method outperforming the transformation method. Compared with standard longitudinal models, our model is shown to be superior in identifying insurers’ unusual behavior.     

Long-tail longitudinal modeling of insurance company expenses

The insurance industry is known to have high operating expenses in the financial services sector. Insurers, investors and regulators are interested in models to understand the behavior of expenses. However, the current practice ignores skewness, occasional negative values as well as their temporal dependence.Addressing these three features, this paper develops a longitudinal model of insurance company expenses that can be used for prediction, to identify unusual behavior, and to measure firm efficiency. Specifically, we use a three-parameter asymmetric Laplace density for the marginal distribution of insurers’ expenses in each year. Copula functions are employed to accommodate their temporal dependence. As a function of explanatory variables, the location parameter allows us to analyze an insurer’s expenses in light of the firm’s characteristics. Our model can be interpreted as a longitudinal quantile regression.The analysis is performed using property–casualty insurance company data from the National Association of Insurance Commissioners of years 2001–2006. Due to the long-tailed nature of insurers’ expenses, two alternative approaches are proposed to improve the performance of the longitudinal quantile regression model: rescaling and transformation. Predictive densities are derived that allow one to compare the predictions for individual insurers in a hold-out-sample. Both predictive models are shown to be reasonable with the rescaling method outperforming the transformation method. Compared with standard longitudinal models, our model is shown to be superior in identifying insurers’ unusual behavior.     

Detecting atypical observations in financial data: the forward search for elliptical copulas

In the last few years, copulas have been widely applied in many field of studies. Concentrating our attention on financial applications, we pursue the goal to detect multivariate atypical observations by extending to elliptical copulas the forward search originally introduced in linear and nonlinear regression by Atkinson and Riani (Robust diagnostic regression analysis. Springer, New York, 2000). Considering that, in the forward search, observations are ranked according to their closeness to the fitted data, we need to define a measure through which to initialize, progress and monitor the search. We achieve this goal building up the forward search for elliptical copulas relying on the squared Mahalanobis distance. Stressing the need to find theoretical boundaries for the inference on outliers, we introduce a procedure for computing envelopes as in Riani and Atkinson (Adv Data Anal Classif 1:123–141, 2007). Once defined our framework, we apply the forward search to a simulated environment where contaminations are exogenously introduced then, we carry out the analysis on n equity log-return real time series.     

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.     

Tail negative dependence and its applications for aggregate loss modeling

Tail order of copulas can be used to describe the strength of dependence in the tails of a joint distribution. When the value of tail order is larger than the dimension, it may lead to tail negative dependence. First, we prove results on conditions that lead to tail negative dependence for Archimedean copulas. Using the conditions, we construct new parametric copula families that possess upper tail negative dependence. Among them, a copula based on a scale mixture with a generalized gamma random variable (GGS copula) is useful for modeling asymmetric tail negative dependence. We propose mixed copula regression based on the GGS copula for aggregate loss modeling of a medical expenditure panel survey dataset. For this dataset, we find that there exists upper tail negative dependence between loss frequency and loss severity, and the introduction of tail negative dependence structures significantly improves the aggregate loss modeling.     

从SPVⅡ分布生成的新的多元偏态t分布的有关性质

随机向量的t分布属于椭球等高分布族,然而,它是对称分布.在许多诸如经济学、生理学、社会学等领域中,有时回归模型中的随机误差不再满足对称性,通常表现出高度的偏态性(skewness).于是就有了偏态椭球等高分布族.本文在已有的多元偏态t分布的基础上,着重研究它的分布性质,包括线性组合分布、边缘分布、条件分布及各阶矩.     

From Archimedean to Liouville copulas

We use a recent characterization of the d-dimensional Archimedean copulas as the survival copulas of d-dimensional simplex distributions (McNeil and Nešlehová (2009) [1]) to construct new Archimedean copula families, and to examine the relationship between their dependence properties and the radial parts of the corresponding simplex distributions. In particular, a new formula for Kendall’s tau is derived and a new dependence ordering for non-negative random variables is introduced which generalises the Laplace transform order. We then generalise the Archimedean copulas to obtain Liouville copulas, which are the survival copulas of Liouville distributions and which are non-exchangeable in general. We derive a formula for Kendall’s tau of Liouville copulas in terms of the radial parts of the corresponding Liouville distributions.     

Multivariate negative binomial models for insurance claim counts

It is no longer uncommon these days to find the need in actuarial practice to model claim counts from multiple types of coverage, such as the ratemaking process for bundled insurance contracts. Since different types of claims are conceivably correlated with each other, the multivariate count regression models that emphasize the dependency among claim types are more helpful for inference and prediction purposes. Motivated by the characteristics of an insurance dataset, we investigate alternative approaches to constructing multivariate count models based on the negative binomial distribution. A classical approach to induce correlation is to employ common shock variables. However, this formulation relies on the NB-I distribution which is restrictive for dispersion modeling. To address these issues, we consider two different methods of modeling multivariate claim counts using copulas. The first one works with the discrete count data directly using a mixture of max-id copulas that allows for flexible pair-wise association as well as tail and global dependence. The second one employs elliptical copulas to join continuitized data while preserving the dependence structure of the original counts. The empirical analysis examines a portfolio of auto insurance policies from a Singapore insurer where claim frequency of three types of claims (third party property damage, own damage, and third party bodily injury) are considered. The results demonstrate the superiority of the copula-based approaches over the common shock model. Finally, we implemented the various models in loss predictive applications.     

Efficient MCMC estimation of some elliptical copula regression models through scale mixtures of normals

This paper proposes an efficient estimation method for some elliptical copula regression models by expressing both copula density and marginal density functions as scale mixtures of normals (SMN). Implementing these models using the SMN is novel and allows efficient estimation via Bayesian methods. An innovative algorithm for the case of complex semicontinuous margins is also presented. We utilize the facts that copulas are invariant to the location and scale of the margins; all elliptical distributions have the same correlation structure; and some densities can be represented by the SMN. Two simulation studies, one on continuous margins and the other on semicontinuous margins, highlight the favorable performance of the proposed methods. Two empirical studies, one on the US excess returns and one on the Thai wage earnings, further illustrate the applicability of the proposals.     

过离散次数分布模型的尾部特征

在保险精算和生物统计等领域,离散型次数分布模型的应用十分广泛.当实际数据的尾部较长(即过离散),且零点的概率较大时,许多模型的拟合效果往往欠佳.本文通过计算概率之比的极限和偏度系数,对混合泊松分布和复合泊松分布的右尾特征和零点概率进行了比较,给出了它们的尾部排列顺序,以及尾部长短与零点概率的关系,从而为模型的构造或选择提供了一种指导.本文最后应用一组实际数据说明了在构造或选择次数分布模型时如何考虑尾部特征,从而改善对实际数据的拟合效果.     

Learning Interactions via Hierarchical Group-Lasso Regularization

We introduce a method for learning pairwise interactions in a linear regression or logistic regression model in a manner that satisfies strong hierarchy: whenever an interaction is estimated to be nonzero, both its associated main effects are also included in the model. We motivate our approach by modeling pairwise interactions for categorical variables with arbitrary numbers of levels, and then show how we can accommodate continuous variables as well. Our approach allows us to dispense with explicitly applying constraints on the main effects and interactions for identifiability, which results in interpretable interaction models. We compare our method with existing approaches on both simulated and real data, including a genome-wide association study, all using our R package glinternet.     

Noncentral matrix quadratic forms of the skew elliptical variables

In this paper, the noncentral matrix quadratic forms of the skew elliptical variables are studied. A family of the matrix variate noncentral generalized Dirichlet distributions is introduced as the extension of the noncentral Wishart distributions, the Dirichlet distributions and the noncentral generalized Dirichlet distributions. Main distributional properties are investigated. These include probability density and closure property under linear transformation and marginalization, the joint distribution of the sub-matrices of the matrix quadratic forms in the skew elliptical variables and the moment generating functions and Bartlett's decomposition of the matrix quadratic forms in the skew normal variables. Two versions of the noncentral Cochran's Theorem for the matrix variate skew normal distributions are obtained, providing sufficient and necessary conditions for the quadratic forms in the skew normal variables to have the matrix variate noncentral generalized Dirichlet distributions. Applications include the properties of the least squares estimation in multivariate linear model and the robustness property of the Wilk's likelihood ratio statistic in the family of the matrix variate skew elliptical distributions.     

Distorted Mix Method for constructing copulas with tail dependence

This paper introduces a method for constructing copula functions by combining the ideas of distortion and convex sum, named Distorted Mix Method. The method mixes different copulas with distorted margins to construct new copula functions, and it enables us to model the dependence structure of risks by handling the central and tail parts separately. By applying the method we can modify the tail dependence of a given copula to any desired level measured by tail dependence function and tail dependence coefficients of marginal distributions. As an application, a tight bound for asymptotic Value-at-Risk of order statistics is obtained by using the method. An empirical study shows that copulas constructed by this method fit the empirical data of SPX 500 Index and FTSE 100 Index very well in both central and tail parts.     

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.     

Multivariate geometric stable distributions in financial applications

We describe a class of multivariate geometric stable laws that can be used in modeling multivariate financial portfolios of securities. These heavy tailed distributions are stable with respect to geometric summation and accommodate the possibility of market crashes. We look at bivariate currency exchange rates data and show that its main features, peakedness and heavy tails, are very well captured by the geometric stable model.     

离散型半参数广义线性纵向数据模型的方差成分检验

在回归分析中, 随机误差是否存在方差非齐性是大家十分关心的问题, 本文根据Laplace展开原理针对随机效应的影响研究了基于纵向数据的离散型半参数广义线性模型的方差成分检验,得到了Score检验统计量, 最后通过一个实例和计算机模拟验证了本文所提出的方法的有效性.     

Portfolio separation properties of the skew-elliptical distributions, with generalizations

The two-fund separation property of the elliptical distributions is extended to the skew-elliptical case by adding a number of funds equaling the rank of the skewness matrix. The singular extended skew-elliptical distributions are covered, as is a further generalization to the case where the set conditioned upon is not an orthant.     

Robust portfolio optimization with copulas

Conditional Value at Risk (CVaR) is widely used in portfolio optimization as a measure of risk. CVaR is clearly dependent on the underlying probability distribution of the portfolio. We show how copulas can be introduced to any problem that involves distributions and how they can provide solutions for the modeling of the portfolio. We use this to provide the copula formulation of the CVaR of a portfolio. Given the critical dependence of CVaR on the underlying distribution, we use a robust framework to extend our approach to Worst Case CVaR (WCVaR). WCVaR is achieved through the use of rival copulas. These rival copulas have the advantage of exploiting a variety of dependence structures, symmetric and not. We compare our model against two other models, Gaussian CVaR and Worst Case Markowitz. Our empirical analysis shows that WCVaR can asses the risk more adequately than the two competitive models during periods of crisis.