Paneliperf: agree.ind

data(chocolates)
res<-paneliperf(chocolates, formul = "~ Product + Panelist + Session + Product:Panelist + Product:Session + Panelist:Session", formul.j = "~Product", col.j = 1, firstvar = 5, synthesis = TRUE)
resagree<-magicsort(res$agree, sort.mat = res$r2.ind, method = "median")
coltable(resagree, level.lower = 0.00, level.upper = 0.85, main.title = "Agreement between panelists")