[R] aov or lme effect size calculation

Doran, Harold HDoran at air.org
Tue Sep 2 20:14:28 CEST 2008


Greg

Upgrade your packages to the supported versions (lme4 and Matrix), and
use lmer and not lme.

### Example
> example(lmer)
> anova(fm1)
Analysis of Variance Table
     Df Sum Sq Mean Sq F value
Days  1  30032   30032  45.854

Your method for eta-squared with a mixed model is another story,
however.


> -----Original Message-----
> From: Greg Trafton [mailto:greg.trafton at nrl.navy.mil] 
> Sent: Tuesday, September 02, 2008 1:57 PM
> To: Doran, Harold
> Cc: r-help at r-project.org
> Subject: Re: [R] aov or lme effect size calculation
> 
> Sorry about that.  My problem is computational, not 
> statistical and exactly as you say:  I don't quite know how 
> to get the correct variance component from either aov or lme. 
>  the way to compute partial eta squared is:
> 
> partial-eta-squared = SS(effect) / (SS(effect) + SS(error))
> 
> AOV gives Sum Squares for both effects and the interaction, 
> but lme doesn't even give that in default format.
> 
> thanks,
> greg
> 
> On Sep 2, 2008, at 11:43 AM, Doran, Harold wrote:
> 
> > Greg
> >
> > You haven't really explained what your problem is. If it is 
> conceptual 
> > (i.e., how do I do it) this is not really the right place 
> for in-depth 
> > statistical advice, but it is often given. OTOH, if your problem is 
> > computational, please explain what that is? For example, maybe you 
> > know how to compute eta-squared, but you want to extract 
> the variance 
> > component and you can't figure that out.
> >
> > Without more info, it is hard to help. Now, with that said, 
> with lme 
> > (or mixed models) you have multiple variance components, so 
> how would 
> > you go about computing eta-squared anyhow?
> >
> >> -----Original Message-----
> >> From: r-help-bounces at r-project.org
> >> [mailto:r-help-bounces at r-project.org] On Behalf Of Greg Trafton
> >> Sent: Tuesday, September 02, 2008 10:25 AM
> >> To: r-help at r-project.org
> >> Subject: [R] aov or lme effect size calculation
> >>
> >> (A repost of this request with a bit more detail)
> >>
> >> Hi, All.  I'd like to calculate effect sizes for aov or 
> lme and seem 
> >> to have a bit of a problem.  partial-eta squared would be my first 
> >> choice, but I'm open to suggestions.
> >>
> >> I have a completely within design with 2 conditions (condition and 
> >> palette).
> >>
> >> Here is the aov version:
> >>
> >>> fit.aov <- (aov(correct ~ cond * palette + Error(subject),
> >> data=data))
> >>> summary(fit.aov)
> >>
> >> Error: subject
> >>          Df  Sum Sq Mean Sq F value Pr(>F) Residuals 15
> >> 0.17326 0.01155
> >>
> >> Error: Within
> >>             Df  Sum Sq Mean Sq F value    Pr(>F)
> >> cond          1 0.32890 0.32890  52.047 4.906e-09 ***
> >> palette       1 0.21971 0.21971  34.768 4.447e-07 ***
> >> cond:palette  1 0.50387 0.50387  79.735 1.594e-11 ***
> >> Residuals    45 0.28437 0.00632
> >> ---
> >> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> >>
> >> and here is the lme version:
> >>
> >>> fm1 <- lme(correct ~ cond * palette, random=~1 | subject,
> >> data=data)  > anova(fm1)
> >>             numDF denDF  F-value p-value
> >> (Intercept)      1    45 4031.042  <.0001
> >> cond             1    45   52.047  <.0001
> >> palette          1    45   34.768  <.0001
> >> cond:palette     1    45   79.735  <.0001
> >>
> >> Thanks so much!
> >> Greg
> >>
> >> ______________________________________________
> >> R-help at r-project.org mailing list
> >> https://stat.ethz.ch/mailman/listinfo/r-help
> >> PLEASE do read the posting guide
> >> http://www.R-project.org/posting-guide.html
> >> and provide commented, minimal, self-contained, reproducible code.
> >>
> 
> 



More information about the R-help mailing list