Parametric per trait repeatabilities
Estimates the variance in the sample not due to measurement error
CalcRepeatability(ID, ind.data)
ID
: identity of individualsind.data
: individual measurementsvector of repeatabilities
Requires at least two observations per individual
num.ind = length(iris[,1]) ID = rep(1:num.ind, 2) ind.data = rbind(iris[,1:4], iris[,1:4]+array(rnorm(num.ind*4, 0, 0.1), dim(iris[,1:4]))) CalcRepeatability(ID, ind.data)
Lessels, C. M., and Boag, P. T. (1987). Unrepeatable repeatabilities: a common mistake. The Auk, 2(January), 116-121.
Guilherme Garcia
Useful links