ATAC seq density calculation

Question

This paper Fig2b

b Changes of chromatin accessibility in AMD (n = 14) relative to
normal (n = 11) in all retina samples. Each dot represents one
ATAC-Seq peak. Blue line in the left panel indicates average fold
changes of peaks with the same ATAC-Seq intensity. The percentage of
reduced peaks is shown under the density curve in the right panel.

How is the density calculated or what is being shown there? average ATAC seq signal or Fold change?

StupidWolf · Accepted Answer

They are not referring to density of reads. It is the density of the fold change, the y variable in the plot to left of Fig2b. It tells you the distribution of the foldchange and 92% of the fold change is <0.
See this page for more details on density and histograms
Suppose your data is something like this (purely making this up):
set.seed(555)
fc = rnorm(1000,-0.5,2)
averATAC = -0.5*(fc - 0.5)^2 + 0.2*fc + 3  + rnorm(100,25,10)

You fit a density to it and you get the plot:
par(mfrow=c(1,2))
plot(averATAC,fc)
plot(dens$y,dens$x,xlab = "density",ylab="fc")

Phoenix Mu · Answer

Just count the number of reads overlapping an ATAC-seq peak and divide it by the length of the peak gives you the read density. Each dot is a fold change, but to calculate fold change, average density from all the samples are taken, that is what the "average" is referring to.

Answered by Phoenix Mu on April 15, 2021

ATAC seq density calculation

2 Answers

Add your own answers!

Ask a Question