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

Bug: OverflowError when running hmmratac #649

Open
atj31 opened this issue Jun 19, 2024 · 0 comments
Open

Bug: OverflowError when running hmmratac #649

atj31 opened this issue Jun 19, 2024 · 0 comments

Comments

@atj31
Copy link

atj31 commented Jun 19, 2024

I'm trying to run hmmratac to call peaks from my bampe files the code I used is
hmmratac -i Alignments/file.bam -f BAMPE --outdir Peaks/ -n name

It runs for a long time but I get an 'OverflowError: signed integer is greater than maximum' error

I tried to troubleshoot but unsure how to proceed - I can run the callpeak on the same bampe file without any issues.

Thank you in advance!

The whole running code pasted below:
[...]
INFO @ 19 Jun 2024 15:58:13: [428 MB] 37000000 fragments parsed
INFO @ 19 Jun 2024 15:58:16: [434 MB] 38000000 fragments parsed
INFO @ 19 Jun 2024 15:58:19: [441 MB] 38850552 fragments have been read.
INFO @ 19 Jun 2024 15:58:50: [441 MB] #2 Use EM algorithm to estimate means and stddevs of fragment lengths
INFO @ 19 Jun 2024 15:58:50: [441 MB] # for mono-, di-, and tri-nucleosomal signals...
INFO @ 19 Jun 2024 15:58:50: [441 MB] # A random seed 10151 has been used in the sampling function
INFO @ 19 Jun 2024 15:58:54: [491 MB] # Downsampled 1336693 fragments will be used for EM training...
INFO @ 19 Jun 2024 16:02:48: [491 MB] # Reached convergence after 11 iterations
INFO @ 19 Jun 2024 16:02:48: [491 MB] # The means and stddevs after EM:
INFO @ 19 Jun 2024 16:02:48: [491 MB] # short mono di tri
INFO @ 19 Jun 2024 16:02:48: [491 MB] # means: 50 173.9 381.2 600
INFO @ 19 Jun 2024 16:02:48: [491 MB] # stddevs: 20 53.6 51.4 20
INFO @ 19 Jun 2024 16:02:48: [491 MB] # Compute the weights for each fragment length for each of the four signal types
INFO @ 19 Jun 2024 16:02:55: [491 MB] # Generate short, mono-, di-, and tri-nucleosomal signals
Traceback (most recent call last):
File "/mnt/data/Utils/python_environ/macs3/bin/macs3", line 1028, in
main()
File "/mnt/data/Utils/python_environ/macs3/bin/macs3", line 101, in main
run( args )
File "/mnt/data/Utils/python_environ/macs3/lib/python3.9/site-packages/MACS3/Commands/hmmratac_cmd.py", line 155, in run
digested_atac_signals = generate_digested_signals( petrack, weight_mapping )
File "MACS3/Signal/HMMR_Signal_Processing.pyx", line 106, in MACS3.Signal.HMMR_Signal_Processing.generate_digested_signals
File "MACS3/Signal/HMMR_Signal_Processing.pyx", line 129, in MACS3.Signal.HMMR_Signal_Processing.generate_digested_signals
File "MACS3/Signal/BedGraph.pyx", line 219, in MACS3.Signal.BedGraph.bedGraphTrackI.add_chrom_data_hmmr_PV
File "MACS3/Signal/BedGraph.pyx", line 232, in MACS3.Signal.BedGraph.bedGraphTrackI.add_chrom_data_hmmr_PV
OverflowError: signed integer is greater than maximum

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

No branches or pull requests

1 participant