Algorithmic Difference between MACS and MACS2
1
1
Entering edit mode
7.1 years ago
taraeicher ▴ 50

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 • 6.0k views
ADD COMMENT
0
Entering edit mode

See here for some links.

ADD REPLY
0
Entering edit mode

The change log is really helpful. Thanks!

ADD REPLY
0
Entering edit mode
7.1 years ago
taraeicher ▴ 50

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

ADD COMMENT

Login before adding your answer.

Traffic: 3067 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6