I have an issue when running Step 5 of the pipeline. This is the traceback:
Traceback (most recent call last):
File "/work/project/ladsie_002/mHiC/s5_prior.py", line 1050, in <module>
main()
File "/work/project/ladsie_002/mHiC/s5_prior.py", line 319, in main
(x,y,yerr)= calculateProbabilities(mainDic, binStats,resolution,os.path.join(outputPath,libName+".fithic_pass1"), observedIntraInRangeSum)
File "/work/project/ladsie_002/mHiC/s5_prior.py", line 755, in calculateProbabilities
var = M2/(possPairsInRange-1)
ZeroDivisionError: float division by zero
possPairsInRange is read from a dictionary called binStats. I logged the content of this dictionary:
possPairsInRange is the second element of each value, so in the last entry (#43) it goes down to 1, which should be the cause of the above error, since var = M2/(possPairsInRange-1) yields a division by zero.
Would you mind sharing your input values and maybe a few lines from the input files? This step is modified from Fit-Hi-C. Have you successfully run Fit-Hi-C for significant interaction detection before?
I have an issue when running Step 5 of the pipeline. This is the traceback:
possPairsInRange
is read from a dictionary calledbinStats
. I logged the content of this dictionary:possPairsInRange
is the second element of each value, so in the last entry (#43) it goes down to 1, which should be the cause of the above error, sincevar = M2/(possPairsInRange-1)
yields a division by zero.