[R] [OT] factorial design
Rolf Turner
r.turner at auckland.ac.nz
Tue Nov 4 23:20:05 CET 2008
On 4/11/2008, at 6:29 PM, Edna Bell wrote:
> Dear R Gurus:
>
> I vaguely remember reading that if interaction was present in a
> factorial design, then the main effect results were suspect.
>
> However, I was reading a text which now uses the tests for main
> effects even if interaction is present.
>
> Which is correct, please?
Both, neither, whatever.
You have to specify what you mean by ``main effect'' if interaction
is present. One possible definition is to average the response of
factor A (at each level) over all levels of factor B. There is an
``A main effect'' (by this definition) if these averages are different.
E.g. (1): mu(1,1) = 1, mu(1,2) = -1, mu(2,1) = -1, mu(2,2) = 1
where mu(i,j) is the population mean when factor A is at level i
and factor B is at level j. Factor A averages out to 0 at both levels
so there is no main A effect in the sense defined. (And likewise
there is no main B effect either.) ``Clearly'' there is interaction!
E.g. (2): mu(1,1) = 1, mu(1,2) = 2, mu(2,1) = -1, mu(2,2) = -2.
Here factor A averages to 1.5 at level 1 and to -1.5 at level 2
so there *is* a main A effect ***in the sense defined***. Note
that factor B averages to 0 at both levels so there is no main
B effect in the sense defined.
This is the sort of thing that the much maligned ``Type III'' sums
of squares test for.
But is it *really* what you are interested in? Probably NOT.
Read Bill Venables' ``Exegesis on Linear Models'' for more discussion
of why such things are a bad idea:
http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf
The important thing is to get it clear in your head what you want
to know about, and then conduct your analysis in such a way as to
find about that. This is sometimes not so easy. But there's no
avoiding it if you want to do a sensible analysis. Doing tests
mindlessly without a clear understanding of exactly what it is that
you are testing is a recipe for nonsense if not disaster.
cheers,
Rolf Turner
######################################################################
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
More information about the R-help
mailing list