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

Hichipper Inter-chromosomal Loop Statistics #66

Open
bdecker01 opened this issue Feb 22, 2019 · 2 comments
Open

Hichipper Inter-chromosomal Loop Statistics #66

bdecker01 opened this issue Feb 22, 2019 · 2 comments

Comments

@bdecker01
Copy link

bdecker01 commented Feb 22, 2019

Hello Hichipper Software Programmers,

I am using your software to analyze some Hi-C / Hi-ChIP data for a research project I am working on, and I am wondering if there is a way to output this same "per-loop FDR measures from the loop proximity bias correction algorithm originally implemented in Mango and presented in the same format" that is output for 'interactions.all.mango' intra-chromosomal loops, for the inter-chromosomal loops results, or if some similar statistics can be easily applied post run on the '.inter.loop_counts.bedpe' output files.

Thanks for the help!
Brandon

@caleblareau
Copy link
Contributor

Hi Brandon,

I can see why this would be of interest, but it's hard to imagine how to make a statistic that would be applicable here. Every statistic for computing an FDR value would be reliant on some background model, which for intrachromosomal interactions is a function of genomic proximity. It's hard for me to imagine what an appropriate background model is for interchromosomal interactions.

To these ends, I'd recommend just using the raw counts between loci as I don't see any reason why they would be biased for any pair of points? You could compute an empirical distribution function to attain something resembling a p-value on the raw counts.

Or do you have an idea of how one could create a background model?

@bdecker01
Copy link
Author

Hello Caleb,

Thank you for the advice, I computed the probability for the scores based on their empirical distribution and it looks like that gives me a pretty reasonable statistical output. I am still trying a few other different approaches to determine statistical significance for the interchromosomal interactions, but I will let you know if I find some other appropriate ways to determine it.

Thanks again for the help!
Brandon

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