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

some questions #66

Open
ahmedabbas811 opened this issue Jun 12, 2023 · 3 comments
Open

some questions #66

ahmedabbas811 opened this issue Jun 12, 2023 · 3 comments

Comments

@ahmedabbas811
Copy link

Dear all,
Why the predicted and real contact frequencies in the '.scc' file are too small (~0.0001)?
And how can I get values like the original values in the contact.gz files?

Thanks

@polyaB
Copy link
Contributor

polyaB commented Jun 13, 2023

  1. Contact values are so small because we divided each contact value by normalized coefficient which we calculated for each cell type. You can find them here https://github.com/labdevgen/3Dpredictor/tree/master/input/normalized_coefficients . The calculation method with specific formula you can find in the Methods Hi-C data processing section in our paper "Quantitative prediction of enhancer-promoter interactions". Genome Res. 2020 Jan;30(1):72-84. doi: 10.1101/gr.249367.119.
  2. You can get these files from .hic file using Juicer tools dump

@ahmedabbas811
Copy link
Author

Thank you,

I have one more question,

Suppose I train with odd chr and test with even chr. Suppose I test with chr 8 which has in the contacts.gz file, it has 1013818 loops

the output ".scc" which has original loops vs predicted ones to calculate the SCC value, this file has 324531 loops only. Is it expected to predict this less number of loops? If expected, what is the reason? if not expected, what can be the mistake that I did in running the program?

Thanks a lot

@polyaB
Copy link
Contributor

polyaB commented Jun 20, 2023 via email

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