Skip to content

Bootstrapping Issue with Resolution cr-like-em in USA ModeΒ #157

@caitlinottaway

Description

@caitlinottaway

I am currently attempting to run alevin-fry to quantify some snRNA-seq data in USA mode with bootstraps. From my understanding, USA-mode should be allowed to be done with bootstrapping, as the bootstrap flag is meaningful with the cr-like-em quantification method, and that method uses USA mode (https://alevin-fry.readthedocs.io/en/latest/quant.html).

However, when I am trying to run this below script, i get the following error. I would appreciate any guidance and clarification!

  alevin-fry quant \
  -t 6 \
  -i $a/alevin_fry/alevinfry_unfiltered_pl_gencode_quants_noforcecells/${c}quant \
  -o $a/alevin_fry/alevinfry_unfiltered_pl_gencode_quants_res_noforcecells/${c}quant \
  --tg-map $a/alevin_fry/salmon_index_files/transcriptome_splici_fl_gencode85/transcriptome_splici_fl85_t2g_3col.tsv \
  -b 30 \
  --resolution cr-like-em \
  --use-mtx

ERROR: assertion left == right failed: currently USA-mode (all-in-one unspliced/spliced/ambiguous) analysis cannot be used with bootstrapping.
left: 30
right: 0

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions