On Universal Optimality of Circular Weakly Neighbor Balanced Designs Under an Interference Model

In this article, the optimal design problem in a fixed effects interference model with left-neighbor effects is studied. It is known (Druilhet, 1999 Druilhet , P. ( 1999 ). Optimality of circular neighbor balanced designs . J. Statist. Plann. Infer. 81 : 141 – 152 .[Crossref], [Web of Science ®] , [Google Scholar]) that circular neighbor balanced designs (CNBDs) are universally optimal in such a model. We prove the universal optimality of circular weakly neighbor balanced designs (CWNBDs), which require a smaller number of blocks than CNBDs. CWNBDs with the number of blocks smaller than the number of treatments belong to the class of partially neighbor balanced designs (PNBDs) defined by Wilkinson et al. (1983 Wilkinson , G. N. , Eckert , R. , Hancock , T. W. , Mayo , O. ( 1983 ). Nearest neighbour (NN) analysis of field experiments . J. Roy. Statist. Soc. Ser. B 45 : 151 – 210 . [Google Scholar]). We give a construction method for some CWNBDs, with examples.     

Circular First- and Second-Order Balanced Repeated Measurements Designs

Sharma (1977 Sharma , V. K. ( 1977 ). Change-over designs with complete balance for first and second order residual effect . Canad. J. Statist. 5 : 121 – 132 .[Crossref] , [Google Scholar]) and Aggarwal et al. (2006 Aggarwal , M. L. , Deng , L.-Y. , Jha , M. K. ( 2006 ). Balanced residual treatment effects designs of first and second order . Statist. Probab. Lett. 76 : 597 – 600 .[Crossref], [Web of Science ®] , [Google Scholar]) considered non circular construction of first- and second-order balanced repeated measurements designs. Sharma et al. (2002 Sharma , V. K. , Varghese , C. , Jaggi , S. ( 2002 ). On optimality of change-over designs balanced for first and second order residual effects . Metron 60 : 153 – 162 . [Google Scholar]) constructed circular first- and second-order balanced repeated measurements designs only for a class with parameters (v, p = 3n, n = v ²) and also showed its universal optimality. In this article, we consider circular construction of first- and second-order balanced repeated measurements designs and strongly balanced repeated measurements designs by using the method of cyclic shifts. Some new circular designs with parameters (v, p, n) for cases p = v, p < v and p > v are given.     

Analyzing Supersaturated Designs by Means of an Information Based Criterion

The cost and time consumption of many industrial experimentations can be reduced using the class of supersaturated designs since this can be used for screening out the important factors from a large set of potentially active variables. A supersaturated design is a design for which there are fewer runs than effects to be estimated. Although there exists a wide study of construction methods for supersaturated designs, their analysis methods are yet in an early research stage. In this article, we propose a method for analyzing data using a correlation-based measure, named as symmetrical uncertainty. This method combines measures from the information theory field and is used as the main idea of variable selection algorithms developed in data mining. In this work, the symmetrical uncertainty is used from another viewpoint in order to determine more directly the important factors. The specific method enables us to use supersaturated designs for analyzing data of generalized linear models for a Bernoulli response. We evaluate our method by using some of the existing supersaturated designs, obtained according to methods proposed by Tang and Wu (1997 Tang , B. , Wu , C. F. J. (1997). A method for constructing supersaturated designs and its E(s ²)-optimality. Canadian Journal of Statistics 25:191–201.[Crossref], [Web of Science ®] , [Google Scholar]) as well as by Koukouvinos et al. (2008 Koukouvinos , C. , Mylona , K. , Simos , D. E. ( 2008 ). E(s ²)-optimal and minimax-optimal cyclic supersaturated designs via multi-objective simulated annealing . Journal of Statistical Planning and Inference 138 : 1639 – 1646 .[Crossref], [Web of Science ®] , [Google Scholar]). The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. Additionally, Receiver Operating Characteristics (ROC) curves methodology is applied as an additional statistical tool for performance evaluation.     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

Optimum Mixture Design Using Deficiency Criterion

In a mixture experiment the measured response is assumed to depend only on the relative proportion of ingredients or components present in the mixture. Scheffe (1958 Scheffe , H. ( 1958 ). Experiments with mixtures . J. R. Statist. Soc. B 20 : 344 – 360 . [Google Scholar], 1963 Scheffe , H. ( 1963 ). Simplex – centroid design for experiments with mixtures . J. R. Statist. Soc. B 25 : 235 – 263 . [Google Scholar]) first systematically considered this problem and introduced different models and designs suitable in such situations. Optimum designs for the estimation of parameters of different mixture models are available in the literature. The problem of estimating the optimum proportion of mixture components is of great practical importance. Pal and Mandal (2006 Pal , M. , Mandal , N. K. ( 2006 ). Optimum designs for optimum mixtures . Statist. Probab. Lett. 76 ( 13 ): 1369 – 1379 .[Crossref], [Web of Science ®] , [Google Scholar], 2007 Pal , M. , Mandal , N. K. (2007). Optimum mixture design via equivalence theory. communicated. JCISS 32:107–126. [Google Scholar]) attempted to find a solution to this problem by adopting a pseudo-Bayesian approach and using the trace criterion. Subsequently, Pal and Mandal (2008 Pal , M. , Mandal , N. K. ( 2008 ). Minimax designs for optimum mixtures . To appear in Statist. Probab. Lett.  [Google Scholar]) solved the problem using minimax criterion. In this article, the deficiency criterion due to Chatterjee and Mandal (1981 Chatterjee , S. K. , Mandal , N. K. ( 1981 ). Response surface designs for estimating the optimal point . Calcutta Statist. Assoc. Bull. . 30 : 145 – 169 .[Crossref] , [Google Scholar]) has been used as a measure for comparing the performance of competing designs.     

A Method for Analyzing Supersaturated Designs with a Block Orthogonal Structure

Supersaturated designs is a large class of factorial designs which can be used for screening out the important factors from a large set of potentially active variables. The huge advantage of these designs is that they reduce the experimental cost drastically, but their critical disadvantage is the confounding involved in the statistical analysis. In this article, we propose a method for analyzing data using a specific type of supersaturated designs. This method heavily uses the special block orthogonal structure of the supersaturated designs given by Tang and Wu (1997 Tang , B. , Wu , C. F. J. ( 1997 ). A method for constructing supersaturated designs and its Es ²-optimality . Canadian J. Statist. 25 : 191 – 201 .[Crossref], [Web of Science ®] , [Google Scholar]). Also, we compare our method with several known statistical analysis methods by using some of the existing supersaturated designs. The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. The results are presented in tables and the discussion to follow.     

Independent Quadratic Forms in 3-Variance-Component Models

This article addresses derivation and existence of quadratic forms that were suggested by Burch (2007 Burch , B. D. ( 2007 ). Generalized confidence intervals for proportions of total variance in mixed linear models . J. Statist. Plann. Infer. 137 : 2394 – 2404 .[Crossref], [Web of Science ®] , [Google Scholar]) for procedures for inference on variance components in mixed linear models in combination with generalized fiducial inference. A relatively simple algorithm leading to the required quadratic forms in a general 3-variance-component model is stated and designs for two-way ANOVA models without interactions that permit Burch's procedure are characterized. This complements developments in the original article by Burch.     

Inverse Sampling with Unequal Selection Probabilities

In some real situations the population of interest is divided into two groups, of which one contains only a few units. In other cases, the population may be considered as subdivided into two group', for example, if only a few units display a value of the variable of interest which is highly different from zero, while all the other units show a value equal to or near zero. In both cases, inverse sampling is more efficient than classical fixed sample-size designs to obtain the parameter estimators for the whole population as well as for its groups (e.g., Salehi and Seber, 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]). In fact, in this design the procedure selection continues until a prefixed number of units with the characteristic of interest is sampled. Since it is not known a priori to which group the population units belong, the sample size is a random variable. Christman and Lan (2001 Christman , M. C. , Lan , F. ( 2001 ). Inverse adaptive cluster sampling . Biometrics 57 : 1096 – 1105 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Salehi and Seber (2001 Salehi , M. M. , Seber , G. A. F. ( 2001 ). A new proof of Murthy's estimator which applies to sequential sampling . Austral. NZ J. Statist. 43 : 281 – 286 . [Google Scholar] 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]) considered inverse sampling designs when all the population units have equal selection probabilities. In this article, we consider the general case in which the units may have unequal probabilities of being included in the sample. In fact, in many real situations different units may have different selection probabilities because of some inherent features of the sampling procedure, or in order to obtain better estimates. We derive unbiased estimators of the totals of the two groups, their variance and the corresponding unbiased variance estimators in inverse sampling with replacement. Finally, we derive similar results for more complex designs, where the selection procedure stops before observing the prefixed number of units from the rare group.     

Circular Strongly Balanced Repeated Measurements Designs

Magda (1980 Magda , C. G. ( 1980 ). Circular balanced repeated measurements designs . Commun. Statist. Theor. Meth. 9 : 1901 – 1918 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Hedayat (1981 Hedayat , A. S. ( 1981 ). Repeated measurements designs-IV: recent advances (with discussion) . Bull. Int. Statist. Inst. 1 : 591 – 610 . [Google Scholar]) first considered the construction of circular strongly balanced repeated measurements designs. Sen and Mukerjee (1987 Sen , M. , Mukerjee , R. ( 1987 ). Optimal repeated measurements designs under interaction . J. Statist. Plann. Infer. 17 : 81 – 91 .[Crossref], [Web of Science ®] , [Google Scholar]) and Roy (1988 Roy , B. K. ( 1988 ). Construction of strongly balanced uniform repeated measurements designs . J. Statist. Plann. Infer. 19 : 341 – 348 .[Crossref], [Web of Science ®] , [Google Scholar]) considered the optimality and existence of circular strongly balanced repeated measurements designs based on the method of differences and Hamiltonian decomposition of lexicographic product of two graphs. In this article, we consider the construction of circular strongly balanced repeated measurements designs using the newly proposed method called cyclic shifts, and propose some new designs for p < v.     

Efficient Ratio and Product Estimators in Stratified Random Sampling

This paper suggests an efficient class of ratio and product estimators for estimating the population mean in stratified random sampling using auxiliary information. It is interesting to mention that, in addition to many, Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]), Kadilar and Cingi (2003 Kadilar , C. , Cingi , H. ( 2003 ). Ratio estimator in stratified sampling . Biometr. J. 45 : 218 – 225 .[Crossref], [Web of Science ®] , [Google Scholar], 2005 Kadilar , C. , Cingi , H. ( 2005 ). A new estimator in stratified random sampling . Commun. Statist. Theor. Meth. 34 : 597 – 602 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Singh and Vishwakarma (2007 Singh , H. P. , Vishwakarma , G. K. ( 2007 ). Modified exponential ratio and product estimators for finite population mean in double sampling . Austr. J. Statist. 36 ( 3 ): 217 – 225 . [Google Scholar]) estimators are identified as members of the proposed class of estimators. The expressions of bias and mean square error (MSE) of the proposed estimators are derived under large sample approximation in general form. Asymptotically optimum estimator (AOE) in the class is identified alongwith its MSE formula. It has been shown that the proposed class of estimators is more efficient than combined regression estimator and Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]) estimator. Moreover, theoretical findings are supported through a numerical example.     

Least Median of Squares Solution of Multiple Linear Regression Models Through the Origin

Barreto and Maharry (2006 Barreto , H. , Maharry , D. ( 2006 ). Least median of squares and regression through the origin . Comput. Statist. Data Anal. 50 : 1391 – 1397 .[Crossref], [Web of Science ®] , [Google Scholar]) showed that PROGRESS algorithm fails to find a correct minimum “Least Median of Squares/LMS” estimate for bivariate regression models which have no intercept. Kayhan and Gunay (2008 Kayhan , Y. , Gunay , S. ( 2008 ). A new approach to least median of squares and regression through the origin . Commun. Statist. Theor. Meth. 37 ( 5 ): 773 – 781 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) presented a different approach for the regression models through the origin which includes at most two unknown parameters. However, LMS estimate for multiple linear regression models still remains an open issue. The aim of this study is to show that finding true LMS estimate for zero intercept multiple linear regression models can be treated as a convex optimization problem and to provide a more general algorithm for any dimensional linear regression models.     

On Estimating the Finite Population Mean Using Two Auxiliary Variables in Two Phase Sampling in the Presence of Non Response

This article considers the problem of estimating the population mean using information on two auxiliary variables in the presence of non response under two-phase sampling. Some improved ratio-in-regression type estimators have been proposed in four different situations of non response along with their properties under large sample approximation. Efficiency comparisons of the proposed estimators with the usual unbiased estimator by Hansen and Hurwitz (1946 Hansen , M. H. , Hurwitz , W. N. ( 1946 ). The problem of non response in sample surveys . J. Amer. Statist. Assoc. 41 : 517 – 529 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), conventional ratio and regression estimators using single auxiliary variable and Singh and Kumar (2010b Singh , H. P. , Kumar , S. ( 2010b ). Improved estimation of population mean under double sampling with sub-sampling the non-respondents . J. Statist. Plann. Infer. 140 ( 9 ): 2536 – 2550 .[Crossref], [Web of Science ®] , [Google Scholar]) estimators using two auxiliary variables have been made. Finally, these theoretical findings are illustrated by a numerical example.     

On Rates of Convergence for Posterior Distributions Under Misspecification

We extend the approach of Walker (2003 Walker , S. ( 2003 ). On sufficient conditions for Bayesian consistency . Biometrika 90 ( 2 ): 482 – 488 .[Crossref], [Web of Science ®] , [Google Scholar]); (2004 Walker , S. ( 2004 ). New approaches to Bayesian consistency . Ann. Statist. 32 ( 5 ): 2028 – 2043 .[Crossref], [Web of Science ®] , [Google Scholar]) to the case of misspecified models. A sufficient condition for establishing rates of convergence is given based on a key identity involving martingales, which does not require construction of tests. We also show roughly that the result obtained by using tests can also be obtained by our approach, which demonstrates the potential wider applicability of this method.     

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.     

From Ordinary to Shrinkage Square-Root Estimators

Consider a skewed population. Suppose an intelligent guess could be made about an interval that contains the population mean. There may exist biased estimators with smaller mean squared error than the arithmetic mean within such an interval. This article indicates when it is advisable to shrink the arithmetic mean towards a guessed interval using root estimators. The goal is to obtain an estimator that is better near the average of natural origins. An estimator proposed. This estimator contains the Thompson (1968 Thompson , J. R. ( 1968 ). Accuracy borrowing in the estimation of the mean by shrinkage towards an interval . J. Amer. Statist. Assoc. 63 : 953 – 963 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) ordinary shrinkage estimator, the Jenkins et al. (1973 Jenkins , O. C. , Ringer , L. J. , Hartley , H. O. ( 1973 ). Root estimators . J Amer. Statist. Assoc. 68 : 414 – 419 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) square-root estimator, and the arithmetic sample mean as special cases. The bias and the mean squared error of the proposed more general estimator is compared with the three special cases. Shrinkage coefficients that yield minimum mean squared error estimators are obtained. The proposed estimator is considerably more efficient than the three special cases. This remains true for highly skewed populations. The merits of the proposed shrinkage square-root estimator are supported by the results of numerical and simulation studies.     

The Research on Two Kinds of Restricted Biased Estimators Based on Mean Squared Error Matrix

Sakall?oglu et al. (2001 Sakall?oglu , Kaç?ranlar , Akdeniz ( 2001 ). Mean squared error comparisons of some biased estimators . Commun. Statist. Theor. Meth. 30 : 347 – 361 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) dealt with the comparisons among the ridge estimator, Liu estimator, and iteration estimator. Akdeniz and Erol (2003 Akdeniz , F. , Erol , H. ( 2003 ). Mean squared error matrix comparisons of some biased estimators in linear regression . Commun. Statist. Theor. Meth. 32 : 2389 – 2413 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) have compared the (almost unbiased) generalized ridge regression estimator with the (almost unbiased) generalized Liu estimator in the matrix mean squared error sense. In this article, we study the ridge estimator and Liu estimator with respect to linear equality restriction, and establish some sufficient conditions for the superiority of the restricted ridge estimator over the restricted Liu estimator and the superiority of the restricted Liu estimator over the restricted ridge estimator under mean squared error matrix, respectively. Furthermore, we give a numerical example.     

Designs Robust Against the Presence of an Aberration in a Diallel Cross Design Setup

In this article, we find designs insensitive to the presence of an outlier in a diallel cross design setup for estimating a complete set of orthonormal contrasts among the effects of the general combining abilities of a set of parental lines. The criterion of robustness, suggested by Mandal (1989 Mandal , N. K. ( 1989 ). On robust designs . Cal. Stat. Assoc. Bull. 38 : 115 – 119 . [Google Scholar]) in block design setup and used by Biswas (2012 Biswas , A. ( 2012 ). Block designs robust against the presence of an aberration in a treatment-control setup . Commun. Statist. Theor. Meth. 41 : 920 – 933 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in treatment-control setup, is adapted here. Complete diallel cross designs, suggested by Gupta and Kageyama (1994 Gupta , S. , Kageyama , S. ( 1994 ). Optimal complete diallel crosses . Biometrika 81 : 420 – 424 .[Crossref], [Web of Science ®] , [Google Scholar]), and partial diallel cross designs, suggested by Gupta et al. (1995 Gupta , S. , Das , A. , Kageyama , S. ( 1995 ). Single replicate orthogonal block designs for circulant partial diallel crosses . Commun. Statist. A 24 : 2601 – 2607 .[Taylor & Francis Online] , [Google Scholar]) and Mukerjee (1997 Mukerjee , R. ( 1997 ). Optimal partial diallel crosses . Biometrika 84 : 939 – 948 .[Crossref], [Web of Science ®] , [Google Scholar]), are found to be robust under certain conditions.     

Data-Driven Tests for Trend

This article proposes a new nonparametric test for the ordered alternatives problem in the k-sample setting for null hypothesis of lack of trend. This article further elaborates upon and extends the results of Ledwina and Wy?upek (2012a Ledwina , T. , Wy?upek , G. ( 2012a ). Two-sample test against one-sided alternatives . Scand. J. Statist. 39 : 358 – 381 .[Crossref], [Web of Science ®] , [Google Scholar]) obtained for k = 2. Simulations show that the new test has high and stable power and is able to control the Type I error to satisfactory extent, thus solving the problem posed in Terpstra and Magel (2003 Terpstra , J. T. , Magel , R. C. ( 2003 ). A new nonparametric test for the ordered alternative problem . J. Nonparametr. Statist. 15 : 289 – 301 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Our theoretical results say that asymptotic errors of both kinds do not exceed significance level, thus implying that the test is asymptotically unbiased.     

On the Asymmetric Marcinkiewicz-Zygmund Strong Law of Large Numbers for Linear Random Fields

In this article, the asymmetric Marcinkiewicz-Zygmund strong law of large numbers for linear random field under negative association is obtained. Our result generalizes a result in Gut and Studtmüller (2009 Gut , A. , Studtmüller , U. ( 2009 ) An asymmetric Marcinkiewicz-Zygmund LLN for random fields . Statist. Probab. Lett. 79 : 1016 – 1020 .[Crossref], [Web of Science ®] , [Google Scholar]). An asymmetric Marcinkiewicz-Zygmund LLN for random fields to the linear random field by using the Beverige-Nelson decomposition.     

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.