Description Usage Arguments Value Author(s) References See Also Examples
Calculates repeatability based on inter-class correlations
1 | rpt.corr(y, groups, CI=0.95, nboot=1000, npermut=1000)
|
y |
Vector of measurements. Missing values are not allowed. |
groups |
Vector of group identitities (will be converted to a factor). Note that each group identity has to appear exactly twice. |
CI |
Width of the confidence interval between 0 and 1 (defaults to 0.95). |
nboot |
Number of bootstrapping runs used when calculating an asymptotic confidence interval (defaults to 1000). |
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: "corr"). |
R |
Point estimate for repeatability R. |
se |
Asymptotic standard error for repeatability based on non-parametric bootstrapping. |
CI |
Asymptotic confidence interval for repeatability based on non-parametric bootstrapping. |
P |
Asymptotic P value from a significance test for the intraclass correlation based on permutation. |
R.permut |
Permutation samples for R. |
Holger Schielzeth (holger.schielzeth@ebc.uu.se) & Shinichi Nakagawa (shinichi.nakagawa@otago.ac.nz)
Sokal, R. R. and F. J. Rohlf (1995). Biometry: The principles and practice of statistics in biological research. New York, W.H. Freeman and Company.
Nakagawa, S. and Schielzeth, H. (2010) Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biological Reviews 85: 935-956
rpt.aov, rpt.remlLMM, rpt.mcmcLMM, rpt, print.rpt
1 2 3 4 5 6 | # repeatability for male breeding success on a transformed scale
data(Fledglings)
Fledglings$sqrtFledge <- sqrt(Fledglings$Fledge)
attach(Fledglings)
(rpt.Fledge <- rpt.corr(sqrtFledge, MaleID, nboot=10, npermut=10)) # reduced number of iterations
detach(Fledglings)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.