The Kruskal–Wallis test by ranks, Kruskal–Wallis H test[1] (named after William Kruskal and W. Allen Wallis), or one-way ANOVA on ranks[1] is a non-parametric method for testing whether samples originate from the same distribution.[2][3][4] It is used for comparing two or more independent samples of equal or different sample sizes. It extends the Mann–Whitney U test, which is used for comparing only two groups. The parametric equivalent of the Kruskal–Wallis test is the one-way analysis of variance (ANOVA).
A significant Kruskal–Wallis test indicates that at least one sample stochastically dominates one other sample. The test does not identify where this stochastic dominance occurs or for how many pairs of groups stochastic dominance obtains. For analyzing the specific sample pairs for stochastic dominance, Dunn's test,[5] pairwise Mann-Whitney tests without Bonferroni correction,[6] or the more powerful but less well known Conover–Iman test[6] are sometimes used.
Since it is a non-parametric method, the Kruskal–Wallis test does not assume a normal distribution of the residuals, unlike the analogous one-way analysis of variance. If the researcher can make the assumptions of an identically shaped and scaled distribution for all groups, except for any difference in medians, then the null hypothesis is that the medians of all groups are equal, and the alternative hypothesis is that at least one population median of one group is different from the population median of at least one other group.
Contents
1Method
2Exact probability tables
3Exact distribution of
4See also
5References
6Further reading
7External links
Method
Rank all data from all groups together; i.e., rank the data from 1 to N ignoring group membership. Assign any tied values the average of the ranks they would have received had they not been tied.
The test statistic is given by:
where:
is the number of observations in group
is the rank (among all observations) of observation from group
is the total number of observations across all groups
is the average rank of all observations in group
is the average of all the .
If the data contain no ties the denominator of the expression for is exactly and . Thus
The last formula only contains the squares of the average ranks.
A correction for ties if using the short-cut formula described in the previous point can be made by dividing by , where G is the number of groupings of different tied ranks, and ti is the number of tied values within group i that are tied at a particular value. This correction usually makes little difference in the value of H unless there are a large number of ties.
Finally, the decision to reject or not the null hypothesis is made by comparing to a critical value obtained from a table or a software for a given significance or alpha level. If is bigger than , the null hypothesis is rejected. If possible (no ties, sample not too big) one should compare to the critical value obtained from the exact distribution of . Otherwise, the distribution of H can be approximated by a chi-squared distribution with g-1 degrees of freedom. If some values are small (i.e., less than 5) the exact probability distribution of can be quite different from this chi-squared distribution. If a table of the chi-squared probability distribution is available, the critical value of chi-squared, , can be found by entering the table at g − 1 degrees of freedom and looking under the desired significance or alpha level.
If the statistic is not significant, then there is no evidence of stochastic dominance between the samples. However, if the test is significant then at least one sample stochastically dominates another sample. Therefore, a researcher might use sample contrasts between individual sample pairs, or post hoc tests using Dunn's test, which (1) properly employs the same rankings as the Kruskal–Wallis test, and (2) properly employs the pooled variance implied by the null hypothesis of the Kruskal–Wallis test in order to determine which of the sample pairs are significantly different.[5] When performing multiple sample contrasts or tests, the Type I error rate tends to become inflated, raising concerns about multiple comparisons.
Exact probability tables
A large amount of computing resources is required to compute exact probabilities for the Kruskal–Wallis test. Existing software only provides exact probabilities for sample sizes less than about 30 participants. These software programs rely on asymptotic approximation for larger sample sizes.
Exact probability values for larger sample sizes are available. Spurrier (2003) published exact probability tables for samples as large as 45 participants.[7] Meyer and Seaman (2006) produced exact probability distributions for samples as large as 105 participants.[8]
Exact distribution of
Choi et al.[9]
make a review of two methods that had been developed to compute the exact distribution of , propose a new one, and compare the exact distribution to its chi-squared approximation.
See also
Friedman test
Jonckheere's trend test
References
^ abKruskal–Wallis H Test using SPSS Statistics, Laerd Statistics
^Kruskal; Wallis (1952). "Use of ranks in one-criterion variance analysis". Journal of the American Statistical Association. 47 (260): 583–621. doi:10.1080/01621459.1952.10483441.
^Corder, Gregory W.; Foreman, Dale I. (2009). Nonparametric Statistics for Non-Statisticians. Hoboken: John Wiley & Sons. pp. 99–105. ISBN 9780470454619.
^Siegel; Castellan (1988). Nonparametric Statistics for the Behavioral Sciences (Second ed.). New York: McGraw–Hill. ISBN 0070573573.
^ abDunn, Olive Jean (1964). "Multiple comparisons using rank sums". Technometrics. 6 (3): 241–252. doi:10.2307/1266041.
^ abConover, W. Jay; Iman, Ronald L. (1979). "On multiple-comparisons procedures" (PDF) (Report). Los Alamos Scientific Laboratory. Retrieved 2016-10-28.
^Spurrier, J. D. (2003). "On the null distribution of the Kruskal–Wallis statistic". Journal of Nonparametric Statistics. 15 (6): 685–691. doi:10.1080/10485250310001634719.
^Meyer; Seaman (April 2006). "Expanded tables of critical values for the Kruskal–Wallis H statistic". Paper presented at the annual meeting of the American Educational Research Association, San Francisco. Critical value tables and exact probabilities from Meyer and Seaman are available for download at http://faculty.virginia.edu/kruskal-wallis/. A paper describing their work may also be found there.
^Won Choi, Jae Won Lee, Myung-Hoe Huh, andSeung-Ho Kang (2003). "An Algorithm for Computing the Exact Distribution of the Kruskal–Wallis Test". Communications in Statistics - Simulation and Computation (32, number 4): 1029–1040. doi:10.1081/SAC-120023876.CS1 maint: Multiple names: authors list (link)
Further reading
Daniel, Wayne W. (1990). "Kruskal–Wallis one-way analysis of variance by ranks". Applied Nonparametric Statistics (2nd ed.). Boston: PWS-Kent. pp. 226–234. ISBN 0-534-91976-6.
… is endemic in some parts of Oceania, including Wallis and Fortuna Islands and Tonga. Infection is usually acquired during the butchering of wild pigs . Papua New Guinea has several venomous snakes …
… 2531 paired maternal-umbilical cord blood samples were tested for dengue-specific immunoglobulin (Ig)M to determine… including American Samoa, Cook Islands, French Polynesia, New Caledonia, and Tonga .…
…False-positive tests are also uncommon . Agreement between the two most commonly used commercial test formats, latex agglutination and enzyme-linked immunoassays (EIAs), is excellent (92 to 98 percent) . Test utility… thought to be largely restricted to tropical and subtropical regions, including Australia and Papua New Guinea. It is now clear that sporadic cases occur in various regions around the world. The clinical …
… Fish can be tested using a mouse bioassay and an IgG immunoassay, but these tests are costly and time consuming and are not widely used .… toxicity although it can also be seen in patients with neurotoxic shellfish poisoning. In the South Pacific, paresis may occur in up to 10 percent of patients . Reversible cerebellar dysfunction also can …
… clinical specimen Serologic testing alone is not a reliable method of diagnosis. Indirect hemagglutination test (IHAT) and other serologic tests are available,… Cases are also reported from Papua New Guinea and New Caledonia .…
English Journal
Interleukin-10 responses from acute exercise in healthy subjects: A systematic review.
Cabral-Santos C, de Lima Junior EA, Fernandes IMDC, Pinto RZ, Rosa-Neto JC, Bishop NC, Lira FS.
Journal of cellular physiology. 2019 Jul;234(7)9956-9965.
Interleukin 10 (IL-10) is a cytokine that plays a critical role with potent anti-inflammatory properties when produced during exercise, limiting host immune response to pathogens and preventing tissue damage. The purpose of this systematic review was to assess the response of IL-10 after acute exerc
Load-bearing capacities of ultra-thin occlusal veneers bonded to dentin.
Maeder M, Pasic P, Ender A, Özcan M, Benic GI, Ioannidis A.
Journal of the mechanical behavior of biomedical materials. 2019 Jul;95()165-171.
To test whether the load-bearing capacity of occlusal veneers made of ceramic or hybrid materials bonded to dentin does differ from those of porcelain-fused-to metal or lithium disilicate glass ceramic crowns. In 80 human molars, occlusal tooth substance was removed so that the defects extended into
Graphene Oxide/Silver Nanoparticle Coating Produced by Electrophoretic Deposition Improved the Mechanical and Tribological Properties of NiTi Alloy for Biomedical Applications.
Rokaya D, Srimaneepong V, Qin J, Siraleartmukul K, Siriwongrungson V.
Journal of nanoscience and nanotechnology. 2019 Jul;19(7)3804-3810.
The aim of this study was to evaluate the mechanical and tribological properties of graphene oxide/silver nanoparticle (GO/AgNP) coated medical grade nickel-titanium (NiTi) alloy. The alloy substrates were coated using electrophoretic deposition at 30 V for 1, 5, and 10 min and were characterized by
鹿児島大学教育学部研究紀要. 教育科学編 = Bulletin of the Faculty of Education, Kagoshima University. Studies in education (70), 45-53, 2019-03-11
… The analysis of the survey results compared grades by way of the Kruskal-Wallis one-way analysis of variance test.Consideration of physique and differences in physical fitness were inspected and researched by way of the Chi-squared test and Spearman's correlation between two variables.The results of the study revealed that parents' anxiety over the safety of first grader's' who learn judo was significantly higher than that of second and third grader's' parents. …
Introduction A Kruskal-Wallis test is typically performed when an analyst would like to test for differences between three or more treatments or conditions. However, unlike a one-way ANOVA, the response variable of interest is not ...
Summary Use the Kruskal–Wallis test when you have one nominal variable and one ranked variable. It tests whether the mean ranks are the same in all the groups. When to use it The most common use of the Kruskal–Wallis test is ...