cLoops2 icon indicating copy to clipboard operation
cLoops2 copied to clipboard

cLoops2 montage error

Open YuZhang337 opened this issue 8 months ago • 0 comments

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.

YuZhang337 avatar Jun 20 '24 08:06 YuZhang337