#====================================================鉴定不同定义细胞类群分子标记
MyMarker_celltype <- function(seurat_object,path_out){
marker_list <- list()
for (cell_type in unique(seurat_object@meta.data$cell_type)){
cell_1 <- rownames(seurat_object@meta.data)[seurat_object@meta.data$cell_type == cell_type]
marker_list[[cell_type]] <- FindMarkers(seurat_object,ident.1=cell_1,only.pos=TRUE,assay="RNA")
}
write.xlsx(marker_list, file=paste(path_out,'all_cell_type_markers.xlsx',sep='/'), row.names=TRUE,overwrite=T)
saveRDS(marker_list,file=paste(path_out,"all_cell_type_markers.rds",sep="/"))
cell_type_marker_top <- unique(unlist(lapply(marker_list,function(x) {top <- rownames(x[order(x$avg_log2FC,decreasing=T),])[1:10]}))) #每个top10作为marker
plot_markers <- DotPlot(seurat_object , features = cell_type_marker_top,assay='RNA',group.by="cell_type" ) + coord_flip()
ggsave("dot_plot_cell_type_markers.png",plot_markers,path=path_out,width=10,height=15)
return(marker_list)
}
输出细胞类群的marker并绘图
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 前面给大家介绍过神兵利器——单细胞细胞类群基因marker鉴定新方法：COSG[https://www.jians...
- 在单细胞数据分析当中，当我们聚类分群完成之后，紧接着就是细胞类型注释，细胞类型的注释就离不开 基因marker ，...
- 忙里偷闲，整理了一套绘图代码，代码300行，效果如上图所示：其实就是把单细胞转录组中常见的气泡图根据分组信息展开了...
- 由于基因表达跟细胞表面蛋白的表达具有时间差，所以要根据经验具体情况具体分析，内容仅供参考 哺乳类动物B细胞的分...
- T 细胞（全部）：CD3+ Naive T cell: CD3, CD4, CCR7, CD62L+, IL-7R...
