Question

我正在使用带有3D维数组的R。我试图将它用作不同时刻的一组2D矩阵。

我发现了一种我真的不理解的行为，我想知道为什么正在发生。我试图在这里和其他地方找到解释，但直到现在我仍然有疑问。

我有这样的3D数组：

array3D=array(1:45,c(5,3,3))

正如我所料，我可以访问单个2D矩阵

 array3D[1,,]
     [,1] [,2] [,3]
[1,]    1   16   31
[2,]    6   21   36
[3,]   11   26   41

然而，尝试访问两个2D矩阵我没有得到我期望的结果

 array3D[1:2,,]
, , 1

     [,1] [,2] [,3]
[1,]    1    6   11
[2,]    2    7   12

, , 2

     [,1] [,2] [,3]
[1,]   16   21   26
[2,]   17   22   27

, , 3

     [,1] [,2] [,3]
[1,]   31   36   41
[2,]   32   37   42

我发现我可以使用aperm(array3D[1:2,,])来解决这个问题，但我不了解正在做什么。

另一个问题是，当我尝试做作业时，我不明白为什么这不起作用

array3D[1:2,,]=matrix(9:1,3,3)
array3D[1,,]
     [,1] [,2] [,3]
[1,]    9    3    6
[2,]    7    1    4
[3,]    5    8    2

我认为我可以通过循环解决这个问题，或者可以使用aaply来解决这个问题，但我认为如果我想使用3D数组对于理解正在发生的事情非常重要。如果有人能指出我正确的方向，我将非常高兴。

我试图在这里找到答案并阅读http://adv-r.had.co.nz/，但到目前为止还没有运气。

更新

我发现一切都有效，如果不使用第一个索引我使用最后一个，但我仍然不明白为什么。是R固有的东西吗？是否可以以其他方式使用第一个？

array3D=array(1:45,c(3,3,5))
array3D[,,1:2]=matrix(9:1,3,3)
array3D[,,2]
     [,1] [,2] [,3]
[1,]    9    6    3
[2,]    8    5    2
[3,]    7    4    1

Answer 1

我认为你想要实现的目标并不十分清楚，但这里有一些例子：在第一点上，您可以通过以下方式从z方向选择三个三乘三矩阵中的两个：

array3D[,,1:2]

因此，您可以使用适当大小的数组替换：

array3D[,,1:2] <- array(18:1,c(3,3,2))

关于为什么必须使用第三个索引的问题：将其想象为3D坐标系中的z方向。行将是x方向（垂直），列是y方向（水平）。索引array3D[1:2,,]时，您选择了前两行，同时保持所有内容都在x和z方向。

R子集和多维数组中的分配

1 个答案: