On Testing Equality of Means of Correlated Variates with Missing Data on One Response

A statistic is proposed for testing the hypothesis of equality of the means of a bivariate normal distribution with unknown common variance and correlation coefficient when observations are missing on one of the variates. The distribution of the statistic is approximated by a normal distribution under the null hypothesis. The empirical powers of the statistic are computed and compared with those of the conventional paired t and the other known statistics. The power comparisons support the use of the proposed test.     

Testing Equality of Means in the Presence of Correlation and Missing Data

A statistic, derived from the combination of two dependent tests, is proposed for testing the hypothesis of equality of the means of a bivariate normal distribution with unknown common variance and correlation coefficient when observations are missing on one or both variates. The null distribution of the statistic is approximated by a well-known distribution. The empirical powers of the statistic are computed and compared with some of the known statistics. The comparisons support the use of the proposed test.     

On Testing Equality of Means of Correlated Variates with Incomplete Data

A statistic is proposed for testing the hypothesis of equality of the means of a bivariate normal distribution with unknown common variance and correlation coefficient when observations are missing on one of the variates. Expressions for the second and fourth central moments of the statistic are obtained. These moments are used to approximate the distribution of the statistic by a Student's t distribution under the null hypothesis. The powers of the test are computed and compared with those of the conventional paired t and the other known statistics.     

Testing Equality of Correlated Means in the Presence of Unequal Variances and Missing Values

Several statistics are proposed for testing the hypothesis of equality of the means of bivariate normal distribution with unknown variances and correlation coefficient when observations are missing on both variatea. The null distributions of the statistics are approximated by well-known distributions. The empirical sizes and powers of the statistics are computed and compared with paired t test and some of the known statistics based on available data. The comparisons support the use of two of the statistics proposed in this paper.     

Prediction of protein structural class by discriminant analysis

Protein structural class--alpha, beta, mixed (alpha/beta or alpha + beta), irregular--can be predicted from the amino acid sequence by discriminant analysis. Discrimination is based on distributions, in the classes, of vectors of attributes characterizing the sequences. In this paper, two sets of attributes and two methods of estimating their distributions are compared using more than 100 proteins from the Protein Data Bank. The best results were obtained when canonical variates of the frequencies of occurrence of 20 amino acids and non-parametric estimates of their distributions were used. Three variates are sufficient to allocate proteins to one of four classes with 83% reliability (estimated by cross-validation) and four variates allowed allocation to one of five classes with 78% reliability.     

Discrimination of DNA ploidy patterns by order statistics

The use of order statistics to discriminate and classify DNA ploidy patterns is proposed, especially for the classification of additional observations: whether a given sample is more likely to have come from a normal or an abnormal tissue, and with what probability, based on its ploidy pattern. The method involves the order of observations within each of several samples (e.g., euploid and aneuploid DNA patterns) and the use of subsets of the obtained order statistics as independent variables in a linear discriminant analysis. It thus replaces univariate observations by (some of) their order statistics, which are then used as the variables in the discriminant analysis. The procedure does not require normality of distributions or the transformation of nonnormal distributions, as do many discriminant functions; order statistics are usually distribution-free and thus are particularly useful for nonparametric inference. Preliminary simulation studies verified the potential usefulness of the order statistics discriminant function method as applied to DNA ploidy analysis. Its advantages as compared to the usual methods for hypothesis testing, e.g., the use of the chi-square or Kolmogorov-Smirnov tests to as certain "goodness-of-fit," is discussed. The proposed method is easy to implement and easy to interpret; it is also applicable to the study of distributions of other types of measurements.     

Interim Treatment Selection in Drug Development

We consider hypothesis testing in a clinical trial with an interim treatment selection. Recently, unconditional and conditional procedures for selecting one treatment as the winner have been proposed when the mean responses are approximately normal. In this paper, we generalize both procedures to multi-winner cases. The distributions of the test statistics are obtained and step-down approaches are proposed. We prove that both unconditional and conditional procedures strongly control the family-wise error rate. We give a brief discussion on power comparisons.     

The Doubly Noncentral Distribution of Wilks' Statistic in MANOVA

When plot effects are modeled in a randomized block design in multivariate analysis of variance the error and hypothesis matrices have independent noncentral Wishart distributions. This gives rise to doubly noncentral distribution of Wilks' statistic. The doubly noncentral distribution of Wilks' statistics, when the noncentrality matrices for the noncentral Wishart distributions in the above setting are of rank one or two, is investigated.     

Multipoint linkage detection in the presence of heterogeneity

Linkage heterogeneity is common for complex diseases. It is well known that loss of statistical power for detecting linkage will result if one assumes complete homogeneity in the presence of linkage heterogeneity. To this end, Smith (1963, Annals of Human Genetics 27, 175-182) proposed an admixture model to account for linkage heterogeneity. It is well known that for this model, the conventional chi-squared approximation to the likelihood ratio test for no linkage does not apply even when the sample size is large. By dealing with nuclear families and one marker at a time for genetic diseases with simple modes of inheritance, score-based test statistics (Liang and Rathouz, 1999, Biometrics 55, 65-74) and likelihood-ratio-based test statistics (Lemdani and Pons, 1995, Biometrics 51, 1033-1041) have been proposed which have a simple large-sample distribution under the null hypothesis of linkage. In this paper, we extend their work to more practical situations that include information from multiple markers and multi-generational pedigrees while allowing for a class of general genetic models. Three different approaches are proposed to eliminate the nuisance parameters in these test statistics. We show that all three approaches lead to the same asymptotic distribution under the null hypothesis of no linkage. Simulation results show that the proposed test statistics have adequate power to detect linkage and that the performances of these two classes of test statistics are quite comparable. We have applied the proposed method to a family study of asthma (Barnes et al., 1996), in which the score-based test shows evidence of linkage with p-value <0.0001 in the region of interest on chromosome 12. Additionally, we have implemented this score-based test within the frequently used computer package GENEHUNTER.     

Omnibus Two-Sample Test Procedures with Randomly Censored Data

Versions of the KOLMOGOROV -SMIRNOV , KUIPER , and cRAMéR -VON MISES statistics are introduced for the problem of comparing the underlying survival distributions of two populations when observations are subject to random censorship, and a two-sided alternative is of interest. These statistics are shown to have readily amenable asymptotic distributions; furthermore, the statistics reduce to their classical counterparts in the absence of censorship. Data pertaining to a cancer chemotherapy experiment are examined with these statistics.     

Artificial cranial deformation in Kow Swamp 1 and 5: A response to

Recently, Curnoe (2007) tested the predictions of competing models of modern human origins using three crania from Australia: Kow Swamp 1 and 5 and Keilor. The Kow Swamp specimens have long been suspected of having been altered through artificial deformation of the skull. Though Curnoe (2007) provided assurances that no evidence of deformation is present in those specimens, the current study retests the hypothesis that these Australian specimens are artificially deformed. The Australian crania are compared to known deformed individuals from New Britian through canonical variates analysis, and the resulting Mahalanobis distances are examined for statistical significance with random expectation statistics. The results show that Kow Swamp 1 and 5 have strong shape similarities to known deformed individuals, and both crania are very different in shape from Keilor. Keilor is statistically significantly different in shape from both Kow Swamp specimens and all of the known deformed specimens. These findings cast doubt on Curnoe's (2007) conclusions of a shared Australian cranial morphology as well as the retention of an archaic suite of morphologies in the Australians.     

New Ranked Set Sampling for One‐Parameter Family of Distributions

Bhoj (1997c) proposed a new ranked set sampling (NRSS) procedure for a specific two‐parameter family of distributions when the sample size is even. This NRSS procedure can be applied to one‐parameter family of distributions when the sample size is even. However, this procedure cannot be used if the sample size is odd. Therefore, in this paper, we propose a modified version of the NRSS procedure which can be used for one‐parameter distributions when the sample size is odd. Simple estimator for the parameter based on proposed NRSS is derived. The relative precisions of this estimator are higher than those of other estimators which are based on other ranked set sampling procedures and the best linear unbiased estimator using all order statistics.     

To permute or not to permute

Permutation test is a popular technique for testing a hypothesis of no effect, when the distribution of the test statistic is unknown. To test the equality of two means, a permutation test might use a test statistic which is the difference of the two sample means in the univariate case. In the multivariate case, it might use a test statistic which is the maximum of the univariate test statistics. A permutation test then estimates the null distribution of the test statistic by permuting the observations between the two samples. We will show that, for such tests, if the two distributions are not identical (as for example when they have unequal variances, correlations or skewness), then a permutation test for equality of means based on difference of sample means can have an inflated Type I error rate even when the means are equal. Our results illustrate permutation testing should be confined to testing for non-identical distributions. CONTACT: calian@raunvis.hi.is.     

Bayesian Prediction Regions for the Extreme Order Statistics

In this paper, prediction bounds for the order statistics are dealt with. For this purpose, predictive distributions derived by Bayesian approach, are utilized. In particular, bounds for the smallest and the largest order statistics are set when a series of samples are drawn from exponential, Pareto and power function populations. These bounds assist in knowing the nature of these predicted statistics even without actually having the sample observations.     

Testing some hypotheses on human evolution and sexual dimorphism

Research on human evolution and sexual dimorphism motivates an interesting test problem. In studying hominid phylogeny it is of interest to test whether parallel evolution plays a role. With regard to sexual dimorphism it is of interest to known whether the directions of sexual dimorphism in the populations being compared are the same. We show that testing these two problems gives rise to the same type of hypothesis testing, viz. the problem of testing the hypothesis that the means of independent, normally distributed random vectors with unit covariance matrices are situated on a straight line through the origin. A test is proposed and applied to study the sexual dimorphism of 20 recent skull populations. In this example the hypothesis of equal directions of sexual dimorphism is rejected. The classical theory of constructing multiple discriminant functions (canonical variates) is adapted to the problem of comparing sexual dimorphisms.     

CONDUCTING PHYLOGENETIC COMPARATIVE STUDIES WHEN THE PHYLOGENY IS NOT KNOWN

A method is proposed to conduct phylogenetic analyses of comparative or interspecific data when the true phylogeny is not known. Standard models of speciation and/or extinction or other methods are used to generate a sample from the set of all possible phylogenies for the measured species. The comparative data are then analyzed on each of the possible trees to obtain a distribution of possible evolutionary statistics for these data. The mean of this distribution is proposed as a reasonable estimate of the true evolutionary statistic of interest. Ways of obtaining confidence intervals and of developing hypothesis tests for this mean statistic are also proposed. The method can be used with any comparative method or phylogenetic analysis technique when phylogenetic relationships among species are not known or when branch lengths for a phylogeny in units of expected character change (as required by most methods) are not available. Computer programs to conduct the analyses are available on request.     

Estimating the accuracy of psychological scales using longitudinal data

This article presents a method for estimating the accuracy of psychological screening scales using receiver operating characteristic curves and associated statistics. Screening scales are typically semicontinuous within a known range with distributions that are nearly symmetric when the target condition is present and highly skewed when the condition is absent. We model screening scale outcomes using truncated normal distributions that accommodate these different distributional shapes and use subject-specific random effects to adjust for multiple assessments within individuals. Using the proposed model, we estimate the accuracy of the Symptom Checklist as a measure of major depression from a repeatedly screened sample of patients.     

A New Statistical Goodness-of-Fit Test Based on Graphical Representation

A statistical goodness-of-fit test, based on representing the sample observations by linked vectors, is developed. The direction and the length of the linked vectors are defined as functions of the expected values of the order statistics and sample order statistics, respectively. The underlying method can be used to test distributional assumptions for any location-scale family. A test statistic Q_n is introduced and some of its properties are studied. It is shown that the proposed test can be generalized to test if two or more independent samples come from the same distribution. The test procedure provides a graphical method of identifying the true distribution when the null hypothesis is rejected.     

A polar transport model of planarian regeneration

Abstract. A mathematical model based on the hypothesis of polar transport is presented for pattern formation and regeneration in planarian worms. The planarian body pattern is determined by two perpendicularly oriented morphogen distributions, one of which is bilaterally symmetrical. The mathematically simulated regeneration responses of this model are in good agreement with experimental observations. Further experimental work is needed to test the polar transport hypothesis.     

Testing equality of survival functions based on both paired and unpaired censored data

We introduce two test procedures for comparing two survival distributions on the basis of randomly right-censored data consisting of both paired and unpaired observations. Our procedures are based on generalizations of a pooled rank test statistic previously proposed for uncensored data. One generalization adapts the Prentice-Wilcoxon score, while the other adapts the Akritas score. The use of these particular scoring systems in pooled rank tests with randomly right-censored paired data has been advocated by several researchers. Our test procedures utilize the permutation distributions of the test statistics based on a novel manner of permuting the scores. Permutation versions of tests for right-censored paired data and for two independent right-censored samples that use the proposed scoring systems are obtained as special cases of our test procedures. Simulation results show that our test procedures have high power for detecting scale and location shifts in exponential and log-logistic distributions for the survival times. We also demonstrate the advantages of our test procedures in terms of utilizing randomly occurring unpaired observations that are discarded in test procedures for paired data. The tests are applied to skin graft data previously reported elsewhere.