# NOT RUN {
data(feems)
cube <- feemscale(
feemscatter(feemcube(feems, FALSE), rep(24, 4)), na.rm = TRUE
)
# will take time
sh <- feemsplithalf(cube, 2:4, random = 5, const = rep('nonneg', 3))
plot(sh)
plot(sh, 'factors')
# }
Run the code above in your browser using DataLab