Question: Macs2 failed for pooled samples
0
gravatar for RT
2.7 years ago by
RT340
European Union
RT340 wrote:

Dear All,

I am new to ChIP analysis. I have 3 replicates for ChIP samples and 2 for control samples. I am using macs2 for peak calling. It was all going well until I called peaks on pooled replicates.

MACS2 command:

macs2 callpeak -t chipdata_pooled.tagAlign.gz -c mocksamples.tagAlign.gz -f BED -n chipdata_pooled_VS_mocksamples -g 1.1e8 -p 1e-3

I get the following error and the program terminates:

Too few paired peaks (91) so I can not build the model! Broader your MFOLD range parameter may erase this error. If it still can't build the model, we suggest to use --nomodel and --extsize 147 or other fixed number instead.

I did not get the above issue with the individual replicates. Can anyone please help.

Many thanks,

chip-seq macs2 • 1.1k views
ADD COMMENTlink modified 2.7 years ago by EagleEye6.5k • written 2.7 years ago by RT340
0
gravatar for EagleEye
2.7 years ago by
EagleEye6.5k
Sweden
EagleEye6.5k wrote:

Hi,

I have a suggestion. Instead of pooling the samples, I would recommend calling peaks for individual replicates and use IDR.

ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by EagleEye6.5k

I am very confused now. I actually started following IDR pipeline. Does not it recommend to run macs2 on both individual replicates and pooled replicates? https://sites.google.com/site/anshulkundaje/projects/idr

ADD REPLYlink written 2.7 years ago by RT340
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1464 users visited in the last hour