# \donttest{
library(VCA)
data(CA19_9)
fits.CA19_9 <- anovaVCA(result~site/day, CA19_9, by="sample")
# extract repeatability
mat.CA19_9 <- get_mat(fits.CA19_9, "error")
res.CA19_9 <- fit_vfp(mat.CA19_9, 1:9)
# default is on variance scale
residuals(res.CA19_9)
# residuals on cv scale for model 6
resid(res.CA19_9, model.no=6, type="cv")
# }
Run the code above in your browser using DataLab