Find covariation of different DNA-binding protein binding between two conditions I have ChIP-seq data of RNAPII and two other factors which we think follow RNAPII during transcription in two different conditions, and I'd like to show that the genes that lose RNAPII also lose the other two factors, and same for gains.
Do you know any methods that can calculate some sort of covariation coefficients?
I know there are methods that calculate covariation of genes between two conditions, but none for different measures on the same gene.
Answer from @stupidwolf converted from comment:
I would try something like this, for every gene you have a logFC of RNAPII in control versus treatment. Do this for other 2 proteins, and you correlate the logFC. In very rough terms, correlation is the scaled form of covariance, so this will work.
If you have 1000 genes, do correlation between 1000 logFC from RNAPII against 1000 logFC from other proteins.
Answer from @m converted from comment:
You use a covariation stat and then perform Pearson's correlation to standardise the response between -1 to 1. An answer of &gt;0.7 is considered robust and &gt;0.8 highly robust. ITs an interesting question and you are right covariation is the way to answerr it