Second-Order Covariance Matrix Formula for Heteroskedastic Generalized Linear Models

For the first time, we provide a matrix formula for second-order covariances of maximum likelihood estimates in heteroskedastic generalized linear models, thus generalizing the results of Cordeiro (2004 Cordeiro , G. M. ( 2004 ). Second-order covariance matrix of maximum likelihood estimates in generalized linear models . Statist. Probab. Lett. 66 : 153 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]) and Cordeiro et al. (2006 Cordeiro , G. M. , Barroso , L. P. , Botter , D. A. (2006). Covariance matrix formula for generalized linear models with unknown dispersion. Commun. Statist. Theor. Meth. 35:113–120.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) related to the generalized linear models with known and unknown dispersion parameter, respectively. The covariance matrix formula does not involve cumulants of log-likelihood derivatives and can be easily obtained using simple matrix operations. We apply our main result to a simple model. Some simulations show that the second-order covariances can be quite pronounced in small to moderate samples. The usual covariances of the maximum likelihood estimates can be corrected by these second-order covariances.     

Bias Correction in Generalized Nonlinear Models with Dispersion Covariates

We derive general formulae for the second-order biases of maximum likelihood estimates of the parameters in generalized nonlinear models with dispersion covariates. This result generalizes previous work by Botter and Cordeiro (1998 Botter , D. A. , Cordeiro , G. M. ( 1998 ). Improved estimates for generalized linear models with dispersion covariates . J. Statist. Comput. Simul. 62 : 91 – 104 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Cordeiro and McCullagh (1991 Cordeiro , G. M. , McCullagh , P. ( 1991 ). Bias correction in generalized linear models . J. Roy Statist. Soc. B 53 : 629 – 643 . [Google Scholar]). The practical use of such bias corrections is illustrated in a simulation study.     

A Computer Program to Improve LR Tests for Generalized Linear Models

Generalized linear models enable the fitting of models to a wide range of data types. These models are based on exponential dispersion distributions. Improved likelihood ratio tests for these models were developed by Cordeiro (1983 Cordeiro , G. M. (1983). Improved likelihood ratio statistics for generalized linear models. Journal of the Royal Statistical Society, Series B: Methodological 45:404–413. [Google Scholar])Cordeiro (1987 Cordeiro , G. M. ( 1987 ). On the corrections to the likelihood ratio statistics . Biometrika 74 : 265 – 274 .[Crossref], [Web of Science ®] , [Google Scholar]). We present a simple R program source for calculating Bartlett corrections to improve likelihood ratio tests in these models. The program was tested on some special models, confirming all of the previously reported numerical results for the Bartlett corrections.     

The Identification of Potential Longitudinal Biomarkers and Measurements of Effectiveness for Biomarkers as Surrogates in Multivariate Survival Data

A proposed method based on frailty models is used to identify longitudinal biomarkers or surrogates for a multivariate survival. This method is an extention of earlier models by Wulfsohn and Tsiatis (1997 Wulfsohn , M. S. , Tsiatis , A. A. ( 1997 ). A joint model for survival and longitudinal data measured with error . Biometrics 53 ( 1 ): 330 – 339 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Song et al. (2002 Song , X. , Davidian , M. , Tsiatis , A. A. ( 2002 ). A Semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data . Biometrics 58 ( 4 ): 742 – 753 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, similar to Henderson et al. (2002 Henderson , R. , Diggle , P. J. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), a joint likelihood function combines the likelihood functions of the longitudinal biomarkers and the multivariate survival times. We use simulations to explore how the number of individuals, the number of time points per individual and the functional form of the random effects from the longitudianl biomarkers influence the power to detect the association of a longitudinal biomarker and the multivariate survival time. The proposed method is illustrate by using the gastric cancer data.     

Generalized Scrambling in Quantitative Optional Randomized Response Models

Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an optional randomized response model using a linear combination scrambling which is a generalization of the multiplicative scrambling of Eichhorn and Hayre (1983 Eichhorn , B. H. , Hayre , L. S. ( 1983 ). Scrambled randomized response methods for obtaining sensitive quantitative data . J. Statist. Plann. Infer. 7 : 307 – 316 .[Crossref], [Web of Science ®] , [Google Scholar]) and the additive scrambling of Gupta et al. (2006, 2010). In this article, we discuss two main issues. (1) Can the Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) model be improved further by using a two-stage approach?; (2) Does the linear combination scrambling provide any benefit over the additive scrambling of Gupta et al. (2010 Gupta , S. N. , Shabbir , J. , Sehra , S. ( 2010 ). Mean and sensitivity estimation in optional randomized response models . J. Statist. Plann. Infer. 140 : 2870 – 2874 .[Crossref], [Web of Science ®] , [Google Scholar])? We will note that the answer to the first question is “yes” but the answer to the second question is “no.”     

Bias-Corrected Maximum Likelihood Estimators in Nonlinear Heteroscedastic Models

Nonlinear heteroscedastic models are widely used in econometrics and statistical applications. We derive matrix formulae for the second-order biases of the maximum likelihood estimators of the parameters in the mean and variance response which generalize previous results by Cook et al. (1986 Cook , D. R. , Tsai , C. L. , Wei , B. C. ( 1986 ). Bias in nonlinear regression . Biometrika 73 : 615 – 623 .[Crossref], [Web of Science ®] , [Google Scholar]) and Cordeiro (1993 Cordeiro , G. M. ( 1993 ). Bartlett corrections and bias correction for two heteroscedastic regression models . Commun. Statist. Theor. Meth. 22 : 169 – 188 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The biases of the estimators are easily obtained as vectors of regression coefficients from suitable weighted linear regressions. The practical use of such biases is illustrated in a simulation study and in an application to a real data set.     

Covariance Matrix Formula for Exponential Family Nonlinear Models

This article gives a matrix formula for second-order covariances of maximum likelihood estimators in exponential family nonlinear models, thus generalizing the result of Cordeiro (2004 Cordeiro , G. M. ( 2004 ). Second-order covariance matrix of maximum likelihood estimates in generalized linear models . Statist. Probab. Lett. 66 : 153 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]) valid for generalized linear models with known dispersion parameter. Some simulations show that the second-order covariances for exponential family nonlinear models can be quite pronounced in small to moderate sample sizes.     

A side sensitive modified group runs control chart to detect shifts in the process mean

Gadre and Rattihalli [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] have introduced the Modified Group Runs (MGR) control chart to identify the increases in fraction non-conforming and to detect shifts in the process mean. The MGR chart reduces the out-of-control average time-to-signal (ATS), as compared with most of the well-known control charts. In this article, we develop the Side Sensitive Modified Group Runs (SSMGR) chart to detect shifts in the process mean. With the help of numerical examples, it is illustrated that the SSMGR chart performs better than the Shewhart's X¯ chart, the synthetic chart [12 Wu, Z. and Spedding, T. A. 2000. A synthetic control chart for detecting small shifts in the process mean. J. Qual. Technol., 32: 32–38. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], the Group Runs chart [4 Gadre, M. P. and Rattihalli, R. N. 2004. A group runs control chart for detecting shifts in the process mean. Econ. Qual. Control, 19: 29–43. [Crossref] [Google Scholar]], the Side Sensitive Group Runs chart [6 Gadre, M. P. and Rattihalli, R. N. 2007. A side sensitive group runs control chart for detecting shifts in the process mean. Stat. Methods Appl., 16: 27–37. [Crossref], [Web of Science ®] [Google Scholar]], as well as the MGR chart [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]]. In some situations it is also superior to the Cumulative Sum chart p9 Page, E. S. 1954. Continuous inspection schemes. Biometrika, 41: 100–114. [Crossref], [Web of Science ®] [Google Scholar]] and the exponentially weighed moving average chart [10 Roberts, S. W. 1959. Control chart tests based on geometric moving averages. Technometrics, 1: 239–250. [Taylor & Francis Online] [Google Scholar]]. In the steady state also, its performance is better than the above charts.     

Hazards Regression for Cancer Studies

The complication in analyzing tumor data is that the tumors detected in a screening program tend to be slowly progressive tumors, which is the so-called length-biased sampling that is inherent in screening studies. Under the assumption that all subjects have the same tumor growth function, Ghosh (2008 Ghosh , D. ( 2008 ). Proportional hazards regression for cancer studies . Biometrics 64 : 141 – 148 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) developed estimation procedures for proportional hazards model. In this article, by modeling growth function as a function of covariates, we demonstrate that Ghosh (2008 Ghosh , D. ( 2008 ). Proportional hazards regression for cancer studies . Biometrics 64 : 141 – 148 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])'s approach can be extended to the case when each subject has a specific growth function. A simulation study is conducted to demonstrate the potential usefulness of the proposed estimators for the regression parameters in the proportional and additive hazards model.     

Robustness of the EWMA control chart for individual observations

The traditional exponentially weighted moving average (EWMA) chart is one of the most popular control charts used in practice today. The in-control robustness is the key to the proper design and implementation of any control chart, lack of which can render its out-of-control shift detection capability almost meaningless. To this end, Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] studied the performance of the traditional EWMA chart for the mean for i.i.d. data. We use a more extensive simulation study to further investigate the in-control robustness (to non-normality) of the three different EWMA designs studied by Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]]. Our study includes a much wider collection of non-normal distributions including light- and heavy-tailed and symmetric and asymmetric bi-modal as well as the contaminated normal, which is particularly useful to study the effects of outliers. Also, we consider two separate cases: (i) when the process mean and standard deviation are both known and (ii) when they are both unknown and estimated from an in-control Phase I sample. In addition, unlike in the study done by Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], the average run-length (ARL) is not used as the sole performance measure in our study, we consider the standard deviation of the run-length (SDRL), the median run-length (MDRL), and the first and the third quartiles as well as the first and the 99th percentiles of the in-control run-length distribution for a better overall assessment of the traditional EWMA chart's in-control performance. Our findings sound a cautionary note to the (over) use of the EWMA chart in practice, at least with some types of non-normal data. A summary and recommendations are provided.     

Computing the Exact Fisher Information Matrix of Periodic State-Space Models

The purpose of this article is to develop algorithms for computing the exact Fisher information matrix of periodic time-varying state-space models. We first present a relatively simple recursive algorithm which computes the elements of the exact information matrix without involving numerical differentiation, since all required derivatives are analytically evaluated. The proposed algorithm extends the procedure due to Cavanaugh and Shumway (1996 Cavanaugh , J. E. , Shumway , R. H. ( 1996 ). On computing the expected Fisher information matrix for state-space model parameters . Statist. Probab. Lett. 26 : 347 – 355 .[Crossref], [Web of Science ®] , [Google Scholar]) to the periodic state-space framework. Exploiting the approach used in Klein et al. (2000 Klein , A. , Mélard , G. , Zahaf , T. ( 2000 ). Construction of the exact Fisher information matrix of Gaussian time series models by means of matrix differential rules . Linear Alg. Applic. 321 : 209 – 232 .[Crossref], [Web of Science ®] , [Google Scholar]), a second algorithm is proposed in order to obtain the exact information matrix as a whole instead of element by element. The algorithms are first developed in a general framework and then specialized to the case of a periodic Gaussian vector autoregressive moving-average (PVARMA) model.     

A Note on Importance Resampling for Multi-Dimensional Statistics

Johns (1988 Johns , M. V. (1988). Importance sampling for bootstrap confidence intervals. Journal of the American Statistical Association 83:709–714.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Davison (1988 Davison , A. C. ( 1988 ). Discussion of paper by D. V. Hinkley . Journal of the Royal Statistical Society Series B 50 : 356 – 357 . [Google Scholar]), and Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]) used importance sampling for calculating bootstrap distributions of one-dimensional statistics. Realizing that their methods can not be extended easily to multi-dimensional statistics, Fuh and Hu (2004 Fuh , C. D. , Hu , I. ( 2004 ). Efficient importance sampling for events of moderate deviations with applications . Biometrika 91 : 471 – 490 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an exponential tilting formula for statistics of multi-dimension, which is optimal in the sense that the asymptotic variance is minimized for estimating tail probabilities of asymptotically normal statistics. For one-dimensional statistics, Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multi-step variance minimization approach that can be viewed as a generalization of the two-step variance minimization approach proposed by Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we generalize the approach of Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) to multi-dimensional statistics, which applies to general statistics and does not resort to asymptotics. Empirical results on a real survival data set show that the proposed algorithm provides significant computational efficiency gains.     

Exact Inference for the pth-Quantile and the Reliability of the Two-Parameter Exponential Distribution with Singly Type II Censoring: A Standard Approach

Starting from a standard pivot, exact inference for the pth-quantile and for the reliability of the two-parameter exponential distribution in case of singly Type II censored samples is developed in this article. Fernandez (2007 Fernandez , A. J. ( 2007 ). On calculating generalized confidence intervals for the two-parameter exponential reliability function . Statistics 41 : 129 – 135 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first obtained some of the results proposed in this article, but, differently from what are proposed here, and developed his theory starting from a generalized pivot. An illustrative example shows that, with the expressions proposed in this article, it is also possible to overcome some shortcomings raising from the formulas by Fernandez (2007 Fernandez , A. J. ( 2007 ). On calculating generalized confidence intervals for the two-parameter exponential reliability function . Statistics 41 : 129 – 135 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Finally, a new expression for the moments of the pivot is obtained.     

On Testing Sample Selection Bias Under the Multicollinearity Problem

ABSTRACT

This paper reviews and extends the literature on the finite sample behavior of tests for sample selection bias. Monte Carlo results show that, when the “multicollinearity problem” identified by Nawata (1993 Nawata , K. ( 1993 ). A note on the estimation of models with sample-selection biases . Economics Letters 42 : 15 – 24 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) is severe, (i) the t-test based on the Heckman–Greene variance estimator can be unreliable, (ii) the Likelihood Ratio test remains powerful, and (iii) nonnormality can be interpreted as severe sample selection bias by Maximum Likelihood methods, leading to negative Wald statistics. We also confirm previous findings (Leung and Yu, 1996 Leung , S. F. , Yu , S. ( 1996 ). On the choice between sample selection and two-part models . Journal of Econometrics 72 : 197 – 229 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) that the standard regression-based t-test (Heckman, 1979 Heckman , J. J. ( 1979 ). Sample selection bias as a specification error . Econometrica 47 : 153 – 161 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]) and the asymptotically efficient Lagrange Multiplier test (Melino, 1982 Melino , A. ( 1982 ). Testing for sample selection bias . Review of Economic Studies 49 : 151 – 153 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]), are robust to nonnormality but have very little power.     

Performance Evaluation of Different Tests for Location Parameters

In this article, our objective is to evaluate the performance of different tests which are used to compare the equality of more than two location parameters. We have considered six tests (including some commonly used) in this study, one of which is parametric and the others are nonparametric. These tests include the usual F test (Fisher and Mackenzie, 1923 Fisher , R. A. , Mackenzie , M. A. ( 1923 ). Studies in crop variation. II. The manurial response of different potato . Journal of Agricultural Science 13 : 311 – 320 .[Crossref], [Web of Science ®] , [Google Scholar]), Kruskal–Wallis test (Kruskall and Wallis, 1952 Kruskall , W. H. , Wallis , W. A. ( 1952 ). Use of ranks in one-criterion variance analysis . Journal of American Statistical Association 47 : 583 – 621 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Kolmogorov–Smirnov test (David, 1958 David , H. T. ( 1958 ). Three-sample Kolnogorov–Smirnov test . The Annals of Mathematical Statistics 29 : 842 – 851 .[Crossref] , [Google Scholar]), the g test (Stekler, 1987 Stekler , H. O. ( 1987 ). Who forecasters better? Journal of Business and Economic Statistics 5 : 155 – 158 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), f test (Batchelor, 1990 Batchelor , R. A. ( 1990 ). All forecasters are equal . Journal of Business and Economic Statistics 8 : 143 – 144 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Extension of Median test (as given in Daniel, 1990 Daniel , W. W. ( 1990 ). Applied Nonparametric Statistics. , 2nd ed. Duxbury Classic Series , Boston . [Google Scholar]). Performance of these tests are compared under different symmetric, skewed and contaminated probability distributions that include Normal, Cauchy, Uniform, Laplace, Lognormal, Exponential, Weibull, Gamma, t, Chi-square, Half Normal, Mixed Weibull, and Mixed Normal. Performances of these tests are measured in terms of power. We have suggested appropriate tests which may perform better under different situations. It is expected that researchers will find these results useful in decision making.     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

An Improved Solution for Multiple Comparisons with a Control in Response Surface Methodology with Two Predictor Variables

Sa and Edwards (1993 Sa , P. , Edwards , D. ( 1993 ). Multiple comparisons with a control in response surface methodology . Technometrics 35 ( 4 ): 436 – 445 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first proposed the Multiple Comparisons with a Control problem in Response Surface Methodology. They provided an exact solution for one predictor variable and a conservative solution when number of predictor variables is more than one. Merchant et al. (1998 Merchant , A. , McCann , M. , Edwards , D. ( 1998 ). Improved multiple comparisons with a control in response surface analysis . Technometrics 40 ( 4 ): 291 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]) improved the solution for the latter case. This article improves Merchant et al.'s solution for the case of rotatable designs in two predictor variables.     

GMM in linear regression for longitudinal data with multiple covariates measured with error

Griliches and Hausman 5 Griliches, Z. and Hausman, J. A. 1986. Errors in variables in panel data. J. Econometrics, 32: 93–118. [Crossref], [Web of Science ®] [Google Scholar] and Wansbeek 11 Wansbeek, T. J. 2001. GMM estimation in panel data models with measurement error. J. Econometrics, 104: 259–268. [Crossref], [Web of Science ®] [Google Scholar] proposed using the generalized method of moments (GMM) to obtain consistent estimators in linear regression models for longitudinal data with measurement error in one covariate, without requiring additional validation or replicate data. For usefulness of this methodology, we must extend it to the more realistic situation where more than one covariate are measured with error. Such an extension is not straightforward, since measurement errors across different covariates may be correlated. By a careful construction of the measurement error correlation structure, we are able to extend Wansbeek's GMM and show that the extended Griliches and Hausman's GMM is equivalent to the extended Wansbeek's GMM. For illustration, we apply the extended GMM to data from two medical studies, and compare it with the naive method and the method assuming only one covariate having measurement error.     

A Nonparametric Test for Symmetry Based on Freeman and Halton's Ideas on Contingency Tables

In this article, we propose a nonparametric method to test for symmetry in bivariate data. By using the extension of Fisher's exact treatment for 2 × 2 contingency tables proposed by Freeman and Halton (1951 Freeman , G. H. , Halton , J. H. ( 1951 ). Note on an exact treatment of contingency tables, goodness of fit and other problems of significance . Biometrika 38 : 141 – 149 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), we can test the hypothesis of equal distribution for two samples of integer valued variables. Then, by counting the number of observations belonging to each cell of a symmetric, appropriately built grid, we can produce the two samples of integers required to use this test for equal distribution. The resulting test for symmetry is potentially extendible to higher dimensions. A simulation study is performed to compare with some known tests (Bowker, 1948 Bowker , A. H. ( 1948 ). A test for symmetry in contingency tables . Journal of the American Statistical Association 43 : 572 – 574 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Hollander, 1971 Hollander , M. ( 1971 ). A nonparametric test for bivariate symmetry . Biometrika 58-1 : 203 – 212 .[Crossref], [Web of Science ®] , [Google Scholar]; and its improvement given in Krampe and Kuhnt, 2007 Krampe , A. , Kuhnt , S. ( 2007 ). Bowker's test for symmetry and modifications within the algebraic framework . Computational Statistics and Data Analysis 51 : 4124 – 4142 .[Crossref], [Web of Science ®] , [Google Scholar]). Our proposal represents a competitive option as a test for symmetry.     

D-Optimal Chemical Balance Weighing Designs with N ≡ 0 (mod 4) and 3 Objects

In this article, the problem of estimation of the individual weights of three objects using a chemical balance weighing design is considered. We use the criterion of D-optimality. We assume that the covariance matrix of errors is the matrix of first-order autoregressive process. Such problems were discussed in Li and Yang (2005 Li , C. H. , Yang , S. Y. ( 2005 ). On a conjecture in D-optimal designs with n ≡ 0 (mod 4) . Lin. Alg. Applic. 400 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]) and also in Yeh and Lo Huang (2005 Yeh , H. G. , Lo Huang , M. N. ( 2005 ). On exact D-optimal designs with 2 two-level factors and n autocorrelated observations . Metrika 61 : 261 – 275 .[Crossref], [Web of Science ®] , [Google Scholar]). We present some results of D-optimal designs in certain class of designs with the design matrix X  ∈ M _n×3(±1) such that each column of matrix X has at least one 1 and one ?1.