$ 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~~