Description Usage Arguments Value Author(s) References See Also Examples
Repeatability calculations based on Analysis of Variance (ANOVA).
1 | rpt.aov(y, groups, CI=0.95, npermut=1000)
|
y |
Vector of measurements. |
groups |
Vector of group identities (will be converted to a factor). |
CI |
Width of the confidence interval between 0 and 1 (defaults to 0.95). |
npermut |
Number of permutations used when calculating asymptotic P values (defaults to 1000). |
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 (here: "ANOVA"). |
R |
Point estimate for the repeatability (denoted as R). |
se |
Asymptotic standard error for repeatability (ICC) based Becker (1982). |
CI |
Asymptotic confidence interval for repeatability based on non-parametric bootstrapping. |
P |
Named vector of two P values (significance tests): |
R.permut |
Repeatability R estimates for each permutation run. |
Holger Schielzeth (holger.schielzeth@ebc.uu.se) & Shinichi Nakagawa (shinichi.nakagawa@otago.ac.nz)
Becker, W. A. (1992) A manual of quantitative genetics. 5th edn. Academic Enterprises, Pullman, WA.
Lessells, C. M. and Boag, P. T. (1987) Unrepeatable repeatabilities: a common mistake. Auk 104: 116-121.
Nakagawa, S. and Schielzeth, H. (2010) Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biological Reviews 85: 935-956
rpt.corr, rpt.remlLMM, rpt.mcmcLMM, 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.aov(Tarsus, BirdID, npermut=10)) # reduced number of npermut iterations
detach(BodySize)
# repeatability estimation for weight (body mass) - a lower R than the previous one
data(BodySize)
attach(BodySize)
(rpt.Weight <- rpt.aov(Weight, BirdID, npermut=10)) # reduced number of npermut iterations
detach(BodySize)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.