TODO

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