deltaRank = rangmax - rangmin + 1
Size = deltaRank^k
Rank = matrix(0, Size, k)
deltaRank = rangmax - rangmin + 1
Size = deltaRank^k
Rank = matrix(0, Size, k)
phi = array(0, dim=c(p,m,k,L*Size))
lvraisemblance = matrix(0, L*Size, 2)
phi = array(0, dim=c(p,m,k,L*Size))
lvraisemblance = matrix(0, L*Size, 2)