Skip to content

Got zero division error when doing insulation score analysis #207

@arisu1017

Description

@arisu1017

Hi, FanC is an amazing and convenient tools!

However I encountered an issue when I tried FanC insulation on my juicer .hic files. No matter what resolution I changes, I always got this error message stating that there were zero division errors.

Here are my error message:

2025-05-17 23:19:33,789 INFO FAN-C version: 0.9.28
/home/mm/miniconda3/lib/python3.9/site-packages/fanc/compatibility/juicer.py:353: UserWarning: Support for Juicer .hic v9 is still in beta. Please report any issues to #92
warnings.warn(f"Support for Juicer .hic v{self.version} is still in beta. "
2025-05-17 23:19:34,223 INFO Chosen window sizes: 1000000
2025-05-17 23:19:36,784 INFO 1 (1/25)
Traceback (most recent call last):
File "/home/mm/miniconda3/bin/fanc", line 127, in
Fanc()
File "/home/mm/miniconda3/bin/fanc", line 93, in init
command([sys.argv[0]] + sys.argv[option_ix:], log_level=log_level, verbosity=verbosity)
File "/home/mm/miniconda3/lib/python3.9/site-packages/fanc/commands/fanc_commands.py", line 3801, in insulation
_domain_scores(parser, input_file, output_file, output_format,
File "/home/mm/miniconda3/lib/python3.9/site-packages/fanc/commands/fanc_commands.py", line 3853, in _domain_scores
scores = score_class.from_hic(matrix, file_name=output_file,
File "/home/mm/miniconda3/lib/python3.9/site-packages/fanc/architecture/domains.py", line 375, in from_hic
values_by_chromosome[w_ix][k] /= (bin_window_size ** 2 - na_total)
ZeroDivisionError: division by zero
/home/mm/miniconda3/lib/python3.9/site-packages/tables/file.py:113: UnclosedFileWarning: Closing remaining open file: fanc_example_100kb.insulation
warnings.warn(UnclosedFileWarning(msg))

I think my hic files are fine cuz I can view the contact map through juicerbox. Could you @kaukrise take a look on the error and give me some advice? Thanks in advance!!

Metadata

Metadata

Assignees

No one assigned

    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