This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Algorithmic Difference between MACS and MACS2

I think I have a pretty good understanding of how the MACS algorithm works. In a nutshell:

  1. Remove duplicate reads.
  2. Model a Poisson distribution for the background based on genome size and read count.
  3. Detect candidate peaks based on p-value.
  4. Model Poisson distributions for local regions of 1 kbp, 2 kbp, 5 kbp, and 10 kbp around each candidate peak. Take the max of these for each point.
  5. Calculate the p-value with respect to these local regions for each candidate peaks to filter out false positives.

However, I'm unable to find any information on the algorithmic difference between MACS and MACS2. I'm wondering what changes, if any, MACS2 makes to this algorithm. Also, what is the difference between the algorithms for finding broad peaks and narrow peaks in MACS2? Anybody have an idea where I can find out more about this?

chip-seq peak calling dnase-seq macs2 macs
• 7,165 views • link
written by taraeicher • taraeicher

See here for some links.

• 0 views • link
written by h.mon h.mon

The change log is really helpful. Thanks!

• 0 views • link
written by taraeicher • taraeicher

1 answer

The change log details the changes made in each version of MACS/MACS2. It can be found here: https://github.com/taoliu/MACS/blob/master/ChangeLog

• 0 views • link
written by taraeicher • taraeicher

Log in to answer this question.

More posts like this