Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/73.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
为多个输入文件设置多个目标文件-sortBam_R_Bioinformatics - Fatal编程技术网

为多个输入文件设置多个目标文件-sortBam

为多个输入文件设置多个目标文件-sortBam,r,bioinformatics,R,Bioinformatics,作为前言,我对使用R是相当陌生的 我正在处理Rsamtools包,我正在尝试对bam文件进行排序。我正在使用sortBam()命令,它需要以下信息: sortBam(file, destination, ...). 我有一个R对象,其中包含5个文件: bamfiles <- list.files(fileDir, pattern = "\\.bam", full.names = TRUE) bamfiles假设在三个不同的目录中有3个bam文件 bams <- list.fil

作为前言,我对使用R是相当陌生的

我正在处理Rsamtools包,我正在尝试对bam文件进行排序。我正在使用sortBam()命令,它需要以下信息:

sortBam(file, destination, ...). 
我有一个R对象,其中包含5个文件:

bamfiles <- list.files(fileDir, pattern = "\\.bam", full.names = TRUE)

bamfiles假设在三个不同的目录中有3个bam文件

bams <- list.files(".", patter="*.bam", full=T, recur=T)
bams
# [1] "./d1/100.bam" "./d2/11.bam"  "./d3/28.bam" 
排序部分:

for (i in 1:length(bams)) {
  sortBam(bams[i], des[i])
}
磁盘上的结果文件

ll -htr *
d1:
total 72K
-rw-rw-r-- 1 35K Apr  8  2016 100.bam
-rw-rw-r-- 1 35K Apr  8  2016 100.sorted.bam

d2:
total 64K
-rw-rw-r-- 1 30K Apr  8  2016 11.bam
-rw-rw-r-- 1 31K Apr  8  2016 11.sorted.bam

d3:
total 80K
-rw-rw-r-- 1 40K Apr  8  2016 28.bam
-rw-rw-r-- 1 40K Apr  8  2016 28.sorted.bam

非常感谢你!这解决了我的问题。我很好奇sortBam是否可以同时在多个文件上使用,所以这也是一个很好的发现。
ll -htr *
d1:
total 72K
-rw-rw-r-- 1 35K Apr  8  2016 100.bam
-rw-rw-r-- 1 35K Apr  8  2016 100.sorted.bam

d2:
total 64K
-rw-rw-r-- 1 30K Apr  8  2016 11.bam
-rw-rw-r-- 1 31K Apr  8  2016 11.sorted.bam

d3:
total 80K
-rw-rw-r-- 1 40K Apr  8  2016 28.bam
-rw-rw-r-- 1 40K Apr  8  2016 28.sorted.bam