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

cLoops2 montage error #11

Open
YuZhang337 opened this issue Jun 20, 2024 · 0 comments
Open

cLoops2 montage error #11

YuZhang337 opened this issue Jun 20, 2024 · 0 comments

Comments

@YuZhang337
Copy link

Hi Yaqiang @YaqiangCao,
Sorry to bother you again. I tried using the cLoops2 montage command to show the interaction between some enhancers. These enhancers are all included in the loop file.My command line is as follows: cLoops2 montage -f ../HiC-540-samp/chr2-chr2.ixy -bed test.bed -o HiC540-samp-chr2-top -simple -p 20. But got the following error:
Traceback (most recent call last):
File "/zhangyu/miniconda2/envs/cooltools/bin/cLoops2", line 8, in
sys.exit(main())
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/cLoops2/cLoops2.py", line 4104, in main
montage(
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/cLoops2/montage.py", line 523, in montage
cr.fill_between(rad + scaffoldw,
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/cLoops2/montage.py", line 229, in fill_between
self.pax.fill_between([ets, ete], [rad, rad], [rad + h],
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/matplotlib/init.py", line 1423, in inner
return func(ax, *map(sanitize_sequence, args), **kwargs)
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/matplotlib/axes/_axes.py", line 5367, in fill_between
return self._fill_between_x_or_y(
File "/zhangyu/miniconda2/envs/cooltools/lib/python3.8/site-packages/matplotlib/axes/_axes.py", line 5286, in _fill_between_x_or_y
raise ValueError(f"where size ({where.size}) does not match "
ValueError: where size (1) does not match x size (2)

The contents of the bed file are as follows:
$cat HiC540-541-merge-both-active-enhancer-chr2-top10-sorted.txt
chr2 20174305 20179255 non-promoter-peak26257
chr2 20178648 20182995 non-promoter-peak26260
chr2 20205946 20210208 non-promoter-peak26262
chr2 20494099 20498894 non-promoter-peak26275
chr2 20524068 20529013 non-promoter-peak26278
chr2 20581781 20586128 non-promoter-peak26284
chr2 20592901 20597454 non-promoter-peak26287
chr2 20594000 20598277 non-promoter-peak26288
chr2 20650944 20655617 non-promoter-peak26294
chr2 20652637 20657872 non-promoter-peak26296
chr2 20663430 20667936 non-promoter-peak26297
chr2 20665906 20671983 non-promoter-peak26299
chr2 20725026 20729363 non-promoter-peak26307
chr2 29095408 29100372 non-promoter-peak26475
chr2 29096802 29101155 non-promoter-peak26476
chr2 29118859 29123489 non-promoter-peak26480
chr2 29626571 29630860 non-promoter-peak26512
chr2 29627260 29631883 non-promoter-peak26513
chr2 29628079 29632738 non-promoter-peak26514
chr2 29629344 29633863 non-promoter-peak26515

I tried expanding the enhancer peak region, but I still got the same error.Is there anything I need to pay special attention to in this bed file?I look forward to your reply.

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

1 participant