tcga压缩包提取合并_TCGA数据库下载后文件合并与转换

本文详细介绍了如何从TCGA数据库下载并处理癌症相关数据,包括使用gdc-client下载原始Seq-count数据,通过Perl脚本整理文件,利用R脚本合并数据成seq-count矩阵,以及如何将基因ID转换为符号名,以备后续差异表达分析。
摘要由CSDN通过智能技术生成

1、从TCGA上下载癌症相关数据到rawdata,下载一般用TCGA_gdc-client。以原始的Seq-count为例。

gdc-client download -m gdc_manifest_20200923_032916.txt -d path/to/gz_file

2、加入cart之后下载manifest.json文件。(PS:这个文件真的无比讨厌,处理起来很费经)

3、将数据解压缩之后生成看起来像乱码的文件夹。此时利用perl脚本将这每个文件夹中的文件放入另外一个file文件夹当中,然后将文件名改成data_in_one便于识别。

###将所有压缩文件移到一个名为files的文件里面

###用法为cmd perl 1.move_gzFileToFiles.pl

use strict;

use warnings;

use File::Copy;

my $newDir="files";

unless(-d $newDir)

{

mkdir $newDir or die $!;

}

my @allFiles=glob("*");

foreach my $subDir(@allFiles)

{

if((-d $subDir) && ($subDir ne $newDir))

{

opendir(SUB,"./$subDir") or die $!;

while(my $file=readdir(SUB))

{

if($file=~/\.gz$/)

{

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值