R 如何在数组中复制矩阵

R 如何在数组中复制矩阵,r,R,假设我有一个维数为2x2的矩阵a。我如何创建一个2x2xk数组,它只是矩阵a的k次迭代 我怀疑abind包是我的朋友,但到目前为止,我还没有弄清楚如何使用某种repa,k公式来复制整个矩阵,而不仅仅是元素。在abind公式中列出矩阵a k次,就像下面的2x2x6数组代码一样,我找不到一种方法来做这件事 a <- matrix(c(1,0,0,1), nrow=2) library(abind) axk <- abind(a,a,a,a,a,a,along=3) 我希望输出相当于ax

假设我有一个维数为2x2的矩阵a。我如何创建一个2x2xk数组,它只是矩阵a的k次迭代

我怀疑abind包是我的朋友,但到目前为止,我还没有弄清楚如何使用某种repa,k公式来复制整个矩阵,而不仅仅是元素。在abind公式中列出矩阵a 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次。

您可以这样做

> 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

天啊,我应该能猜出来的。非常感谢。