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

armatus aborted #18

Open
aminakur opened this issue Dec 10, 2022 · 0 comments
Open

armatus aborted #18

aminakur opened this issue Dec 10, 2022 · 0 comments

Comments

@aminakur
Copy link

aminakur commented Dec 10, 2022

I believe that my issue is similar to #12
I have a sparse matrix file of the entire genome in hicpro format.
The resolution of my matrix is 1kb. Genome size is 379.6Mb. So I have a sparse matrix file with 379633 bins.

$ head az.hicpro
1	1	0.298138568943344
1	2	0.11464277349213885
1	3	0.012447999473593192
1	4	0.008286857583595878
$tail az.hicpro
379631	379632	0.049243019444128566
379631	379633	0.006873960422137417
379632	379632	0.1314174872134235
379632	379633	0.0387687115012471
379633	379633	0.08176332692011251

When I run armatus with this command:
armatus -S -g 0.9 -r 1000 -s 0.1 -i az.hicpro -o az_armatus
I get an output with 72 consensus domains, which is too low (I found around 4000 TADs with HiCExplorer). The log file:

Reading input from az.hicpro.
Building matrix for chromosome N/A at resolution 1000bp with 380 rows.
Initializing matrix to zero elements
0.0459572%
0.0919145%

I read the comment to #12 and thought that maybe armatus thought my bin numbers were chromosomal coordinates and divided them by resolution, hence low number of domains. I tried the command without -r:
armatus -S -g 0.9 -s 0.1 -i az.hicpro -o az_armatus1
I got an error:

Reading input from az_mrg.hicpro.
Building matrix for chromosome N/A at resolution 1bp with 379634 rows.
Initializing matrix to zero elements
terminate called after throwing an instance of 'std::bad_alloc'
  what():  std::bad_alloc
/opt/slurm/data/slurmd/job27956122/slurm_script: line 14: 2093009 Aborted                 (core dumped) armatus -S -i az.hicpro -g 0.9 -o az_armatus

How can I resolve this?

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

1 participant