Index Error when running hic_data.find_compartments

[LAlvarezGlez]

Hi!

I am trying to run hic_data.find_compartments and no matter what parameter I change or what input data I use I keep getting this error:

eigenvectors = hic_data.find_compartments(label_compartments="none", perc_zero=99, max_ev=5) File "/home/PYTHON/ENVS/TADbit/lib/python2.7/site-packages/pytadbit/hic_data.py", line 884, in find_compartments if max_ev else len(matrix))] IndexError: index -3 is out of bounds for axis 1 with size 2

I am new with this kind of analysis so I have no clue on what is happening.

Many thanks for your help

Lucía

[david-castillo]

Hi Lucía,

I would first try to run the command with the default parameters.

eigenvectors = hic_data.find_compartments()

If you still get the error I recommend you to update your TADbit version to the master branch, it could be that there was a bug in that spot and it's been corrected in the last version.

Regards

David