Cách chức năng để xếp chồng danh sách ma trận 2d vào ma trận 3d

Question

Sau khi thông minh lapply, tôi đã để lại danh sách ma trận 2 chiều.

Ví dụ:

set.seed(1) 
test <- replicate(5, matrix(runif(25),ncol=5), simplify=FALSE) 
> test 
[[1]] 
      [,1]  [,2]  [,3]  [,4]  [,5] 
[1,] 0.8357088 0.29589546 0.9994045 0.2862853 0.6973738 
[2,] 0.2377494 0.14704832 0.0348748 0.7377974 0.6414624 
[3,] 0.3539861 0.70399206 0.3383913 0.8340543 0.6439229 
[4,] 0.8568854 0.10380669 0.9150638 0.3142708 0.9778534 
[5,] 0.8537634 0.03372777 0.6172353 0.4925665 0.4147353 

[[2]] 
      [,1]  [,2]  [,3]  [,4]  [,5] 
[1,] 0.1194048 0.9833502 0.9674695 0.6687715 0.1928159 
[2,] 0.5260297 0.3883191 0.5150718 0.4189159 0.8967387 
[3,] 0.2250734 0.2292448 0.1630703 0.3233450 0.3081196 
[4,] 0.4864118 0.6232975 0.6219023 0.8352553 0.3633005 
[5,] 0.3702148 0.1365402 0.9859542 0.1438170 0.7839465 

[[3]] 
... 

Tôi muốn biến chúng thành một mảng 3 chiều:

set.seed(1) 
replicate(5, matrix(runif(25),ncol=5))  

Rõ ràng, nếu tôi đang sử dụng sao chép Tôi chỉ có thể bật simplify, nhưng sapply không đơn giản hóa kết quả đúng cách và stack không hoàn toàn. do.call(rbind,mylist) biến nó thành ma trận 2d thay cho mảng 3d.

Tôi có thể làm điều này với một vòng lặp, nhưng tôi đang tìm một cách gọn gàng và chức năng để xử lý nó.

Cách gần nhất tôi đã đưa ra là:

array(do.call(c, test), dim=c(dim(test[[1]]),length(test))) 

Nhưng tôi có cảm giác như đó là không thanh nha (vì nó disassembles và sau đó ráp các thuộc tính mảng của các vectơ, và cần rất nhiều thử nghiệm để làm cho an toàn (ví dụ rằng kích thước của mỗi phần tử đều giống nhau).

Answer 1

Bạn có thể sử dụng gói abind và sau đó sử dụng do.call(abind, c(test, along = 3))

library(abind) 
testArray <- do.call(abind, c(test, along = 3)) 

Hoặc bạn có thể sử dụng simplify = 'array' trong một cuộc gọi đến sapply, (thay vì lapply). simplify = 'array' là không giống như simplify = TRUE, vì nó sẽ thay đổi lập luận higher trong simplify2array

ví dụ

foo <- function(x) matrix(1:10, ncol = 5) 
# the default is simplify = TRUE 
sapply(1:5, foo) 
     [,1] [,2] [,3] [,4] [,5] 
[1,] 1 1 1 1 1 
[2,] 2 2 2 2 2 
[3,] 3 3 3 3 3 
[4,] 4 4 4 4 4 
[5,] 5 5 5 5 5 
[6,] 6 6 6 6 6 
[7,] 7 7 7 7 7 
[8,] 8 8 8 8 8 
[9,] 9 9 9 9 9 
[10,] 10 10 10 10 10 
# which is *not* what you want 
# so set `simplify = 'array' 
sapply(1:5, foo, simplify = 'array') 
, , 1 

    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 5 7 9 
[2,] 2 4 6 8 10 

, , 2 

    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 5 7 9 
[2,] 2 4 6 8 10 

, , 3 

    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 5 7 9 
[2,] 2 4 6 8 10 

, , 4 

    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 5 7 9 
[2,] 2 4 6 8 10 

, , 5 

    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 5 7 9 
[2,] 2 4 6 8 10 

Answer 2

test2 <- unlist(test) 
dim(test2) <- c(dim(test[[1]]),5) 

---

hoặc nếu bạn không biết kích thước dự kiến ​​trước thời hạn:

dim3 <- c(dim(test[[1]]), length(test2)/prod(dim(test[[1]]))) 
dim(test2) <- dim3 

Answer 3

Mảng chỉ đơn giản là vec nguyên tử tor với kích thước. Mỗi thành phần ma trận của test thực sự chỉ là một vectơ có kích thước. Do đó giải pháp đơn giản nhất mà tôi có thể nghĩ đến là để bỏ danh sách test thành một vectơ và chuyển đổi nó thành một mảng bằng cách sử dụng array và các thứ nguyên được cung cấp phù hợp.

set.seed(1) 
foo <- replicate(5, matrix(runif(25),ncol=5)) 
tmp <- array(unlist(test), dim = c(5,5,5)) 

> all.equal(foo, tmp) 
[1] TRUE 
> is.array(tmp) 
[1] TRUE 
> dim(tmp) 
[1] 5 5 5 

Nếu bạn không muốn hardcode kích thước, chúng ta phải thực hiện một số giả định nhưng có thể dễ dàng điền vào các chiều kích từ test, ví dụ

tmp2 <- array(unlist(test), dim = c(dim(test[[1]]), length(test))) 

> all.equal(foo, tmp2) 
[1] TRUE 

Giả định rằng kích thước của mỗi thành phần đều giống nhau, nhưng sau đó tôi không thấy cách bạn có thể đặt ma trận phụ vào mảng 3-d nếu điều kiện đó không giữ.

Điều này có vẻ khó hiểu, để hủy danh sách, nhưng điều này đơn giản là khai thác cách R xử lý ma trận và mảng dưới dạng vectơ có kích thước.

Answer 4

Hãy thử điều này:

simplify2array(test)