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 :

Code:
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.