Displaying 1 result from an estimated 1 matches for "scanbam".
Did you mean:
scanam
2012 Dec 17
1
Code works standalone, yet same code fails when part of package
...pos", "qwidth", "mapq", "cigar", "mrnm", "mpos",
"isize", "seq"), tag=c("NM"), removeN=TRUE) {
which <- RangesList(chr = IRanges(start,end))
names(which) <- chr
param <- ScanBamParam(which = which, what = what, tag=tag)
bam <- scanBam(bamfile, param=param)
lst <- lapply(names(bam[[1]]), function(elt) {do.call(c, unname(lapply
(bam, "[[", elt)))})
names(lst) <- names(bam[[1]])
# convert to data.frame
df <- do...