A Simple Corrected Score for Logistic Regression with Errors-in-Covariates

We develop a simple corrected score for logistic regression with errors-in-covariates. The new method is an extension of the consistent functional methods proposed by Huang and Wang (2001) and is closely related to the corrected score method by Nakamura (1990 Nakamura, T. (1990). Corrected score function for errors-in-variables models: Methodology and application to generalized linear models. Biometrika. 77:127–137.[Crossref], [Web of Science ®] , [Google Scholar]) and Stefanski (1989 Stefanski, L.A. (1989). Unbiased estimation of a nonlinear function a normal mean with application to measurement error models. Commun. Stat. Ser. A - Theory Methods. 18:4335–4358.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The new method requires that the measurement error distribution is known, but does not require normality. The new method yields a consistent and asymptotically normal estimator under regularity conditions. We examine the finite-sample performance of the new estimator through simulation studies. Finally, we illustrate the new method by applying it to an AIDS study.     

Generalized linear regression models incorporating original outcome distributions

ABSTRACT

In this article, we propose an approach for incorporating continuous and discrete original outcome distributions into the usual exponential family regression models. The new approach is an extension of the works of Suissa (1991 Suissa, S. (1991). Binary methods for continuous outcomes: A parametric alternative. J. Clin. Epidemiol. 44:241–248.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Suissa and Blais (1995 Suissa, S., Blais, L. (1995). Binary regression with continuous outcomes. Stat. Med. 14:247–255.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which present methods to estimate the risk of an event defined in a sample subspace of an original continuous outcome variable. Simulation studies are presented in order to illustrate the performance of the developed methodology. Real data sets are analyzed by using the proposed models.     

More Efficient Tests Robust to Heteroskedasticity of Unknown Form

ABSTRACT

In the presence of heteroskedasticity of unknown form, the Ordinary Least Squares parameter estimator becomes inefficient, and its covariance matrix estimator inconsistent. Eicker (1963 Eicker , B. ( 1963 ). Limit theorems for regression with unequal and dependant errors . Ann. Math. Statist. 34 : 447 – 456 .[Crossref] , [Google Scholar]) and White (1980 White , H. ( 1980 ). A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity . Econometrica 48 : 817 – 838 .[Crossref], [Web of Science ®] , [Google Scholar]) were the first to propose a robust consistent covariance matrix estimator, that permits asymptotically correct inference. This estimator is widely used in practice. Cragg (1983 Cragg , J. G. ( 1983 ). More efficient estimation in the presence of heteroskedasticity of unknown form . Econometrica 51 : 751 – 63 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a more efficient estimator, but concluded that tests basd on it are unreliable. Thus, this last estimator has not been used in practice. This article is concerned with finite sample properties of tests robust to heteroskedasticity of unknown form. Our results suggest that reliable and more efficient tests can be obtained with the Cragg estimators in small samples.     

Robustness of the Affine Equivariant Scatter Estimator Based on the Spatial Rank Covariance Matrix

Visuri et al. (2000 Visuri, S., Koivunen, V., Oja, H. (2000). Sign and rank covariance matrices. J. Stat. Plann. Inference 91:557–575.[Crossref], [Web of Science ®] , [Google Scholar]) proposed a technique for robust covariance matrix estimation based on different notions of multivariate sign and rank. Among them, the spatial rank based covariance matrix estimator that utilizes a robust scale estimator is especially appealing due to its high robustness, computational ease, and good efficiency. Also, it is orthogonally equivariant under any distribution and affinely equivariant under elliptically symmetric distributions. In this paper, we study robustness properties of the estimator with respective to two measures: breakdown point and influence function. More specifically, the upper bound of the finite sample breakdown point can be achieved by a proper choice of univariate robust scale estimator. The influence functions for eigenvalues and eigenvectors of the estimator are derived. They are found to be bounded under some assumptions. Moreover, finite sample efficiency comparisons to popular robust MCD, M, and S estimators are reported.     

On goodness-of-fit tests for Aalen's additive model based on left-truncated right-censored or doubly censored data

ABSTRACT

Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) proposed goodness-of-fit tests for Aalen's additive risk model. In this article, we demonstrate that the approach of Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) can be applied to left-truncated right-censored (LTRC) data and doubly censored data. A simulation study is conducted to investigate the performance of the proposed tests. The proposed tests are illustrated using heart transplant data.     

MULTI-LEVEL FACTORIAL DESIGNS WITH MINIMUM NUMBERS OF LEVEL CHANGES

The order of experimental runs in a fractional factorial experiment is essential when the cost of level changes in factors is considered. The generalized foldover scheme given by [1] Coster, D. C. and Cheng, C. S. 1988. Minimum cost trend free run orders of fractional factorial designs. The Annals of Statistics, 16: 1188–1205. [Crossref], [Web of Science ®] , [Google Scholar]gives an optimal order to experimental runs in an experiment with specified defining contrasts. An experiment can be specified by a design requirement such as resolution or estimation of some interactions. To meet such a requirement, we can find several sets of defining contrasts. Applying the generalized foldover scheme to these sets of defining contrasts, we obtain designs with different numbers of level changes and then the design with minimum number of level changes. The difficulty is to find all the sets of defining contrasts. An alternative approach is investigated by [2] Cheng, C. S., Martin, R. J. and Tang, B. 1998. Two-level factorial designs with extreme numbers of level changes. The Annals of Statistics, 26: 1522–1539. [Crossref], [Web of Science ®] , [Google Scholar]for two-level fractional factorial experiments. In this paper, we investigate experiments with all factors in slevels.     

On Augmented Franses Tests for Seasonal Unit Roots

This article extends the results reported in del Barrio Castro, Osborn and Taylor (2012 del Barrio Castro, T., Osborn, D.R., Taylor, A. M.R. (2012). On augmented HEGY tests for seasonal unit roots. Econometric Theor. 18:1121–1143.[Crossref], [Web of Science ®] , [Google Scholar]) to the approach followed by Franses (1991a Franses, P. H. (1991a). Model selection and seasonality in time series. Tibergen Institute Series, 18. [Google Scholar],b Franses, P.H. (1991b). Seasonality, non-stationarity and the forecasting of monthly time series. Int. J. Forecast. 7:199–208.[Crossref], [Web of Science ®] , [Google Scholar]) to test for seasonal unit roots, providing the asymptotic representation to the seasonal unit roots tests proposed by Franses for a general number of seasons S.     

Restricted estimation and testing of hypothesis in linear measurement errors models

ABSTRACT

In this article, the linear models with measurement error both in the response and in the covariates are considered. Following Shalabh et al. (2007 Shalabh, Garg, G., Misra, N. (2007). Restricted regression estimation in measurement error models. Comput. Stat. Data Anal. 52:1149–1166.[Crossref], [Web of Science ®] , [Google Scholar], 2009 Shalabh, Garg, G., Misra, N. (2009). Use of prior information in the consistent estimation of regression coefficients in measurement error models. J. Multivariate Anal. 100:1498–1520.[Crossref], [Web of Science ®] , [Google Scholar]), we propose several restricted estimators for the regression coefficients. The consistency and asymptotic normality of the restricted estimators are established. Furthermore, we also discuss the superiority of the restricted estimators to unrestricted estimators under Pitman closeness criterion. We also develop several variance estimators and establish their asymptotic distributions. Wald-type statistics are constructed for testing the linear restrictions. Finally, Monte Carlo simulations are conducted to illustrate the finite-sample properties of the proposed estimators.     

A New Estimator of Population Mean in Stratified Sampling

Kadilar and Cingi (2005 Kadilar , C. , Cingi , H. ( 2005 ). A new ratio estimator in stratified sampling . Comm. Statist. Theory Meth. 34 : 1 – 6 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) have suggested a new ratio estimator in stratified sampling. The efficiency of this estimator is compared with the traditional combined ratio estimator on the basis of mean square error (MSE). We propose another estimator by utilizing a simple transformation introduced by Bedi (1996 Bedi , P. K. ( 1996 ). Efficient utilization of auxiliary information at estimation stage . Biomet. J. 38 ( 8 ): 973 – 976 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]). The proposed estimator is found to be more efficient than the traditional combined ratio estimator as well as the Kadilar and Cingi (2005 Kadilar , C. , Cingi , H. ( 2005 ). A new ratio estimator in stratified sampling . Comm. Statist. Theory Meth. 34 : 1 – 6 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) ratio estimator.     

COMBINING THE LIU ESTIMATOR AND THE PRINCIPAL COMPONENT REGRESSION ESTIMATOR

In this paper we introduce a class of estimators which includes the ordinary least squares (OLS), the principal components regression (PCR) and the Liu estimator [1] Liu, K. 1993. A new class of biased estimate in linear regression. Communications in Statistics – Theory and Methods, 22(2): 393–402. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]. In particular, we show that our new estimator is superior, in the scalar mean-squared error (mse) sense, to the Liu estimator, to the OLS estimator and to the PCR estimator.     

Autocovariance Function of the Fractionally Integrated Separable Spatial ARMA (FISSARMA) Models

Spatial modeling is important in many fields and there are various kinds of spatial models. One of such models is known as the fractionally integrated separable spatial ARMA (FISSARMA) model. In the area of time series analysis, Sowell (1992 Sowell, F. (1992). Maximum likelihood estimation of stationary univariate fractionally integrated time series models. J. Econ. 53:165–188.[Crossref], [Web of Science ®] , [Google Scholar]) has established the autocovariance function of the long-memory models using hypergeometric function. In this paper we will extend Sowell’s work for FISSARMA models.     

EVALUATION OF POWER OF TESTS INVOLVING COVARIATES UNDER NONCONSTANT VARIANCE

In many experiments where pre-treatment and post-treatment measurements are taken, investigators wish to determine if there is a difference between two treatment groups. For this type of data, the post-treatment variable is used as the primary comparison variable and the pre-treatment variable is used as a covariate. Although most of the discussion in this paper is written with the pre-treatment variable as the covariate the results are applicable to other choices of a covariate. Tests based on residuals have been proposed as alternatives to the usual covariance methods. Our objective is to investigate how the powers of these tests are affected when the conditional variance of the post-treatment variable depends on the magnitude of the pre-treatment variable. In particular, we investigate two cases. [1] Crager, Michael R. 1987. Analysis of Covariance in Parallel-Group Clinical Trials With Pretreatment Baselines. Biometrics, 43: 895–901. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] The conditional variance of the post-treatment variable gradually increases as the magnitude of the pre-treatment variable increases. (In many biological models this is the case.) [2] Knoke, James D. 1991. Nonparametric Analysis of Covariance for Comparing Change in Randomized Studies with Baseline Values Subject to Error. Biometrics, 47: 523–533. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] The conditional variance of the post-treatment variable is dependent upon natural or imposed subgroups contained within the pre-treatment variable. Power comparisons are made using Monte Carlo techniques.     

Asymptotics of improved generalized moment estimators for spatial autoregressive error models

Abstract

This article considers linear models with a spatial autoregressive error structure. Extending Arnold and Wied (2010) Arnold, M., Wied, D. (2010). Improved GMM estimation of the spatial autoregressive error model. Econ. Lett. 108:65–68.[Crossref], [Web of Science ®] , [Google Scholar], who develop an improved generalized method of moment (GMM) estimator for the parameters of the disturbance process to reduce the bias of existing estimation approaches, we establish the asymptotic normality of a new weighted version of this improved estimator and derive the efficient weighting matrix. We also show that this efficiently weighted GMM estimator is feasible as long as the regression matrix of the underlying linear model is non stochastic and illustrate the performance of the new estimator by a Monte Carlo simulation and an application to real data.     

The Exact General Formulas for the Moments of a Ridge Regression Estimator when the Regression Error Terms Follow a Multivariate t Distribution

Huang (1999 Huang , J. C. ( 1999 ). Improving the estimation precision for a selected parameter in multiple regression analysis: an algebraic approach . Econ. Lett. 62 : 261 – 264 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a feasible ridge regression (FRR) estimator to estimate a specific regression coefficient. Assuming that the error terms follow a normal distribution, Huang (1999 Huang , J. C. ( 1999 ). Improving the estimation precision for a selected parameter in multiple regression analysis: an algebraic approach . Econ. Lett. 62 : 261 – 264 .[Crossref], [Web of Science ®] , [Google Scholar]) examined the small sample properties of the FRR estimator. In this article, assuming that the error terms follow a multivariate t distribution, we derive an exact general formula for the moments of the FRR estimator to estimate a specific regression coefficient. Using the exact general formula, we obtain exact formulas for the bias, mean squared error (MSE), skewness, and kurtosis of the FRR estimator. Since these formulas are very complex, we compare the bias, MSE, skewness, and kurtosis of the FRR estimator with those of ordinary least square (OLS) estimator by numerical evaluations. Our numerical results show that the range of MSE dominance of the FRR estimator over the OLS estimator is widen under a fat tail distributional assumption.     

A Selection Procedure for Selecting the Least Increasing Failure Rate Average Distribution

For k independent absolutely continuous increasing failure rate average (IFRA) life distributions F_i, i = 1, 2, …, k, Link (1989 Link, W.A. (1989). Testing for exponentiality against monotone failure rate average alternatives. Commun. Statist. Theor. Meth. 18(8): 3009–3017.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered a measure of departure of against monotone failure rate average alternatives. In this paper, we use the measure defined by Link for detection of IFRA-ness of life distribution F_i. A two-stage selection procedure to select the least IFRA distribution is proposed. This selection procedure is based on a U-statistic which is an estimator of the measure and can be implemented even when the IFRA life distributions belong to different families. The applications of this procedure are discussed for some well known distributions.     

On the Additive Weibull Distribution

We derive explicit expressions for the moments, incomplete moments, quantile function and generating function of the additive Weibull model pioneered by Xie and Lai (1995 Xie, M., Lai, C.D. (1995). Reliability analysis using an additive Weibull model with bathtub-shaped failure rate function. Reliab. Eng. Syst. Safety 52:87–93.[Crossref], [Web of Science ®] , [Google Scholar]), which is a quite flexible distribution for fitting lifetime data with bathtub-shaped failure rate function. In addition, we estimate the model parameters by maximum likelihood and determine the observed information matrix. The flexibility of the additive Weibull distribution is illustrated by means of one application to real data.     

Nonparametric Estimators of the Distribution Function for One Modified Model of Current Status Data

In this article, we consider the estimation of distribution function for one modified form of current status data. An inverse-probability-weighted (IPW) estimator and a self-consistent estimator (SCE) are proposed. The asymptotic properties of the IPW estimator are derived. A simulation study is conducted to compare the performances among the IPW estimator, SCE, and the product-limit estimator proposed by Patilea and Rolin (2006 Patilea , V. , Rolin , J.-M. (2006). Product-limit estimators of the survival function for two modified forms of current-status data. Bernoulli 12(5):801–819.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results indicate that when right censoring is light and left censoring is heavy, both IPW estimator and SCE can outperform the product-limit estimator. The performances of the IPW estimator and SCE are close to each other.     

Assessing One-Step-Ahead Prediction Error Based on the Median for First-Order Autoregressive Models in the Presence Of Outliers

The prediction of the one-step-ahead observation of the first-order autoregressive process in the presence of outliers is considered. The mean square of the prediction error is obtained based on the median estimator of the model parameter for a stationary process. Monte Carlo simulation methods are employed to investigate the performance of the proposed estimator as well as the conventional ordinary least squares estimators proposed by Zhang and Shaman (1995 Zhang , P. , Shaman , P. ( 1995 ). Assessing prediction error in autoregressive models . Trans. Amer. Mathemat. Soc. 347 : 627 – 637 .[Crossref], [Web of Science ®] , [Google Scholar]) and Kabaila and He (1999 Kabaila , P. , He , Z. ( 1999 ). On assessing prediction error in autoregressive models . J. Time Ser. Anal. 20 : 663 – 670 .[Crossref] , [Google Scholar]) for a process without outliers. The results show that the proposed method outperforms the conventional method. These conclusions are substantiated with results from actual datasets.     

A two-stage design for multivariate estimation of proportions

ABSTRACT

In this article, the problem of estimation of the several proportions of two inter-dependent sensitive characteristics prevailing in a given population is considered. A new two-stage randomized response model has been proposed by extending Lee et al.’s (2013 Lee, C., Sedory, S., Singh, S (2013). Estimation at least seven measures of qualitative variables from a single sample using randomized response technique. Stat. Probab. Lett. 83(1):399–409.[Crossref], [Web of Science ®] , [Google Scholar]) model. The expressions for biases, variances, and mean square errors of different estimators have been provided. The relative efficiency comparisons of the proposed estimators has been made, numerically by considering the different practicable choices of the design parameters. It was observed that the proposed extended version performs efficiently than simple and crossed models of Lee et al. (2013 Lee, C., Sedory, S., Singh, S (2013). Estimation at least seven measures of qualitative variables from a single sample using randomized response technique. Stat. Probab. Lett. 83(1):399–409.[Crossref], [Web of Science ®] , [Google Scholar]).     

Semiparametric analysis of transformation models with doubly censored data

Double censoring arises when T represents an outcome variable that can only be accurately measured within a certain range, [L, U], where L and U are the left- and right-censoring variables, respectively. In this note, using Martingale arguments of Chen et al. [3 Chen, K., Jin, Z. and Ying, Z. 2002. Semiparametric analysis of transformation models with censored data. Biometrika, 89: 659–668. [Crossref], [Web of Science ®] [Google Scholar]], we propose an estimator (denoted by ?β) for estimating regression coefficients of transformation model when L is always observed. Under Cox proportional hazards model, the proposed estimator is equivalent to the partial likelihood estimator for left-truncated and right-censored data if the left-censoring variables L were regarded as left-truncated variables. In this case, the estimator ?β can be obtained by the standard software. A simulation study is conducted to investigate the performance of ?β. For the purpose of comparison, the simulation study also includes the estimator proposed by Cai and Cheng [2 Cai, T. and Cheng, S. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] [Google Scholar]] for the case when L and U are always observed.