[R-sig-eco] ICC confidence intervals and power analysis for random effects in lmer?
Bob O'Hara
bohara at senckenberg.de
Thu Apr 12 08:50:20 CEST 2012
On 04/12/2012 03:50 AM, Bradley Carlson wrote:
> I'm performing an analysis of behavioral variation among individual
> tadpoles, using individual ID as a random effect and time as a continuous
> fixed covariate in the lmer function in lmer4 package. I'm really
> interested in making inferences about the random effect (i.e. the extent of
> variation among individuals). I'd like to do two things that I can't seem
> to find straightforward answers about and I'm hoping someone can help or
> point me to a good resource.
>
> 1) The intraclass correlation coefficient is of particular interest to me,
> as it describes the proportion of variation that occurs among individuals.
> Ideally I'd like to report a confidence interval of the ICC but I can't
> find any way to calculate one, other than a function in the psychometric
> package that appears to only work when there are no covariates in the model
> (random effect only).
MCMC has already been mentioned and lme4 still has its mcmcsamp()
function. Failing that, you could try a parametric bootstrap, which
requires a little bit of coding but simulate() makes it much easier.
> 2) A reviewer requested a power analysis of the ability to detect a
> significant random effect. Any tips on how to approach that?
Report the random effect and confidence intervals. Retrospective power
analyses are pretty pointless (e.g. see
http://beheco.oxfordjournals.org/content/14/3/446.full), unless you're
planning to repeat the experiment.
Bob
--
Bob O'Hara
Biodiversity and Climate Research Centre
Senckenberganlage 25
D-60325 Frankfurt am Main,
Germany
Tel: +49 69 798 40226
Mobile: +49 1515 888 5440
WWW: http://www.bik-f.de/root/index.php?page_id=219
Blog: http://blogs.nature.com/boboh
Journal of Negative Results - EEB: www.jnr-eeb.org
More information about the R-sig-ecology
mailing list