Question

我在3个不同的向量中有来自物理对象的三个轴的数据。

outer(x, y, FUN="*")
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
 [1,]   12    8   20   52   88   80   92   76  124    56
 [2,]   18   12   30   78  132  120  138  114  186    84
 [3,]   12    8   20   52   88   80   92   76  124    56
 [4,]   18   12   30   78  132  120  138  114  186    84
 [5,]   51   34   85  221  374  340  391  323  527   238
 [6,]   30   20   50  130  220  200  230  190  310   140
 [7,]   39   26   65  169  286  260  299  247  403   182
 [8,]   81   54  135  351  594  540  621  513  837   378
 [9,]   99   66  165  429  726  660  759  627 1023   462
[10,]   57   38   95  247  418  380  437  361  589   266

我想通过简单地将它们的全部产品组合成三维矩阵（如果存在于R中）。如果我只有两个轴，说“x”和“y”我可以做以下

{{1}}

但是，因为我有三个轴，我不知道该怎么做。在“x”上使用外部＆amp; “y”然后将其乘以“z”与使用“x”＆amp;上的外部不同“z”并乘以“y”。生成3D矩阵后的目标是绘制3D矩阵。

我已经查看了其他类似的帖子，但发现它令人困惑Creating a 3D matrix with R?

Answer 1

将outer放入Reduce后，我们可以vector与list一起使用。

res1 <- Reduce(outer, list(x,y,z))

或另一个选项expand.grid rowProds来自library(matrixStats)，然后将产品向量（'v1'）转换为array。

library(matrixStats)
v1 <- rowProds(as.matrix(expand.grid(x,y,z)))
res2 <- array(v1, dim=c(length(x), length(y), length(z)))
identical(res1, res2)
#[1] TRUE
dim(res2)
#[1] 10 10 10

数据

x <- c(4, 6, 4, 6, 17, 10, 13, 27, 33, 19)
y <- c(3, 2, 5, 13, 22, 20, 23, 19, 31, 14)
z <- c(21, 6, 12, 35, 19, 44, 34, 23, 26, 33)

将三个向量组合成R中的3d矩阵对象

1 个答案:

数据