rptR: repeatability estimation and variance decomposition by generalized linear mixed‐effects models

Summary Intra-class correlations (ICC) and repeatabilities (R) are fundamental statistics for quantifying the reproducibility of measurements and for understanding the structure of biological variation. Linear mixed effects models offer a versatile framework for estimating ICC and R. However, while point estimation and significance testing by likelihood ratio tests is straightforward, the quantification of uncertainty is not as easily achieved. A further complication arises when the analysis is conducted on data with non-Gaussian distributions because the separation of the mean and the variance is less clear-cut for non-Gaussian than for Gaussian models. Nonetheless, there are solutions to approximate repeatability for the most widely used families of generalized linear mixed models (GLMMs). Here, we introduce the R package rptR for the estimation of ICC and R for Gaussian, binomial and Poisson-distributed data. Uncertainty in estimators is quantified by parametric bootstrapping and significance testing is implemented by likelihood ratio tests and through permutation of residuals. The package allows control for fixed effects and thus the estimation of adjusted repeatabilities (that remove fixed effect variance from the estimate) and enhanced agreement repeatabilities (that add fixed effect variance to the denominator). Furthermore, repeatability can be estimated from random-slope models. The package features convenient summary and plotting functions. Besides repeatabilities, the package also allows the quantification of coefficients of determination R2 as well as of raw variance components. We present an example analysis to demonstrate the core features and discuss some of the limitations of rptR.

[1]  K. McGraw,et al.  Forming inferences about some intraclass correlation coefficients. , 1996 .

[2]  Kate L. Laskowski,et al.  The repeatability of behaviour: a meta-analysis , 2009, Animal Behaviour.

[3]  Shinichi Nakagawa,et al.  General Methods for Evolutionary Quantitative Genetic Inference from Generalized Mixed Models , 2016, Genetics.

[4]  J. Fleiss,et al.  Intraclass correlations: uses in assessing rater reliability. , 1979, Psychological bulletin.

[5]  Paul C. Johnson Extension of Nakagawa & Schielzeth's R2GLMM to random slopes models , 2014, Methods in ecology and evolution.

[6]  Shinichi Nakagawa,et al.  Repeatability for Gaussian and non‐Gaussian data: a practical guide for biologists , 2010, Biological reviews of the Cambridge Philosophical Society.

[7]  J. Stamps,et al.  Using repeatability to study physiological and behavioural traits: ignore time-related change at your peril , 2015, Animal Behaviour.

[8]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[9]  W. Revelle psych: Procedures for Personality and Psychological Research , 2017 .

[10]  Matthew E. Wolak,et al.  Guidelines for estimating repeatability , 2012 .

[11]  D. Bates,et al.  Fitting Linear Mixed-Effects Models Using lme4 , 2014, 1406.5823.

[12]  Hadley Wickham,et al.  Tools to Make Developing R Packages Easier , 2016 .

[13]  Shinichi Nakagawa,et al.  A general and simple method for obtaining R2 from generalized linear mixed‐effects models , 2013 .

[14]  C. Boake Repeatability: Its role in evolutionary studies of mating behavior , 1989, Evolutionary Ecology.