【代码】美化 | 单细胞转录组多组差异基因火山图/环形火山图

文摘   2024-10-28 10:30   江苏  

#

合集说明

您现在看到的是单细胞转录组高级分析的部分内容(持更中)...在阅读全文之前,推荐您先掌握单细胞转录组标准分析流程。

单细胞转录组标准分析

单细胞转录组高级分析

推荐阅读顺序

Cellranger上游分析(植物适用)

去批次-liger

美化-featureplot
美化-dotplot美化-featureplot
美化-降维图

数据库-植物干细胞

细胞分化-基本概念1
细胞分化-基本概念2
细胞分化-monocle1原理
细胞分化-monocle2原理
细胞分化-monocle2实战
细胞分化-velocyto.R

细胞通讯-CellPhoneDBv5简介
细胞通讯-CellPhoneDBv5安装及使用
细胞通讯-可视化-热图
细胞通讯-可视化-网络图

联合bulk-Scissor
联合单细胞atac-Liger

文献复现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计算差异基因后的结果文件,包括上下调的结果。如果只想展示上/下调,可以换成第三种形式。先介绍第一种。

#

以下是正文部分

你好我是一只羊
个人号,内容主要涉及种质资源、分子标记开发及遗传多样性分析,表观遗传、编程语言在生物信息学中的应用、转录组、基因组、单细胞测序多组学数据分析等;其它更新平台:B站&小红书-一只羊做生信/捡羊毛的咩/生信小羊🐑
 最新文章