Hi everyone, I'm trying to call peaks with macs2 but run into the next warning when doing so:
"WARNING @ Tue, 14 Jun 2016 20:52:53: Too few paired peaks (0) 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.
WARNING @ Tue, 14 Jun 2016 20:52:53: Process for pairing-model is terminated! "
I am working with treatment and control samples, the number of reads for each is: 18448586 for control and 26581254 for treatment. All the reads are mapped and uniquely paired, as I already filtered them.
Could anyone explain the cause for the warning? I've tried to find information but haven't find many...Is the error related to the difference in number of reads? I thought macs scaled the samples automatically.
Thanks in advance!
Update: So I'm guessing this is the reason:
"DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of unique tags on + strand: 627266
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of peaks in + strand: 19221
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of unique tags on - strand: 0
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of peaks in - strand: 0
DEBUG @ Thu, 16 Jun 2016 16:00:24: Chrom chr16 is discarded!"
And it happens for all the chromosomes. Any ideas on how to fix it?
"WARNING @ Tue, 14 Jun 2016 20:52:53: Too few paired peaks (0) 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.
WARNING @ Tue, 14 Jun 2016 20:52:53: Process for pairing-model is terminated! "
I am working with treatment and control samples, the number of reads for each is: 18448586 for control and 26581254 for treatment. All the reads are mapped and uniquely paired, as I already filtered them.
Could anyone explain the cause for the warning? I've tried to find information but haven't find many...Is the error related to the difference in number of reads? I thought macs scaled the samples automatically.
Thanks in advance!
Update: So I'm guessing this is the reason:
"DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of unique tags on + strand: 627266
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of peaks in + strand: 19221
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of unique tags on - strand: 0
DEBUG @ Thu, 16 Jun 2016 16:00:24: Number of peaks in - strand: 0
DEBUG @ Thu, 16 Jun 2016 16:00:24: Chrom chr16 is discarded!"
And it happens for all the chromosomes. Any ideas on how to fix it?
Comment