for(r in 1:k){
b[r] = sum(sum(abs(phi[,,r])))
}
- gam2 = sum(gam[1,]) #BIG DOUTE
+ gam2 = colSums(gam)
a = sum(gam*t(log(Pi)))
#tant que les props sont negatives
Pi = transpose(Pi)
return(list(phi=phi, rho=rho, Pi=Pi, LLF=LLF, S=S))
-}
\ No newline at end of file
+}