Average signal for multiple bw over one bed

usage: signal_plot_bed.py [-h] [-j JID] [-u U] [-d D]
                          [--computeMatrix_addon_parameters COMPUTEMATRIX_ADDON_PARAMETERS]
                          [-c COLORS] --bw_files BW_FILES [BW_FILES ...] -f
                          BED

signal plot 2

optional arguments:
  -h, --help            show this help message and exit
  -j JID, --jid JID     enter a job ID, which is used to make a new directory.
                        Every output will be moved into this folder. (default:
                        signal_plot_bed_yli11_2020-11-13)
  -u U                  upstream flanking length (default: 5000)
  -d D                  downstream flanking length (default: 5000)
  --computeMatrix_addon_parameters COMPUTEMATRIX_ADDON_PARAMETERS
                        add user-defined parameters to computeMatrix (default:
                        )
  -c COLORS, --colors COLORS
                        colors, seperated by comma, hex color is OK (default:
                        red,green,blue,yellow,grey,purple,darkgreen,darkred,pi
                        nk,orange)
  --bw_files BW_FILES [BW_FILES ...]
                        bw_file_list (default: None)
  -f BED, --bed BED     one bed file (default: None)
  --plotProfile_addon_parameters PLOTPROFILE_ADDON_PARAMETERS
                        plotProfile_addon_parameters (default: )

Summary

Given one bed file and multiple bw files, plot the average signal over the bed file. Users can assign colors.

Example

../../_images/signal_plot_on_bed.png

Input

Input bw file names, separated by space.

Output

You will receive the figure through email.

Usage

Copy your bw files and bed files and run the following for each bed file.

You can use Hex color code for colors, e.g., --colors #03fc24,#7522d4, separated by comma.

hpcf_interactive

module load python/2.7.13

signal_plot_bed.py -f loci.bed -u 1000 -d 1000 --bw_files CADD.bdg.bw HbFBase.bdg.bw DeepSEA.bdg.bw --colors red,green,blue

To fix Y-axis, use:

signal_plot_bed.py -f loci.bed -u 1000 -d 1000 --bw_files CADD.bdg.bw HbFBase.bdg.bw DeepSEA.bdg.bw --colors red,green,blue --plotProfile_addon_parameters " --yMin 0 --yMax 10"

code @ github.