# \donttest{
library(VCA)
data(VCAdata1)
lst <- anovaVCA(y~(device+lot)/day/run, VCAdata1, by="sample")
mat <- get_mat(lst) # automatically selects "total"
res <- fit_vfp(model.no=1:9, Data=mat)
coef(res)
# }
Run the code above in your browser using DataLab