通过 row.names 对矩阵进行子集化
Posted
技术标签:
【中文标题】通过 row.names 对矩阵进行子集化【英文标题】：Subsetting a matrix by 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)
但是，这个特定的代码似乎无法完成这项工作。如何以这种方式执行子集？
【问题讨论】：
【参考方案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), ]
【讨论】：
以上是关于通过 row.names 对矩阵进行子集化的主要内容，如果未能解决你的问题，请参考以下文章