R语言：如何移除基因表达矩阵中指定基因向量对应的行？

嘿，这个需求其实用基础R的索引功能就能轻松搞定！你已经知道怎么选中要排除的基因，反向操作只需要取反筛选条件或者明确保留的基因集合就行，下面给你两种常用的方法：

方法一：直接用逻辑取反筛选 #

这是最直接的写法，一步到位：

# 移除excluded_genes中的基因，保留剩余所有基因
filtered_expression_matrix <- expression_matrix[!rownames(expression_matrix) %in% excluded_genes, ]

拆解一下这个代码的逻辑：

• rownames(expression_matrix)：获取矩阵所有行（也就是所有基因）的名称
• %in% excluded_genes：判断每个基因是否在你的排除列表里，返回一个和行数一致的逻辑向量（TRUE表示在排除列表，FALSE表示不在）
• 前面加!：把逻辑向量取反，变成“不在排除列表里的基因”才为TRUE
• 最后用这个取反后的逻辑向量筛选矩阵的行，就得到了移除指定基因后的矩阵

方法二：先确定保留的基因集合再筛选 #

如果觉得逻辑取反有点绕，也可以先明确要保留的基因，再进行筛选，可读性更强：

# 找出所有不在排除列表里的基因
keep_genes <- setdiff(rownames(expression_matrix), excluded_genes)
# 用保留的基因筛选矩阵
filtered_expression_matrix <- expression_matrix[keep_genes, ]

这里setdiff()函数的作用是返回第一个集合（矩阵行名）中存在，但第二个集合（排除基因）中不存在的元素，正好就是你需要保留的基因列表，之后直接用这个列表筛选行就可以了。

小提醒 #

• 一定要确保excluded_genes里的基因名称和矩阵行名完全匹配（包括大小写、拼写、空格等），不然会出现漏删或者误删的情况
• 如果你的矩阵没有设置行名，那可以用行号来操作：假设排除的基因对应的行号是excluded_row_numbers，那么代码就是expression_matrix[-excluded_row_numbers, ]（负号表示移除这些行）

内容的提问来源于stack exchange，提问作者Peter Nilsson