tests/testKEEPeffs.R

In SensMixed: Analysis of Sensory and Consumer Data in a Mixed Model Framework

require(SensMixed)

#convert some variables to factors in TVbo
TVbo <- convertToFactors(TVbo, c("Assessor", "Repeat", "Picture"))

res <- sensmixed(c("Coloursaturation", "Colourbalance"),
                 prod_effects = c("TVset", "Picture"), replication="Repeat", 
                 assessor="Assessor", data=TVbo, 
                 control = list(keep.effs = "Assessor"), MAM = FALSE)

stopifnot(res$step_res[[1]]$rand.table["Assessor","elim.num"] == "kept")
stopifnot(res$step_res[[2]]$rand.table["Assessor","elim.num"] == "kept")

res <- sensmixed(c("Coloursaturation", "Colourbalance"),
                 prod_effects = c("TVset", "Picture"), replication="Repeat", 
                 assessor="Assessor", data=TVbo, MAM = FALSE)

stopifnot(res$step_res[[1]]$rand.table["Assessor","elim.num"] == "kept")
stopifnot(res$step_res[[2]]$rand.table["Assessor","elim.num"] == "kept")
stopifnot(res$step_res[[1]]$rand.table["TVset:Picture:Assessor","elim.num"] == "kept")
stopifnot(res$step_res[[2]]$rand.table["TVset:Picture:Assessor","elim.num"] == "kept")