arrays - R 应用数组输出维度

Question

与此问题相关:Multiply multidimensional array with same-sized matrix

如果我创建以下数组和矩阵:

a <- array(1, dim=c(2,2,3))
b <- matrix(c(1,1,1,1), nrow=2)

并且想要执行 b 的元素乘法，并且 a 的每个切片沿第三维移动，我会使用 apply 。但是，我得到以下奇怪的结果:

> dim(apply(a, 3, `*`, b))
[1] 4 3
> newa <- array(0, dim=c(2,2,3))
> newa[] <- apply(a, 3, `*`, b)
> dim(newa)
[1] 2 2 3

为什么这两件事给出不同的答案？

Answer 1

所得尺寸为 apply(a, 3, * , b)是 4 3因为根据文档，请参阅 ?apply :

If each call to FUN returns a vector of length n, then apply returns an array of dimension c(n, dim(X)[MARGIN]) if n > 1.

在这种情况下，每次调用FUN返回长度为 4 的向量尺寸为 MARGIN=3是 3 。注意每次调用FUN的结果在apply被强制为 vector不是matrix .

现在，与 newa[] <- array(... ，您将替换 newa 中的值对象的值来自 array计算。 newa 的属性没有改变，所以它的尺寸仍然是 2 2 3 .

希望这有帮助。