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.     

Non-parametric Estimation of a Survival Function with Two-stage Design Studies

Abstract.  The two-stage design is popular in epidemiology studies and clinical trials due to its cost effectiveness. Typically, the first stage sample contains cheaper and possibly biased information, while the second stage validation sample consists of a subset of subjects with accurate and complete information. In this paper, we study estimation of a survival function with right-censored survival data from a two-stage design. A non-parametric estimator is derived by combining data from both stages. We also study its large sample properties and derive pointwise and simultaneous confidence intervals for the survival function. The proposed estimator effectively reduces the variance and finite-sample bias of the Kaplan–Meier estimator solely based on the second stage validation sample. Finally, we apply our method to a real data set from a medical device postmarketing surveillance study.     

Log-density Deconvolution by Wavelet Thresholding

Abstract.  This paper proposes a new wavelet-based method for deconvolving a density. The estimator combines the ideas of non-linear wavelet thresholding with periodized Meyer wavelets and estimation by information projection. It is guaranteed to be in the class of density functions, in particular it is positive everywhere by construction. The asymptotic optimality of the estimator is established in terms of the rate of convergence of the Kullback–Leibler discrepancy over Besov classes. Finite sample properties are investigated in detail, and show the excellent empirical performance of the estimator, compared with other recently introduced estimators.     

An unbiased estimator for the survival function of censored data

A monotonic. pointwise unbiased and uniformly consistent estimator for the survival function of failure time under the random censorship model is proposed. This estimator is closely related to the Kaplan-Meier. the Nelson-Aalen. and the reduced sample estimator. Large sample properties of the new estimator are discussed.     

Unbiased Estimation of the Distribution Function of an Exponential Population Using Order Statistics with Application in Ranked Set Sampling

In this paper we consider the problem of unbiased estimation of the distribution function of an exponential population using order statistics based on a random sample. We present a (unique) unbiased estimator based on a single, say ith, order statistic and study some properties of the estimator for i = 2. We also indicate how this estimator can be utilized to obtain unbiased estimators when a few selected order statistics are available as well as when the sample is selected following an alternative sampling procedure known as ranked set sampling. It is further proved that for a ranked set sample of size two, the proposed estimator is uniformly better than the conventional nonparametric unbiased estimator, further, for a general sample size, a modified ranked set sampling procedure provides an unbiased estimator uniformly better than the conventional nonparametric unbiased estimator based on the usual ranked set sampling procedure.     

Quantile regression and variable selection for the single-index model

In this paper, we propose a new full iteration estimation method for quantile regression (QR) of the single-index model (SIM). The asymptotic properties of the proposed estimator are derived. Furthermore, we propose a variable selection procedure for the QR of SIM by combining the estimation method with the adaptive LASSO penalized method to get sparse estimation of the index parameter. The oracle properties of the variable selection method are established. Simulations with various non-normal errors are conducted to demonstrate the finite sample performance of the estimation method and the variable selection procedure. Furthermore, we illustrate the proposed method by analyzing a real data set.     

A Mixed Model-assisted Regression Estimator that Uses Variables Employed at the Design Stage

The Generalized regression estimator (GREG) of a finite population mean or total has been shown to be asymptotically optimal when the working linear regression model upon which it is based includes variables related to the sampling design. In this paper a regression estimator assisted by a linear mixed superpopulation model is proposed. It accounts for the extra information coming from the design in the random component of the model and saves degrees of freedom in finite sample estimation. This procedure combines the larger asymptotic efficiency of the optimal estimator and the greater finite sample stability of the GREG. Design based properties of the proposed estimator are discussed and a small simulation study is conducted to explore its finite sample performance.     

Estimators for a Poisson parameter using ranked set sampling

Using ranked set sampling, a viable BLUE estimator is obtained for estimating the mean of a Poisson distribution. Its properties, such as efficiency relative to the ranked set sample mean and to the maximum likelihood estimator, have been calculated for different sample sizes and values of the Poisson parameter. The estimator (termed the normal modified r.s.s. estimator is more efficient than both the ranked set sample mean and the MLE. It is recommended as a reasonable estimator of the Poisson mean ( λ) to be used in a ranked set sampling environment.     

Estimation of non-occurrence probability in periodic model

Dimitrov and Khalil (1992) introduced a class of new probability distributions for modeling environmental evolution with periodic behavior. One of the key parameters in these distributions is α, the probability that the event being studied does not occur. In that article the authors derive an estimator for this parameter assuming a series of conditions. In this article it is shown that the estimator is valid under more general conditions, i.e. same of the assumptions are not necessary. It is shown that under the assumption that the elapsed time measured from the starting point of a period until the first occurrence time of the event given that the event occurred in this cycle is related to α, an approximate maximum likelihood estimator of a is proposed. The large sample properties of the estimator are discussed. Monte Carlo study is done for supporting the theoretical results.     

Simple non-parametric estimators for unemployment duration analysis

Summary.  We consider an extension of conventional univariate Kaplan–Meier-type estimators for the hazard rate and the survivor function to multivariate censored data with a censored random regressor. It is an Akritas-type estimator which adapts the non-parametric conditional hazard rate estimator of Beran to more typical data situations in applied analysis. We show with simulations that the estimator has nice finite sample properties and our implementation appears to be fast. As an application we estimate non-parametric conditional quantile functions with German administrative unemployment duration data.     

Double AR model without intercept: An alternative to modeling nonstationarity and heteroscedasticity

This paper presents a double AR model without intercept (DARWIN model) and provides us a new way to study the nonstationary heteroscedastic time series. It is shown that the DARWIN model is always nonstationary and heteroscedastic, and its sample properties depend on the Lyapunov exponent. An easy-to-implement estimator is proposed for the Lyapunov exponent, and it is unbiased, strongly consistent, and asymptotically normal. Based on this estimator, a powerful test is constructed for testing the ordinary oscillation of the model. Moreover, this paper proposes the quasi-maximum likelihood estimator (QMLE) for the DARWIN model, which has an explicit form. The strong consistency and asymptotic normality of the QMLE are established regardless of the sign of the Lyapunov exponent. Simulation studies are conducted to assess the performance of the estimation and testing, and an empirical example is given for illustrating the usefulness of the DARWIN model.     

Spectral Density Based Goodness-of-Fit Tests for Time Series Models

A new goodness-of-fit test for time series models is proposed. The test statistic is based on the distance between a kernel estimator of the ratio between the true and the hypothesized spectral density and the expected value of the estimator under the null. It provides a quantification of how well a parametric spectral density model fits the sample spectral density (periodogram). The asymptotic distribution of the statistic proposed is derived and its power properties are discussed. To improve upon the large sample (Gaussian) approximation of the distribution of the test statistic under the null, a bootstrap procedure is presented and justified theoretically. The finite sample performance of the test is investigated through a simulation experiment and applications to real data sets are given.     

Ratio estimators based on a ranked set sample in a finite population setting

This paper considers the ratio estimator in a finite population setting in a ranked set sampling (RSS) design, where the sample is constructed either with or without replacement policies. It is shown that the proposed ratio estimator is slightly biased, but the amount of bias is smaller than the amount of bias of a simple random sample (SRS) ratio estimator. For the proposed ratio estimator, the paper provides explicit expressions for its mean square error and precision relative to the other competing estimators. It is shown that the new estimator has a substantial amount of improvement in efficiency with respect to SRS estimator. The proposed estimator is applied to two different finite population settings to estimate population mean.     

On Small Sample Properties of the Wald,LR and LM Tests in a Linear Model with AR(1) Errors

When the error terms are autocorrelated, the conventional t-tests for individual regression coefficients mislead us to over-rejection of the null hypothesis. We examine, by Monte Carlo experiments, the small sample properties of the unrestricted estimator of ρ and of the estimator of ρ restricted by the null hypothesis. We compare the small sample properties of the Wald, likelihood ratio and Lagrange multiplier test statistics for individual regression coefficients. It is shown that when the null hypothesis is true, the unrestricted estimator of ρ is biased. It is also shown that the Lagrange multiplier test using the maximum likelihood estimator of ρ performs better than the Wald and likelihood ratio tests.     

Monotone Smoothing with Application to Dose-Response Curve

Motivated by problems that arise in dose-response curve estimation, we developed a new method to estimate a monotone curve. The resulting monotone estimator is obtained by combining techniques from smoothing splines with nonnegativity properties of cubic B-splines. Numerical experiments are given to exemplify the method.     

A Seemingly Unrelated Nonparametric Additive Model with Autoregressive Errors

This article considers a nonparametric additive seemingly unrelated regression model with autoregressive errors, and develops estimation and inference procedures for this model. Our proposed method first estimates the unknown functions by combining polynomial spline series approximations with least squares, and then uses the fitted residuals together with the smoothly clipped absolute deviation (SCAD) penalty to identify the error structure and estimate the unknown autoregressive coefficients. Based on the polynomial spline series estimator and the fitted error structure, a two-stage local polynomial improved estimator for the unknown functions of the mean is further developed. Our procedure applies a prewhitening transformation of the dependent variable, and also takes into account the contemporaneous correlations across equations. We show that the resulting estimator possesses an oracle property, and is asymptotically more efficient than estimators that neglect the autocorrelation and/or contemporaneous correlations of errors. We investigate the small sample properties of the proposed procedure in a simulation study.     

Pooling reliability functions

In this article large sample pooling procedures for reliability functions of an exponential life testing model is considered. Asymptotic properties of shrinkage estimation procedure subsequent to preliminary tests are developed. It is shown that the proposed estimator possesses substantially snakker asymptotic mean squared error than the usual estimator in a region of the lparameter space. Relative efficiencies of the purposed estimators to the usual estimators are obtained and recommendations of the level of the preliminary tests are provided. Relative dominance picture of the estimators is presented. It is shown that the proposed estimator provides a wider dominance range over usual estimator than the usual preliminary test estimator. More importantly, the size of the preliminary test is meaningful. Simulation studies is also carried out to appraise the performance of the estimators when samples are small.     

The effects of error magnitude and bandwidth selection for deconvolution with unknown error distribution

The error distribution is generally unknown in deconvolution problems with real applications. A separate independent experiment is thus often conducted to collect the additional noise data in those studies. In this paper, we study the nonparametric deconvolution estimation from a contaminated sample coupled with an additional noise sample. A ridge-based kernel deconvolution estimator is proposed and its asymptotic properties are investigated depending on the error magnitude. We then present a data-driven bandwidth selection algorithm with combining the bootstrap method and the idea of simulation extrapolation. The finite sample performance of the proposed methods and the effects of error magnitude are evaluated through simulation studies. A real data analysis for a gene Illumina BeadArray study is performed to illustrate the use of the proposed methods.     

Two Stochastic Restricted Principal Components Regression Estimator in Linear Regression

In this article, we propose two stochastic restricted principal components regression estimator by combining the approach followed in obtaining the ordinary mixed estimator and the principal components regression estimator in linear regression model. The performance of the two new estimators in terms of matrix MSE criterion is studied. We also give an example and a Monte Carlo simulation to show the theoretical results.     

Recursive location and scale estimators

It is shown that a recursive estimator with the same asymptotic properties as the median has convergence properties in finite samples which depend heavily on the scale of the data. A simple modification which adjusts for the scale is suggested and its application illustrated on simulated data. The modified estimator has much improved properties which are similar to those of the sample (non-recursive) median.