Pipeline didn't provide "_model.r " outputs from MEDIP-seq data

[chantisakee]

Hi, I've analysed Medipseq data using MACS2

here is my script:

macs2 callpeak -t Tumor_rmdup.bam -c Control_rmdup.bam \ -f BAMPE -g hs \ -p 0.01 \ -n RCC10_L4 \ --bdg \ --outdir ${OUTPUT_DIR}

here is all the outputs i've got : RCC10_L4_peaks.narrowPeak RCC10_L4_summits.bed RCC10_L4_peaks.xls ;

Just wondering Why didn't I get another "_model.r" output from this script? and How do I get them?

Thank you in advance, Chantisa

[taoliu]

@chantisakee With -f BAMPE, MACS entered the Paired-end mode. In this mode, the 'fragment size' MACS uses to build pileup depends on the actual insertion size detected from read alignment. Therefore, it won't build the read shifting model and there won't be any model.r file output.