通过 row.names 对矩阵进行子集化答案

【问题标题】：Subsetting a matrix by row.names通过 row.names 对矩阵进行子集化
【发布时间】：2014-05-05 10:02:47
【问题描述】：

我有一个包含以下 row.names 的矩阵：

"X1"   "X5"   "X33"  "X37"  "X52"  "X566"

现在我只想选择与列表条目匹配的行，例如：

include_list <- c("X1", "X5", "X33")

我想我会做这样的事情：

data.subset <- subset(data, row.names == include_list)

但是，这个特定的代码似乎无法完成这项工作。如何以这种方式执行子集？

【问题讨论】：

标签： r data-manipulation

【解决方案1】：

设置一些假数据：

m <- matrix(1:30, 6, 5)
rownames(m) <- c("X1", "X5", "X33", "X37", "X52", "X566")
m
#      [,1] [,2] [,3] [,4] [,5]
# X1      1    7   13   19   25
# X5      2    8   14   20   26
# X33     3    9   15   21   27
# X37     4   10   16   22   28
# X52     5   11   17   23   29
# X566    6   12   18   24   30

这里可能最容易使用矩阵索引进行子集化 ([)：

include_list <- c("X1", "X5", "X33")
m[include_list, ]
#     [,1] [,2] [,3] [,4] [,5]
# X1     1    7   13   19   25
# X5     2    8   14   20   26
# X33    3    9   15   21   27

subset() 函数的替代方案：

subset(m, rownames(m) %in% include_list)
#      [,1] [,2] [,3] [,4] [,5]
# X1     1    7   13   19   25
# X5     2    8   14   20   26
# X33    3    9   15   21   27

【讨论】：

【解决方案2】：

这似乎也有效：

include_list <- head(read.csv("/Users/histelheim/include_list.csv", header = FALSE))
include_list <- c(do.call("cbind", include_list)) 
data[row.names(data) %in% include_list, ]

【讨论】：

【解决方案3】：

以防万一某人的 include_list 是类 list（而不是向量），那么这是对 Rich Scriven 的答案的修改，该答案基于this 在这种情况下有效。应添加unlist()：

m <- matrix(1:30, 6, 5)
rownames(m) <- c("X1", "X5", "X33", "X37", "X52", "X566")
include_list <- list("X1", "X5", "X33")


m[unlist(include_list), ]

【讨论】：