Description Usage Arguments Details Value Author(s) References See Also Examples

Calculates repeatability from a linear mixed-effects models fitted by MCMC

1 | ```
rpt.mcmcLMM(y, groups, CI=0.95, prior=NULL, verbose=FALSE, ...)
``` |

`y` |
Vector of a response values |

`groups` |
Vector of group identities |

`CI` |
Width of the Bayesian credible interval (defaults to 0.95) |

`prior` |
List of prior values passed to the MCMCglmm function in MCMCglmm (see there for more details). Default priors will be used if prior is |

`verbose` |
Whether or not MCMCglmm should print MH diagnostics are printed to screen. Defaults to FALSE. |

`...` |
Additonal arguements that are passed on to MCMCglmm (e.g. length of chain, thinning interval). |

Models are fitted using the MCMCglmm function in MCMCglmm.
Models are fitted with `prior=list(R=list(V=1,n=10e-2), G=list(G1=list(V=1,n=10e-2)))`

unless other priors are specified in the call.

Returns an object of class rpt that is a a list with the following elements:

`datatype` |
Response distribution (here: "Gaussian"). |

`method` |
Method used to calculate repeatability (intra-class correlation, ICC) (here: "MCMC"). |

`CI` |
Width of the Bayesian credibility interval. |

`R` |
Point estimate for repeatability (intra-class correlation, ICC), i.e. the mode of the posterior distribution. |

`se` |
Standard error ( |

`CI.R` |
Bayesian credibility interval for the repeatability (ICC) based on the posterior distribution of |

`P` |
Significace test, returned as |

`R.post` |
MCMC samples form the posterior distributions of |

Holger Schielzeth ([email protected]) & Shinichi Nakagawa ([email protected])

Carrasco, J. L. and Jover, L. (2003). *Estimating the generalized concordance correlation coefficient through variance components*. Biometrics 59: 849-858.

Nakagawa, S. and Schielzeth, H. (2010) *Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists*. Biological Reviews 85: 935-956

rpt.remlLMM, rpt.aov, rpt.corr, rpt, print.rpt

1 2 3 4 5 6 7 8 9 10 11 | ```
# repeatability estimation for tarsus length - a very high R
data(BodySize)
attach(BodySize)
(rpt.BS <- rpt.mcmcLMM(Tarsus, BirdID))
detach(BodySize)
# repeatability estimation for weight (body mass) - a lower R than the previous one
data(BodySize)
attach(BodySize)
(rpt.Weight <- rpt.mcmcLMM(Weight, BirdID))
detach(BodySize)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.