r - 基于 Seurat 对象的读取计数的子集数据（FetchData 中的错误）

Question

我正在尝试根据感兴趣的基因计数对 Seurat 对象（称为 dNSC_cells）进行子集化。具体来说，我有一个基因列表，我计划遍历它们以对我的数据进行子集化并进行一些 Wilcox 测试。我到目前为止看起来像这样：

pro_genes_list <- c("Bcl2", "Bid")

for (p in pro_genes_list) { 
  
  median_prog <- median(GetAssayData(dNSC_cells, slot = 'counts')[p,])
  with_p <- colnames(subset(dNSC_cells, subset = as.name(p) > median_prog))

但是，它在最后一行阻塞，出现以下错误：

Error in FetchData(object = object, vars = unique(x = expr.char[vars.use]),  : 
  None of the requested variables were found: 

我也尝试过使用

subset = p > median_prog

但它给出了同样的错误。将非常感谢任何指针！

Answer 1

该subset函数不支持变量中的基因符号，您需要先提取一个数据框。

for (p in pro_genes_list) { 
  median_prog <- median(GetAssayData(dNSC_cells, slot = 'counts')[p,])
  expr <- FetchData(object = dNSC_cells, vars = p)
  with_p <- colnames(dNSC_cells[, which(expr > median_prog)])
}