metaseq.plotutils.input_ip_plots

metaseq.plotutils.input_ip_plots(iparr, inputarr, diffed, x, sort_ind, prefix=None, limits1=(None, None), limits2=(None, None), hlines=None, vlines=None)[source]

All-in-one plotting function to make a 5-panel figure.

Panels are IP, input, and diffed; plus 2 line plots showing averages.

Parameters:
  • inputarr (iparr,) – NumPy arrays constructed by a genomic_signal object
  • diffed – Difference of iparr and inputarr, but can be some other transformation.
  • x – Extent to use – for TSSs, maybe something like np.linspace(-1000, 1000, bins), or for just bin IDs, something like np.arange(bins).
  • sort_ind – row order for each of the 3 panels – usually interesting to use clustered_sortind or tip_zscores
  • prefix – Used to prefix plot titles with ‘%(prefix)s IP”, etc
  • limits1 – Tuple passed to the Normalize function for IP and input.
  • limits2 – Tuple passed tot he Normalize function for the diffed array
  • hlines – List of (position, kwarg) tuples for plotting horizontal lines. Kwargs are passed directly to axhline. Useful for delimiting clusters, if you used clustered_sortind and have both row_order and breaks.
  • vlines – List of (position, kwargs) tuples. A vertical line will be plotted at each position using kwargs.