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