Paneliperf: r2.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)
resr2<-magicsort(res$r2.ind, method = "median", ascending = FALSE)
coltable(resr2, level.lower = 0.00, level.upper = 0.85, main.title = "Adjusted R-square (by panelist)")