Hazard rate estimation for left truncated and right censored data

Left truncation and right censoring (LTRC) presents a unique challenge for nonparametric estimation of the hazard rate of a continuous lifetime because consistent estimation over the support of the lifetime is impossible. To understand the problem and make practical recommendations, the paper explores how the LTRC affects a minimal (called sharp) constant of a minimax MISE convergence over a fixed interval. The corresponding theory of sharp minimax estimation of the hazard rate is presented, and it shows how right censoring, left truncation and interval of estimation affect the MISE. Obtained results are also new for classical cases of censoring or truncation and some even for the case of direct observations of the lifetime of interest. The theory allows us to propose a relatively simple data-driven estimator for small samples as well as the methodology of choosing an interval of estimation. The estimation methodology is tested numerically and on real data.     

Nonparametric estimation of interval-censored failure time data in the presence of informative censoring

Nonparametric estimation of a survival function is one of the most commonly asked questions in the analysis of failure time data and for this, a number of procedures have been developed under various types of censoring structures (Kalbfleisch and Prentice, 2002). In particular, several algorithms are available for interval-censored failure time data with independent censoring mechanism (Sun, 2006; Turnbull, 1976). In this paper, we consider the interval-censored data where the censoring mechanism may be related to the failure time of interest, for which there does not seem to exist a nonparametric estimation procedure. It is well-known that with informative censoring, the estimation is possible only under some assumptions. To attack the problem, we take a copula model approach to model the relationship between the failure time of interest and censoring variables and present a simple nonparametric estimation procedure. The method allows one to conduct a sensitivity analysis among others.     

Estimation and confidence bands of a conditional survival function with censoring indicators missing at random

The nonparametric estimator of the conditional survival function proposed by Beran is a useful tool to evaluate the effects of covariates in the presence of random right censoring. However, censoring indicators of right censored data may be missing for different reasons in many applications. We propose some estimators of the conditional cumulative hazard and survival functions which allow to handle this situation. We also construct the likelihood ratio confidence bands for them and obtain their asymptotic properties. Simulation studies are used to evaluate the performances of the estimators and their confidence bands.     

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

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

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.     

Robust estimation of the Pickands dependence function under random right censoring

We consider robust nonparametric estimation of the Pickands dependence function under random right censoring. The estimator is obtained by applying the minimum density power divergence criterion to properly transformed bivariate observations. The asymptotic properties are investigated by making use of results for Kaplan–Meier integrals. We investigate the finite sample properties of the proposed estimator with a simulation experiment and illustrate its practical applicability on a dataset of insurance indemnity losses.     

Estimation of a Support Curve via Order Statistics

This paper deals with nonparametric estimation of the boundary curve of the support of a bivariate density function. This estimation problem arises in various contexts, such as for example scatterpoint image analysis and frontier estimation in econometrics. The setup in this paper is a general one, allowing the bivariate density function to be infinite, bounded away from zero or zero at the boundary. Two estimators for the boundary curve are introduced, both based on order statistics. The asymptotic distribution of the estimators and their rate of convergence are established. Via a comparison of the rates of convergence we recommend which estimator to use in a particular situation. Both estimators can be used as an initial estimator in a two-stage procedure, designed for getting a better estimation. Simulation studies demonstrate the finite-sample behavior of the estimators and the proposed two-stage procedure. We illustrate the procedure on a data set on American electric utility companies.     

Hazard Rate Estimation in Nonparametric Regression with Censored Data

Consider a regression model in which the responses are subject to random right censoring. In this model, Beran studied the nonparametric estimation of the conditional cumulative hazard function and the corresponding cumulative distribution function. The main idea is to use smoothing in the covariates. Here we study asymptotic properties of the corresponding hazard function estimator obtained by convolution smoothing of Beran's cumulative hazard estimator. We establish asymptotic expressions for the bias and the variance of the estimator, which together with an asymptotic representation lead to a weak convergence result. Also, the uniform strong consistency of the estimator is obtained.     

����Copula������Ӧ��-ǿ��ģ�Ϳɿ��ȵķDz�������

In this paper, a nonparametric method for reliability of the stress-strength model is proposed when the dependent stress variable and strength variable are subject to right censoring. The dependence between variables is measured by the common Farlie-Gumbel-Morgenstern copula function and Clayton copula function. Using the empirical process theory, consistency and asymptotic normality of the proposed estimator is established in this paper. The results of numerical simulation show that the proposed method performs well in the case of finite sample. The method proposed in this paper has a wide application prospect in practice.     

左截断双删失数据下可加可乘风险率模型的估计

在生存分析中,可加可乘风险率模型常用来研究协变量对初始事件和终止事件之间持续时间的影响效应。在本文中,我们考虑了在初始事件存在部分区间删失,同时终止事件存在左截断右删失的情形下,持续时间的可加可乘风险率模型的估计问题。我们提出了一个两阶段估计过程来估计模型的回归参数。并通过模拟分析验证了估计的大样本性质。最后利用该方法分析了恶性黑色素瘤手术治疗数据。     

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.     

Bivariate option pricing with copulas

The adoption of copula functions is suggested in order to price bivariate contingent claims. Copulas enable the marginal distributions extracted from vertical spreads in the options markets to be imbedded in a multivariate pricing kernel. It is proved that such a kernel is a copula function, and that its super-replication strategy is represented by the Fréchet bounds. Applications provided include prices for binary digital options, options on the minimum and options to exchange one asset for another. For each of these products, no-arbitrage pricing bounds, as well as values consistent with the independence of the underlying assets are provided. As a final reference value, a copula function calibrated on historical data is used.     

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.     

Estimation in additive cox models by marginal integration

Assuming an additive model on the covariate effect in proportional hazards regression, we consider the estimation of the component functions. The estimator is based on the marginal integration method. Then we use a new kind of nonparametric estimator as the pilot estimator of the marginal integration. The pilot estimator is constructed by an analogy to the two-sample problems and by appealing to the principles of local partial likelihood and local linear fitting. We derive the asymptotic distribution of the marginal integration estimator of the component functions. The result of a simulation study is also given.     

纵向数据部分线性测量误差模型的二次推断函数估计

基于纵向数据部分线性测量误差模型, 研究了模型中兴趣参数部分回归系数的估计问题. 首先采用B样条方法逼近模型中的非参数函数, 然后提出修正的二次推断函数(QIF)方法对模型中参数部分的回归系数进行估计, 所提方法可以提高估计的效率. 在一定的正则条件下, 证明了所得到的估计量具有相合性和渐近正态性. 最后, 通过模拟研究和实例分析验证了所提出估计方法的有限大样本性质.     

On testing equality of pairwise rank correlations in a multivariate random vector

Spearman’s rank-correlation coefficient (also called Spearman’s rho) represents one of the best-known measures to quantify the degree of dependence between two random variables. As a copula-based dependence measure, it is invariant with respect to the distribution’s univariate marginal distribution functions. In this paper, we consider statistical tests for the hypothesis that all pairwise Spearman’s rank correlation coefficients in a multivariate random vector are equal. The tests are nonparametric and their asymptotic distributions are derived based on the asymptotic behavior of the empirical copula process. Only weak assumptions on the distribution function, such as continuity of the marginal distributions and continuous partial differentiability of the copula, are required for obtaining the results. A nonparametric bootstrap method is suggested for either estimating unknown parameters of the test statistics or for determining the associated critical values. We present a simulation study in order to investigate the power of the proposed tests. The results are compared to a classical parametric test for equal pairwise Pearson’s correlation coefficients in a multivariate random vector. The general setting also allows the derivation of a test for stochastic independence based on Spearman’s rho.     

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.     

On the Estimation of a Support Curve of Indeterminate Sharpness

We propose nonparametric methods for estimating the support curve of a bivariate density, when the density decreases at a rate which might vary along the curve. Attention is focused on cases where the rate of decrease is relatively fast, this being the most difficult setting. It demands the use of a relatively large number of bivariate order statistics. By way of comparison, support curve estimation in the context of slow rates of decrease of the density may be addressed using methods that employ only a relatively small number of order statistics at the extremities of the point cloud. In this paper we suggest a new type of estimator, based on projecting onto an axis those data values lying within a thin rectangular strip. Adaptive univariate methods are then applied to the problem of estimating an endpoint of the distribution on the axis. The new method is shown to have theoretically optimal performance in a range of settings. Its numerical properties are explored in a simulation study.     

Parametric estimation of a bivariate stable Lévy process

We propose a parametric model for a bivariate stable Lévy process based on a Lévy copula as a dependence model. We estimate the parameters of the full bivariate model by maximum likelihood estimation. As an observation scheme we assume that we observe all jumps larger than some ε>0 and base our statistical analysis on the resulting compound Poisson process. We derive the Fisher information matrix and prove asymptotic normality of all estimates when the truncation point ε→0. A simulation study investigates the loss of efficiency because of the truncation.     

Hazard function estimation with cause-of-death data missing at random

Hazard function estimation is an important part of survival analysis. Interest often centers on estimating the hazard function associated with a particular cause of death. We propose three nonparametric kernel estimators for the hazard function, all of which are appropriate when death times are subject to random censorship and censoring indicators can be missing at random. Specifically, we present a regression surrogate estimator, an imputation estimator, and an inverse probability weighted estimator. All three estimators are uniformly strongly consistent and asymptotically normal. We derive asymptotic representations of the mean squared error and the mean integrated squared error for these estimators and we discuss a data-driven bandwidth selection method. A simulation study, conducted to assess finite sample behavior, demonstrates that the proposed hazard estimators perform relatively well. We illustrate our methods with an analysis of some vascular disease data.