R语言移除低/零方差基因时rowVars报错，寻求替代方案

错误原因是matrixStats::rowVars()仅支持矩阵（matrix）或类矩阵对象，不直接兼容data.frame类型的raw_counts。以下是几种可行的替代方案：

方案1：转换为矩阵后使用rowVars #

先把data.frame转为矩阵，再调用rowVars：

library(matrixStats)
# 将data.frame转换为数值矩阵
raw_counts_mat <- as.matrix(raw_counts)
# 计算行方差
read_vars <- rowVars(raw_counts_mat)
# 筛选非零方差基因
raw_counts <- raw_counts[which(read_vars > 1e-10), ]
# 查看处理后的维度
dim(raw_counts)

方案2：用base R的apply手动计算行方差 #

无需额外依赖包，直接用base R的apply函数逐行计算方差：

# 计算每行的方差，na.rm=TRUE用于处理可能存在的NA值
read_vars <- apply(raw_counts, 1, var, na.rm = TRUE)
# 筛选方差大于1e-10的行
raw_counts <- raw_counts[which(read_vars > 1e-10), ]
dim(raw_counts)

方案3：用tidyverse工具链实现 #

如果习惯使用dplyr等tidyverse包，可以用行操作实现：

library(dplyr)
raw_counts <- raw_counts %>%
  rowwise() %>%
  # 计算每行所有列的方差
  mutate(gene_var = var(c_across(everything()), na.rm = TRUE)) %>%
  ungroup() %>%
  # 过滤掉方差过小的基因
  filter(gene_var > 1e-10) %>%
  # 移除临时生成的方差列
  select(-gene_var)
dim(raw_counts)

内容的提问来源于stack exchange，提问作者sid more