Penalized and Shrinkage Estimation in the Cox Proportional Hazards Model

This article considers the shrinkage estimation procedure in the Cox's proportional hazards regression model when it is suspected that some of the parameters may be restricted to a subspace. We have developed the statistical properties of the shrinkage estimators including asymptotic distributional biases and risks. The shrinkage estimators have much higher relative efficiency than the classical estimator, furthermore, we consider two penalty estimators—the LASSO and adaptive LASSO—and compare their relative performance with that of the shrinkage estimators numerically. A Monte Carlo simulation experiment is conducted for different combinations of irrelevant predictors and the performance of each estimator is evaluated in terms of simulated mean squared error. Simulation study shows that the shrinkage estimators are comparable to the penalty estimators when the number of irrelevant predictors in the model is relatively large. The shrinkage and penalty methods are applied to two real data sets to illustrate the usefulness of the procedures in practice.     

Revisiting Beal's Confidence Intervals for the Difference of Two Binomial Proportions

Confidence interval construction for the difference of two independent binomial proportions is a well-known problem with a full panoply of proposed solutions. In this paper, we focus largely on the family of intervals proposed by Beal (1987 Beal , S. ( 1987 ). Asymptotic confidence intervals for the difference between two binomial parameters for use with small samples . Biometrics 43 : 941 – 950 . [CSA] [CROSSREF] [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). This family, which includes the Haldane and Jeffreys–Perks intervals as special cases, assumes a symmetric prior distribution for the population proportions p ₁ and p ₂. We propose new methods that allow the currently observed data to set the prior distribution by taking a parametric empirical-Bayes approach; in addition, we also provide an investigation of the new interval' behaviors in small-sample situations. Unlike other solutions, our intervals can be used adaptively for experiments conducted in multiple stages over time. We illustrate this notion using data from an Argentinean study involving the Mal Rio Cuarto virus and its transmission to susceptible maize crops.     

On Multi-Treatment Adaptive Allocation Design for Dichotomous Response

The use of ridit, as a probability score, is a very common practice to compare discrete random variables in discrete data analysis. In the present work we formulate ridit reliability functionals for some comparison of K independent binary random variables. We use such functionals to provide a generalized response-adaptive design (GRAD) on K(≥ +2) treatment-arms for dichotomous response variables. We exhibit some properties of the proposed design and compare it with some of the existing competitors by computing its various performance measures. We also provide a discussion towards a possible modification of the GRAD in the presence of covariates.     

Inference under Heteroscedasticity of Unknown Form Using an Adaptive Estimator

The heteroscedasticity consistent covariance matrix estimators are commonly used for the testing of regression coefficients when error terms of regression model are heteroscedastic. These estimators are based on the residuals obtained from the method of ordinary least squares and this method yields inefficient estimators in the presence of heteroscedasticity. It is usual practice to use estimated weighted least squares method or some adaptive methods to find efficient estimates of the regression parameters when the form of heteroscedasticity is unknown. But HCCM estimators are seldom derived from such efficient estimators for testing purposes in the available literature. The current article addresses the same concern and presents the weighted versions of HCCM estimators. Our numerical work uncovers the performance of these estimators and their finite sample properties in terms of interval estimation and null rejection rate.     

Estimation of rare and clustered population variance in adaptive cluster sampling

Abstract

Many researchers used auxiliary information together with survey variable to improve the efficiency of population parameters like mean, variance, total and proportion. Ratio and regression estimation are the most commonly used methods that utilized auxiliary information in different ways to get the maximum benefits in the form of high precision of the estimators. Thompson first introduced the concept of Adaptive cluster sampling, which is an appropriate technique for collecting the samples from rare and clustered populations. In this article, a generalized exponential type estimator is proposed and its properties have been studied for the estimation of rare and highly clustered population variance using single auxiliary information. A numerical study is carried out on a real and artificial population to judge the performance of the proposed estimator over the competing estimators. It is shown that the proposed generalized exponential type estimator is more efficient than the adaptive and non adaptive estimators under conventional sampling design.     

On the Use of Adaptive Nearest Neighbors for Missing Value Imputation

A popular nonparametric treatment of missing value imputation uses methods based on k-nearest neighbors, where the number k of nearest neighbors is fixed without any consideration of the local features of missing values. This article proposes an alternative imputation method based on adaptive nearest neighbors, which takes into account the local features of the data. The proposed method adapts the number of neighbors in imputing the missing values according to the location of the missing values. Efficiency evaluation is then gauged through simulation studies using both simulated and real data. It is shown that the proposed method has distinct advantages over the imputation method based on k-nearest neighbors.     

Algorithmic Inference of Two-Parameter Gamma Distribution

We provide an estimation procedure of the two-parameter Gamma distribution based on the Algorithmic Inference approach. As a key feature of this approach, we compute the joint probability distribution of these parameters without assuming any prior. To this end, we propose a numerical algorithm which is often beneficial of a highly efficient speed up based on an approximate analytical expression of the probability distribution. We contrast our interval and point estimates with those recently obtained in Son and Oh (2006 Son , Y. S. , Oh , M. ( 2006 ). Bayesian estimation of the two-parameter gamma distribution . Communications in Statistics – Simulation and Computation 35 : 285 – 293 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) for the same problem. From this benchmark we realize that our estimates are both unbiased and more accurate, albeit more dispersed, in some cases, than the competitors' methods, where the dispersion drawback is notably mitigated w.r.t. Bayesian methods by a greater estimate decorrelation. We also briefly discuss the theoretical novelty of the adopted inference paradigm which actually represents a brush up on a Fisher perspective dating to almost a century, made feasible today by the available computational tools.     

Adaptive Test for Periodicity in Autoregressive Conditional Heteroskedastic Processes

This article is concerned with the periodicity testing problem in Autoregressive Conditional Heteroskedastic (ARCH) process. Adaptive locally asymptotically optimal test is derived, when the innovation density is unspecified but symmetric satisfying only some general technical assumptions, for the null hypothesis of classical ARCH process against an alternative of periodically correlated ARCH dependence. The main technical tool is LeCam's (1960 LeCam , L. ( 1960 ). Locally Asymptotically Normal Families of Distributions . University California Publ. Statistics 3:27–98 . [Google Scholar]) Local Asymptotic Normality (LAN) property. The LAN property of the central sequence is shown via the adapted sufficient Swensen's conditions (1985 Swensen , A. R. ( 1985 ). The asymptotic distribution of the likelihood ratio for autoregressive time series with a regression trend . Journal of Multivariate Analysis 16 : 54 – 70 .[Crossref], [Web of Science ®] , [Google Scholar]). The performance of the established test is shown via simulation studies.