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