Generate '.hcs' file from '.cool' raw data #3
Description
Hi,
I've tried to generate .hcs file from .cool file, but without any success. Error traceback as following:
Traceback (most recent call last):
File "/home/diablo/analysis/igm-master/demo/a.py", line 12, in
m.matrix = alabtools.matrix.sss_matrix((spm.data, spm.indices, spm.indptr))
File "/home/diablo/mamba/envs/igm/lib/python3.6/site-packages/alabtools/matrix.py", line 95, in init
self._pop_diag()
File "/home/diablo/mamba/envs/igm/lib/python3.6/site-packages/alabtools/matrix.py", line 151, in _pop_diag
self.diagonal += td
ValueError: operands could not be broadcast together with shapes (26350,) (26348,) (26350,)
I've noticed other file formats are also supported in alabtools.Contactmatrix function, however, the '_load_hic' method seems not implemented when using .hic file as input.
The genome version that I used is mm10, Is this kind of problem raised due to this genome is not supported?
So, is there any way to solve this problem? Thanks!