[R-lang] Odd summary(aov()) result

Francisco Torreira ftorrei2 at uiuc.edu
Mon May 14 10:52:25 PDT 2007


Dear Adam,

> I'm just surprised
> that commutativity doesn't hold.

A*B is is not the equivalent notation of any combination of A, B and
A:B. It is only the equivalent of A+B+A:B. When your design is
completely balanced, the output should be the same regardless of the
order. However, with unbalanced/incomplete designs you will get
different outputs. This is because of the famous Type I / Type III
sums of squares issue. To be honest I don't think this is the answer
to your problem, since repeated-measures designs require only 1 data
point per cell. Is any data point missing?

Francisco

On 5/14/07, Adam Baker <adamb924 at gmail.com> wrote:
> I wonder if somebody could comment on why I get different summaries
> with these two commands, which differ only in the order of the
> variables.
>
> summary(aov(MidAvg ~ Ccon*Rpres + Error(Subject/Ccon*Rpres),data=nr))
> summary(aov(MidAvg ~ Rpres*Ccon + Error(Subject/Ccon*Rpres),data=nr))
>
> With the first command, R doesn't report the main effect of Rpres.
> With the second command, it does (code & output below). It's not very
> important since the interaction is significant; I'm just surprised
> that commutativity doesn't hold.
>
> Thanks,
> Adam
>
>
> > summary(aov(MidAvg ~ Ccon*Rpres + Error(Subject/Ccon*Rpres),data=nr))
>
> Error: Subject
>           Df   Sum Sq  Mean Sq F value   Pr(>F)
> Ccon       1 13245765 13245765  17.767 0.002256 **
> Residuals  9  6709794   745533
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Rpres
>      Df   Sum Sq  Mean Sq
> Ccon  1 38923433 38923433
>
> Error: Subject:Ccon
>            Df   Sum Sq  Mean Sq F value    Pr(>F)
> Ccon        2 16646448  8323224 239.589 3.334e-14 ***
> Ccon:Rpres  1   906946   906946  26.107 6.229e-05 ***
> Residuals  19   660053    34740
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Subject:Rpres
>            Df Sum Sq Mean Sq F value  Pr(>F)
> Ccon:Rpres  1 645835  645835  9.5416 0.01295 *
> Residuals   9 609176   67686
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Subject:Ccon:Rpres
>            Df  Sum Sq Mean Sq F value    Pr(>F)
> Ccon:Rpres  2 7636019 3818009  676.48 < 2.2e-16 ***
> Residuals  20  112879    5644
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Within
>            Df  Sum Sq Mean Sq F value Pr(>F)
> Residuals 499 9862093   19764
>
>
>
>
>
>
> > summary(aov(MidAvg ~ Rpres*Ccon + Error(Subject/Ccon*Rpres),data=nr))
>
> Error: Subject
>           Df   Sum Sq  Mean Sq F value   Pr(>F)
> Rpres      1 13245765 13245765  17.767 0.002256 **
> Residuals  9  6709794   745533
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Rpres
>       Df   Sum Sq  Mean Sq
> Rpres  1 38923433 38923433
>
> Error: Subject:Ccon
>            Df   Sum Sq  Mean Sq F value    Pr(>F)
> Ccon        2 16646448  8323224 239.589 3.334e-14 ***
> Rpres:Ccon  1   906946   906946  26.107 6.229e-05 ***
> Residuals  19   660053    34740
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Subject:Rpres
>            Df Sum Sq Mean Sq F value  Pr(>F)
> Rpres:Ccon  1 645835  645835  9.5416 0.01295 *
> Residuals   9 609176   67686
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Subject:Ccon:Rpres
>            Df  Sum Sq Mean Sq F value    Pr(>F)
> Rpres:Ccon  2 7636019 3818009  676.48 < 2.2e-16 ***
> Residuals  20  112879    5644
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Error: Within
>            Df  Sum Sq Mean Sq F value Pr(>F)
> Residuals 499 9862093   19764
> _______________________________________________
> R-lang mailing list
> R-lang at ling.ucsd.edu
> https://ling.ucsd.edu/mailman/listinfo.cgi/r-lang
>


-- 
Francisco Torreira
PhD Candidate in Romance Linguistics
University of Illinois at Urbana-Champaign

https://netfiles.uiuc.edu/ftorrei2/www/index.html
tel: (+1) 217 - 778 8510


More information about the R-lang mailing list