我目前有一个数组:
vector1 <- c(5,9,3)
vector2 <- c(10,11,12,13,14,15)
new.array <- array(c(vector1,vector2),dim = c(3,3,100))
print(new.array)
, , 1
[,1] [,2] [,3]
[1,] 5 10 13
[2,] 9 11 14
[3,] 3 12 15
, , 2
[,1] [,2] [,3]
[1,] 5 10 13
[2,] 9 11 14
[3,] 3 12 15
...
我想知道如何从每个矩阵的第三行中减去第二行,同时保留数组结构。例如,我想获得:
, , 1
[,1] [,2] [,3]
[1,] -6 1 1
, , 2
[,1] [,2] [,3]
[1,] -6 1 1
...
感谢。
答案 0 :(得分:8)
这样做:
new.array[3, , , drop = FALSE] - new.array[2, , , drop = FALSE]
drop = FALSE
是为了确保数组的子集具有与原始数组相同的维度。
答案 1 :(得分:5)
你可以这样做:
如果你运行dim(new.array)
,它将按行,列和最后一维的顺序返回三个维数,表示你的矩阵数。按顺序访问任何特定的矩阵。你必须这样做:array [,, 1]将获取该特定数组的第一个矩阵,依此类推。在这种情况下,形成的矩阵总数为100.运行dim(new.array)
并亲眼看看。在这100个小矩阵上使用lapply
。你可以从第二行中减去第三行。如下。
> lapply(1:100, function(x)new.array[,,x][2,] - new.array[,,x][3,])