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

How to draw a beautiful HIC map #14

Open
songxh1996 opened this issue Oct 24, 2023 · 5 comments
Open

How to draw a beautiful HIC map #14

songxh1996 opened this issue Oct 24, 2023 · 5 comments

Comments

@songxh1996
Copy link

songxh1996 commented Oct 24, 2023

Hello, Professor Wang, thank you for providing so much useful software. I tried a lot of software, such as fanc,hicexplorer, but did not draw the hic map I wanted. I was so lucky to see this beautiful HIC map in you. But I didn't find out how you drew it, whether I missed what it was.
How do you draw a beautiful HIC map? Which script is used

fig1 svg
p4jkei4a

@XiaoTaoWang
Copy link
Owner

Hi, you may refer to this page as a reference https://xiaotaowang.github.io/TADLib/visualize.html

@songxh1996
Copy link
Author

songxh1996 commented Oct 26, 2023

Thank you very much for your help. When I finished running the test data according to your guidance, it was successful. But when I used my own data, it failed. I can't figure out why.

my cool file
##########.hic2.cool
hicConvertFormat -m LD_nomt.allValidPairs.hic --inputFormat hic --outputFormat cool -o KR_matrix.cool -r 10000 --correction_name KR
###########hic-pro result covert cool
cooler load -f coo LD_nomt_10000_abs.bed LD_nomt_10000_iced.matrix 10000.cool --count-as-float

vis = Triangle('GM12878-MboI.hg19.cool::10000', '21', 45000000, 46500000)
vis = Triangle('10000.cool::10000', 'NC_059157.1', 4500000, 4650000)
Traceback (most recent call last):
File "", line 1, in
File "/home/songxianhan/miniconda3/envs/test/lib/python3.7/site-packages/tadlib/visualize/heatmaps.py", line 18, in init
self.clr = cooler.Cooler(uri)
File "/home/songxianhan/miniconda3/envs/test/lib/python3.7/site-packages/cooler/api.py", line 85, in init
self._refresh()
File "/home/songxianhan/miniconda3/envs/test/lib/python3.7/site-packages/cooler/api.py", line 106, in _refresh
raise KeyError(err_msg) from None
KeyError: "No cooler found at: 10000.cool. Coolers found in ['/']. Use '::' to specify a group path"

@XiaoTaoWang
Copy link
Owner

Can you try the following command instead:

vis = Triangle('10000.cool', 'NC_059157.1', 4500000, 4650000)

@songxh1996
Copy link
Author

When I use this code, the error persists.

vis = Triangle('40000.cool', 'NC_059157.1', 4500000, 4650000)
Traceback (most recent call last):
File "", line 1, in
File "/home/songxianhan/miniconda3/envs/tset/lib/python3.7/site-packages/tadlib/visualize/heatmaps.py", line 18, in init
self.clr = cooler.Cooler(uri)
File "/home/songxianhan/miniconda3/envs/tset/lib/python3.7/site-packages/cooler/api.py", line 85, in init
self._refresh()
File "/home/songxianhan/miniconda3/envs/tset/lib/python3.7/site-packages/cooler/api.py", line 106, in _refresh
raise KeyError(err_msg) from None
KeyError: "No cooler found at: 40000.cool. Coolers found in ['/40000']. Use '::' to specify a group path"

@XiaoTaoWang
Copy link
Owner

Can you show the result after you execute the following command?

$ cooler ls 40000.cool

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