Paneliperf: prob.ind

res<-paneliperf(chocolates, formul = "~Product+Panelist+Session+
Product:Panelist+Product:Session+Panelist:Session",
  formul.j = "~Product", col.j = 1, firstvar = 5, synthesis = TRUE)

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)")