[R] ANOVA: Does a Between-Subjects Factor belong in the Error Term?
Peter Dalgaard
p.dalgaard at biostat.ku.dk
Mon Jul 9 22:23:09 CEST 2007
Alex Baugh wrote:
> I am executing a Repeated Measures Analysis of Variance with 1 DV (LOCOMOTOR
> RESPONSE), 2 Within-Subjects Factors (AGE, ACOUSTIC CONDITION), and 1
> Between-Subjects Factor (SEX).
>
> Does anyone know whether the between-subjects factor (SEX) belongs in the
> Error Term of the aov or not? And if it does belong, where in the Error Term
> does it go? The 3 possible scenarios are listed below:
>
>
>
> e.g.,
>
> 1. Omit Sex from the Error Term:
>
>
>> My.aov = aov(Locomotor.Response~(Age*AcousticCond*Sex) + Error
>>
> (Subject/(Timepoint*Acx.Cond)), data=locomotor.tab)
>
> note: Placing SEX outside the double paretheses of the Error Term has the
> same statistical outcome effect as omitting it all together from the Error
> Term (as shown above in #1).
>
>
>
> 2. Include SEX inside the Error Term (inside Double parentheses):
>
>
>> My.aov = aov(Locomotor.Response~(Age*AcousticCond*Sex) + Error
>>
> (Subject/(Timepoint*Acx.Cond+Sex)), data=locomotor.tab)
>
>
>
> 3. Include SEX inside the Error Term (inside Single parentheses):
>
>
>
>> My.aov = aov(Locomotor.Response~(Age*AcousticCond*Sex) + Error
>>
> (Subject/(Timepoint*Acx.Cond)+Sex), data=locomotor.tab)
>
> note: Placing SEX inside the single parentheses (as shown above in #3)
> generates no main effect of Sex. Thus, I'm fairly confident that option #3
> is incorrect.
>
>
>
> Scenarios 1,2, and 3 yield different results in the aov summary.
>
>
You don't generally want terms with systematic effects to appear as
error terms also, so 3 is wrong.
In 2 you basically have a random effect of sex within subject, which is
nonsensical since the subjects presumably have only one sex each. This
presumably generates an error stratum with 0 DF, which may well be harmless.
That leaves 1 as the likely solution.
You'll probably do yourself a favour if you learn to expand error terms,
a/b == a + a:b, etc.; that's considerably more constructive than trying
to think in terms of whether things are inside or outside parentheses.
>
> Thanks for your help!
>
> Alex
>
>
>
>
>
>
More information about the R-help
mailing list