SCENIC-Regulon的热图怎么实现？

愿武艺晴小朋友一定得每天都开心！ 

---

当按流程，算出SCENIC的结果后；

基于scRNA-data,运用pySCENIC寻找细胞群里面活跃的调节子_pyscenic 流程-优快云博客

像可视化出一张具有生物学意义的热图，该怎么做呢？

其中，用什么数据来作为热图的输入呢？ 有两个选择1是用RSS值；二是用AUC值

这里我们选AUC函数，并将AUC值按组间归一化后进行热图展示。

第一步：得到AUC值矩阵

AUC = getAUC(sub_regulonAUC)

第二步：归一化矩阵

# Split the cells by cluster:
cellsPerGroup <- split(rownames(metadata),metadata[,selectedResolution]) 

# 按分组 Calculate average expression:
regulonActivity_byGroup <- sapply(cellsPerGroup,
                                          function(cells) 
                                                      rowMeans(getAUC(sub_regulonAUC)[,cells]))
# Scale expression:
regulonActivity