search

VCCRI / Sierra

Discover differential transcript usage from polyA-captured single cell RNA-seq data
GNU General Public License v3.0
49 stars 17 forks source link

Error in SplitBam (extra argument to Rsamtools::ScanBamParam) #28

Closed borshago closed 3 years ago

borshago commented 3 years ago

Hi,

I'm following the vignette in the wiki and running the code in chunks. The following chunk raised an error:

outdir = "bam_subsets/"
dir.create(outdir)
SplitBam(bamfile[1], cells.df, outdir)

and the error reads:

Error in Rsamtools::ScanBamParam(tag = bamTags, what = what, tag = bamTags) : formal argument "tag" matched by multiple actual arguments

Looking at the code in split_bams.R, all the function calls to Rsamtools::ScanBamParam contain the tag argument twice. Shouldn't the argument be passed just once?

Best, Daniel

davhum commented 3 years ago

Hi Daniel,

You are correct. Thanks for bringing to our attention. Issue should be resolved now. Let us know if you have further issues. Regards,