An empirical comparison of several methods for testing the equality of dependent proportions

Three methods for testing the equality of nonindependent proportions were compared with, the use of Monte Carlo techniques. The three methods included Cochran's test, an ANOVA F test, and Hotelling's T² test. With respect to empirical significance levels, the ANOVA F test is recommended as the preferred method of analysis.

Oftentimes an experimenter is interested in testing the equality of several proportions. When the proportions are independent Kemp and Butcher (1972) and Butcher and Kemp (1974) compared several methods for analysing large sample binomial data for the case of a 3 x 3 factorial design without replication. In addition, Levy and Narula (1977) compared many of the same methods for analyzing binomial data; however, Levy and Narula investigated the relative utility of the methods for small sample sizes.     

Tests for homogeneity of variance

Six procedures which convert tests of homogeneity of variance into tests for mean equality for independent groups are compared. The tests are the analysis of variance (ANOVA) and Welch F statistics. The Welch statistics are included since it was anticipated that ANOVA would not provide a robust test when samples of unequal sizes are obtained from non-normal populations. However, the Welch tests are not found to be uniformly preferrable. In addition, a prior recommendation for Miller's jackknife procedure is not supported for the unequal sample size case. The data indicates that the current tests for variance heterogeneity are either sensitive to non-normality or, if robust, lacking in power. Therefore, these tests cannot be recommended for the purpose of testing the validity of the ANOVA homogeneity assumption.     

A Comparison of Significance Testing Procedures for the Intraclass Correlation from Family Data

Different procedures for testing problems concerning intraclass correlation from familial data are considered in the case of varying number of siblings per family. Under the assumption of multivariate normality, the hypotheses that the intraclass correlation is equal to a specified value are tested. To assess the performance of the tests, Monte Carlo simulations are designed to compare their powers. The Neyman's (1959) C(α) test and the test based on the modified ANOVA F statistic are shown to be consistently more powerful than other procedures.     

The Effects of Misconceptions on the Properties of Friedman's Test

Friedman's test is a widely used rank-based alternative to the analysis of variance (ANOVA) F-test for identifying treatment differences in a randomized complete block design. Many texts provide incomplete or misleading information about when Friedman's test may be appropriately applied. We discuss the assumptions needed for the test and common misconceptions. We show via simulation that when the variance or skew of the treatment distributions differ, application of Friedman's test to detect differences in treatment location can result in Type I error probabilities larger than the nominal α, and even when α is unaffected, the power of the test can be less than expected.     

Robust estimation and testing in one-way ANOVA for Type II censored samples: skew normal error terms

Censoring can be occurred in many statistical analyses in the framework of experimental design. In this study, we estimate the model parameters in one-way ANOVA under Type II censoring. We assume that the distribution of the error terms is Azzalini's skew normal. We use Tiku's modified maximum likelihood (MML) methodology which is a modified version of the well-known maximum likelihood (ML) in the estimation procedure. Unlike ML methodology, MML methodology is non-iterative and gives explicit estimators of the model parameters. We also propose new test statistics based on the proposed estimators. The performances of the proposed estimators and the test statistics based on them are compared with the corresponding normal theory results via Monte Carlo simulation study. A real life data is analysed to show the implementation of the methodology presented in this paper at the end of the study.     

Some comparative studies on testing parallelism of several straight lines under heteroscedastic variances

This paper examines the robustness of the Welch test, the James test as well as Tan's ANOVA test (to be referred as F_β test) for testing parallelism in k straight lines under heteroscedasticity and nonnormality. Results of Monte Carlo studies demonstrate the robustness of all tests with respect to departure from normality. Further, there is hardly any difference between these methods with respect to both power and size of the test.     

The bias of the mle of box's degrees of freedom correction factor for correlation in anova

Approximate expressions for the mean and variance of the MLE of Box's 2-way ANOVA degrees of freedom factor are given for the case when the usual F test for column effects is appropriate even though there is correlation across columns. A simulation study is performed showing the approximations are good for a variety of numbers of rows and columns. The results confirm that the estimated degrees of freedom factor is badly biased when no factor at all should be used.     

Modification of Tukey's additivity test

In this paper we discuss testing for an interaction in the two-way ANOVA with just one observation per cell. The known results are reviewed and a simulation study is performed to evaluate type I and type II risks of the tests. It is shown that the Tukey and Mandel additivity tests have very low power in case of more general interaction scheme. A modification of Tukey's test is developed to resolve this issue. All tests mentioned in the paper have been implemented in R package Additivity Tests.     

A weighted test of internal symmetry

This study examines extensions of McNemar's Test with multinomial responses, and proposes a linear weighting scheme, based on the distance of the response change, that is applied to one of these extensions (Bowker's test). This weighted version of Bowker's test is then appropriate for ordinal response variables. A Monte Carlo simulation was conducted to examine the Type I error rate of the weighted Bowker's test for a cross-classification table based on a five-category ordinal response scale. The weighted Bowker's test was also applied to a data set involving change in student attitudes towards mathematics. The results of the weighted Bowker's test were compared with the results of Bowker's test applied to the same set of data.     

An ANOVA-type test for multiple change points

We consider the problem of testing the null hypothesis of no change against the alternative of multiple change points in a series of independent observations. We propose an ANOVA-type test statistic and obtain its asymptotic null distribution. We also give approximations of its limiting critical values. We report the results of Monte Carlo studies conducted to compare the power of the proposed test against a number of its competitors. As illustrations we analyzed three real data sets.     

A fiducial-based approach to the one-way ANOVA in the presence of nonnormality and heterogeneous error variances

In this study, we propose a new test for testing the equality of the treatment means in one-way ANOVA when the usual normality and the homogeneity of variances assumptions are not met. In developing the proposed test, we benefit from the Fisher's fiducial inference [1–3]. Distribution of the error terms is assumed to be long-tailed symmetric (LTS) which includes the normal distribution as a limiting case. Modified maximum likelihood (MML) estimators are used in the test statistics rather than the traditional least squares (LS) estimators, since LS estimators have very low efficiencies under nonnormal distributions, see Tiku [4] for the details of MML methodology. An extensive Monte Carlo simulation study is done to compare the efficiency of the proposed test with the corresponding test based on normal theory, see Li et al. [5]. Finally, we give a real life example to show the applicability of the proposed methodology.     

Powers of Discrete Goodness-of-Fit Test Statistics for a Uniform Null Against a Selection of Alternative Distributions

The comparative powers of six discrete goodness-of-fit test statistics for a uniform null distribution against a variety of fully specified alternative distributions are discussed. The results suggest that the test statistics based on the empirical distribution function for ordinal data (Kolmogorov–Smirnov, Cramér–von Mises, and Anderson–Darling) are generally more powerful for trend alternative distributions. The test statistics for nominal (Pearson's chi-square and the nominal Kolmogorov–Smirnov) and circular data (Watson's test statistic) are shown to be generally more powerful for the investigated triangular (∨), flat (or platykurtic type), sharp (or leptokurtic type), and bimodal alternative distributions.     

Robust tests for the equality of variances for clustered data

Tests for the equality of variances are often needed in applications. In genetic studies the assumption of equal variances of continuous traits, measured in identical and fraternal twins, is crucial for heritability analysis. To test the equality of variances of traits, which are non-normally distributed, Levene [H. Levene, Robust tests for equality of variances, in Contributions to Probability and Statistics, I. Olkin, ed. Stanford University Press, Palo Alto, California, 1960, pp. 278–292] suggested a method that was surprisingly robust under non-normality, and the procedure was further improved by Brown and Forsythe [M.B. Brown and A.B. Forsythe, Robust tests for the equality of variances, J. Amer. Statis. Assoc. 69 (1974), pp. 364–367]. These tests assumed independence of observations. However, twin data are clustered – observations within a twin pair may be dependent due to shared genes and environmental factors. Uncritical application of the tests of Brown and Forsythe to clustered data may result in much higher than nominal Type I error probabilities. To deal with clustering we developed an extended version of Levene's test, where the ANOVA step is replaced with a regression analysis followed by a Wald-type test based on a clustered version of the robust Huber–White sandwich estimator of the covariance matrix. We studied the properties of our procedure using simulated non-normal clustered data and obtained Type I error rates close to nominal as well as reasonable powers. We also applied our method to oral glucose tolerance test data obtained from a twin study of the metabolic syndrome and related components and compared the results with those produced by the traditional approaches.     

Generalisation of Mood's Test to Factorial Designs to Detect Heterogeneity of Variances

Mood's test, which is a relatively old test (and the oldest non‐parametric test among those tests in its class) for determining heterogeneity of variance, is still being widely used in different areas such as biometry, biostatistics and medicine. Although it is a popular test, it is not suitable for use on a two‐way factorial design. In this paper, Mood's test is generalised to the 2 × 2 factorial design setting and its performance is compared with that of Klotz's test. The power and robustness of these tests are examined in detail by means of a simulation study with 10,000 replications. Based on the simulation results, the generalised Mood's and Klotz's tests can especially be recommended in settings in which the parent distribution is symmetric. As an example application we analyse data from a multi‐factor agricultural system that involves chilli peppers, nematodes and yellow nutsedge. This example dataset suggests that the performance of the generalised Mood test is in agreement with that of the generalised Klotz's test.     

The analysis of multicentre clinical trials when there is heterogeneity between centres

Two-treatment multicentre clinical trials are very common in practice. In cases where a non-parametric analysis is appropriate, a rank-sum test for grouped data called the van Elteren test can be applied. As an alternative approach, one may apply a combination test such as Fisher's combination test or the inverse normal combination test (also called Liptak's method) in order to combine centre-specific P-values. If there are no ties and no differences between centres with regard to the groups’ sample sizes, the inverse normal combination test using centre-specific Wilcoxon rank-sum tests is equivalent to the van Elteren test. In this paper, the van Elteren test is compared with Fisher's combination test based on Wilcoxon rank-sum tests. Data from two multicentre trials as well as simulated data indicate that Fisher's combination of P-values is more powerful than the van Elteren test in realistic scenarios, i.e. when there are large differences between the centres’ P-values, some quantitative interaction between treatment and centre, and/or heterogeneity in variability. The combination approach opens the possibility of using statistics other than the rank sum, and it is also a suitable method for more complicated designs, e.g. when covariates such as age or gender are included in the analysis.     

Analysis of ranked data in randomized blocks when there are missing values

Data consisting of ranks within blocks are considered for randomized block designs when there are missing values. Tied ranks are possible. Such data can be analysed using the Skillings–Mack test. Here we suggest a new approach based on carrying out an ANOVA on the ranks using the general linear model platform available in many statistical packages. Such a platform allows an ANOVA to be calculated when there are missing values. Indicative sizes and powers show the ANOVA approach performs better than the Skillings–Mack test.     

On a goodness-of-fit test for normality with unknown parameters and type-II censored data

We propose a new goodness-of-fit test for normal and lognormal distributions with unknown parameters and type-II censored data. This test is a generalization of Michael's test for censored samples, which is based on the empirical distribution and a variance stabilizing transformation. We estimate the parameters of the model by using maximum likelihood and Gupta's methods. The quantiles of the distribution of the test statistic under the null hypothesis are obtained through Monte Carlo simulations. The power of the proposed test is estimated and compared to that of the Kolmogorov–Smirnov test also using simulations. The new test is more powerful than the Kolmogorov–Smirnov test in most of the studied cases. Acceptance regions for the PP, QQ and Michael's stabilized probability plots are derived, making it possible to visualize which data contribute to the decision of rejecting the null hypothesis. Finally, an illustrative example is presented.     

Empirical study of some non-parametric tests for dispersion of correlated data

Results of a computer simulation study of power and robustness of three competitor tests for comparing scales, for use with correlated data: Rothstein, Richardson and Bell (RRB), Arvesen, and Pitman, are presented. It is found that unless one could ímprove the approximate null distributions for Arvesen's and Pitman's test, RRB's procedure is best, having simulated probabilities of Type I error closest to the test's nominal α and being reasonably robust and powerful, for all distributions considered.     

Bayes factor consistency for unbalanced ANOVA models

In practical situations, most experimental designs often yield unbalanced data which have different numbers of observations per unit because of cost constraints, missing data, etc. In this paper, we consider the Bayesian approach to hypothesis testing or model selection under the one-way unbalanced fixed-effects analysis-of-variance (ANOVA) model. We adopt Zellner's g-prior with the beta-prime distribution for g, which results in an explicit closed-form expression of the Bayes factor without integral representation. Furthermore, we investigate the model selection consistency of the Bayes factor under three different asymptotic scenarios: either the number of units goes to infinity, the number of observations per unit goes to infinity, or both go to infinity. The results presented extend some existing ones of the Bayes factor for the balanced ANOVA models in the literature.     

Multivariate analysis of variance with additional data from an unknown subset of the populations

In this paper we consider the problem of testing the means of k multivariate normal populations with additional data from an unknown subset of the k populations. The purpose of this research is to offer test procedures utilizing all the available data for the multivariate analysis of variance problem because the additional data may contain valuable information about the parameters of the k populations. The standard procedure uses only the data from identified populations. We provide a test using all available data based upon Hotelling' s generalized T²statistic. The power of this test is computed using Betz's approximation of Hotelling' s generalized T²statistic by an F-distribution. A comparison of the power of the test and the standard test procedure is also given.