I was wondering, if I have a Hi-C matrices and the average ICE-normalized interaction vs. genomic distance, how can I create a plot like the following one?
I created something similar:
but it's not contact probability. I want to draw as well the fractal lineal regression.
I'd like to refer you to Organization of the Mitotic Chromosome, Naumova et al 2013 Science. There is a paragraph about calculating contact probability in supplemental, as well as some references.
I cant see the first image, but from the second, it looks like you would be better off by doing logarithmic binning. I made that plot (and the log-binned one) using code from this tutorial:
Hi,
can you solve this problems? I have trouble in this figure. I read the paper about Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome, But I am more confused. So, if you could solve this, could you add a comment? Or how to calculate the frequency is helpful.