tests/testRep.R

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

require(SensMixed)
load(system.file("testdata","senscherry.RData",package="SensMixed"))

testRep <- FALSE

if(testRep){
#convert some variables to factors in TVbo
TVbo <- convertToFactors(TVbo, c("Assessor", "Repeat", "Picture"))
                         
res_paral <- sensmixed(c("Coloursaturation", "Colourbalance"),
                       prod_effects = c("TVset", "Picture"), replication="Repeat", 
                       assessor="Assessor", data=TVbo, error_structure="ASS-REP",
                       MAM = FALSE)

stopifnot(nrow(res_paral$random$Chi)==11)

senscherry$Replicate <- as.factor(senscherry$Replicate)
res_cherry<- sensmixed(c("Apple.aroma", "Cherry.aroma"),
                       prod_effects = c("Fla", "Fib"), replication="Replicate", 
                       assessor="Assessor", data = senscherry, 
                       error_structure="ASS-REP", MAM = FALSE)

stopifnot(nrow(res_cherry$random$Chi)==11)

res_cherryMAMmult<- sensmixed(c("Apple.aroma", "Cherry.aroma"),
                       prod_effects = c("Fla", "Fib"), replication="Replicate", 
                       assessor="Assessor", data = senscherry, 
                       error_structure="ASS-REP", MAM=TRUE)

res_cherryMAMmult$scaling

## with one sample
res_MAM <- sensmixed(c("Apple.aroma", "Cherry.aroma"),
                     prod_effects=c("Sample"), replication="Replicate", 
                     assessor="Assessor", data=senscherry, MAM=TRUE, 
                     control=list(reduce.random = FALSE))

res_MAM_BAL <- sensmixed(c("Apple.aroma", "Cherry.aroma"),
                     prod_effects=c("Sample"), replication="Replicate", 
                     assessor="Assessor", data=senscherry, 
                     control=list(MAM_balanced=TRUE))


stopifnot(all.equal(res_MAM$scaling$FScaling[,"Apple.aroma"],
                    res_MAM_BAL[[3]][,,1]["Scaling","F"], tol=1e-2))

stopifnot(all.equal(res_MAM$scaling$FScaling[,"Cherry.aroma"],
                    res_MAM_BAL[[3]][,,2]["Scaling","F"], tol=1e-2))

}