Question: Macs2 failed for pooled samples
0
gravatar for RT
19 months ago by
RT330
European Union
RT330 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 • 809 views
ADD COMMENTlink modified 19 months ago by EagleEye5.7k • written 19 months ago by RT330
0
gravatar for EagleEye
19 months ago by
EagleEye5.7k
Sweden
EagleEye5.7k 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 19 months ago • written 19 months ago by EagleEye5.7k

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 19 months ago by RT330
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: 1011 users visited in the last hour