dotProduct = tcrossprod(Y[i,]%*%Rho[,,r]-X[i,]%*%phi[,,r])
logGamIR = log(Pi[r]) + log(det(Rho[,,r])) - 0.5*dotProduct
#Z[i] = index of max (gam[i,])
dotProduct = tcrossprod(Y[i,]%*%Rho[,,r]-X[i,]%*%phi[,,r])
logGamIR = log(Pi[r]) + log(det(Rho[,,r])) - 0.5*dotProduct
#Z[i] = index of max (gam[i,])