A to czterowymiarowa tablica ze ściemniaczem 100 * 100 * 100 * 100. Chcę wybrać 10000 sub macierzy z ostatnich dwóch wymiarów A. B i C są wektorami o długości 10000. Są to kryteria wyboru. B określa numer wiersza A, a C określa numer kolumny.R, wybierz sekwencję macierzy podrzędnej z wysokiej tablicy wymiarowej
A <- array(rnorm(100^4), dim=c(100,100,100,100))
B <- sample(nrow(A) , 10000 , repl = TRUE)
C <- sample(ncol(A) , 10000 , repl = TRUE)
D <- array(0, dim=c(10000,100,100))
pomocą pętli:
system.time(
for (i in 1:10000){
D[i,,] <- A[B[i],C[i],,]
})
user system elapsed
10.20 0.14 10.34
z mapply:
sub_array <- function(b,c) return(A[b,c,,])
system.time(D <- mapply(FUN = sub_array, B, C, SIMPLIFY='array'))
user system elapsed
9.77 3.75 29.17
która jest nawet dłuższy. Czy jest to szybszy sposób? Dzięki.