ON THE ESTIMATION OF SURVIVAL FUNCTION UNDER RANDOM CENSORSHIP

ABSTRACT

The paper deals with an improvement of the well-known Kaplan–Meier estimator of survival function when the censoring mechanism is random and independent of the failure times. Small sample size properties of the new estimator, as well as the original Kaplan–Meier estimator are inspected by means of Monte Carlo simulations. It follows from the simulations that the proposed estimator prevails with respect to some basic statistical characteristics.     

Maximum likelihood estimation of bivariate survival probabilities and its existence and uniqueness

Nonparametric maximum likelihood estimation of bivariate survival probabilities is developed for interval censored survival data. We restrict our attention to the situation where response times within pairs are not distinguishable, and the univariate survival distribution is the same for any individual within any pair. Campbell's (1981) model is modified to incorporate this restriction. Existence and uniqueness of maximum likelihood estimators are discussed. This methodology is illustrated with a bivariate life table analysis of an angioplasty study where each patient undergoes two procedures.     

Profile likelihood approaches for semiparametric copula and frailty models for clustered survival data

ABSTRACT

In clustered survival data, the dependence among individual survival times within a cluster has usually been described using copula models and frailty models. In this paper we propose a profile likelihood approach for semiparametric copula models with different cluster sizes. We also propose a likelihood ratio method based on profile likelihood for testing the absence of association parameter (i.e. test of independence) under the copula models, leading to the boundary problem of the parameter space. For this purpose, we show via simulation study that the proposed likelihood ratio method using an asymptotic chi-square mixture distribution performs well as sample size increases. We compare the behaviors of the two models using the profile likelihood approach under a semiparametric setting. The proposed method is demonstrated using two well-known data sets.     

Nonparametric analysis of recurrent gap time data

Abstract

Recurrent event data are frequently encountered in longitudinal studies. In many applications, the times between successive recurrent events (gap times) are often of interest and lead to problems that have received much attention recently. In this article, using the approach of inverse probability-of-censoring weights (IPCW), we propose nonparametric estimators for the estimation of the bivariate distribution and survival functions for gap times of recurrent event data. We also consider the estimation of Kendall’s tau for two gap times by expressing it as an integral functional of the bivariate survival function. The asymptotic properties of the proposed estimators are established. Simulation studies are conducted to investigate their finite sample performance.     

A two-stage estimation in the Clayton–Oakes model with marginal linear transformation models for multivariate failure time data

This paper considers the analysis of multivariate survival data where the marginal distributions are specified by semiparametric transformation models, a general class including the Cox model and the proportional odds model as special cases. First, consideration is given to the situation where the joint distribution of all failure times within the same cluster is specified by the Clayton–Oakes model (Clayton, Biometrika 65:141–151, l978; Oakes, J R Stat Soc B 44:412–422, 1982). A two-stage estimation procedure is adopted by first estimating the marginal parameters under the independence working assumption, and then the association parameter is estimated from the maximization of the full likelihood function with the estimators of the marginal parameters plugged in. The asymptotic properties of all estimators in the semiparametric model are derived. For the second situation, the third and higher order dependency structures are left unspecified, and interest focuses on the pairwise correlation between any two failure times. Thus, the pairwise association estimate can be obtained in the second stage by maximizing the pairwise likelihood function. Large sample properties for the pairwise association are also derived. Simulation studies show that the proposed approach is appropriate for practical use. To illustrate, a subset of the data from the Diabetic Retinopathy Study is used.     

Mimicking Cox-Regression

Abstract

A class of objective functions, related to the Cox partial likelihood, that generates unbiased estimating equations is proposed. These equations allow for estimation of interest parameters when nuisance parameters are proportional to expectations. Examples of the objective functions are applied to binary data with a log-link in three situations: independent observations, independent groups of observations with common random intercept and discrete survival data. It is pointed out that the Peto–Breslow approximation to the partial likelihood with discrete failure times fits a conditional model with a log-link.     

Survival Estimation with Missing Censoring Times for the Generalized Koziol–Green Model

ABSTRACT

When analyzing time-to-event data, there are various situations in which right censoring times for unfailed units are missing. In that context, by taking a supplementary sample of a convenient percentage of unfailed units, we propose a semi-parametric method for estimating a survival function under the natural extension of the Koziol–Green model to double random censoring. Some large sample properties of this estimator are derived. We prove uniform strong consistency and asymptotic weak convergence to a Gaussian process. A simulation study is also presented in order to analyze the behavior of the proposed estimator.     

On exact strong laws of large numbers for ratios of random variables and their applications

Abstract

We study the almost sure convergence of weighted sums of ratios of independent random variables satisfying some general, mild conditions. The obtained results are applied to exact laws for order statistics. An exact law for independent random variables which are nonidentically distributed is also proved and applied to ratios of adjacent order statistics for a sample of uniformly distributed random variables.     

Admissible and optimal sampling strategy for estimating finite population mean in randomized response surveys with multiple responses

ABSTRACT

We consider the problem of estimation of a finite population mean (or proportion) related to a sensitive character under a randomized response model when independent responses are obtained from each sampled individual as many times as he/she is selected in the sample and prove the admissibility of a sampling strategy in a class of comparable linear unbiased strategies. We prove that the admissible strategy is also optimal in this class under a super-population model.     

Conditional Quantile Estimation for Truncated and Associated Data

Abstract

In survival or reliability studies, it is common to have data which are not only incomplete but weakly dependent too. Random truncation and censoring are two common forms of such data when they are neither independent nor strongly mixing but rather associated. The focus of this paper is on estimating conditional distribution and conditional quantile functions for randomly left truncated data satisfying association condition. We aim at deriving strong uniform consistency rates and asymptotic normality for the estimators and thereby, extend to association case some results stated under iid and α-mixing hypotheses. The performance of the quantile function estimator is evaluated on simulated data sets.     

Clustering and classification problems in genetics through U-statistics

ABSTRACT

Genetic data are frequently categorical and have complex dependence structures that are not always well understood. For this reason, clustering and classification based on genetic data, while highly relevant, are challenging statistical problems. Here we consider a versatile U-statistics-based approach for non-parametric clustering that allows for an unconventional way of solving these problems. In this paper we propose a statistical test to assess group homogeneity taking into account multiple testing issues and a clustering algorithm based on dissimilarities within and between groups that highly speeds up the homogeneity test. We also propose a test to verify classification significance of a sample in one of two groups. We present Monte Carlo simulations that evaluate size and power of the proposed tests under different scenarios. Finally, the methodology is applied to three different genetic data sets: global human genetic diversity, breast tumour gene expression and Dengue virus serotypes. These applications showcase this statistical framework's ability to answer diverse biological questions in the high dimension low sample size scenario while adapting to the specificities of the different datatypes.     

Estimation of the joint survival function for successive duration times under double-truncation and right-censoring

ABSTRACT

In incident cohort studies, survival data often include subjects who have had an initiate event at recruitment and may potentially experience two successive events (first and second) during the follow-up period. When disease registries or surveillance systems collect data based on incidence occurring within a specific calendar time interval, the initial event is usually subject to double truncation. Furthermore, since the second duration process is observable only if the first event has occurred, double truncation and dependent censoring arise. In this article, under the two sampling biases with an unspecified distribution of truncation variables, we propose a nonparametric estimator of the joint survival function of two successive duration times using the inverse-probability-weighted (IPW) approach. The consistency of the proposed estimator is established. Based on the estimated marginal survival functions, we also propose a two-stage estimation procedure for estimating the parameters of copula model. The bootstrap method is used to construct confidence interval. Numerical studies demonstrate that the proposed estimation approaches perform well with moderate sample sizes.     

A METHOD FOR ESTIMATING PARAMETER IN NONNEGATIVE MA(1) MODELS

ABSTRACT

A method for estimating parameter in nonnegative MA(1) models is proposed and investigated in the paper. The method also gives nontrivial confidence sets on confidence level 1. Small sample properties of new estimator are demonstrated in a simulation study.     

Variance Components Testing in ANOVA‐Type Mixed Models

The purpose of this article is threefold. First, variance components testing for ANOVA ‐type mixed models is considered, in which response may not be divided into independent sub‐vectors, whereas most of existing methods are for models where response can be divided into independent sub‐vectors. Second, testing that a certain subset of variance components is zero. Third, as normality is often violated in practice, it is desirable to construct tests under very mild assumptions. To achieve these goals, an adaptive difference‐based test and an adaptive trace‐based test are constructed. The test statistics are asymptotically normal under the null hypothesis, are consistent against all global alternatives and can detect local alternatives distinct from the null at a rate as close to n ^{? 1 ∕ 2} as possible with n being the sample size. Moreover, when the dimensions of variance components in different sets are bounded, we develop a test with chi‐square as its limiting null distribution. The finite sample performance of the tests is examined via simulations, and a real data set is analysed for illustration.     

Estimating the association parameter for copula models under dependent censoring

Summary. Many biomedical studies involve the analysis of multiple events. The dependence between the times to these end points is often of scientific interest. We investigate a situation when one end point is subject to censoring by the other. The model assumptions of Day and co-workers and Fine and co-workers are extended to more general structures where the level of association may vary with time. Two types of estimating function are proposed. Asymptotic properties of the proposed estimators are derived. Their finite sample performance is studied via simulations. The inference procedures are applied to two real data sets for illustration.     

Testing for the partial Koziol–Green model with covariates

In this paper we consider some non-parametric goodness-of-fit statistics for testing the partial Koziol–Green regression model. In this model, the response at a given covariate value is subject to random right censoring by two independent censoring times. One of these censoring times is informative in the sense that its survival function is some power of the survival function of the response. The goodness-of-fit statistics are based on an underlying empirical process for which large sample theory is obtained.     

Exact inference on multiple exponential populations under a joint type-II progressive censoring scheme

ABSTRACT

Recently Mondal and Kundu [Mondal S, Kundu D. A new two sample type-II progressive censoring scheme. Commun Stat Theory Methods. 2018. doi:10.1080/03610926.2018.1472781] introduced a Type-II progressive censoring scheme for two populations. In this article, we extend the above scheme for more than two populations. The aim of this paper is to study the statistical inference under the multi-sample Type-II progressive censoring scheme, when the underlying distributions are exponential. We derive the maximum likelihood estimators (MLEs) of the unknown parameters when they exist and find out their exact distributions. The stochastic monotonicity of the MLEs has been established and this property can be used to construct exact confidence intervals of the parameters via pivoting the cumulative distribution functions of the MLEs. The distributional properties of the ordered failure times are also obtained. The Bayesian analysis of the unknown model parameters has been provided. The performances of the different methods have been examined by extensive Monte Carlo simulations. We analyse two data sets for illustrative purposes.     

Joint modeling of longitudinal data with a dependent terminal event

ABSTRACT

Longitudinal data often arise in longitudinal follow-up studies, and there may exist a dependent terminal event such as death that stops the follow-up. In this article, we propose a new joint modeling for the analysis of longitudinal data with informative observation times via a dependent terminal event and two latent variables. Estimating equations are developed for parameter estimation, and asymptotic properties of the resulting estimators are established. In addition, a generalization of the joint model with time-varying coefficients for the longitudinal response variable is considered, and goodness-of-fit methods for assessing the adequacy of the model are also provided. The proposed method works well in our simulation studies, and is applied to a data set from a bladder cancer study.     

Testing for abrupt breaks in variance structures with smooth changes

Abstract

We investigate the problem of testing for variance breaks in the case where the variance structure is assumed to be smoothly time-varying under the null. Since the classical tests are aimed to detect any change in the variance, they are not able to distinguish between smooth non constant variance and abrupt breaks. In this paper a new procedure for detecting variance breaks taking into account for smooth changes in the variance under the null is proposed. The finite sample properties of the test we introduce are investigated by Monte Carlo experiments. The theoretical outputs are illustrated using U.S. macroeconomic data.     

Stochastic comparisons of two-unit repairable systems

Abstract

We consider two models of two-unit repairable systems: cold standby system and warm standby system. We suppose that the lifetimes and repair times of the units are all independent exponentially distributed random variables. Using stochastic orders we compare the lifetimes of systems under different assumptions on the parameters of exponential distributions. We also consider a cold standby system where the lifetimes and repair times of its units are not necessarily exponentially distributed.