君子生非异也,善假于物也
——《荀子·劝学》
上游分析请看生信技能树的历史教程 假设你已经全面学习了肿瘤外显子测序数据的分析,得到了你所有配对样本的somatic mutation信息,而且制作成了还算标准的maf格式文件,那么就可以跟着这篇教程进行可视化啦。
全部流程在R里面运行,maftools提供了直接读取maf文件的接口,而且存储为S4对象,非常方便进行一系列可视化操作。
而且大部分都是一行代码出图,基本上不需要太高深的R语言知识,现学现用。
而且得到zuguan大神的complexHeatmap的oncoplot也被封装成了一行代码:
oncoplot(maf = laml, top = 10, removeNonMutated = TRUE)
即可出图如下:
还有以前我在生信菜鸟团博客介绍过的,MutationMapper on cBioPortal 的lollipopPlot
最近比较火的TMB也囊括进来了:
下载TCGA所有癌症的maf文件计算TMB
华而不实的词云也少不了咯:
还有探索MAF分布:
其实图并不重要,重要的是理解图背后的生物学医院。
生信技能树GATK4系列教程 GATK4的gvcf流程 你以为的可能不是你以为的 新鲜出炉的GATK4培训教材全套PPT,赶快下载学习吧 曾老湿最新私已:GATK4实战教程 GATK4的CNV流程-hg38