data(chocolates)
resperf<-paneliperf(chocolates, formul = "~Product+Panelist+Product:Panelist", formul.j = "~Product",
col.j = 1, firstvar = 5, lastvar = 12, synthesis = FALSE, graph = FALSE)
resperfprob<-magicsort(resperf$prob.ind, method = "median")
coltable(resperfprob, level.lower = 0.05, level.upper = 1, main.title = "P-value of the F-test (by panelist)")