The new version of csaw
introduced the very helpful dumpPE
function to improve efficiency of downstream methods.
However, is one expected to use this method one file at a time or can multiple files be supplied.
For example, csaw
suggests this:
outBam <- dumpPE(bamFile, "whee", param=xparam, overwrite=TRUE)
out2 <- windowCounts(outBam, param=reform(xparam, fast.pe=TRUE), filter=1)
However, I want to do something like this
outBams <- dumpPE(bamFiles, prefix=list_of_prefixes, param=xparam, overwrite=TRUE) out2 <- windowCounts(outBams, param=reform(xparam, fast.pe=TRUE), filter=1)
Is this a possibility ?
Currently the specs say that bam.file should be of type character(1).
Any help would be appreciated.
Thanks !
I should point out that
dumpPE
is going to be retired in the next release. I made the defaultpe="both"
option faster on the original BAM files, which avoids the need to construct a temporary file. It also makes the downstream processing more consistent, which was difficult to do with the dumped alignments.Thank you, Aaron. This helps!