R 遍历目录下所有子目录中的某类文件
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 可以输出文本,还可以存储文本/向量
本文来自博客园,作者:BioinformaticsMaster,转载请注明原文链接:https://www.cnblogs.com/koujiaodahan/p/15908367.html
posted on 2022-02-18 14:00 BioinformaticsMaster 阅读(123) 评论(0) 收藏 举报
浙公网安备 33010602011771号