data(croissant)
# Use a subset of 40 subjects for a faster example
subjects <- unique(croissant[[1]])[1:40]
croissant40 <- croissant[croissant[[1]] %in% subjects, ]
prepro <- preprocess_JAR_liking(
Data = croissant40,
nprod = 6,
nsub = 40,
liking_col = "OVL_Overall Liking",
levelsJAR = 5,
scale = FALSE
)
CATA <- prepro$CATA
liking <- prepro$liking
Data <- prepro$Datafinal
res <- cluscata_liking(
Data,
nblo = 40,
NameBlocks = prepro$NameSub,
printlevel = FALSE
)
summary(res)
plot(
res,
prepro$CATA,
prepro$liking,
scale = FALSE
)
Run the code above in your browser using DataLab