复现Cell图表：pyscenic分析之转录因子二项值热图

接上一节视频教程的分析结果:
来源：
https://www.bilibili.com/video/BV1Yj411r7dE/?spm_id_from=333.999.0.0&vd_source=05b5479545ba945a8f5d7b2e7160ea34

今天我们复现一篇cell子刊的图表，这篇文章有一副关于转录因子的图表，观察这个图有什么特点呢？第一是热图是二项值热图，只有0，1两个值，我们知道，在R语言版本的SCENIC分析中，最后可以得到二项值热图，那么pyscenic的分析结果中也是可以进行二项值分析并做热图的。第二是热图左侧注释有auc值的注释，热图是按照分组split的，而且两组的celltype分布是对称的。第三则是左侧行名的展示。

image.png

（Single-cell RNA-seq unravels alterations of the human spermatogonial stem cell compartment in patients with impaired spermatogenesis）

我们的复现结果如下，因为数据是随意挑选的，所以结果看上去没有原图那么明显，但是所有的元素我们都完成了。

image.png

学习这个帖子，我们可以学习到： 1、单细胞seurat转python对象，并进行python数据处理****2、pyscenic分析结果分析得到TF二项值 3、复杂热图的制作以及行列注释，这里我们提供了两种方法****4、文字注释以及与heatmap的结合****5、一些基本数据操作处理的学习接下来我们具体看看做法。首先在pyscenic分析完成后，我们最后一步得到一个loom文件，我们在python中将其转化，并计算转录因子二项值并保存用于后续作图。