data(chocolates)
resperf=panelperf(chocolates,formul="~Product+Panelist+Session+Product:Panelist+Session:Product+Panelist:Session",firstvar=5)
## Sort results by product p.values.
coltable(magicsort(resperf$p.value,sort.mat=resperf$p.value[,1],bycol=FALSE),main.title="Panel performance (sorted by product P-value)")
Run the code above in your browser using DataLab