Empirical likelihood inference for censored median regression model via nonparametric kernel estimation

An alternative to the accelerated failure time model is to regress the median of the failure time on the covariates. In the recent years, censored median regression models have been shown to be useful for analyzing a variety of censored survival data with the robustness property. Based on missing information principle, a semiparametric inference procedure for regression parameter has been developed when censoring variable depends on continuous covariate. In order to improve the low coverage accuracy of such procedure, we apply an empirical likelihood ratio method (EL) to the model and derive the limiting distributions of the estimated and adjusted empirical likelihood ratios for the vector of regression parameter. Two kinds of EL confidence regions for the unknown vector of regression parameters are obtained accordingly. We conduct an extensive simulation study to compare the performance of the proposed methods with that normal approximation based method. The simulation results suggest that the EL methods outperform the normal approximation based method in terms of coverage probability. Finally, we make some discussions about our methods.     

Semiparametric analysis in double-sampling designs via empirical likelihood

Double-sampling designs are commonly used in real applications when it is infeasible to collect exact measurements on all variables of interest. Two samples, a primary sample on proxy measures and a validation subsample on exact measures, are available in these designs. We assume that the validation sample is drawn from the primary sample by the Bernoulli sampling with equal selection probability. An empirical likelihood based approach is proposed to estimate the parameters of interest. By allowing the number of constraints to grow as the sample size goes to infinity, the resulting maximum empirical likelihood estimator is asymptotically normal and its limiting variance-covariance matrix reaches the semiparametric efficiency bound. Moreover, the Wilks-type result of convergence to chi-squared distribution for the empirical likelihood ratio based test is established. Some simulation studies are carried out to assess the finite sample performances of the new approach.     

Empirical likelihood inference for linear transformation models

Empirical likelihood inference is developed for censored survival data under the linear transformation models, which generalize Cox's [Regression models and life tables (with Discussion), J. Roy. Statist. Soc. Ser. B 34 (1972) 187-220] proportional hazards model. We show that the limiting distribution of the empirical likelihood ratio is a weighted sum of standard chi-squared distribution. Empirical likelihood ratio tests for the regression parameters with and without covariate adjustments are also derived. Simulation studies suggest that the empirical likelihood ratio tests are more accurate (under the null hypothesis) and powerful (under the alternative hypothesis) than the normal approximation based tests of Chen et al. [Semiparametric of transformation models with censored data, Biometrika 89 (2002) 659-668] when the model is different from the proportional hazards model and the proportion of censoring is high.     

Corrected empirical likelihood inference for right-censored partially linear single-index model

This article deals with the inference on a right-censored partially linear single-index model (RCPLSIM). The main focus is the local empirical likelihood-based inference on the nonparametric part in RCPLSIM. With a synthetic data approach, an empirical log-likelihood ratio statistic for the nonparametric part is defined and it is shown that its limiting distribution is not a central chi-squared distribution. To increase the accuracy of the confidence interval, we also propose a corrected empirical log-likelihood ratio statistic for the nonparametric function. The resulting statistic is proved to follow a standard chi-squared limiting distribution. Simulation studies are undertaken to assess the finite sample performance of the proposed confidence intervals. A real example is also considered.     

Semiparametric analysis based on weighted estimating equations for transformation models with missing covariates

Missing covariate data are very common in regression analysis. In this paper, the weighted estimating equation method (Qi et al., 2005) [25] is used to extend the so-called unified estimation procedure (Chen et al., 2002) [4] for linear transformation models to the case of missing covariates. The non-missingness probability is estimated nonparametrically by the kernel smoothing technique. Under missing at random, the proposed estimators are shown to be consistent and asymptotically normal, with the asymptotic variance estimated consistently by the usual plug-in method. Moreover, the proposed estimators are more efficient than the weighted estimators with the inverse of true non-missingness probability as weight. Finite sample performance of the estimators is examined via simulation and a real dataset is analyzed to illustrate the proposed methods.     

Estimation in generalized linear models for functional data via penalized likelihood

We analyze in a regression setting the link between a scalar response and a functional predictor by means of a Functional Generalized Linear Model. We first give a theoretical framework and then discuss identifiability of the model. The functional coefficient of the model is estimated via penalized likelihood with spline approximation. The L² rate of convergence of this estimator is given under smoothness assumption on the functional coefficient. Heuristic arguments show how these rates may be improved for some particular frameworks.     

Empirical likelihood inference in partially linear single-index models for longitudinal data

The empirical likelihood method is especially useful for constructing confidence intervals or regions of parameters of interest. Yet, the technique cannot be directly applied to partially linear single-index models for longitudinal data due to the within-subject correlation. In this paper, a bias-corrected block empirical likelihood (BCBEL) method is suggested to study the models by accounting for the within-subject correlation. BCBEL shares some desired features: unlike any normal approximation based method for confidence region, the estimation of parameters with the iterative algorithm is avoided and a consistent estimator of the asymptotic covariance matrix is not needed. Because of bias correction, the BCBEL ratio is asymptotically chi-squared, and hence it can be directly used to construct confidence regions of the parameters without any extra Monte Carlo approximation that is needed when bias correction is not applied. The proposed method can naturally be applied to deal with pure single-index models and partially linear models for longitudinal data. Some simulation studies are carried out and an example in epidemiology is given for illustration.     

Bartlett correctable two-sample adjusted empirical likelihood

We propose a two-sample adjusted empirical likelihood (AEL) to construct confidence regions for the difference of two d-dimensional population means. This method eliminates the non-definition of the usual two-sample empirical likelihood (EL) and is shown to be Bartlett correctable. We further show that when the adjustment level is half the Bartlett correction factor for the usual two-sample EL, the two-sample AEL has the same high-order precision as the EL with Bartlett correction. To enhance the performance of the two-sample AEL with adjustment level being half the Bartlett correction factor, we propose a less biased estimate of the Bartlett correction factor. The efficiency of the proposed method is illustrated by simulations and a real data example.     

Nonparametric inference under competing risks and selection-biased sampling

The aim of this paper is to carry out statistical inference in a competing risks setup when only selection-biased observation of the data of interest is available. We introduce estimators of the cumulative incidence functions and study their joint large sample behavior.     

Empirical likelihood for single-index models

The empirical likelihood method is especially useful for constructing confidence intervals or regions of the parameter of interest. This method has been extensively applied to linear regression and generalized linear regression models. In this paper, the empirical likelihood method for single-index regression models is studied. An estimated empirical log-likelihood approach to construct the confidence region of the regression parameter is developed. An adjusted empirical log-likelihood ratio is proved to be asymptotically standard chi-square. A simulation study indicates that compared with a normal approximation-based approach, the proposed method described herein works better in terms of coverage probabilities and areas (lengths) of confidence regions (intervals).     

Asymptotic properties of the maximum likelihood estimator for the proportional hazards model with doubly censored data

Doubly censored data, which include left as well as right censored observations, are frequently met in practice. Though estimation of the distribution function with doubly censored data has seen much study, relatively little is known about the inference of regression coefficients in the proportional hazards model for doubly censored data. In particular, theoretical properties of the maximum likelihood estimator of the regression coefficients in the proportional hazards model have not been proved yet. In this paper, we show the consistency and asymptotic normality of the maximum likelihood estimator and prove its semiparametric efficiency. The proposed methods are illustrated with simulation studies and analysis of an application from a medical study.     

Uniform convergence of nonparametric regressions in competing risk models with right censoring

We consider, in the presence of covariates, non-independent competing risks that are subject to right censoring. We define a nonparametric estimator of the incident regression function through the generalized product-limit estimator of the conditional censorship distribution function. Under suitable conditions, we establish the almost sure uniform convergence of those estimators with an appropriate rate.     

Sieve maximum likelihood estimation for doubly semiparametric zero-inflated Poisson models

For nonnegative measurements such as income or sick days, zero counts often have special status. Furthermore, the incidence of zero counts is often greater than expected for the Poisson model. This article considers a doubly semiparametric zero-inflated Poisson model to fit data of this type, which assumes two partially linear link functions in both the mean of the Poisson component and the probability of zero. We study a sieve maximum likelihood estimator for both the regression parameters and the nonparametric functions. We show, under routine conditions, that the estimators are strongly consistent. Moreover, the parameter estimators are asymptotically normal and first order efficient, while the nonparametric components achieve the optimal convergence rates. Simulation studies suggest that the extra flexibility inherent from the doubly semiparametric model is gained with little loss in statistical efficiency. We also illustrate our approach with a dataset from a public health study.     

Statistical inference for the index parameter in single-index models

In this paper, we are concerned with statistical inference for the index parameter in the single-index model . Based on the estimates obtained by the local linear method, we extend the generalized likelihood ratio test to the single-index model. We investigate the asymptotic behaviour of the proposed test and demonstrate that its limiting null distribution follows a χ²-distribution, with the scale constant and the number of degrees of freedom being independent of nuisance parameters or functions, which is called the Wilks phenomenon. A simulated example is used to illustrate the performance of the testing approach.     

Empirical likelihood for linear models under m -dependent errors

In this paper, the empirical likelihood confidence regions for the regression coefficient in a linear model are constructed under m-dependent errors. It is shown that the blockwise empirical likelihood is a good way to deal with dependent samples. Partly supported by the National Natural Science Foundation of China and the SF of Guangxi Normal University.     

Blockwise empirical likelihood for time series of counts

Time series of counts have a wide variety of applications in real life. Analyzing time series of counts requires accommodations for serial dependence, discreteness, and overdispersion of data. In this paper, we extend blockwise empirical likelihood (Kitamura, 1997 [15]) to the analysis of time series of counts under a regression setting. In particular, our contribution is the extension of Kitamura’s (1997) [15] method to the analysis of nonstationary time series. Serial dependence among observations is treated nonparametrically using a blocking technique; and overdispersion in count data is accommodated by the specification of a variance-mean relationship. We establish consistency and asymptotic normality of the maximum blockwise empirical likelihood estimator. Simulation studies show that our method has a good finite sample performance. The method is also illustrated by analyzing two real data sets: monthly counts of poliomyelitis cases in the USA and daily counts of non-accidental deaths in Toronto, Canada.     

An empirical likelihood approach to data analysis under two-stage sampling designs

A new empirical likelihood approach is developed to analyze data from two-stage sampling designs, in which a primary sample of rough or proxy measures for the variables of interest and a validation subsample of exact information are available. The validation sample is assumed to be a simple random subsample from the primary one. The proposed empirical likelihood approach is capable of utilizing all the information from both the specific models and the two available samples flexibly. It maintains some nice features of the empirical likelihood method and improves the asymptotic efficiency of the existing inferential procedures. The asymptotic properties are derived for the new approach. Some numerical studies are carried out to assess the finite sample performance.     

Likelihood and conditional likelihood inference for generalized additive mixed models for clustered data

Lin and Zhang (J. Roy. Statist. Soc. Ser. B 61 (1999) 381) proposed the generalized additive mixed model (GAMM) as a framework for analysis of correlated data, where normally distributed random effects are used to account for correlation in the data, and proposed to use double penalized quasi-likelihood (DPQL) to estimate the nonparametric functions in the model and marginal likelihood to estimate the smoothing parameters and variance components simultaneously. However, the normal distributional assumption for the random effects may not be realistic in many applications, and it is unclear how violation of this assumption affects ensuing inferences for GAMMs. For a particular class of GAMMs, we propose a conditional estimation procedure built on a conditional likelihood for the response given a sufficient statistic for the random effect, treating the random effect as a nuisance parameter, which thus should be robust to its distribution. In extensive simulation studies, we assess performance of this estimator under a range of conditions and use it as a basis for comparison to DPQL to evaluate the impact of violation of the normality assumption. The procedure is illustrated with application to data from the Multicenter AIDS Cohort Study (MACS).     

Empirical likelihood for linear models under negatively associated errors

In this paper, we discuss the construction of the confidence intervals for the regression vector β in a linear model under negatively associated errors. It is shown that the blockwise empirical likelihood (EL) ratio statistic for β is asymptotically χ²-type distributed. The result is used to obtain an EL based confidence region for β.     

Bias-corrected empirical likelihood in a multi-link semiparametric model

In this paper, we investigate the empirical likelihood for constructing a confidence region of the parameter of interest in a multi-link semiparametric model when an infinite-dimensional nuisance parameter exists. The new model covers the commonly used varying coefficient, generalized linear, single-index, multi-index, hazard regression models and their generalizations, as its special cases. Because of the existence of the infinite-dimensional nuisance parameter, the classical empirical likelihood with plug-in estimation cannot be asymptotically distribution-free, and the existing bias correction is not extendable to handle such a general model. We then propose a link-based correction approach to solve this problem. This approach gives a general rule of bias correction via an inner link, and consists of two parts. For the model whose estimating equation contains the score functions that are easy to estimate, we use a centering for the scores to correct the bias; for the model of which the score functions are of complex structure, a bias-correction procedure using simpler functions instead of the scores is given without loss of asymptotic efficiency. The resulting empirical likelihood shares the desired features: it has a chi-square limit and, under-smoothing technique, high order kernel and parameter estimation are not needed. Simulation studies are carried out to examine the performance of the new method.