- Gam = matrix(0, ncol = K, nrow = n)
- gam = Gam
- for (i in 1:n){
- for (r in 1:K){
- sqNorm2 = sum( (Y[i,]%*%model$rho[,,r]-X[i,]%*%model$phi[,,r])^2 )
- Gam[i,r] = model$pi[r] * exp(-0.5*sqNorm2)* det(model$rho[,,r])
- }
- gam[i,] = Gam[i,] / sum(Gam[i,])
- }
- affec = apply(gam, 1,which.max)