矩阵列表上的R冒号运算符

Question

我在R中创建了一个矩阵列表。在列表的所有矩阵中，我想“拉出”特定索引的矩阵元素集合。 我当时认为冒号运算符可能允许我在一行中实现它。 例如，这是尝试访问列表中所有矩阵的[1,1]元素：

myList = list() #list of matrices
myList[[1]] = matrix(1:9, nrow=3, ncol=3, byrow=TRUE) #arbitrary data
myList[[2]] = matrix(2:10, nrow=3, ncol=3, byrow=TRUE)

#I expected the following line to output myList[[1]][1,1], myList[[2]][1,1]
slice = myList[[1:2]][1,1] #prints error: "incorrect number of dimensions"

上面代码的最后一行会引发错误“维数不正确”。

作为参考，这是我正在尝试做的工作（但不太优雅）的实现：

#assume myList has already been created (see the code snippet above)
slice = c()
for(x in 1:2) {
    slice = c(slice, myList[[x]][1,1])
}
#this works. slice = [1 2]

有谁知道如何在一行中完成上述操作？

请注意，我的“矩阵列表”可以替换为其他内容。 如果有人可以建议一个允许我执行上述操作的替代“矩阵集合”数据结构，那么这将得到解决。

也许这个问题很愚蠢......我真的希望有一个干净的单线实现。

Answer 1

两件事情。 首先， [和[[ 。 来自?'['的相关句子：

[，[[和$]之间最重要的区别是[可以选择多个元素，而另外两个选择单个元素。

所以你可能想做myList[1:2] 。 其次，您无法以您描述的方式组合子集化操作。 完成myList[1:2]您将获得两个矩阵的列表。 列表通常只有一个维度，因此在您的情况下执行myList[1:2][1,1]是没有意义的。 （参见有关例外情况的评论。）

你可以尝试lapply ： lapply(myList,'[',1,1) 。

Answer 2

如果您的矩阵都具有相同的尺寸，您可以将它们存储在三维数组中。 这肯定会使索引和提取元素更容易......

## One way to get your data into an array
a <- array(c(myList[[1]], myList[[2]]), dim=c(3,3,2))

## Extract the slice containing the upper left element of each matrix
a[1,1,]
# [1] 1 2

Answer 3

这有效：

> sapply(myList,"[",1,1)
[1] 1 2

编辑：哦，对不起，我在早期答案结束时看到几乎相同的想法。 但无论如何，sapply可能更接近你想要的东西