don <- read.table("QDA_saxo_2021.csv",header=TRUE, sep=";",stringsAsFactors = TRUE)
for (i in 1:3) don[,i] <- as.factor(don[,i])
summary(don)
library(SensoMineR)

## Caractérisation des produits
resdecat<-decat(don, formul="~Morceau+Nom+Morceau:Nom+Session+Nom:Session+Morceau:Session", firstvar = 4)
barrow(resdecat$tabT)
barrow(t(resdecat$tabT), numr = 3, numc = 3)

## Construction espace produit 
res <- panellipse(don, col.p = 2, col.j = 1, firstvar = 4)
coltable(res$hotelling, main.title = "P-values for the Hotelling's T2 tests")

## Performance du jury
res=panelperf(don, firstvar = 4, formul = "~Morceau+Nom+Session+Morceau:Nom+Morceau:Session+Nom:Session")
coltable(magicsort(res$p.value, sort.mat = res$p.value[,1], bycol = FALSE,
                   method = "median"), main.title = "Panel performance (sorted by product P-value)")


## Performances individuelles
res=paneliperf(don, firstvar = 4, formul = "~Morceau+Nom+Session+Morceau:Nom+Morceau:Session+Nom:Session", formul.j="~Morceau", col.j=1)
resprob<-magicsort(res$prob.ind, method = "median")
coltable(resprob, level.lower = 0.05, level.upper = 1,
         main.title = "P-value of the F-test (by panelist)")