Adding loghist for logarithmic histogram

Author: mnshgl0110

hometools/hometools.py

@@ -2151,7 +2151,7 @@ def reg_str_to_list(regstr):
 
 def mapbp(args):
     """
-        Outputs mapping positions for the given reference genome coordinate
+    Outputs mapping positions for the given reference genome coordinate
     """
     import pysam
     from collections import defaultdict, deque

hometools/plot.py

@@ -533,4 +533,17 @@ def inlineplotsr(ax, args):
     #     sys.exit('Error in saving the figure. Try using a different backend.' + '\n' + e.with_traceback())
     # logger.info('Finished')
 
-    return ax
+    return ax
+
+
+def loghist(x, ax, bins=10):
+    """
+    Generate the plot on a logarithmic histogram
+    """
+    import numpy as np
+    hist, bins = np.histogram(x, bins=bins)
+    logbins = np.logspace(np.log10(bins[0]), np.log10(bins[-1]), len(bins))
+    ax.hist(x, bins=logbins)
+    ax.set_xscale('log')
+    return ax
+# END