Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

plot multiple bed files as bar plots below the matrix #36

Open
girdhar3 opened this issue Feb 11, 2019 · 2 comments
Open

plot multiple bed files as bar plots below the matrix #36

girdhar3 opened this issue Feb 11, 2019 · 2 comments

Comments

@girdhar3
Copy link

I am trying to plot H3K27ac peaks, H3K4me3 peaks and TAD regions below the matrix
python HiCPlotter.py -f $file -chr chr22 -o Example -r 10000 -n Neuron_PFC -s 1800 -e 4000 -o Example -ptd 1 -pcd 3 -pcdf $TAD,$H3K4me3_peaks,$H3K27ac_peaks -pdb 1 -ptr 1

I am getting this error

Traceback (most recent call last):
File "HiCPlotter.py", line 2171, in
HiCplotter(**args)
File "HiCPlotter.py", line 663, in HiCplotter
matrix,nums,tricks=read_HiCdata(files[exp],fileHeader,fileFooter,cleanNANs,smoothNoise,window,tadRange,plotInsulation,plotTadDomains,randomBins)
File "HiCPlotter.py", line 61, in read_HiCdata
if plotInsulation or plotTadDomains and not randomBins: nums,tricks=insulation(matrix,ins_window,rel_window)
File "HiCPlotter.py", line 503, in insulation
if current[0] not in pBorders and current[0]+1 not in regions: pBorders.append(current[0])
IndexError: index 0 is out of bounds for axis 0 with size 0

@girdhar3
Copy link
Author

when I do this

python HiCPlotter.py -f $file -chr chr22 -o Example -r 10000 -n Neuron_PFC -s 1800 -e 4000 -o Example -pcd 1 -pcdf $TAD $H3K4me3_peaks $H3K27ac_peaks

It gives Upps!! Please provide equal number of HiC matrix and custom domains

@akdemirlab
Copy link
Owner

Hi,
If you'd like to plot your data as bar plots then you should use "--barPlots"
python HiCPlotter.py -f $file -chr chr22 -o Example -r 10000 -n Neuron_PFC -s 1800 -e 4000 -o Example -ptd 1 -bl TAD,K4me3,K27Ac -b $TAD,$H3K4me3_peaks,$H3K27ac_peaks -pdb 1 -ptr 1

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants