# NOT RUN {
plot(x <- feem(matrix(1:42, 7), 400:406, 350:355))
data(feems)
dataset <- feemcube(feems, FALSE)[1:30*6, 1:9*6,]
dataset <- feemscatter(dataset, rep(24, 4), 'pchip')
dataset <- feemife(dataset, absorp)
plot(dataset <- feemscale(dataset, na.rm = TRUE))
# }
# NOT RUN {
# takes a long time
(sh <- feemsplithalf(
dataset, nfac = 2:5, const = rep('nonneg', 3), splits = 4)
)
plot(sh)
jk <- feemjackknife(dataset, nfac = 3, const = rep('nonneg', 3))
plot(jk)
# }
# NOT RUN {
pf <- feemparafac(dataset, nfac = 2, const = rep('nonneg', 3))
plot(pf)
# }
Run the code above in your browser using DataLab