首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >在R中创建循环

在R中创建循环
EN

Stack Overflow用户
提问于 2019-11-18 07:12:40
回答 1查看 40关注 0票数 0

我想在R中为下面的函数写一个循环:

代码语言:javascript
运行
复制
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“。

代码语言:javascript
运行
复制
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)
}

但它不断给出错误,它不能识别函数。

EN

回答 1

Stack Overflow用户

发布于 2019-11-18 12:59:48

我得到它是为了保存单独的文件:

解决方案:

代码语言:javascript
运行
复制
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"))
}
票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/58906097

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档