Subtract Matrix von n, k Dimensionen von Array von Matrizen von n, k Dimensionen

Question

Wenn ich ein Array haben eine

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

und wenn ich Matrix B

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

Der Code B zu subtrahieren von der jede Matrix der Matrix A wäre:

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

jedoch, was der Code wäre, die gleiche Berechnung unter Verwendung einer Funktion von „a auszuführen pled "Familie?

Answer 1

Dies ist, was sweep für ist.

sweep(A, 1:2, B) 

Answer 2

Vielleicht nicht sehr intuitiv:

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

Answer 3

Weil Sie auf der letzten Array-Dimension sind Looping, können Sie einfach tun:

d <- A - as.vector(B) 

und es wird viel schneller sein. Es ist die gleiche Idee wie wenn Sie einen Vektor von einer Matrix subtrahieren: Der Vektor wird recycelt, so dass er von jeder Spalte subtrahiert wird.