Invalid divide value in brie_MH_Heuristic

[trevortknguyen]

Dearest Yuanhua,

Running brie with a specific bam file generates this error. Other bam files that I have tried work perfectly fine.

Traceback (most recent call last):                                                                      
  File "/home/user/anaconda3/envs/brie-env/bin/brie", line 11, in <module>
    sys.exit(main())
  File "/home/user/anaconda3/envs/brie-env/lib/python2.7/site-packages/brie/brie.py", line 268, in main
    M=M,  Mmin=initial, gap=gap, nproc=nproc)
  File "/home/user/anaconda3/envs/brie-env/lib/python2.7/site-packages/brie/models/model_brie.py", line 451, in brie_MH_Heuristic
    FPKM_all = Cnt_all / tranLen.reshape(-1, 1) / total_count * 10**9
RuntimeWarning: invalid value encountered in divide 

I believe that this is some sort of divide by zero error that prevents the rest of the analysis from running. No output files are created.

I have yet to fully investigate your code myself, but would there be any conditions that could cause one of the divisors to be zero?

Regards,
Trevor