# \donttest{
data(eem_list)
dim_min <- 3 # minimum number of components
dim_max <- 7 # maximum number of components
nstart <- 25 # random starts for PARAFAC analysis, models built simulanuously, best selected
# cores <- parallel::detectCores(logical=FALSE) # use all cores but do not use all threads
cores <- 2 # package checks only run with 2 cores
maxit = 2500
ctol <- 10^-7 # tolerance for parafac
pfres_comps <- eem_parafac(eem_list, comps = seq(dim_min, dim_max),
normalise = TRUE, maxit = maxit, nstart = nstart, ctol = ctol, cores = cores)
## with a defined number of converging models
#pfres_comps <- eem_parafac(eem_list, comps = seq(dim_min, dim_max),
# normalise = TRUE, maxit = maxit, nstart = nstart, ctol = ctol,
# output = "all", strictly_converging = TRUE, cores = cores, verbose = TRUE)
pfres_comps2 <- eem_parafac(eem_list, comps = seq(dim_min, dim_max),
normalise = TRUE, maxit = maxit, nstart = nstart, ctol = ctol, cores = cores, output = "all")
# }
Run the code above in your browser using DataLab