J'ai un tableau Sigma à 3 dimensions m x n x p ainsi qu'une matrice A de dimensions l x m.
Je souhaite alors créer le tableau Omega à 3 dimensions l x n x p tel que
Code : Tout sélectionner
Omega[,,i] = A%*%Sigma[,,i]
Pour cela je peux faire la boucle:
Code : Tout sélectionner
Omega <- array(NA, dim=c(l,n,p))
for(i in 1:p){
Omega[,,i] = A%*%Sigma[,,i]
}
mais y aurait-il un moyen plus direct ? Je n'y parviens pas avec apply().