Sottrazione di n, dimensioni k dalla matrice di matrici di n, k dimensioni

Question

Se ho un array A

A <- array(0, c(4, 3, 5)) 
for(i in 1:5) { 
    set.seed(i) 
    A[, , i] <- matrix(rnorm(12), 4, 3) 
} 

e se devo matrice B

set.seed(6) 
B <- matrix(rnorm(12), 4, 3) 

Il codice di sottrarre B dal ciascuna matrice della matrice a sarebbe:

d<-array(0, c(4,3,5)) 
for(i in 1:5){ 
    d[,,i]<-A[,,i]-B 
} 

Tuttavia, quale sarebbe il codice per eseguire lo stesso calcolo utilizzando una funzione da un" pply "famiglia?

Answer 1

Questo è ciò che è per sweep.

sweep(A, 1:2, B) 

Answer 2

Forse non molto intuitivo:

A[] <- apply(A, 3, `-`, B) 

Answer 3

Poiché state loop sull'ultima dimensione dell'array, si può semplicemente fare:

d <- A - as.vector(B) 

e sarà molto più veloce. È la stessa idea di quando sottrai un vettore da una matrice: il vettore viene riciclato, quindi viene sottratto a ogni colonna.