单细胞转录组数据处理之细胞亚群比例比较
前面我们一起学习到了单细胞转录组数据的降维聚类分群，而且拿到的亚群，也可以解释其生物学意义，见：单细胞转录组数据处理之细胞亚群注释，接下来可以把拿到的亚群进行更细致的分群，或者看不同样本不同组别的内部的细胞亚群的比例变化。
2个10X样本可以看细胞亚群比例
可以看到如下是2个单细胞转录组样本，它们的批次效应去除的挺好的，tSNE图如下：
作者就根据聚类分群结果，对细胞亚群命名后，重新拆分开来看两个样本的不同细胞亚群的比例，如下：
d, Bar chart of the relative frequency of immune cell types derived from aggregated MC data in blood and plaque.
很明显，有细胞亚群的比例变化：
有趣的是，还可以使用火山图的形式来展现，这两个分组的细胞亚群比例变化情况。不过，需要
e, Volcano plot of the fold change of MC frequency in plaque and blood from 15 patients
另外一个例子也是2个样本
同样的是2个样本，所以理论上分群后，注释了，就可以同样的看不同亚群在这两个样本的比例差异情况。
只需要简单的一个table函数，就可以计算出不同样本的不同细胞亚群的细胞数量，进而计算比例后绘制图表。
如果是4个样本
这个情况下，比例条形图仍然是可以的，但是火山图就没有可能啦，火山图只能是两个对比。
可以看到，上图的4个样本的去除批次效应也挺好的，有趣的是分群得到的13个亚群，在4个样本里面的比例差异不大。
文末友情宣传
全国巡讲全球听（买一得五）第3期（4月6日开始） ，你的生物信息学入门课。 数据挖掘线上班来袭（两天变三周，实力加量），医学生/医生首选技能提高课。 生信技能树的2019年终总结 ，你的生物信息学成长宝藏 2020学习主旋律，B站74小时免费教学视频为你领路
赞 (0)