# \donttest{
library(VCA)
data(VCAdata1)
lst <- anovaVCA(y~(device+lot)/day/run, VCAdata1, by="sample")
get_mat(lst) # automatically selects 'total'
# pooled version of intermediate precision (error+run+day)
get_mat(lst, 4:6)
# only repeatability ('error')
get_mat(lst, "error")
# use remlVCA instead
lst2 <- remlVCA(y~(device+lot)/day/run, VCAdata1, by="sample")
get_mat(lst2)
# }
Run the code above in your browser using DataLab