我想在R中为下面的函数写一个循环:
Idents(immune) = "seurat_clusters"
immune.0 = subset(immune, idents = ("0"))
Idents(immune.0) = "stim"
immune.0_markers = FindAllMarkers(object=immune.0, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.0_markers,"immune.0markers.csv")
save(file='immune.0.Rdata',immune.0)
Idents(immune) = "seurat_clusters"
immune.1 = subset(immune, idents = ("1"))
Idents(immune.1) = "stim"
immune.1_markers = FindAllMarkers(object=immune.1, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.1_markers,"immune.1markers.csv")
save(file='immune.1.Rdata',immune.1)
Idents(immune) = "seurat_clusters"
immune.2 = subset(immune, idents = ("2"))
Idents(immune.2) = "stim"
immune.2_markers = FindAllMarkers(object=immune.2, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.2_markers,"immune.2markers.csv")
save(file='immune.2.Rdata',immune.2)
Idents(immune) = "seurat_clusters"
immune.3 = subset(immune, idents = ("3"))
Idents(immune.3) = "stim"
immune.3_markers = FindAllMarkers(object=immune.3, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.3_markers,"immune.3markers.csv")
save(file='immune.3.Rdata',immune.3)
我尝试了以下方法,集群被称为从0到20的"seurat_clusters“。
for (i in 0:20){
Idents(immune) = "seurat_clusters"
immune.i = subset(immune, idents = ("i"))
Idents(immune.i) = "stim"
immune.i_markers = FindAllMarkers(object=immune.i, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.i_markers,"immune.imarkers.csv")
save(file='immune.i.Rdata',immune.i)
}
但它不断给出错误,它不能识别函数。
发布于 2019-11-18 12:59:48
我得到它是为了保存单独的文件:
解决方案:
for (i in 0:20){
Idents(immune) = "seurat_clusters"
immune.i = subset(immune, idents = as.character(i))
Idents(immune.i) = "stim"
immune.i_markers = FindAllMarkers(object=immune.i, only.pos=TRUE, logfc.threshold=0.25)
write.csv(immune.i_markers, paste0(i,".csv"))
}
https://stackoverflow.com/questions/58906097
复制相似问题