[R-sig-ME] Calculating an Intraclass Correlation Coefficient or Repeatability estimate from lmer() output

Jarrod Hadfield j.hadfield at ed.ac.uk
Wed Oct 14 16:38:44 CEST 2009


Hi,

You have the components the wrong way round:  the ID variance is the  
among groups, and the Residual variance is the within groups, but  
other than that it looks fine.

Jarrod
On 14 Oct 2009, at 13:24, amelie lescroel wrote:

> Dear all,
>
> I am interested in calculating an Intraclass Correlation Coefficient  
> (ICC) or Repeatability estimate from a mixed model output (lmer()).  
> Lessels & Boag (1987) defined repeatability as the intraclass  
> correlation coefficient based on variance components derived from a  
> one-way anova:
> r = among-groups variance components / (within-group variance  
> components + among-group variance components)
> What I would like to know is where/how to find these variance  
> components in a lmer() output.
>
> Specifically, we measured foraging efficiency (CPUE) of birds during  
> 10 consecutive years. We have several measures per bird for each  
> year and the same birds were measured over multiple years. What I  
> would like to get is an estimate of the intra-individual consistency  
> of foraging efficiency over time. And I thought that the ICC or  
> repeatability would be this estimate (recently used in the same  
> manner in Lang et al. 2009 Ecology 90(9): 2513-2523).
>
> Let’s run a model with the ID of the birds as a random factor:
>
>> (fm1 <- lmer(log.CPUE~(1|ID)))
> Linear mixed model fit by REML
> Formula: log.CPUE ~ (1 | ID)
>   AIC   BIC logLik deviance REMLdev
> -1144 -1126  575.2    -1158   -1150
> Random effects:
> Groups   Name        Variance Std.Dev.
> ID       (Intercept) 0.010292 0.10145
> Residual             0.036602 0.19132
> Number of obs: 3320, groups: ID, 341
>
> Fixed effects:
>            Estimate Std. Error t value
> (Intercept)   0.3484     0.0068   51.24
>
> In this case, is the within-individual variance = 0.010 and the  
> among-individuals variance = 0.037?
>
> Thanks for your help,
>
> Amelie
>
> Amélie Lescroël
> Seabird ecologist
> URU 420, Université de Rennes I - Service du Patrimoine Naturel,  
> Museum National d’Histoire Naturelle
> 263 Av. du Gal Leclerc
> 35042 Rennes Cedex
> France
>
>
>
>
> __________ Information from ESET Mail Security, version of virus  
> signature database 4506 (20091014) __________
>
> The message was checked by ESET Mail Security.
> http://www.eset.com
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.




More information about the R-sig-mixed-models mailing list