按集群导出Seurat对象数据

Question

我正在使用Seurat执行单个单元分析，并且对导出每个集群中所有单元的数据很感兴趣。我尝试使用下面的代码，但没有成功。

我的Seurat对象名为Patients。我还附上了我的Seurat对象的截图。我正在寻找提取所有集群(即Ductal1、Macrophage1、Macrophage2等)。

meta.data.cluster <- unique(x = Patients@meta.data$active.ident)

for(group in meta.data.cluster) {
  group.cells <- WhichCells(object = Patients, subset.name = "active.ident" , accept.value = group)
  data_to_write_out <- as.data.frame(x = as.matrix(x = Patients@raw.data[, group.cells]))
  write.csv(x = data_to_write_out, row.names = TRUE, file = paste0(save_dir,"/",group, "_cluster_outfile.csv"))
}

​

​

我是R和编码的新手，所以非常感谢您的帮助！:)

Answer 1

它不起作用，因为您的元数据下没有active.ident列。例如，如果我们使用与您的数据集类似的示例数据集并设置ident：

library(Seurat)

M = matrix(rnbinom(5000,mu=20,size=1),ncol=50)
colnames(M) = paste0("P",1:50)
rownames(M) = paste0("gene",1:100)

Patients = CreateSeuratObject(M)
Patients$grp = sample(c("Ductal1","Macrophage1","Macrophage2"),50,replace=TRUE)
Idents(Patients) = Patients$grp

您可以看到这行代码没有给您提供任何值：

meta.data.cluster <- unique(x = Patients@meta.data$active.ident)
meta.data.cluster
NULL

您可以执行以下操作：

meta.data.cluster <- unique(Idents(Patients))

for(group in meta.data.cluster) {
  group.cells <- WhichCells(object = Patients, idents = group)
  data_to_write_out <- as.data.frame(GetAssayData(Patients,slot = 'counts')[,group.cells])
  write.csv(data_to_write_out, row.names = TRUE, file = paste0(save_dir,"/",group, "_cluster_outfile.csv"))
}

另请注意，您可以使用GetAssayData获得计数。你可以设置一个组的子集，然后这样写出来：

wh <- which(Idents(Patients) =="Macrophage1" )
da = as.data.frame(GetAssayData(Patients,slot = 'counts')[,wh])
write.csv(da,...)