A remark on a paper by misra

MISRA (1978) sets confidence intervals for a double linear compound of multivariate normal regression coefficients by using ROY'S maximum root test criterion. The exact test statistic to be used is STUDENT'S t. The t statistic gives narrower confidence bounds than those given by ROY's maximum root statistic. A result given by MORRISON (1975, p. 18, equation 10) for profile analysis is also obtained by using the STUDENT'S t test.     

Nonparametric All‐Pairs Multiple Comparisons

Nonparametric all‐pairs multiple comparisons based on pairwise rankings can be performed in the one‐way design with the Steel‐Dwass procedure. To apply this test, Wilcoxon's rank sum statistic is calculated for all pairs of groups; the maximum of the rank sums is the test statistic. We provide exact calculations of the asymptotic critical values (and P‐values, respectively) even for unbalanced designs. We recommend this asymptotic method whenever large sample sizes are present. For small sample sizes we recommend the use of the new statistic according to Baumgartner , Weiss , and Schindler (1998, Biometrics 54 , 1129–1135) instead of Wilcoxon's rank sum for the multiple comparisons. We show that the resultant procedure can be less conservative and, according to simulation results, more powerful than the original Steel‐Dwass procedure. We illustrate the methods with a practical data set.     

Network algorithm for the Exact Test of Hardy‐Weinberg Proportion for Multiple Alleles

We propose a new technique for the exact test of Hardy‐Weinberg proportion that considerably extends the bounds of computational feasibility. Our algorithm is constructed analogously to the network algorithm for Freeman‐Halton exact test in two‐way contingency tables. In this algorithm, the smallest and the largest values for the statistic are important and some interesting new theorems are proved for computing these values. Numerical examples are given to illustrate the practicality of our algorithm.     

Multiple Comparisons of Polynomial Distributions

Asymptotically correct 90 and 95 percentage points are given for multiple comparisons with control and for all pair comparisons of several independent samples of equal size from polynomial distributions. Test statistics are the maxima of the X²-statistics for single comparisons. For only two categories the asymptotic distributions of these test statistics result from DUNNETT'S many-one tests and TUKEY'S range test (cf. MILLER, 1981). The percentage points for comparisons with control are computed from the limit distribution of the test statistic under the overall hypothesis H₀. To some extent the applicability of these bounds is investigated by simulation. The bounds can also be used to improve Holm's sequentially rejective Bonferroni test procedure (cf. HOLM, 1979). The percentage points for all pair comparisons are obtained by large simulations. Especially for 3×3-tables the limit distribution of the test statistic under H₀ is derived also for samples of unequal size. Also these bounds can improve the corresponding Bonferroni-Holm procedure. Finally from SKIDÁK's probability inequality for normal random vectors (cf. SKIDÁK, 1967) a similar inequality is derived for dependent X²-variables applicable to simultaneous X²-tests.     

An Exact Unconditional test for the Hardy-Weinberg Equilibrium

An exact test based on the unconditional distribution of a test statistic for the Hardy-Weinberg equilibrium is introduced. This test is usually more powerful and requires less extensive tabulations compared with the (ordinary) exact conditional test. Tables of critical values are given for N = 5(1)100 and α =·10, ·05, ·01. Sample sizes required to attain a given power with the unconditional and conditional exact tests are provided.     

Using Discrete Distributions to Compute Powers for Goodness-of-Fit Test Statistics in a One-Way Multinomial Setting

Both the Bionomial and Poisson distributions are employed in this study to compute approximate powers for goodness-of-fit test statistics. The procedure adopted involves simulating 1000 samples from each of these distributions. These samples are then employed to compute both the randomized nominal critical values and estimated powers. The type I error rates returned from the use of the randomized critical levels Cα fall within the acceptable regions. We illustrate the use of the procedure with the Pearson's X² test statistic and show that this can readily be extended to any of the other well known goodness-of-fit test statistics.     

Median Tests for Censored Survival Data; a Contingency Table Approach

Summary The median failure time is often utilized to summarize survival data because it has a more straightforward interpretation for investigators in practice than the popular hazard function. However, existing methods for comparing median failure times for censored survival data either require estimation of the probability density function or involve complicated formulas to calculate the variance of the estimates. In this article, we modify a K ‐sample median test for censored survival data ( Brookmeyer and Crowley, 1982 , Journal of the American Statistical Association 77, 433–440) through a simple contingency table approach where each cell counts the number of observations in each sample that are greater than the pooled median or vice versa. Under censoring, this approach would generate noninteger entries for the cells in the contingency table. We propose to construct a weighted asymptotic test statistic that aggregates dependent χ² ‐statistics formed at the nearest integer points to the original noninteger entries. We show that this statistic follows approximately a χ² ‐distribution with k? 1 degrees of freedom. For a small sample case, we propose a test statistic based on combined p ‐values from Fisher’s exact tests, which follows a χ² ‐distribution with 2 degrees of freedom. Simulation studies are performed to show that the proposed method provides reasonable type I error probabilities and powers. The proposed method is illustrated with two real datasets from phase III breast cancer clinical trials.     

Signed-Rank Tests in the Presence of Ties (with extended tables)

We begin with a review of the areas of application of the signed-rank tests (SRTs) and we conclude that the results are exact only if no ties of non-null differences exist. In order to apply the SRTs according to WILCOXON and according to PRATT also in the presence of ties, by assigning midranks, we derive their null distributions. As special cases the null distributions for the problem without ties are obtained. In order to save the practising statistician the time-consuming calculations of the distribution functions, we compute tables of critical values (for reasons of volume they will be published as part of the reprints only). For N₀ = 0 (1) 5 null differences and M = = 1(1) 10 non-null differences the critical values of all distributions with all possible tie vectors are calculated. Instructions are provided and an example serves to illustrate the use of the table. The extension of the tables are obtained by means of counting formulas given in the text. Approximations are provided in order to make the application of tests possible for larger samples as well. It is shown that the approximation of the null distribution in the presence of ties by the null distributions under the assumption of no ties in some cases overstates and sometimes understates the exact rejection probability. For N₀ = 0 (1) 10 and M = 1 (1) 10 all distributions with all possible tie vectors for the SRTs with WILCOXON and PRATT ranking are examined with respect to the lattice type of the test statistic. The result is given in table 6. It is evident that the portion of PRATT -distributions with lattice character decreases as the number of null differences increases. Continuity corrections are obtained for the asymptotic normal distribution which take into account the lattice character of the distribution of the test statistic.     

An Adaptive Location‐Scale Test

Increasing locations are often accompanied by an increase in variability. In this case apparent heteroscedasticity can indicate that there are treatment effects and it is appropriate to consider an alternative involving differences in location as well as in scale. As a location‐scale test the sum of a location and a scale test statistic can be used. However, the power can be raised through weighting the sum. In order to select values for this weighting an adaptive design with an interim analysis is proposed: The data of the first stage are used to calculate the weights and with the second stage's data a weighted location‐scale test is carried out. The p‐values of the two stages are combined through Fisher's combination test. With a Lepage‐type location‐scale test it is illustrated that the resultant adaptive test can be more powerful than the ‘optimum’ test with no interim analysis. The principle to calculate weights, which cannot be reasonably chosen a priori, with the data of the first stage may be useful for other tests which utilize weighted statistics, too. Furthermore, the proposed test is illustrated with an example from experimental ecology.     

Testing the Equality of Dependent Variances

The likelihood ratio test for testing equality of vgE;2 correlated variables is developed. In general, evaluation of the test statistic involves the iterative optimization of a likelihood function with 1 + v(v – 1)/2 parameters. The explicit form of the test statistic is derived in the bivariate case, and an iterative algorithm for determining the maximum likelihood estimates is suggested. A limited Monte Carlo study determines the behavior of the proposed procedure under the null hypothesis and variety of parameter values.     

Approximating the Non-Null Distribution of the Likelihood Ratio Test Statistic

This paper is concerned with investigation into the behavior of the likelihood ratio test statistic G² when the alternative hypothesis M (QQ) is the true model. Exact moments of G² are computed empirically and three approximations are considered for approximating the non-null distribution of G². Our results show that the two parameter gamma distribution provides a closer approximation to the exact powers of G². A randomized procedure was employed to obtain critical values based on 1000 simulated samples.     

A Confidence‐Limit‐Based Approach to the Assessment of Hardy–Weinberg Equilibrium

The classical χ²‐procedure for the assessment of Hardy–Weinberg equilibrium (HWE) is tailored for detecting violations of HWE. However, many applications in genetic epidemiology require approximate compatibility with HWE. In a previous contribution to the field (Wellek, S. (2004). Biometrics, 60 , 694–703), the methodology of statistical equivalence testing was exploited for the construction of tests for problems in which the assumption of approximate compatibility of a given genotype distribution with HWE plays the role of the alternative hypothesis one aims to establish. In this article, we propose a procedure serving the same purpose but relying on confidence limits rather than critical bounds of a significance test. Interval estimation relates to essentially the same parametric function that was previously chosen as the target parameter for constructing an exact conditional UMPU test for equivalence with a HWE conforming genotype distribution. This population parameter is shown to have a direct genetic interpretation as a measure of relative excess heterozygosity. Confidence limits are constructed using both asymptotic and exact methods. The new approach is illustrated by reanalyzing genotype distributions obtained from published genetic association studies, and detailed guidance for choosing the equivalence margin is provided. The methods have been implemented in freely available SAS macros.     

On the exact null distribution of the generalised likelihood ratio test for analysing unreplicated factorial designs

Al-Shiha and Yang (1999) proposed a multistage procedure for analysing unreplicated factorial experiments, which is based on the statistic that is derived from the generalised likelihood ratio test statistic under the assumption of normality. It was shown by their simulation study that the method is quite competitive with Lenth's (1989) method. In their paper, because of the difficulty of determining the null distribution analytically, the quantiles of the null distribution were empirically simulated. In this paper, we give the exact null distribution of their test statistic, which makes it possible to calculate the critical values of the test.     

Some Simple Statistics to Test Against Ordered Alternatives in Complete Block Designs

Consider testing the hypothesis of no treatment effects against a postulated ranking of the m treatments, given data from n Complete Blocks. A suitable test statistic is the weighted average rank correlation w = σb_QiC_i where C_i is the correlation between the postulated ranking and the ranking observed within the ith block, Q_i is the rank of the ith block with respect to credibility, and the b_i's are weights such that 0 ≦ b¹ ≦ … ≦ bn. In this paper we introduce some simple statistics: the first extends the signed-rank statistic to m ≦ 3, the second uses a simple measure of correlation based on the antirank, and the third a statistic based on Spearman's footrule. Tables for critical values are provided and the normal approximation is investigated.     

gener: a server‐based analysis of pollen pool structure

The server‐based program gener performs the two‐generation analysis of pollen flow for data consisting of mother/offspring arrays using genetic markers. The gener program decomposes the genetic variance sampled by maternal individuals within and among pollen pool components of genetic variance and is accessible from http://dyerlab.bio.vcu.edu . These estimates are used to construct the test statistic, Φ_ft , whose significance is tested via permutation. The Φ_ft statistic can subsequently be used to quantify genetic effective pollen donor population size (N_ep), effective mating area and dispersal distance. Furthermore, the gener program can calculate Φ_ft values for all pairs of substrata within the data set.     

On Exact Unconditional Test for Linear Trend in Dose‐Response Studies

The one‐degree‐of‐freedom Cochran‐Armitage (C‐A) test statistic for linear trend has been widely applied in various dose‐response studies (e.g., anti‐ulcer medications and short‐term antibiotics, animal carcinogenicity bioassays and occupational toxicant studies). This approximate statistic relies, however, on asymptotic theory that is reliable only when the sample sizes are reasonably large and well balanced across dose levels. For small, sparse, or skewed data, the asymptotic theory is suspect and exact conditional method (based on the C‐A statistic) seems to provide a dependable alternative. Unfortunately, the exact conditional method is only practical for the linear logistic model from which the sufficient statistics for the regression coefficients can be obtained explicitly. In this article, a simple and efficient recursive polynomial multiplication algorithm for exact unconditional test (based on the C‐A statistic) for detecting a linear trend in proportions is derived. The method is applicable for all choices of the model with monotone trend including logistic, probit, arcsine, extreme value and one hit. We also show that this algorithm can be easily extended to exact unconditional power calculation for studies with up to a moderately large sample size. A real example is given to illustrate the applicability of the proposed method.     

A Robust Genome‐Wide Scan Statistic of the Wellcome Trust Case–Control Consortium

Summary In genome‐wide association (GWA) studies, test statistics that are efficient and robust across various genetic models are preferable, particularly for studying multiple diseases in the Wellcome Trust Case–Control Consortium ( WTCCC, 2007 , Nature 447 , 661–678). A new test statistic, the minimum of the p‐values of the trend test and Pearson's test, was considered by the WTCCC. It is referred to here as MIN2. Because the minimum of two p‐values is no longer a valid p‐value itself, the WTCCC only used it to rank single nucleotide polymorphisms (SNPs) but did not report the p‐values of the associated SNPs when MIN2 was used for ranking. Given its importance in practice, we derive the asymptotic null distribution of MIN2, study some of its analytical properties related to GWA studies, and compare it with existing methods (the trend test, Pearson's test, MAX3, and the constrained likelihood ratio test [CLRT]) by simulations across a wide range of possible genetic models: the recessive (REC), additive (ADD), multiplicative (MUL), dominant (DOM), and overdominant models. The results show that MAX3 and CLRT have greater efficiency robustness than other tests when the REC, ADD/MUL, and DOM models are possible, whereas Pearson's test and MIN2 have greater efficiency robustness if the possible genetic models also include the overdominant model. We conclude that robust tests (MAX3, MIN2, CLRT, and Pearson's test) are preferable to a single trend test for initial GWA studies. The four robust tests are applied to more than 100 SNPs associated with 11 common diseases identified by the two WTCCC GWA studies.     

Tests for Homogeneity of Variances Using Robust Weighted Likelihood Estimates

The classical normal-theory tests for testing the null hypothesis of common variance and the classical estimates of scale have long been known to be quite nonrobust to even mild deviations from normality assumptions for moderate sample sizes. Levene (1960) suggested a one-way ANOVA type statistic as a robust test. Brown and Forsythe (1974) considered a modified version of Levene's test by replacing the sample means with sample medians as estimates of population locations, and their test is computationally the simplest among the three tests recommended by Conover , Johnson , and Johnson (1981) in terms of robustness and power. In this paper a new robust and powerful test for homogeneity of variances is proposed based on a modification of Levene's test using the weighted likelihood estimates (Markatou , Basu , and Lindsay , 1996) of the population means. For two and three populations the proposed test using the Hellinger distance based weighted likelihood estimates is observed to achieve better empirical level and power than Brown-Forsythe's test in symmetric distributions having a thicker tail than the normal, and higher empirical power in skew distributions under the use of F distribution critical values.     

TESTING FOR PHYLOGENETIC SIGNAL IN BIOLOGICAL TRAITS: THE UBIQUITY OF CROSS‐PRODUCT STATISTICS

To evaluate rates of evolution, to establish tests of correlation between two traits, or to investigate to what degree the phylogeny of a species assemblage is predictive of a trait value so‐called tests for phylogenetic signal are used. Being based on different approaches, these tests are generally thought to possess quite different statistical performances. In this article, we show that the Blomberg et al. K and K*, the Abouheif index, the Moran's I, and the Mantel correlation are all based on a cross‐product statistic, and are thus all related to each other when they are associated to a permutation test of phylogenetic signal. What changes is only the way phylogenetic and trait similarities (or dissimilarities) among the tips of a phylogeny are computed. The definitions of the phylogenetic and trait‐based (dis)similarities among tips thus determines the performance of the tests. We shortly discuss the biological and statistical consequences (in terms of power and type I error of the tests) of the observed relatedness among the statistics that allow tests for phylogenetic signal. Blomberg et al. K* statistic appears as one on the most efficient approaches to test for phylogenetic signal. When branch lengths are not available or not accurate, Abouheif's C_mean statistic is a powerful alternative to K*.     

Testing for Homogeneity in Meta‐Analysis I. The One‐Parameter Case: Standardized Mean Difference

Summary Meta‐analysis seeks to combine the results of several experiments in order to improve the accuracy of decisions. It is common to use a test for homogeneity to determine if the results of the several experiments are sufficiently similar to warrant their combination into an overall result. Cochran’s Q statistic is frequently used for this homogeneity test. It is often assumed that Q follows a chi‐square distribution under the null hypothesis of homogeneity, but it has long been known that this asymptotic distribution for Q is not accurate for moderate sample sizes. Here, we present an expansion for the mean of Q under the null hypothesis that is valid when the effect and the weight for each study depend on a single parameter, but for which neither normality nor independence of the effect and weight estimators is needed. This expansion represents an order O(1/n) correction to the usual chi‐square moment in the one‐parameter case. We apply the result to the homogeneity test for meta‐analyses in which the effects are measured by the standardized mean difference (Cohen’s d‐statistic). In this situation, we recommend approximating the null distribution of Q by a chi‐square distribution with fractional degrees of freedom that are estimated from the data using our expansion for the mean of Q. The resulting homogeneity test is substantially more accurate than the currently used test. We provide a program available at the Paper Information link at the Biometrics website http://www.biometrics.tibs.org for making the necessary calculations.