cata (version 0.0.10.9)

bcluster: Wrapper function for b-cluster analysis

Description

By default, bcluster calls a function to perform b-cluster analysis by a non-hierarchical iterative ascent algorithm, then inspects results if there are multiple runs.

Usage

bcluster(X, inspect = TRUE, inspect.plot = TRUE, algorithm = "n", 
measure = "b", G = NULL, M = NULL, max.iter = 500, tol = exp(-32), 
runs = 1, seed = 2021, verbose = FALSE)

Value

list with elements:

• runs : b-cluster analysis results from bcluster.n or bcluster.h (in a list if runs>1)

• inspect : result from inspect (the plot from this function is rendered if inspect.plot is TRUE)

Arguments

X

three-way array with I assessors, J products, M attributes where CATA data have values 0 (not checked) and 1 (checked)

inspect

default (TRUE) calls the inspect function to evaluate all solutions (when runs>1)

inspect.plot

default (TRUE) plots results from the inspect function

algorithm

default is n for non-hierarchical; h for hierarchical

measure

default is b for the b-measure; Q for Cochran's Q test

G

number of clusters (required for non-hierarchical algorithm)

M

initial cluster memberships

max.iter

maximum number of iteration allowed (default 500)

tol

non-hierarchical algorithm stops if variance over 5 iterations is less than tol (default: exp(-32))

runs

number of runs (defaults to 1)

seed

for reproducibility (default is 2021)

verbose

maximum number of iterations

References

Castura, J.C., Meyners, M., Varela, P., & Næs, T. (2022). Clustering consumers based on product discrimination in check-all-that-apply (CATA) data. Food Quality and Preference, 104564. tools:::Rd_expr_doi("10.1016/j.foodqual.2022.104564").

Examples

data(bread)

# b-cluster analysis on the first 14 consumers and the first 6 attributes
(b1 <- bcluster(bread$cata[1:14,,1:6], G=2))
# identical to: 
# (b2 <- bcluster.n(bread$cata[1:10,,1:6], G=2))