用html5画瀑布图,maftools | 从头开始绘制发表级oncoplot(瀑布图)

本文详细介绍了如何使用R语言的maftools包来绘制oncoplot(瀑布图),包括载入数据、MAF对象汇总、绘制基础oncoplots以及添加不同信息,如SCNA、significance values和临床信息注释。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

对于组学数据的分析和展示来说,maftools算是一个宝藏“R包”,可用于MAF格式的组学数据的汇总,分析和可视化展示。随着癌症基因组学的进步, [突变注释格式](https://wiki.nci.nih.gov/display/TCGA/Mutation+Annotation+Format+(MAF)+Specification) (MAF) 被广泛用于存储检测到的somatic variants。[The Cancer Genome Atlas](http://cancergenome.nih.gov/) 项目对30多种不同的癌症进行了测序,每种癌症类型的样本量超过200种。maftools-R包能够有效的汇总,分析和可视化MAF格式的文件。

maftools函数主要分为可视化和分析两个模块,其主要功能及简短的描述如下所示,使用时只需读取MAF文件然后降MAF对象传递给所需要的绘图或分析功能即可。

170069543_1_2019090310395778.jpg

本次主要使用R-maftools包绘制组学突变结果(MAF)的oncoplot或者叫“瀑布图”。

一、 载入R包,数据

1) 载入maftools包if (!require("BiocManager"))

install.packages("BiocManager")

BiocManager::install("maftools")

2) 载入数据

通过read.maf函数读入MAF文件,将各种数据(组学基因突变,拷贝数变异,临床数据,表达数据等)汇总并将其存储为MAF对象。#TCGA-LAML MAF file (gz)laml.maf = system.file('extdata', 'tcga_laml.maf.gz', package = 'maftools')

#clinical information (optional)laml.clin = system.file('extdata', 'tcga_laml_annot.tsv', package = 'maftools')

laml = read.maf(maf = laml.maf, clinicalData = laml.clin)

MAF对象中除了上面基因突变数据和对应的临床数据外,还可以加入拷贝数变异,表达数据等其他数据类型,后面需要的时候会添加。

二、 MAF对象汇总

1) 展示MAF重点变量的summary信息#Shows sample summry.getSampleSummary(laml)#Shows gene summary.getGeneSummary(laml)#Shows all fields in MAFgetFields(laml)#shows clinical data associated with samplesgetClinicalData(laml)#Writes maf su

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值