如何在数组中复制矩阵

Question

说我有一个尺寸为2x2的矩阵a。我如何创建一个仅2个矩阵a的k次迭代的2x2xk数组？

我怀疑abind包是我的朋友，但是到目前为止，我还没有弄清楚如何执行某种rep（a，k）公式来复制整个矩阵，而不仅仅是元素。除了在abind公式中k次列出矩阵之外，我找不到一种方法，就像下面的2x2x6数组的代码一样。

a <- matrix(c(1,0,0,1), nrow=2)
library(abind)
axk <- abind(a,a,a,a,a,a,along=3)

我希望输出等于axk，但无需列出矩阵k次。

Answer 1

你可以

> replicate(6, a)
, , 1

     [,1] [,2]
[1,]    1    0
[2,]    0    1

, , 2

     [,1] [,2]
[1,]    1    0
[2,]    0    1

, , 3

     [,1] [,2]
[1,]    1    0
[2,]    0    1

, , 4

     [,1] [,2]
[1,]    1    0
[2,]    0    1

, , 5

     [,1] [,2]
[1,]    1    0
[2,]    0    1

, , 6

     [,1] [,2]
[1,]    1    0
[2,]    0    1