打开APP
userphoto
未登录

开通VIP,畅享免费电子书等14项超值服

开通VIP
批量将多个文件内的压缩文件复制至一个文件内

有感于 果子学生信

应用:从TCGA GDC中下载下来的文件进行汇总

library(TCGAbiolinks)

query <- GDCquery(project = "TCGA-LUSC",

                    data.category = "Transcriptome Profiling",

                    data.type = "Gene Expression Quantification",

                    workflow.type = "HTSeq - Counts" )

  GDCdownload(query)

  

# 所有下载的counts数据存在“./GDCdata/TCGA-LUSC/harmonized/Transcriptome_Profiling”;

# 并且被保存在一个个文件夹内

# 接下来首先将每个文件夹内的文件转自一个文件内"./expr_data/raw_data_in_one"

dirname <- "001a2541-a9d3-4811-97bf-d8fcd89b4c66"

for (dirname in dir("./GDCdata/TCGA-LUSC/harmonized/Transcriptome_Profiling/Gene_Expression_Quantification/")) {

  ## use list.files function to find single file within “./GDCdata/TCGA-LUSC/harmonized/Transcriptome_Profiling”

  file <- list.files(path = paste0(getwd(),

                                   "/GDCdata/TCGA-LUSC/harmonized/Transcriptome_Profiling/Gene_Expression_Quantification/",

                                   dirname),pattern = "*.counts.gz") #找到文件中的内容,pattern可以是正则表达式

  ## 使用file.copy函数复制黏贴压缩文件到raw_data_in_one

  file.copy(paste0(getwd(),"/GDCdata/TCGA-LUSC/harmonized/Transcriptome_Profiling/Gene_Expression_Quantification/",

                   dirname,"/",file),to = "./expr_data/raw_data_in_one") #复制内容到新的文件夹

}

本站仅提供存储服务,所有内容均由用户发布,如发现有害或侵权内容,请点击举报
打开APP,阅读全文并永久保存 查看更多类似文章
猜你喜欢
类似文章
1行代码提取6种TCGA表达矩阵2.0版
超友好的TCGA数据下载方式!
TCGAbiolinks下载TCGA数据总结
TCGA数据库悄咪咪更新了—RNAseq没有HTSeq-Counts了
把GDC下载的多个TCGA文件批量读入R
批量读入TCGA的miRNA数据
更多类似文章 >>
生活服务
热点新闻
分享 收藏 导长图 关注 下载文章
绑定账号成功
后续可登录账号畅享VIP特权!
如果VIP功能使用有故障,
可点击这里联系客服!

联系客服