On testing homogeneity of t normal means against ordered alternatives in r groups

Given k=rt normal populations with a common but unknown variance consisting of t treatments applied to r different groups of units, and supposing that in each group the means are monoto-nically non-decreasing (or non-increasing), then the likelihood ratio test of homogeneity of the means in each group against the simple order alternative is considered. Critical values are provided when one observation is drawn from each of the k populations.     

Testing homoosenetty of t normal means against ordered alternatives in 2 groups having uncommon and unknown variance

Given two series of normal populations with variances unknown and common within the series but possibly differ between the series, the likelihood ratio test of simultaneous homogeneity of means against sinpie order alternative is considered. Sole approxirations Par the null distribution of thetest statistic are also proposed.     

The power functions of the likelihood ratio tests for a simple tree ordering in normal means: unequal weights

Likelihood ratio tests are considered for two testing situations; testing for the homogeneity of k normal means against the alternative restricted by a simple tree ordering trend and testing the null hypothesis that the means satisfy the trend against all alternatives. Exact expressions are given for the power functions for k = 3 and 4 and unequal sample sizes, both for the case of known and unknown population variances, and approximations are discussed for larger k. Also, Bartholomew’s conjectures concerning minimal and maximal powers are investigated for the case of equal and unequal sample sizes. The power formulas are used to compute powers for a numerical example.     

Approximations to the powers of the likelihood ratio tests: the loop ordering and slippage alternatives

Approximations to the power functions of the likelihood ratio tests of homogeneity of normal means against the simple loop ordering at slippage alternatives are considered. If a researcher knows which mean is smallest and which is largest, but does not know how the other means are ordered, then a simple loop ordering is appropriate. The accuracy of the several moment approximations are studied for the case of known variances and it is found that for powers in the range typically of interest, the two-moment approximation seems quite adequate. Approximations based on mixtures of noncentral F variables are developed for the case of unknown variances. The critical values of the test statistics are also tabulated for selected levels of significance.     

Inference on measurements of agreement using marginal association

A statistical test concerning the comparison of two agreements for dependent observations is studied. The concept of stochastic ordering plays an important role in defining order of agreements. A one-sided likelihood ratio test for equality of two agreements is proposed. This test is closely related to the test of marginal homogeneity against marginal stochastic ordering. A real example is analyzed for illustration purposes.     

One-sided range test for testing against an ordered alternative under heteroscedasticity

In a one-way fixed effects analysis of variance model, when normal variances are unknown and possibly unequal, a one-sided range test for testing the null hypothesis H ₀ : μ ₁ = … = μ_k against an ordered alternative H_a : μ ₁ ≤ … ≤ μ_k by a single-stage and a two-stage procedure, respectively, is proposed. The critical values under H ₀ and the power under a specific alternative are calculated. Relation between the one-stage and the two-stage test procedures is discussed. A numerical example to illustrate these procedures is given.     

Optimal Allocation of Observations for Inference on k Ordered Normal Population Means

The optimal allocation of observations when there is a natural ordering in the k normal population means is discussed. It is shown that the design which minimizes the total mean square error of the maximum likelihood estimators in the null case allocates half the observations to each of the two extreme populations. The design is obviously optimal for testing the homogeneity of means against the simple ordered alternative. It is, however, hardly acceptable for the estimation in the nonnull case. It is, therefore, shown that the observations could be allocated to the non-extreme populations according to weights which are proportional to the absolute values of the Abelson and Tukey scores at the same time keeping the minimum local power for testing the simple ordered alternative to be maximal. The design gives also the maximum minimum power, not local, for the alternative in the class of linear tests. It, of course, suffers from a small loss of efficiency for the estimation under the null case but is much better under the nonnull case than the extreme design which allocates half the observations to each of the two extreme populations. Some numerical comparisons of the mean square errors are given.     

The one-way analysis of variance with ordered alternatives: A modification of bartholomew's Ē2 Test

Tests of homogeneity of normal means with the alternative restricted by an ordering on the means are considered. The simply ordered case, μ₁ ≤ μ₂ ≤ ··· ≤ μ_k, and the simple tree ordering, μ₁ ≤ μ_j, for; j= 2, 3,…, k, are emphasized. A modification of the likelihood-ratio test is proposed which is asymptotically equivalent to it but is more robust to violations of the hypothesized orderings. The new test has power at the points satisfying the hypothesized ordering which is similar to that of the likelihood-ratio test provided the degrees of freedom are not too small. The modified test is shown to be unbiased and consistent.     

Preliminary Goodness-of-Fit Tests for Normality do not Validate the One-Sample Student t

One of the most basic topics in many introductory statistical methods texts is inference for a population mean, μ. The primary tool for confidence intervals and tests is the Student t sampling distribution. Although the derivation requires independent identically distributed normal random variables with constant variance, σ², most authors reassure the readers about some robustness to the normality and constant variance assumptions. Some point out that if one is concerned about assumptions, one may statistically test these prior to reliance on the Student t. Most software packages provide optional test results for both (a) the Gaussian assumption and (b) homogeneity of variance. Many textbooks advise only informal graphical assessments, such as certain scatterplots for independence, others for constant variance, and normal quantile–quantile plots for the adequacy of the Gaussian model. We concur with this recommendation. As convincing evidence against formal tests of (a), such as the Shapiro–Wilk, we offer a simulation study of the tails of the resulting conditional sampling distributions of the Studentized mean. We analyze the results of systematically screening all samples from normal, uniform, exponential, and Cauchy populations. This pretest does not correct the erroneous significance levels and makes matters worse for the exponential. In practice, we conclude that graphical diagnostics are better than a formal pretest. Furthermore, rank or permutation methods are recommended for exact validity in the symmetric case.     

A robust modification of the ordered-heterogeneity test

An ordered heterogeneity (OH) test is a test for a trend that combines a non-directional heterogeneity test with the rank-order information specified under the alternative. We propose two modifications of the OH test procedure: (1) to use the mean ranks of the groups rather than the sample means to determine the observed ordering of the groups, and (2) to use the maximum correlation out of the 2^k???1 – 1 possibilities under the alternative rather than the single ordering (1, 2, … , k), where k is the number of independent groups. A simulation study indicates that these two changes increase the power of the ordered heterogeneity test when, as common in practice, the underlying distribution may deviate from a normal distribution and the trend pattern is a priori unknown. In contrast to the original OH test, the modified OH test can detect all possible patterns under the alternative with a relatively high power.     

Improved estimation in Lognormal regression models

Lognormal regression model with unknown error variance is considered. We give a class of estimators of the regression coefficients vector improving upon traditional estimator when the number of independent variables is at least three. The relationship between these estimators on one hand and James-Stein type estimators of the normal mean and improved estimators of the normal variance on another hand is discussed.     

Test for homogeneity in Hardy–Weinberg normal mixture model

The phenotype of a quantitative trait locus (QTL) is often modeled by a finite mixture of normal distributions. If the QTL effect depends on the number of copies of a specific allele one carries, then the mixture model has three components. In this case, the mixing proportions have a binomial structure according to the Hardy–Weinberg equilibrium. In the search for QTL, a significance test of homogeneity against the Hardy–Weinberg normal mixture model alternative is an important first step. The LOD score method, a likelihood ratio test used in genetics, is a favored choice. However, there is not yet a general theory for the limiting distribution of the likelihood ratio statistic in the presence of unknown variance. This paper derives the limiting distribution of the likelihood ratio statistic, which can be described by the supremum of a quadratic form of a Gaussian process. Further, the result implies that the distribution of the modified likelihood ratio statistic is well approximated by a chi-squared distribution. Simulation results show that the approximation has satisfactory precision for the cases considered. We also give a real-data example.     

Asymptotically valid single-stage multiple-comparison procedures

We establish general conditions for the asymptotic validity of single-stage multiple-comparison procedures (MCPs) under the following general framework. There is a finite number of independent alternatives to compare, where each alternative can represent, e.g., a population, treatment, system or stochastic process. Associated with each alternative is an unknown parameter to be estimated, and the goal is to compare the alternatives in terms of the parameters. We establish the MCPs’ asymptotic validity, which occurs as the sample size of each alternative grows large, under two assumptions. First, for each alternative, the estimator of its parameter satisfies a central limit theorem (CLT). Second, we have a consistent estimator of the variance parameter appearing in the CLT. Our framework encompasses comparing means (or other moments) of independent (not necessarily normal) populations, functions of means, quantiles, steady-state means of stochastic processes, and optimal solutions of stochastic approximation by the Kiefer–Wolfowitz algorithm. The MCPs we consider are multiple comparisons with the best, all pairwise comparisons, all contrasts, and all linear combinations, and they allow for unknown and unequal variance parameters and unequal sample sizes across alternatives.     

Conditional interval estimation of the mean following rejection of a two sided test

The effect of rejecting a two-sided preliminary test of significance for the mean of a normal distribution upon subsequent interval estimation of the mean is examined. For the case where the variance is known, conditional confidence intervals may be shorter than unconditional intervals, in contrast to the one-sided preliminary test case examined by Meeks and D’Agostino (1983, The American Statistician, 7, 134-136) . For the case where the variance is unknown and must be estimated by the sample variance, it is shown that customary intervals do not offer uniformly greater or lesser coverage than the nominal level.     

Comparison of ANOVA-F and ANOM tests with regard to type I error rate and test power

A Monte Carlo simulation was conducted to compare the type I error rate and test power of the analysis of means (ANOM) test to the one-way analysis of variance F-test (ANOVA-F). Simulation results showed that as long as the homogeneity of the variance assumption was satisfied, regardless of the shape of the distribution, number of group and the combination of observations, both ANOVA-F and ANOM test have displayed similar type I error rates. However, both tests have been negatively affected from the heterogeneity of the variances. This case became more obvious when the variance ratios increased. The test power values of both tests changed with respect to the effect size (Δ), variance ratio and sample size combinations. As long as the variances are homogeneous, ANOVA-F and ANOM test have similar powers except unbalanced cases. Under unbalanced conditions, the ANOVA-F was observed to be powerful than the ANOM-test. On the other hand, an increase in total number of observations caused the power values of ANOVA-F and ANOM test approach to each other. The relations between effect size (Δ) and the variance ratios affected the test power, especially when the sample sizes are not equal. As ANOVA-F has become to be superior in some of the experimental conditions being considered, ANOM is superior in the others. However, generally, when the populations with large mean have larger variances as well, ANOM test has been seen to be superior. On the other hand, when the populations with large mean have small variances, generally, ANOVA-F has observed to be superior. The situation became clearer when the number of the groups is 4 or 5.     

On the level probabilities for useful partially ordered alternatives in the analysis of variance

In the analysis of variance, we often encounter situations in which we want to test the null hypothesis of homogeneity of the normal means against various partially ordered alternative hypotheses. We study likelihood ratio tests for three useful types of alternatives: d-star, bipartite and broom tree. Especially, we give computational formulas for the level probabilities of the alternative types. The results permit us to obtain critical values for practical use.     

The null distribution of the likelihood ratio test for a mixture of two normals after a restricted box-cox transformation

Maclean et al. (1976) applied a specific Box-Cox transformation to test for mixtures of distributions against a single distribution. Their null hypothesis is that a sample of n observations is from a normal distribution with unknown mean and variance after a restricted Box-Cox transformation. The alternative is that the sample is from a mixture of two normal distributions, each with unknown mean and unknown, but equal, variance after another restricted Box-Cox transformation. We developed a computer program that calculated the maximum likelihood estimates (MLEs) and likelihood ratio test (LRT) statistic for the above. Our algorithm for the calculation of the MLEs of the unknown parameters used multiple starting points to protect against convergence to a local rather than global maximum. We then simulated the distribution of the LRT for samples drawn from a normal distribution and five Box-Cox transformations of a normal distribution. The null distribution appeared to be the same for the Box-Cox transformations studied and appeared to be distributed as a chi-square random variable for samples of 25 or more. The degrees of freedom parameter appeared to be a monotonically decreasing function of the sample size. The null distribution of this LRT appeared to converge to a chi-square distribution with 2.5 degrees of freedom. We estimated the critical values for the 0.10, 0.05, and 0.01 levels of significance.     

On Testing the Mean Equivalence of Treatments from Correlated Normal Populations

Let there be k equally correlated treatment populations under consideration, a Studentized range test is proposed to test the hypothesis of average mean equivalence against the alternative hypothesis of inequivalence. The maximum level and minimum power at some least favorable configurations of means are used to calculate the critical value and the required sample size simultaneously when testing a null against an alternative hypothesis. The range test is applied to a real world problem to find out if the stress levels among children at four time periods due to a newly built nearby airport are equivalent.     

On testing the bioequivalence of several treatments using the measure of distance

A studentized range test is proposed to test the hypothesis of bioequivalence of normal means in terms of a standardized distance among means. A least favourable configuration (LFC) of means to guarantee the maximum level at a null hypothesis and an LFC of means to guarantee the minimum power at an alternative hypothesis are obtained. This level and power of the test are fully independent of the unknown means and variances. For a given level, the critical value of the test under a null hypothesis can be determined. Furthermore, if the power under an alternative is also required at a given level, then both the critical value and the required sample size for an experiment can be simultaneously determined. In situations where the common population variance is unknown and the bioequivalence is the actual distance between means without standardization, a two-stage sampling procedure can be employed to find these solutions.     

Change Detection in INAR(p) Processes Against Various Alternative Hypotheses

Change in the coefficients or the mean of the innovation of an INAR(p) process is a sign of disturbance that is important to detect. The proposed methods can test for change in any one of these quantities separately, or in any collection of them. They make both one-sided and two-sided tests possible, furthermore, they can be used to test against the “epidemic” alternative. The tests are based on a CUSUM process using CLS estimators of the parameters. Under the one-sided and two-sided alternatives, consistency of the tests is proved and the properties of the change-point estimator are also explored.