rpt.aov: ANOVA-based Repeatability

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

View source: R/rpt.aov.R

Description

Repeatability calculations based on Analysis of Variance (ANOVA).

Usage

1
	rpt.aov(y, groups, CI=0.95, npermut=1000)

Arguments

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).

Value

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): P.aov is the P value for the ANOVA F test, P.permut is the permutation based P value.

R.permut

Repeatability R estimates for each permutation run.

Author(s)

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

References

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

See Also

rpt.corr, rpt.remlLMM, rpt.mcmcLMM, rpt, print.rpt

Examples

 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)

rptR documentation built on May 31, 2017, 2:53 a.m.