Question: macs2 callpeak error
0
gravatar for jolin0701-dy
3.4 years ago by
jolin0701-dy70
jolin0701-dy70 wrote:

$ macs2 callpeak -t RIP.bam -c Control.bam -f BAM -g mm -n test -B -q 0.01

INFO @ Wed, 04 May 2016 14:56:39: start X-correlation...

Traceback (most recent call last):

File "/usr/local/bin/macs2", line 617, in <module> main()

File "/usr/local/bin/macs2", line 57, in main run( args )

File "/Library/Python/2.7/site-packages/MACS2/callpeak_cmd.py", line 178, in run opt = options

File "MACS2/PeakModel.pyx", line 108, in MACS2.PeakModel.PeakModel.__init__ (MACS2/PeakModel.c:2438)

File "MACS2/PeakModel.pyx", line 150, in MACS2.PeakModel.PeakModel.build (MACS2/PeakModel.c:3010)

File "MACS2/PeakModel.pyx", line 230, in MACS2.PeakModel.PeakModel.__paired_peak_model (MACS2/PeakModel.c:4266)

ValueError: max() arg is an empty sequence

But when I run the command below, it works fine.

$macs2 callpeak -B -t RIP.bam -c Control.bam --broad -g mm --broad-cutoff 0.1 -n CZ3RIP_vs_CZ1_Input --nomodel --extsize=25 --gsize mm

Any suggestions?

Thanks so much~~

chip-seq • 1.8k views
ADD COMMENTlink written 3.4 years ago by jolin0701-dy70
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: 1064 users visited in the last hour