data(chocolates)
resaverage<-averagetable(chocolates, formul = "~Product+Panelist",
firstvar = 5)
coltable(magicsort(resaverage), level.upper = 6,level.lower = 4,
main.title = "Average by chocolate")
resdecat<-decat(chocolates, formul = "~Product+Panelist", firstvar = 5,
graph = FALSE)
resaverage.sort = resaverage[rownames(magicsort(resdecat$tabT)),
colnames(magicsort(resdecat$tabT))]
coltable(resaverage.sort, magicsort(resdecat$tabT),
level.lower = -1.96, level.upper = 1.96,
main.title = "Average by chocolate")
get(getOption("device"))()
pca(resaverage, scale.unit = TRUE)
Run the code above in your browser using DataLab