[R] Robust ANOVA with variance heterogeneity
kjetilbrinchmannhalvorsen at gmail.com
Tue Oct 6 19:51:18 CEST 2009
Hola! see below.
On Sat, Oct 3, 2009 at 2:45 PM, David Winsemius <dwinsemius at comcast.net> wrote:
> Do you have a citation for that statement? I cannot convince myself that it
> should be true.
OK. that took some time, since I have no nonparametrics book with me,
but it is a fairly
standard assumption the friedman.test shares with wilcox.test and
others. One online reference giving this is:
"Friedman's test makes the following assumptions about the data in X:
All data come from populations having the same continuous
distribution, apart from possibly different locations due to column
and row effects.
All observations are mutually independent. "
This is also easy to investigate by simulation in R:
> A[, 1] <- rnorm(100, 0, 1)
> A[, 2] <- rnorm(100, 0, 5)
> A[, 3] <- rnorm(100, 0, 500)
Friedman rank sum test
Friedman chi-squared = 2.96, df = 2, p-value = 0.2276
which surprised me! This test seems to be somewhat robust against
variance heterogeneity ???, but that case is not included in the
> After looking at the CRAN Task View, I would suggest the OP look at
> rlm(MASS) or lmrob(robustbase).
> On Oct 2, 2009, at 11:05 AM, Kjetil Halvorsen wrote:
>> On Fri, Oct 2, 2009 at 8:45 AM, David Winsemius <dwinsemius at comcast.net>
>>> There are multiple routes to "robust" statistics, but the quick answer to
>>> this question is probably friedman.test
>> I don't think friedman.test is robust to variance heterogeneity. It is
>> only robust to
>>> I seem to remember a CRAN Task View on the area of Robust Statistics.
>>> David Winsemius
>>> On Oct 2, 2009, at 3:05 AM, Maike Luhmann wrote:
>>>> Dear list members,
>>>> I am looking for an alternative function for a two-way ANOVA in the case
>>>> variance heterogeneity. For one-way ANOVA, I found oneway.test(), but I
>>>> didn't find anything alike for two-way ANOVA. Does anyone have a
>>>> Thank you!
>>>> Maike Luhmann
>>>> Freie Universität Berlin
>>> R-help at r-project.org mailing list
>>> PLEASE do read the posting guide
>>> and provide commented, minimal, self-contained, reproducible code.
> David Winsemius, MD
> Heritage Laboratories
> West Hartford, CT
More information about the R-help