Dependent competing risks: Cause elimination and its impact on survival

The dependent competing risks model of human mortality is considered, assuming that the dependence between lifetimes is modelled by a multivariate copula function. The effect on the overall survival of removing one or more causes of death is explored under two alternative definitions of removal, ignoring the causes and eliminating them. Under the two definitions of removal, expressions for the overall survival functions in terms of the specified copula (density) and the net (marginal) survival functions are given. The net survival functions are obtained as a solution to a system of non-linear differential equations, which relates them through the specified copula (derivatives) to the crude (sub-) survival functions, estimated from data. The overall survival functions in a model with four competing risks, cancer, cardiovascular diseases, respiratory diseases and all other causes grouped together, have been implemented and evaluated, based on cause-specific mortality data for England and Wales published by the Office for National Statistics, for the year 2007. We show that the two alternative definitions of removal of a cause of death have different effects on the overall survival and in particular on the life expectancy at birth and at age 65, when one, two or three of the competing causes are removed. An important conclusion is that the eliminating definition is better suited for practical use in competing risks’ applications, since it is more intuitive, and it suffices to consider only positive dependence between the lifetimes which is not the case under the alternative ignoring definition.     

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.     

Modelling Operational Risk Losses with Graphical Models and Copula Functions

The management of Operational Risk has been a difficult task due to the lack of data and the high number of variables. In this project, we treat operational risks as multivariate variables. In order to model them, copula functions are employed, which are a widely used tool in finance and engineering for building flexible joint distributions. The purpose of this research is to propose a new methodology for modelling Operational Risks and estimating the required capital. It combines the use of graphical models and the use of copula functions along with hyper-Markov law. Historical loss data of an Italian bank is used, in order to explore the methodology’s behaviour and its potential benefits.      

Order statistics from mixed exchangeable random variables

Two different exchangeable samples are considered and these two samples are assumed to be independent of each other. From these two samples a new sample is combined and treated as a single set of observations. The distribution of a single order statistic and the joint distribution of two order statistics for a new mixed sample are derived and expressed in terms of joint distribution functions. As a special case the distribution of a single order statistic and the joint distribution of two nonadjacent order statistics from exchangeable random variables are obtained. The results presented in this paper allows widespread applications in modelling of various lifetime data, biomedical sciences, reliability and survival analysis, actuarial sciences etc., where the assumption of independence of data cannot be accepted and the exchangeability is a more realistic assumption.     

Stochastic comparisons of lifetimes of series and parallel systems with dependent and heterogeneous components

This work considers stochastic comparisons of lifetimes of series and parallel systems with dependent and heterogeneous components having lifetimes following the proportional odds (PO) model. The joint distribution of component lifetimes is modeled by Archimedean survival copula. We discuss some potential applications of our findings in system reliability and actuarial science.     

Estimating the joint survival probabilities of married individuals

We estimate the joint survival probability of spouses using a large random sample drawn from a Dutch census. As benchmarks we use two bivariate Weibull models. We consider more flexible models, using a semi-nonparametric approach, by extending the independent Weibull distribution using squared polynomials. Also based on a nonparametric comparison, we find that extending the independent Weibull distribution by a squared third order polynomial shows the best performance. We illustrate our model by calculating remaining life expectancies and annuity values. We find that the husbands life expectancy at birth is generally increasing with his wifes age of death and the wifes life expectancy at birth is generally increasing with her husbands age of death. Ignoring the dependence between the remaining lifetimes of spouses may lead to an underestimation of the value of a joint annuity and an overestimation of the value of a single-life annuity, but less than suggested on the basis of the previous literature.     

Nonparametric estimation of competing risks models with covariates

In competing risks model, several failure times arise potentially. The smallest failure time and its index only are observed. Without specific assumptions, the joint or even the marginal distribution functions of the underlying failure times are not identifiable (A. Tsiatis, Proc. Natl. Acad. Sci. USA 72 (1975) 20). Nonetheless, if each individual is characterized by a “sufficiently informative” set of covariates, these distributions are identifiable under some conditions of regularity (J.J. Heckman and B. Honoré, Biometrika 76 (1989) 325). In this paper, nonparametric kernel estimators of the joint distribution function of failure times conditional on the covariates are proposed. Their weak and strong consistency are discussed.     

递增年金的双随机模型

The dual random models about the life insurance and social pension insurance have received considerable attention in the recent articles on actuarial theory and applications. This paper discusses a general kind of increasing annuity based on its force of interest accumulationfunction as a general random process. The dual random model of the present value of the benefits of the increasing annuity has been set, and their moments have been calculated under certainconditions.     

An application of comonotonicity theory in a stochastic life annuity framework

A life annuity contract is an insurance instrument which pays pre-scheduled living benefits conditional on the survival of the annuitant. In order to manage the risk borne by annuity providers, one needs to take into account all sources of uncertainty that affect the value of future obligations under the contract. In this paper, we define the concept of annuity rate as the conditional expected present value random variable of future payments of the annuity, given the future dynamics of its risk factors. The annuity rate deals with the non-diversifiable systematic risk contained in the life annuity contract, and it involves mortality risk as well as investment risk. While it is plausible to assume that there is no correlation between the two risks, each affects the annuity rate through a combination of dependent random variables. In order to understand the probabilistic profile of the annuity rate, we apply comonotonicity theory to approximate its quantile function. We also derive accurate upper and lower bounds for prediction intervals for annuity rates. We use the Lee-Carter model for mortality risk and the Vasicek model for the term structure of interest rates with an annually renewable fixed-income investment policy. Different investment strategies can be handled using this framework.     

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.     

The impact of the determinants of mortality on life insurance and annuities

Extended risk classification has become an important issue recently in life insurance and annuity markets. Various risk factors have been explored and identified by past research. Using those risk factors, one can construct various risk classes. This enables insurers to provide more equitable life insurance and annuity benefits for individuals in different risk classes and to manage mortality/longevity risk more efficiently. The challenge of modeling mortality using various risk factors is to reflect complicated mortality dynamics in a model while maintaining statistical significance. This paper discusses the development of a mortality model that reflects the impact of various risk factors on mortality. Longitudinal survey data from the Canadian National Population Health Survey was used to determine the significant risk factors and quantify their effect on mortality. The model is used to illustrate how the various risk factors influence actuarial present values of life insurance and annuity benefits.     

����ɾʧ�»������Ӻ����IJ���ģ�͵�ͳ�Ʒ���

In collecting clinical data, data would be censored due to competing risks or patient withdrawal. The statistical inference for censoring data is always based on the assumption that the failure time and censoring time is independent. But in practice the failure time and censoring time are often dependent. Dependent censoring make the job to deal with censoring data more complicated. In this paper, we assume that the joint distribution of the failure time variable and censoring time variable is a function of their marginal distributions. This function is called a copula. Under prespecified copulas, the maximum likelihood estimators for cox proportional hazards models are worked out. Statistical analysis results are carried by simulations. When dependent censoring happens, the proposed method will do better than the traditional method used in independent situations. Simulation results show that the proposed method can get efficient estimations.     

A note on multiple life premiums for dependent lifetimes

We study the properties of multiple life annuity and insurance premiums for general symmetric and survival statuses in the case when the joint distribution of future lifetimes has a dependence structure belonging to some nonparametric neighbourhood of independence. The size of the neighbourhood is controlled by a single parameter, which enables us to model really weak as well as stronger dependencies. We provide bounds on the difference of multiple life premiums for vectors of dependent and independent future lifetimes with the same univariate marginal distributions. Each such upper bound can be treated as a premium loading related to the strength of lifetimes’ dependence.     

Likelihood-based inference for bivariate latent failure time models with competing risks under the generalized FGM copula

Many existing latent failure time models for competing risks do not provide closed form expressions of sub-distribution functions. This paper suggests a generalized FGM copula models with the Burr III failure time distribution such that the sub-distribution functions have closed form expressions. Under the suggested model, we develop a likelihood-based inference method along with its computational tools and asymptotic theory. Based on the expressions of the sub-distribution functions, we propose goodness-of-fit tests. Simulations are conducted to examine the performance of the proposed methods. A real data from the reliability analysis of the radio transmitter-receivers are analyzed to illustrate the proposed methods. The computational programs are made available in the R package GFGM.copula.     

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.     

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

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

Dependence structure between LIBOR rates by copula method

This paper discusses the correlation structure between London Interbank Offered Rates (LIBOR) by using the copula function. We start from one simplified model of A. Brace, D. Gatarek, and M. Musiela (1997) and find out that the copula function between two LIBOR rates can be expressed as a sum of an infinite series, where the main term is a distribution function with Gaussian copula. Partial differential equation method is used for deriving the copula expansion. Numerical results show that the copula of the LIBOR rates and Gaussian copula are very close in the central region and differ in the tail, and the Gaussian copula approximation to the copula function between the LIBOR rates provides satisfying results in the normal situation.     

Upper comonotonicity

In this article, we study a new notion called upper comonotonicity, which is a generalization of the classical notion of comonotonicity. A random vector is upper-comonotonic if its components are moving in the same direction simultaneously when their values are greater than some thresholds. We provide a characterization of this new notion in terms of both the joint distribution function and the underlying copula. The copula characterization allows us to study the coefficient of upper tail dependence as well as the distributional representation of an upper-comonotonic random vector. As an application to financial economics, we show that the several commonly used risk measures, like the Value-at-Risk, the Tail Value-at-Risk, and the expected shortfall, are additive, not only for sum of comonotonic risks, but also for sum of upper-comonotonic risks, provided that the level of probability is greater than a certain threshold.     

Longevity risk management for life and variable annuities: The effectiveness of static hedging using longevity bonds and derivatives

For many years, the longevity risk of individuals has been underestimated, as survival probabilities have improved across the developed world. The uncertainty and volatility of future longevity has posed significant risk issues for both individuals and product providers of annuities and pensions. This paper investigates the effectiveness of static hedging strategies for longevity risk management using longevity bonds and derivatives (q-forwards) for the retail products: life annuity, deferred life annuity, indexed life annuity, and variable annuity with guaranteed lifetime benefits. Improved market and mortality models are developed for the underlying risks in annuities. The market model is a regime-switching vector error correction model for GDP, inflation, interest rates, and share prices. The mortality model is a discrete-time logit model for mortality rates with age dependence. Models were estimated using Australian data. The basis risk between annuitant portfolios and population mortality was based on UK experience. Results show that static hedging using q-forwards or longevity bonds reduces the longevity risk substantially for life annuities, but significantly less for deferred annuities. For inflation-indexed annuities, static hedging of longevity is less effective because of the inflation risk. Variable annuities provide limited longevity protection compared to life annuities and indexed annuities, and as a result longevity risk hedging adds little value for these products.     

��I�ͻ�Ͻ�β��ָ��——�������ֲ�����ʧЧģ�͵�ͳ�Ʒ���

In this paper, we investigate a competing risks model based on exponentiated Weibull distribution under Type-I progressively hybrid censoring scheme. To estimate the unknown parameters and reliability function, the maximum likelihood estimators and asymptotic confidence intervals are derived. Since Bayesian posterior density functions cannot be given in closed forms, we adopt Markov chain Monte Carlo method to calculate approximate Bayes estimators and highest posterior density credible intervals. To illustrate the estimation methods, a simulation study is carried out with numerical results. It is concluded that the maximum likelihood estimation and Bayesian estimation can be used for statistical inference in competing risks model under Type-I progressively hybrid censoring scheme.