n = 50; m = 10; p = 5 X = matrix(runif(n*p, -10, 10), nrow=n) Y = matrix(runif(n*m, -5, 15), nrow=n) beta = array(0, dim=c(p,m,2)) beta[,,1] = 1 beta[,,2] = 2 data = generateXY(n, c(0.4,0.6), rep(0,p), beta, diag(0.5, p), diag(0.5, m)) X = data$X Y = data$Y class = data$class V1 = runValse(X, Y, fast=FALSE) Error in while (!pi2AllPositive) { : missing value where TRUE/FALSE needed V2 = runValse(X, Y, fast=TRUE) list() Error in out[[ind_uniq[l]]] : attempt to select less than one element in get1index ==> Error identified: line 61 in initSmallEM.R, division by 0 It occurs also for smallers values of n and m, e.g.: n = 20; m = 20; p = 3 ===== Also: X <- matrix(runif(100), nrow=50) Y <- matrix(runif(100), nrow=50) (...) Error: cannot allocate vector of size 16.0 Gb