Empirical Likelihood for the Additive Hazards Model with Current Status Data

This paper considers joint analysis of current status and marker data using a threshold model based on first hitting times. A failure time is defined as the time at which a subject's latent health status process first decreases to zero. We extend the bivariate Wiener process model in Whitmore et al. (1998) to the case when only current status data are available. We develop maximum likelihood estimation procedures and provide simulation studies. We apply our methods to a motivating example involving liver tumors in mice.

For regression analysis of doubly truncated data, we propose two empirical likelihood (EL) inference approaches, called non-smooth EL and non-smooth Jackknife EL (JEL), to make inference about regression parameters based on the generalized estimating equations of existing weighted rank estimators. The limiting distributions of non-smooth log-EL and log-JEL ratios statistics are derived and non-smooth EL, and JEL confidence intervals for any specified component of regression parameters are obtained. We carry out extensive simulation studies to compare the proposed approaches with the random weighting (RW) approach. The simulation results demonstrate that the non-smooth EL and JEL confidence intervals have better performances than the RW confidence intervals based on coverage probability and average length of confidence intervals of regression parameters when the dependent variable is subject to the double truncation. A real data example is provided to illustrate the proposed approaches.

In this article, we consider the application of the empirical likelihood method to a partially linear single-index model. We focus on the case where some covariates are measured with additive errors. It is shown that the empirical log-likelihood ratio at the true parameter converges to the standard chi-square distribution. Simulations show that the proposed confidence region has coverage probability which is closer to the nominal level, as well as narrower than those of normal approximation method. A real data example is given.

SUMMARY A great deal of recent attention has focused on the estimation of survival distributions based on current status data, an extreme form of interval censored data. This particular data structure arises in a wide variety of applications where cross-sectional observation either naturally occurs or is preferred to more traditional forms of follow-up. Here we consider current status data in the context of competing risks. We briefly consider simple parametric models as a backdrop to nonparametric procedures. We make some brief comparisons and remarks regarding the nonparametric maximum likelihood estimator. The ideas are illustrated on the data of Krailo & Pike (1983) which considers estimation of the age distribution at both natural and operative menopause. We also consider the case where there is exact observation of failure times due to one of the competing risks when failure occurs prior to the monitoring time. Basic techniques of survival analysis focus on estimation of the survival distribution based on various forms of censored data. An extreme form of censoring arises where the only information on studied individuals is their survival status at a single monitoring time. This particular data structure is known as current status data. Nonparametric estimation of the survival function and semiparametric techniques for related regression models, based on current status data, have been much studied of late. Here, we consider estimation of sub-survival functions based on current status observations in the presence of competing risks. For simplicity, we describe only two competing risks, but the ideas are easily extended when there are more than two possible causes of failure. Let )j denote the cause-specific hazard function for cause j = 1, 2 (Kalbfleisch & Prentice, 2002, p. 252). If T and J are the random variables that measure time to failure and cause of failure, respectively, the two sub-distribution functions of interest are

In this paper, the authors consider the empirical likelihood method for a first-order generalized random coefficient integer-valued autoregressive process. The authors establish the log empirical likelihood ratio statistic and obtain its limiting distribution. Furthermore, the authors investigate the point estimation, confidence regions and hypothesis testing for the parameters of interest. The performance of empirical likelihood method is illustrated by a simulation study and a real data example.

This article develops empirical likelihood methodology for a class of long range dependent processes driven by a stationary Gaussian process. We consider population parameters that are defined by estimating equations in the time domain. It is shown that the standard block empirical likelihood (BEL) method, with a suitable scaling, has a non‐standard limit distribution based on a multiple Wiener–Itô integral. Unlike the short memory time series case, the scaling constant involves unknown population quantities that may be difficult to estimate. Alternative versions of the empirical likelihood method, involving the expansive BEL (EBEL) methods are considered. It is shown that the EBEL renditions do not require an explicit scaling and, therefore, remove this undesirable feature of the standard BEL. However, the limit law involves the long memory parameter, which may be estimated from the data. Results from a moderately large simulation study on finite sample properties of tests and confidence intervals based on different empirical likelihood methods are also reported.

This paper discusses regression analysis of current status data with dependent censoring, a problem that often occurs in many areas such as cross-sectional studies, epidemiological investigations and tumorigenicity experiments. Copula model-based methods are commonly employed to tackle this issue. However, these methods often face challenges in terms of model and parameter identification. The primary aim of this paper is to propose a copula-based analysis for dependent current status data, where the association parameter is left unspecified. Our method is based on a general class of semiparametric linear transformation models and parametric copulas. We demonstrate that the proposed semiparametric model is identifiable under certain regularity conditions from the distribution of the observed data. For inference, we develop a sieve maximum likelihood estimation method, using Bernstein polynomials to approximate the nonparametric functions involved. The asymptotic consistency and normality of the proposed estimators are established. Finally, to demonstrate the effectiveness and practical applicability of our method, we conduct an extensive simulation study and apply the proposed method to a real data example.

ABSTRACTIn this paper, we introduce a partially linear single-index additive hazards model with current status data. Both the unknown link function of the single-index term and the cumulative baseline hazard function are approximated by B-splines under a monotonicity constraint on the latter. The sieve method is applied to estimate the nonparametric and parametric components simultaneously. We show that, when the nonparametric link function is an exact B-spline, the resultant estimator of regression parameter vector is asymptotically normal and achieves the semiparametric information bound and the rate of convergence of the estimator for the cumulative baseline hazard function is optimal. Simulation studies are presented to examine the finite sample performance of the proposed estimation method. For illustration, we apply the method to a clinical dataset with current status outcome.

Empirical likelihood based on synthetic right censored data

The current status and panel count data frequently arise from cancer and tumorigenicity studies when events currently occur. A common and widely used class of two sample tests, for current status and panel count data, is the permutation class. We manipulate the double saddlepoint method to calculate the exact mid-p-values of the underlying permutation distributions of this class of tests. Permutation simulations are replaced by analytical saddlepoint computations which provide extremely accurate mid-p-values that are exact for most practical purposes and almost always more accurate than normal approximations. The method is illustrated using two real tumorigenicity panel count data. To compare the saddlepoint approximation with the normal asymptotic approximation, a simulation study is conducted. The speed and accuracy of the saddlepoint method facilitate the calculation of the confidence interval for the treatment effect. The inversion of the mid-p-values to calculate the confidence interval for the mean rate of development of the recurrent event is discussed.

In this paper, we consider the empirical likelihood confidence intervals for mean difference under the assumption of missing at random (MAR). We prove that the empirical log-likelihood ratio of mean difference converges to a standard chi-squared distribution. We conduct simulation study to compare the proposed empirical likelihood with the normal approximation method in terms of the coverage accuracies and average lengths of confidence intervals. By Monte Carlo simulations, we verify that the proposed method performs well.

Negative binomial group distribution was proposed in the literature which was motivated by inverse sampling when considering group inspection: products are inspected group by group, and the number of non-conforming items of a group is recorded only until the inspection of the whole group is finished. The non-conforming probability p of the population is thus the parameter of interest. In this paper, the confidence interval construction for this parameter is investigated. The common normal approximation and exact method are applied. To overcome the drawbacks of these commonly used methods, a composite method that is based on the confidence intervals of the negative binomial distribution is proposed, which benefits from the relationship between negative binomial distribution and negative binomial group distribution. Simulation studies are carried out to examine the performances of our methods. A real data example is also presented to illustrate the application of our method.

Kendall and Gehan estimating functions are commonly used to estimate the regression parameter in accelerated failure time model with censored observations in survival analysis. In this paper, we apply the jackknife empirical likelihood method to overcome the computation difficulty about interval estimation. A Wilks’ theorem of jackknife empirical likelihood for U-statistic type estimating equations is established, which is used to construct the confidence intervals for the regression parameter. We carry out an extensive simulation study to compare the Wald-type procedure, the empirical likelihood method, and the jackknife empirical likelihood method. The proposed jackknife empirical likelihood method has a better performance than the existing methods. We also use a real data set to compare the proposed methods.

By introducing the notion of “empirical likelihood function of observing sums”, unequally-spaced time series data and longitudinal data generated from Levy processes can be analyzed. Characteristic function is further incorporated to handle the situations where the density function of the increments is difficult to obtain. In the situations where both characteristic function and the density function are available, it is shown through the simulation examples that the proposed empirical maximum likelihood method does not suffer from significant information loss comparing to the maximum likelihood estimation method. The performances of the proposed method is tested for both equally-spaced and unequally-spaced observations.

We consider a varying-coefficient partially linear proportional odds model with current status data. This model enables one to examine the extent to which some covariates interact nonlinearly with an exposure variable, while other covariates present linear effects. B-spline approach and sieve maximum likelihood estimation method are used to get an integrated estimate for the linear coefficients, the varying-coefficient functions and the baseline function. The proposed parameter estimators are proved to be semiparametrically efficient and asymptotically normal, and the estimators for the nonparametric functions achieve the optimal rate of convergence. Simulation studies and a real data analysis are used for assessment and illustration.