Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Alignments error while using PairedGuideSet object as input #12

Open
abearab opened this issue Feb 15, 2023 · 2 comments
Open

Alignments error while using PairedGuideSet object as input #12

abearab opened this issue Feb 15, 2023 · 2 comments
Assignees
Labels
enhancement New feature or request medium-priority

Comments

@abearab
Copy link

abearab commented Feb 15, 2023

Is it possible to run addSpacerAlignmentsIterative on PairedGuideSet object?! I an error running this code:

pairsGuideSet <- addSpacerAlignmentsIterative(
    pairsGuideSet,
    txObject=txdb,
    aligner_index=bowtie_index,
    bsgenome=bsgenome,
    n_mismatches=2
)
R[write to console]: Error in METHOD(x, i) : 
  Subsetting operation on CompressedGRangesList object 'x' produces a
  result that is too big to be represented as a CompressedList object.
  Please try to coerce 'x' to a SimpleList object first (with 'as(x,
  "SimpleList")').


Error in METHOD(x, i) : 
  Subsetting operation on CompressedGRangesList object 'x' produces a
  result that is too big to be represented as a CompressedList object.
  Please try to coerce 'x' to a SimpleList object first (with 'as(x,
  "SimpleList")').

This is the size of my pairsGuideSet:

PairedGuideSet object with 110196 pairs and 4 metadata columns:
@abearab abearab changed the title addSpacerAlignmentsIterative addSpacerAlignmentsIterative error while using PairedGuideSet object as input Feb 15, 2023
@abearab abearab changed the title addSpacerAlignmentsIterative error while using PairedGuideSet object as input Alignments error while using PairedGuideSet object as input Feb 15, 2023
@Jfortin1
Copy link
Member

Hi @abearab,

Yes, addSpacerAlignmentsIterative should work on PairedGuideSet objects, but it seems the data are too big to be handled for the underlying CompressedGRangesList class from the GRanges package; I'll work on a fix, but in the meantime, I'd suggest to do some filtering on your individual gRNAs before pairing them to reduce the size of the object, or alternatively, split your PairedGuideSet object into a list of PairedGuideSet objects and apply the functional sequentially.

@Jfortin1 Jfortin1 self-assigned this Feb 17, 2023
@Jfortin1 Jfortin1 added enhancement New feature or request medium-priority labels Feb 17, 2023
@abearab
Copy link
Author

abearab commented Apr 19, 2023

Got it, thank you for your response. I'll get back to you if I ran to another related issue.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request medium-priority
Projects
None yet
Development

No branches or pull requests

2 participants