New tests for stochastic order with application to case control studies

It is shown that if a binary regression function is increasing then retrospective sampling induces a stochastic ordering of the covariate distributions among the responders, which we call cases, and the non-responders, which we call controls. We also show that if the covariate distributions are stochastically ordered then the regression function must be increasing. This means that testing whether the regression function is monotone is equivalent to testing whether the covariate distributions are stochastically ordered. Capitalizing on these new probabilistic observations we proceed to develop two new non-parametric tests for stochastic order. The new tests are based on either the maximally selected, or integrated, chi-bar statistic of order one. The tests are easy to compute and interpret and their large sampling distributions are easily found. Numerical comparisons show that they compare favorably with existing methods in both small and large samples. We emphasize that the new tests are applicable to any testing problem involving two stochastically ordered distributions.     

A new jackknife test for homogeneity of variances

A new jackknife test is proposed to test the equality of variances in several populations. The new test is based on jackknifing one group of observations at a time, instead of one observation in each group as recommended by Miller for a two sample case, and by Layard for several samples. The proposed test is examined, and compared with other tests, in terms of power and robustness with respect to a wide variety of non-normal distributions. It is found that the new test is robust and has reasonably high power for normal as well as for non-normal observations, irrespective of the sample size. Furthermore, the proposed test is certainly superior to all other tests considered here in small to moderate size samples, and is as good as or better than the other tests in large samples, irrespective of the distribution of sampling observations.     

Small Sample Tests for Shape Parameters of Gamma Distributions

The introduction of shape parameters into statistical distributions provided flexible models that produced better fit to experimental data. The Weibull and gamma families are prime examples wherein shape parameters produce more reliable statistical models than standard exponential models in lifetime studies. In the presence of many independent gamma populations, one may test equality (or homogeneity) of shape parameters. In this article, we develop two tests for testing shape parameters of gamma distributions using chi-square distributions, stochastic majorization, and Schur convexity. The first one tests hypotheses on the shape parameter of a single gamma distribution. We numerically examine the performance of this test and find that it controls Type I error rate for small samples. To compare shape parameters of a set of independent gamma populations, we develop a test that is unbiased in the sense of Schur convexity. These tests are motivated by the need to have simple, easy to use tests and accurate procedures in case of small samples. We illustrate the new tests using three real datasets taken from engineering and environmental science. In addition, we investigate the Bayes’ factor in this context and conclude that for small samples, the frequentist approach performs better than the Bayesian approach.     

Testing stochastic orders in tails of contingency tables

Testing for the difference in the strength of bivariate association in two independent contingency tables is an important issue that finds applications in various disciplines. Currently, many of the commonly used tests are based on single-index measures of association. More specifically, one obtains single-index measurements of association from two tables and compares them based on asymptotic theory. Although they are usually easy to understand and use, often much of the information contained in the data is lost with single-index measures. Accordingly, they fail to fully capture the association in the data. To remedy this shortcoming, we introduce a new summary statistic measuring various types of association in a contingency table. Based on this new summary statistic, we propose a likelihood ratio test comparing the strength of association in two independent contingency tables. The proposed test examines the stochastic order between summary statistics. We derive its asymptotic null distribution and demonstrate that the least favorable distributions are chi-bar distributions. We numerically compare the power of the proposed test to that of the tests based on single-index measures. Finally, we provide two examples illustrating the new summary statistics and the related tests.     

Some linear and nonlinear rank tests for competing risks models

Using the methods of asymptotic decision theory asymptotically optimal for translation and scale families as well as for certian nonparmetric families. Moreover, two new classes of nonlinear rank tests are introduced. These tests are designed for detecting either “ omnibus alternatives ” or “ one sided alternatives of trend ”. Under the null hypothesis of randomness all tests are distribution - free. The asymptotic distributions of the test statistics are derived under contiguous alternatives.     

Robust analysis of variance

SUMMARY For the c -sample location problem with equal and unequal variances, we compare the classical F -test and its robustified version-the Welch test-with some nonparametric counterparts defined for two-sided and one-sided ordered alternatives, such as trend and umbrella alternatives. A new rank test for long-tailed distributions is proposed. The comparison is referred to level alpha and power beta of the tests, and is carried out via Monte Carlo simulation, assuming short-, medium- and long-tailed as well as asymmetric distributions. It turns out that the Welch test is the best one in the case of unequal variances but in the case of equal variances special non-parametric tests are to prefer.     

Customizing Generalizations of the Kolmogorov-Smirnov Goodness-of-fit Test

This paper presents a method of customizing goodness-of-fit tests that transforms the empirical distribution function in such a way as to create tests for certain alternatives. Using the @ , g transform described in Blom(1958), one can create non-parametric tests for an assortment of alternative distributions. As examples, three new ( f , g )-corrected Kolmogorov-Smirnov tests for goodness-of-fit are discussed. One of these tests is powerful for testing whether or not the data come from an alternative that is heavier in the tails. Another test identifies whether or not the data come from an alternative which is heavier in the middle of the distribution. The last test identifies if the data come from an alternative in which the first or third quartile is far from the corresponding quartile of the hypothesized distribution. The behavior of the three new tests is investigated through a power study.     

Conditional tests in a competing risks model

Testing for equality of competing risks based on their cumulative incidence functions (CIFs) or their cause specific hazard rates (CSHRs) has been considered by many authors. The finite sample distributions of the existing test statistics are in general complicated and the use of their asymptotic distributions can lead to conservative tests. In this paper we show how to perform some of these tests using the conditional distributions of their corresponding test statistics instead (conditional on the observed data). The resulting conditional tests are initially developed for the case of k = 2 and are then extended to k > 2 by performing a sequence of two sample tests and by combining several risks into one. A simulation study to compare the powers of several tests based on their conditional and asymptotic distributions shows that using conditional tests leads to a gain in power. A real life example is also discussed to show how to implement such conditional tests.     

NONPARAMETRIC TWO SAMPLE TESTS BASED ON AREA STATISTICS

ABSTRACT

Area statistics are sample versions of areas occurring in a probability plot of two distribution functions F and G. This paper presents a unified basis for five statistics of this type. They can be used for various testing problems in the framework of the two sample problem for independent observations, such as testing equality of distributions against inequality or testing stochastic dominance of distributions in one or either direction against nondominance. Though three of the statistics considered have already been suggested in literature, two of them are new and deserve our interest. The finite sample distributions of the statistics (under F=G) can be calculated via recursion formulae. Two tables with critical values of the new statistics are included. The asymptotic distribution of the properly normalized versions of the area statistics are functionals of the Brownian bridge. The distribution functions and quantiles thereof are obtained by Monte Carlo simulation. Finally, the power functions of the two new tests based on area statistics are compared to the power functions of the tests based on the corresponding supremum statistics, i.e., statistics of the Kolmogorov–Smirnov type.     

Inference for comparing a multinomial distribution with a known standard

We propose a new type of stochastic ordering which imposes a monotone tendency in differences between one multinomial probability and a known standard one. An estimation procedure is proposed for the constrained maximum likelihood estimate, and then the asymptotic null distribution is derived for the likelihood ratio test statistic for testing equality of two multinomial distributions against the new stochastic ordering. An alternative test is also discussed based on Neyman modified minimum chi-square estimator. These tests are illustrated with a set of heart disease data.     

An exact distribution-free test comparing two multivariate distributions based on adjacency

Summary.  A new test is proposed comparing two multivariate distributions by using distances between observations. Unlike earlier tests using interpoint distances, the new test statistic has a known exact distribution and is exactly distribution free. The interpoint distances are used to construct an optimal non-bipartite matching, i.e. a matching of the observations into disjoint pairs to minimize the total distance within pairs. The cross-match statistic is the number of pairs containing one observation from the first distribution and one from the second. Distributions that are very different will exhibit few cross-matches. When comparing two discrete distributions with finite support, the test is consistent against all alternatives. The test is applied to a study of brain activation measured by functional magnetic resonance imaging during two linguistic tasks, comparing brains that are impaired by arteriovenous abnormalities with normal controls. A second exact distribution-free test is also discussed: it ranks the pairs and sums the ranks of the cross-matched pairs.     

Some tests for uniformity of circular distributions powerful against multimodal alternatives

The author introduces new statistics suited for testing uniformity of circular distributions and powerful against multimodal alternatives. One of them has a simple expression in terms of the geometric mean of the sample of chord lengths. The others belong to a family indexed by a continuous parameter. The asymptotic distributions under the null hypothesis are derived. We compare the power of the new tests against Stephens's alternatives with those of Ajne, Watson, and Hermans‐Rasson's tests. Some of the new tests are the most powerful when the alternative has three or four modes. A heuristic justification of this feature is given. An application to the analysis of archaeological data is provided. The Canadian Journal of Statistics 38:80–96; 2010 © 2010 Statistical Society of Canada     

A wide review on exponentiality tests and two competitive proposals with application on reliability

In this paper, two new statistics based on comparison of the theoretical and empirical distribution functions are proposed to test exponentiality. Critical values are determined by means of Monte Carlo simulations for various sample sizes and different significance levels. Through an extensive simulation study, 50 selected exponentiality tests are studied for a wide collection of alternative distributions. From the empirical power study, it is concluded that, firstly, one of our proposals is preferable for IFR (increasing failure rate) and UFR (unimodal failure rate) alternatives, whereas the other one is preferable for DFR (decreasing failure rate) and BFR (bathtub failure rate) alternatives and, secondly, the new tests can be considered serious and powerful competitors to other existing proposals, since they have the same (or higher) level of performance than the best tests in the statistical literature.     

An outlier problem in the determination of ore grade

Data from recordings of ore assays from the Western Australian goldfields provide motivation to devise new tests for outliers when observations are distributed with the same mean but diff ering variances. In the case of equal variances, tests for a single outlier reduce to well-known tests of discordancy. A block discordancy test for k outliers is also described. The question of whether or not one should omit any observation(s) in the calculation of the mean recoverable gold content is addressed in the context of whether or not the data contain outliers, as judged by a normal model for the 'logged' ore assay values. The given data suggest that models with 'logged' values that follow long-tailed approximately normal distributions may be appropriate.     

A review of tests for exponentiality with Monte Carlo comparisons

In this paper, 91 different tests for exponentiality are reviewed. Some of the tests are universally consistent while others are against some special classes of life distributions. Power performances of 40 of these different tests for exponentiality of datasets are compared through extensive Monte Carlo simulations. The comparisons are conducted for different sample sizes of 10, 25, 50 and 100 for different groups of distributions according to the shape of their hazard functions at 5 percent level of significance. Also, the techniques are applied to two real-world datasets and a measure of power is employed for the comparison of the tests. The results show that some tests which are very good under one group of alternative distributions are not so under another group. Also, some tests maintained relatively high power over all the groups of alternative distributions studied while some others maintained poor power performances over all the groups of alternative distributions. Again, the result obtained from real-world datasets agree completely with those of the simulation studies.KEYWORDS: Classes of life distributions, empirical power of a test, exponentiality, goodness-of-fit test, Monte Carlo simulationSubject Classifications: 62E10, 62E20, 62F03     

Performance of bartlett adjustment for certain likelihood ratio tests

The effectiveness of Bartlett adjustment, using one of several methods of deriving a Bartlett factor, in improving the chi-squared approximation to the distribution of the log likelihood ratio statistic is investigated by computer simulation in three situations of practical interest:tests of equality of exponential distributions, equality of normal distributions and equality of coefficients of variation of normal distributions.     

Stringency-based ranking of normality tests

Many parametric statistical inferential procedures in finite samples depend crucially on the underlying normal distribution assumption. Dozens of normality tests are available in the literature to test the hypothesis of normality. Availability of such a large number of normality tests has generated a large number of simulation studies to find a best test but no one arrived at a definite answer as all depends critically on the alternative distributions which cannot be specified. A new framework, based on stringency concept, is devised to evaluate the performance of the existing normality tests. Mixture of t-distributions is used to generate the alternative space. The LR-tests, based on Neyman–Pearson Lemma, have been computed to construct a power envelope for calculating the stringencies of the selected normality tests. While evaluating the stringencies, Anderson–Darling (AD) statistic turns out to be the best normality test.     

Linear signed rank tests for hultivariate location

Previously proposed linear signed rank tests for multivariate location are not invariant under linear transformations of the observations, The asymptotic relative efficiencies of the tests 2 with respect to Hotelling's T²test depend on the direction of shift and the covariance matrix of the alternative distributions. For distributions with highly correlated components, the efficiencies of some of these tests can be arbitrarily low; they approach zero for certain multivariate normal alternatives, This article proposes a transformation of the data to be performed prior to standard linear signed rank tests, The resulting procedures have attractive power and efficiency properties compared to the original tests, In particular, for elliptically symmetric contiguous alternafives, the efficiencies of the new tests equal those of corresponding univariate linear signed rank tests with respect to the t test.     

Tests of fit for exponentiality based on a characterization via the mean residual life function

We study two new omnibus goodness of fit tests for exponentiality, each based on a characterization of the exponential distribution via the mean residual life function. The limiting null distributions of the tests statistics are the same as the limiting null distributions of the Kolmogorov-Smirnov and Cramér-von Mises statistics proposed when testing the simple hypothesis that the distribution of the sample variables is uniform on the interval [0, 1]. Work supported by the Deutsche Forschungsgemeinschaft     

Tests for Equality of Survival Distributions Against Non-Location Alternatives

We propose new two andk-sample tests for evaluating the equality of survival distributions against alternatives that include crossing of survival functions, and proportional and monotone hazard ratios. The tests allow for right censored data. The asymptotic power against local alternatives is investigated. Simulation results demonstrate that the new tests are more powerful than known tests when survival functions cross. We apply the tests to a well known study of chemo- and radio-therapy conducted by the Gastrointestinal Tumor Study Group. TheP-values for both proposed tests are much smaller than for other known tests.