Bounds on total economic capital: the DNB case study

Most banks use the top-down approach to aggregate their risk types when computing total economic capital. Following this approach, marginal distributions for each risk type are first independently estimated and then merged into a joint model using a copula function. Due to lack of reliable data, banks tend to manually select the copula as well as its parameters. In this paper we assess the model risk related to the choice of a specific copula function. The aim is to compute upper and lower bounds on the total economic capital for the aggregate loss distribution of DNB, the largest Norwegian bank, and the key tool for computing these bounds is the Rearrangement Algorithm introduced in Embrechts et al. (J. Bank. Financ. 37(8):2750–2764 2013). The application of this algorithm to a real situation poses a series of numerical challenges and raises a number of warnings which we illustrate and discuss.     

Distribution modeling for reliability analysis: Impact of multiple dependences and probability model selection

Reliability analysis requires modeling of joint probability distribution of uncertain parameters, which can be a challenge since the random variables representing the parameter uncertainties may be correlated. For convenience, a Gaussian data dependence is commonly assumed for correlated random variables. This paper first investigates the effect of multidimensional non-Gaussian data dependences underlying the multivariate probability distribution on reliability results. Using different bivariate copulas in a vine structure, various data dependences can be modeled. The associated copula parameters are identified from available statistical information by moment matching techniques. After the development of the vine copula model for representing the multivariate probability distribution, the reliability involving correlated random variables is evaluated based on the Rosenblatt transformation. The impact of data dependence is significant because a large deviation in failure probability is observed, which emphasizes the need for accurate dependence characterization. A practical method for dependence modeling based on limited data is thus provided. The result demonstrates that the non-Gaussian data dependences can be real in practice, and the reliability can be biased if the Gaussian dependence is used inappropriately. Moreover, the effect of conditioning order on reliability should not be overlooked except that the vine structure contains only one type of copula.     

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

有偏抽样下带信息观察和删失的面板数据的统计分析

面板数据经常出现在许多研究领域, 比如纵向跟踪研究. 在很多情况下, 纵向反应变量与观察 时间和删失时间都有关系. 本文在有偏抽样下, 针对这些相关性存在的情况, 利用一个不能观察的潜在 变量, 提出了一个联合建模方法来刻画纵向反应变量与观察时间和删失时间的相关性, 获得了模型中 回归参数的估计方程以及估计的渐近性质, 并通过数值模拟验证了这些估计在小样本下也是有效的, 同时把该估计方法用于一组实际的膀胱癌数据分析中.     

Quantile regression based on counting process approach under semi-competing risks data

In this paper, we investigate the quantile regression analysis for semi-competing risks data in which a non-terminal event may be dependently censored by a terminal event. The estimation of quantile regression parameters for the non-terminal event is complicated. We cannot make inference on the non-terminal event without extra assumptions. Thus, we handle this problem by assuming that the joint distribution of the terminal event and the non-terminal event follows a parametric copula model with unspecified marginal distributions. We use the stochastic property of the martingale method to estimate the quantile regression parameters under semi-competing risks data. We also prove the large sample properties of the proposed estimator, and introduce a model diagnostic approach to check model adequacy. From simulation results, it shows that the proposed estimator performs well. For illustration, we apply our proposed approach to analyze a real data.     

Bayesian Nonparametric Inference for a Multivariate Copula Function

The paper presents a general Bayesian nonparametric approach for estimating a high dimensional copula. We first introduce the skew–normal copula, which we then extend to an infinite mixture model. The skew–normal copula fixes some limitations in the Gaussian copula. An MCMC algorithm is developed to draw samples from the correct posterior distribution and the model is investigated using both simulated and real applications.     

Parameter uncertainty, sensitivity analysis and prediction error in a water-balance hydrological model

Analysis of uncertainty is often neglected in the evaluation of complex systems models, such as computational models used in hydrology or ecology. Prediction uncertainty arises from a variety of sources, such as input error, calibration accuracy, parameter sensitivity and parameter uncertainty. In this study, various computational approaches were investigated for analysing the impact of parameter uncertainty on predictions of streamflow for a water-balance hydrological model used in eastern Australia. The parameters and associated equations which had greatest impact on model output were determined by combining differential error analysis and Monte Carlo simulation with stochastic and deterministic sensitivity analysis. This integrated approach aids in the identification of insignificant or redundant parameters and provides support for further simplifications in the mathematical structure underlying the model. Parameter uncertainty was represented by a probability distribution and simulation experiments revealed that the shape (skewness) of the distribution had a significant effect on model output uncertainty. More specifically, increasing negative skewness of the parameter distribution correlated with decreasing width of the model output confidence interval (i.e. resulting in less uncertainty). For skewed distributions, characterisation of uncertainty is more accurate using the confidence interval from the cumulative distribution rather than using variance. The analytic approach also identified the key parameters and the non-linear flux equation most influential in affecting model output uncertainty.     

A multi-criteria decision making approach for location planning for urban distribution centers under uncertainty

Location planning for urban distribution centers is vital in saving distribution costs and minimizing traffic congestion arising from goods movement in urban areas. In this paper, we present a multi-criteria decision making approach for location planning for urban distribution centers under uncertainty. The proposed approach involves identification of potential locations, selection of evaluation criteria, use of fuzzy theory to quantify criteria values under uncertainty and application of fuzzy TOPSIS to evaluate and select the best location for implementing an urban distribution center. Sensitivity analysis is performed to determine the influence of criteria weights on location planning decisions for urban distribution centers.The strength of the proposed work is the ability to deal with uncertainty arising due to a lack of real data in location planning for new urban distribution centers. The proposed approach can be practically applied by logistics operators in deciding on the location of new distribution centers considering the sustainable freight regulations proposed by municipal administrations. A numerical application is provided to illustrate the approach.     

Quantifying the risk using copulae with nonparametric marginals

We show that copulae and kernel estimation can be mixed to estimate the risk of an economic loss. We analyze the properties of the Sarmanov copula. We find that the maximum pseudo-likelihood estimation of the dependence parameter associated with the copula with double transformed kernel estimation to estimate marginal cumulative distribution functions is a useful method for approximating the risk of extreme dependent losses when we have large data sets. We use a bivariate sample of losses from a real database of auto insurance claims.     

Sampling Schemes for Bayesian Variable Selection in Generalized Linear Models

Bayesian approaches to prediction and the assessment of predictive uncertainty in generalized linear models are often based on averaging predictions over different models, and this requires methods for accounting for model uncertainty. When there are linear dependencies among potential predictor variables in a generalized linear model, existing Markov chain Monte Carlo algorithms for sampling from the posterior distribution on the model and parameter space in Bayesian variable selection problems may not work well. This article describes a sampling algorithm based on the Swendsen-Wang algorithm for the Ising model, and which works well when the predictors are far from orthogonality. In problems of variable selection for generalized linear models we can index different models by a binary parameter vector, where each binary variable indicates whether or not a given predictor variable is included in the model. The posterior distribution on the model is a distribution on this collection of binary strings, and by thinking of this posterior distribution as a binary spatial field we apply a sampling scheme inspired by the Swendsen-Wang algorithm for the Ising model in order to sample from the model posterior distribution. The algorithm we describe extends a similar algorithm for variable selection problems in linear models. The benefits of the algorithm are demonstrated for both real and simulated data.     

Improved probabilistic point estimation schemes for uncertainty analysis

The probabilistic point estimation (PPE) methods replace the probability distribution of the random parameters of a model with a finite number of discrete points in sample space selected in such a way to preserve limit probabilistic information of involved random parameters. Most PPE methods developed thus far match the distribution of random parameters up to the third statistical moment and, in general, could provide reasonable accurate estimation only for the first two statistical moments of model output. This study proposes two optimization-based point selection schemes for the PPE methods to enhance the accuracy of higher-order statistical moments estimation for model output. Several test models of varying degrees of complexity and nonlinearity are used to examine the performance of the proposed point selection schemes. The results indicate that the proposed point selection schemes provide significantly more accurate estimation of model output uncertainty features than the existing schemes.     

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.     

A hierarchical model for the joint mortality analysis of pension scheme data with missing covariates

A hierarchical model is developed for the joint mortality analysis of pension scheme datasets. The proposed model allows for a rigorous statistical treatment of missing data. While our approach works for any missing data pattern, we are particularly interested in a scenario where some covariates are observed for members of one pension scheme but not the other. Therefore, our approach allows for the joint modelling of datasets which contain different information about individual lives. The proposed model generalizes the specification of parametric models when accounting for covariates. We consider parameter uncertainty using Bayesian techniques. Model parametrization is analysed in order to obtain an efficient MCMC sampler, and address model selection. The inferential framework described here accommodates any missing-data pattern, and turns out to be useful to analyse statistical relationships among covariates. Finally, we assess the financial impact of using the covariates, and of the optimal use of the whole available sample when combining data from different mortality experiences.     

Model selection of copulas: AIC versus a cross validation copula information criterion

Akaike Information Criterion (AIC) is frequently employed in the semiparametric setting of selection of copula models, even though as a model selection tool it was developed in a parametric setting. Recently a Copula Information Criterion (CIC) has been especially designed for copula model selection. In this paper we examine the two approaches and present a simulation study where the performance of a cross-validated version of CIC is compared with the AIC criterion. Only minor differences are observed.     

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.     

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

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

单指标众数模型的统计诊断及在波士顿房价分析中的应用

在金融、经济、社会科学、气候科学、环境科学、工程技术和生物医学等领域,数据分布常常呈现出尖峰厚尾的特征,且密度分布是不对称的有偏分布。此时,单指标众数模型是刻画这些特征的一个重要方法。为此,非常有必要研究该模型下的统计诊断。本文将考虑单指标众数模型基于数据删除模型和众数漂移模型的统计诊断与局部影响分析。模拟研究和波士顿房价数据的结果表明所提出的方法是有效和可行的。     

联合订货下基于ø-散度的多产品库存鲁棒优化模型

针对多产品联合库存决策问题,在市场需求不确定条件下,建立了考虑联合订货成本的多产品库存鲁棒优化模型。针对不确定市场需求,采用一系列未知概率的离散情景进行描述,给出了基于最小最大准则的鲁棒对应模型,并证明了(s,S)库存策略的最优性。进一步,在仅知多产品市场需求历史数据基础上,采用基于ø-散度的数据驱动方法构建了满足一定置信度要求的关于未知需求概率分布的不确定集。在此基础上,为获得(s,S)库存策略的相关参数,运用拉格朗日对偶方法将所建模型等价转化为易于求解的数学规划问题。最后,通过数值计算分析了Kullback-Leibler散度和Cressie-Read散度以及不同的置信水平下的多产品库存绩效,并将其与真实分布下应用鲁棒库存策略得到的库存绩效进行对比。结果表明,需求分布信息的缺失虽然会导致一定的库存绩效损失,但损失值很小,表明基于文中方法得到的库存策略能够有效抑制需求不确定性扰动,具有良好的鲁棒性。     

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.     

Locally averaged Bayesian Dirichlet metrics for learning the structure and the parameters of Bayesian networks

The marginal likelihood of the data computed using Bayesian score metrics is at the core of score+search methods when learning Bayesian networks from data. However, common formulations of those Bayesian score metrics rely on free parameters which are hard to assess. Recent theoretical and experimental works have also shown that the commonly employed BDe score metric is strongly biased by the particular assignments of its free parameter known as the equivalent sample size. This sensitivity means that poor choices of this parameter lead to inferred BN models whose structure and parameters do not properly represent the distribution generating the data even for large sample sizes. In this paper we argue that the problem is that the BDe metric is based on assumptions about the BN model parameters distribution assumed to generate the data which are too strict and do not hold in real settings. To overcome this issue we introduce here an approach that tries to marginalize the meta-parameter locally, aiming to embrace a wider set of assumptions about these parameters. It is shown experimentally that this approach offers a robust performance, as good as that of the standard BDe metric with an optimum selection of its free parameter and, in consequence, this method prevents the choice of wrong settings for this widely applied Bayesian score metric.