手把手教你用R语言下载TCGA数据库:GDCRNAtools

2019-05-23 23:12:09 浏览数 (1)

各位科研芝士的朋友,大家好,今天我们继续分享关于TCGA数据下载的专题,如果你看完了前面所有的关于TCGA推文,那么你对TCGA将不再陌生,这个时候的你不仅学会了九阳神功,还学会了乾坤大挪移了,哈哈,看到这希望你可以继续看下去,毕竟TCGA的生活也应该充满快乐。

或许你会问,之前的网页版工具都有相应的工具包,那么针对官网下载有没有相应的R语言工具包呢?如果,你可以体提出这种疑问,说明你不再是小白了,至少你已经入门了,那今天就给大家推荐,支持官网数据下载的神器GDCRNAtools工具包。

甩出网址链接:

https://bioconductor.org/packages/devel/bioc/vignettes/GDCRNATools/inst/doc/GDCRNATools.html

这个工具包主要是用于之前大火的ceRNA网络构建的工具,但是同时提供了数据下载通道。

GitHub:https://github.com/Jialab-UCR/GDCRNATools

下面开启你的R界面,学习该包:

1. GDCRNATools安装,借助BiocManager安装,前提也是你要安装好BiocManager,命令如下:

2. 加载该包:

Ok,可以看到没有任何问题,这也表明,我们安装并成功加载该工具包

3. 下载数据,这里以胃癌为例子,首先要定义胃癌在TCGA里面的ID信息,为TCGA-STAD,用该包下载就十分简单了,直接用gdcRNADownload命令即可

首先定义项目的id,本例为STAD,接着是数据的类型,我们可以下载RNAseq也可以下载miRNAseq,这里我们选择RNAseq数据下载,在method参数我们可以看到我们选择的是gdc-client,代表的就是借助该工具进行的数据下载,这个时候,如果你还记得我们的第一个官网下载的工具,那就明白了该工具的作用了。

4. 下载个miRNA的数据试试水,由于该数据数据量小,适合演示操作,如下:

结果如下:

你会发现,他是先下载gdc-client工具,接着自动调用该工具下载,中间可能会由于网络等原因报错,不用怕,在运营一次即可,其会在上次下载的基础上,继续接着下载。

如下:

最终下载的结果文件如下:

一些列的文件夹,每个文件夹代表一个样本

5. 解析元文件metadata,进行下一步数据整理

接着还需要下面两个操作进行剔除重复的样本和冗余样本,如下:

6. 文件合并,将所有的表达信息整理成一个矩阵:

这里你需要把path参数换成你的存储路径,可以看到使用gdcRNAMerge我们可以成功的将数据合并成一个表达矩阵,不过该表达矩阵实际上是count。

如下:

获得了这个矩阵,就表示你的数据下载成功了!

0 人点赞