Logiciel R : structure array

[invitec5eb4b89]

Bonjour,

Depuis quelque jour, j'essaie sans succès de trouver une commande simple permettant de faire (par exemple) la somme de toutes les "matrices" d'un objet array (à trois dimensions) à travers la "troisième dimension"...

Par exemple, pour l'array M suivant :

> M <- array(1:8,c(2,2,2))
> M
, , 1

     [,1] [,2]
[1,]    1    3
[2,]    2    4

, , 2

     [,1] [,2]
[1,]    5    7
[2,]    6    8

Je voudrais calculer la matrice somme des deux sous matrices ! J'ai essayé les fonctions "apply", "sapply", "rapply" sans succès...
Pour l'instant je fais une boucle for :

> temp <- matrix(0,dim(M)[1],dim(M)[2])
> for (i in 1:dim(M)[3]) temp <- temp+M[,,i]
> temp
     [,1] [,2]
[1,]    6   10
[2,]    8   12

Mais ça ne me satisfait pas vraiment... Qu'en pensez-vous ?
V.
-----

[invitec5eb4b89]

Mon problème n'a pas déchaîné les foules !!

Je vais quand même poster une nouvelle solution, proposée par mon collègue de bureau, sous la forme d'une fonction prenant comme argument l'objet de type array en question et la fonction que l'on veut appliquer à tous les éléments de la troisième dimension :

fonc_array <- function(m,fonction){
  n <- dim(m)[1]; p <- dim(m)[2]; k <- dim(m)[3]
  temp <- matrix(as.matrix(m), nrow = n*p, ncol = k)
  res  <- matrix(apply(temp, 1, fonction), n,p)
}

V.