1 / 12

One-way nonparametric ANOVA with trigonometric scores

Introducing a new rank test of the χ2 type that allows working with data violating normality assumption, commonly found in biometrical applications and log-transformed Cauchy data.

timothywillis
Download Presentation

One-way nonparametric ANOVA with trigonometric scores

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. One-way nonparametric ANOVA with trigonometric scores by Kravchuk, O.Y. School of Land and Food Sciences, University of Queensland

  2. Inspired by the simplicity of the Kruskal-Wallis k-sample procedure, we introduce a new rank test of the χ2 type that allows one to work with data that violates the normality assumption, being unimodal and symmetric but more heavier tailed than the normal. This type of non-normality is common in biometrical applications and also describes the distribution of the log-transformed Cauchy data. The distribution of the test statistic corresponds to the distribution of the first component of the well-known Cramer-von Mises test statistic.

  3. The test is asymptotically most efficient for the hyperbolic secant distribution that is compared to the normal and logistic distributions in the diagram below. Fig1: Standardised normal, hyperbolic secant and logistic densities

  4. The test is a one way rank based ANOVA, where we assume that within the k treatments the populations are continuous, belong to the same location family and may differ in the location parameter only. There are N experimental units, where the jth treatment accumulates nj units.The test statistic is built on k “bridges” corresponding to k linear contrasts of type T1=(A1<A2,A3). The asymptotic distribution of the test statistic is the χ2 with k-1 degrees of freedom. Computationally, the exact distribution is easy to construct on the basis of k-1 orthogonal contrasts (for example, for k = 3, T1, 2U3 and T2,3).

  5. For small samples (max(nj)<6), the chi-square approximation is more conservative than the exact null distribution. The diagram and table below provide the exact distribution for n1=3, n2=n3=2

  6. We illustrate the method by an artificial example of three normal populations different in location only. The populations are, correspondingly, N(0,1), N(-1,1), N(2,1). Random samples of size 8 are taken from these populations. One-way ANOVA F = 18.07, p = 0.000 Kruskal-Wallis KW = 14.11, p = 0.001 Trigonometric ANOVA Q=14.32, p = 0.001

  7. Illustrating the procedure… When there is a certain linear trend among the treatments, the corresponding bridge tends to have the U-shape. We measure the strength of such a tendency by the first coefficient of the Fourier sine-decomposition of the bridge. S1=1.17 S2=2.58 S3=-3.75

  8. The larger the sample sizes, the smoother the bridge. The actual shape of the bridge depends on the difference in location as well as on the distributions of the underlying populations. If the difference is large, the shape is strictly triangular regardless of the underlying distribution and the median k-sample test works well. If the difference in location is small, for symmetric, unimodal distributions, the shape of the bridges is determined by the tails of the distributions.

  9. The difference in scale among several Cauchy distributions may be analysed by means of the current test. To illustrate such an application, we perform the following ANOVA on the log-transformed Cauchy populations: Cauchy(0,1), Cauchy(0,5) and Cauchy(0,2). The log-transformation of the absolute values of the data makes it more normal-like. However, the analysis of the residuals of one-way ANOVA shows a departure from normality.The test allows us to perform the formal analysis and detect the difference in scale. The Kruskal-Wallis test gives a similar conclusion.

  10. The trigonometric ANOVA on log-transformed Cauchy… Random samples of size 8 were taken from the parent populations. One-way ANOVA Kruskal-Wallis F = 5.78, p = 0.01 KW = 11.26, p = 0.004 Trigonometric ANOVA Q=11.33, p = 0.003

  11. The multiple comparisons and contrasts are to be further developed for this test. The two-way test with trigonometric scores is to be investigated.The test performances are to be compared to the k-sample Cramer-von Mises test. Olena Kravchuk, LAFS, UQo.kravchuk@uq.edu.au(07) 33652171

More Related