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!
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!