您现在看到的是单细胞转录组高级分析的部分内容(持更中)...在阅读全文之前,推荐您先掌握单细胞转录组标准分析流程。
单细胞转录组标准分析
单细胞转录组高级分析
推荐阅读顺序
美化-featureplot
美化-dotplot美化-featureplot
美化-降维图
细胞分化-基本概念1
细胞分化-基本概念2
细胞分化-monocle1原理
细胞分化-monocle2原理
细胞分化-monocle2实战
细胞分化-velocyto.R
细胞通讯-CellPhoneDBv5简介
细胞通讯-CellPhoneDBv5安装及使用
细胞通讯-可视化-热图
细胞通讯-可视化-网络图
文献复现1-公共数据读取
文献复现2-Dotplot组图
文献复现3-GSE192742人肝细胞图谱1
文献复现4-GSE192742人肝细胞图谱2
合集完整内容查看方式:
1. 合集付费:在公众号首页找到“单细胞转录组(高级分析)”合集标签,点击进入;或在CCtalk平台搜索《单细胞转录组分析流程-高级分析》,即可找到本合集全部已更新内容。
2. 单篇付费:在公众号首页找到“单细胞测序”合集标签,点击进入,即可找到单篇阅读入口。
本篇文章将介绍适用于单细胞的多组差异基因火山图代码,结果如下:
安装
如果要绘制环形,为了让clsuter及标注的基因也“变弯”,需要安装geomtextpath包,当前安装这个包需要将textshaping>= 0.4.0,以及ggplot2升到最新,解决完这两个包再安装geomtextpath。不过问题是,ggplot2升到最新后将会影响Seurat的一些绘图功能,谨慎操作吧。
> DimPlot(PRO,group.by = "seurat_clusters",reduction = "umap")
Error in Ops.data.frame(guide_loc, panel_loc) :
‘==’ only defined for equally-sized data frames
textshaping安装问题
如果升级textshaping遇到报错
Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) :
载入了名字空间‘textshaping’ 0.3.6,但需要的是>= 0.4.0
报错原因
<stdin>:1:10: fatal error: 'hb-ft.h' file not found
解决方法:
https://blog.csdn.net/weixin_30783611/article/details/140677544
以Mac为例:
brew install harfbuzz fribidi
重新安装,ggplot2需要升到最新
remotes::install_github("AllanCameron/geomtextpath")
多组火山图
diffData为seurat计算差异基因后的结果文件,包括上下调的结果。如果只想展示上/下调,可以换成第三种形式。先介绍第一种。