#====================================================鑒定不同定義細(xì)胞類群分子標(biāo)記
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)
}
輸出細(xì)胞類群的marker并繪圖
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進(jìn)店門于颖，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人厉萝，你說我怎么就攤上這事。” “怎么了谴垫？”我有些...
- 文/不壞的土叔 我叫張陵章母，是天一觀的道長。 經(jīng)常有香客問我翩剪，道長乳怎，這世上最難降的妖魔是什么？ 我笑而不...
- 正文 為了忘掉前任前弯，我火速辦了婚禮蚪缀，結(jié)果婚禮上，老公的妹妹穿的比我還像新娘恕出。我一直安慰自己询枚，他們只是感情好，可當(dāng)我...
- 文/花漫 我一把揭開白布浙巫。 她就那樣靜靜地躺著金蜀，像睡著了一般。 火紅的嫁衣襯著肌膚如雪的畴。 梳的紋絲不亂的頭發(fā)上渊抄，一...
- 文/蒼蘭香墨 我猛地睜開眼逊桦，長吁一口氣：“原來是場噩夢啊……” “哼眨猎！你這毒婦竟也來了？” 一聲冷哼從身側(cè)響起强经，我...
- 正文 年R本政府宣布橱野，位于F島的核電站，受9級特大地震影響善玫，放射性物質(zhì)發(fā)生泄漏水援。R本人自食惡果不足惜，卻給世界環(huán)境...
- 文/蒙蒙 一茅郎、第九天 我趴在偏房一處隱蔽的房頂上張望蜗元。 院中可真熱鬧，春花似錦系冗、人聲如沸奕扣。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽成畦。三九已至，卻和暖如春涝开，著一層夾襖步出監(jiān)牢的瞬間循帐，已是汗流浹背。 一陣腳步聲響...
推薦閱讀更多精彩內(nèi)容
- 前面給大家介紹過神兵利器——單細(xì)胞細(xì)胞類群基因marker鑒定新方法：COSG[https://www.jians...
- 在單細(xì)胞數(shù)據(jù)分析當(dāng)中棋弥，當(dāng)我們聚類分群完成之后，緊接著就是細(xì)胞類型注釋诚欠，細(xì)胞類型的注釋就離不開 基因marker 顽染，...
- 忙里偷閑，整理了一套繪圖代碼轰绵，代碼300行粉寞，效果如上圖所示：其實就是把單細(xì)胞轉(zhuǎn)錄組中常見的氣泡圖根據(jù)分組信息展開了...
- 由于基因表達(dá)跟細(xì)胞表面蛋白的表達(dá)具有時間差，所以要根據(jù)經(jīng)驗具體情況具體分析左腔，內(nèi)容僅供參考 ??哺乳類動物B細(xì)胞的分...
- T 細(xì)胞（全部）：CD3+ Naive T cell: CD3, CD4, CCR7, CD62L+, IL-7R...
