R 遍历目录下所有子目录中的某类文件

1
2
3
4
5
6
7
8
9
10
11
testvcfdir="/jdfssz1/ST_HEALTH/P18Z10200N0124/qdNB/"
testvcfList <- list.files(testvcfdir, pattern = "*.genotype.vcf.gz$",
                             include.dirs = FALSE, full.names = TRUE, recursive = TRUE)
test_vcf_sh <- c()
content<-""
for (vcfFile in testvcfList) {
  sampleId <- unlist(str_split(basename(vcfFile),pattern = "\\*"))[[1]]
  content<-append(test_vcf_sh,paste("sh /jdfssz1/ST_HEALTH/P21Z10200N0047/lizhichao/zbolt_test/process/script/concordance.sh",vcfFile,sampleId,sep=" "))
 
}
cat(test_vcf_sh,file="./run_conincidence,sh",sep="\n")

basename() 保留文件名  

cat 可以输出文本,还可以存储文本/向量

 

posted on   BioinformaticsMaster  阅读(115)  评论(0编辑  收藏  举报

相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5
点击右上角即可分享
微信分享提示