intervene icon indicating copy to clipboard operation
intervene copied to clipboard

Support for zero interval files

Open jma1991 opened this issue 6 years ago • 0 comments

I use Intervene in a computational pipeline to check overlaps from multiple BED files. The command I use is:

intervene pairwise --input <INPUT> --compute jaccard --output <OUTPUT>

Sometimes one of the BED files contains zero intervals. When this occurs Intervene throws the following error:

Performing a pairwise intersection analysis. Please wait...

Traceback (most recent call last):
  File "/home/s1437643/conda/bin/intervene", line 606, in <module>
    main()
  File "/home/s1437643/conda/bin/intervene", line 426, in main
    pairwise.pairwise_intersection(label_names, options)
  File "/home/s1437643/conda/lib/python3.6/site-packages/intervene/modules/pairwise/pairwise.py", line 480, in pairwise_intersection
    barplot(series, matrix, outfile, options, max_size=max(bed_sizes))
  File "/home/s1437643/conda/lib/python3.6/site-packages/intervene/modules/pairwise/pairwise.py", line 135, in barplot
    cax, order = heatmap_triangle(matrix, ax, options)
  File "/home/s1437643/conda/lib/python3.6/site-packages/intervene/modules/pairwise/pairwise.py", line 211, in heatmap_triangle
    Z = sch.linkage(D, method='average')
  File "/home/s1437643/conda/lib/python3.6/site-packages/scipy/cluster/hierarchy.py", line 713, in linkage
    raise ValueError("The condensed distance matrix must contain only "
ValueError: The condensed distance matrix must contain only finite values.

Would you be able to provide support for zero interval files? Preferably, the file is still included in all the analysis results rather than just not using the file.

jma1991 avatar Jul 28 '18 14:07 jma1991