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

Cannot compute TSSes or TESes without strand info #6

Open
EmilVorsteveld opened this issue May 7, 2024 · 3 comments
Open

Cannot compute TSSes or TESes without strand info #6

EmilVorsteveld opened this issue May 7, 2024 · 3 comments

Comments

@EmilVorsteveld
Copy link

Using GTF files from GENCODE, I get the following error:

cerberus gtf_to_bed --gtf gencode.v45.annotation.gtf -o gencode.v45.annotation.out --mode tss

Exception: Cannot compute TSSes or TESes without strand info. Perhaps use extend() or subsequence() or spliced_subsequence() instead?

@fairliereese
Copy link
Member

Without digging into the issue too much, I have a feeling this might be caused by strange PyRanges stranded vs. unstranded behavior that I've been encountering lately (see this issue). Either way, I recommend installing my version of PyRanges for use with Cerberus. Try this and let me know if it solves your issue!

@EmilVorsteveld
Copy link
Author

This indeed solves this issue, thanks.
Unfortunately, for me, the agg_ics function throws the following error:
File "/Users/emilvorsteveld/miniforge3/envs/cerberus/lib/python3.9/site-packages/pandas/core/internals/base.py", line 57, in _validate_set_axis
raise ValueError(
ValueError: Length mismatch: Expected axis has 4 elements, new values have 3 elements

I have tried older pandas versions (1.3.0 and 1.5.3), but not luck...

@fairliereese
Copy link
Member

Hi, sorry for the delay. Can you please attach the complete error message stack trace so I can see what line in Cerberus threw the error?

Thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants