Question: MACS2 peaks filtering
0
gravatar for GK1610
6 months ago by
GK161040
United States
GK161040 wrote:

I have created a merged_BAM file by merging 100 individual bam files.

I am calling chip-seq peaks on merged_BAM file using MACS2

I want to obtain p-value threshold to filter out the MACS2 peaks

One of the ways I am doing is using IDR on 3 copies of merged_BAM that are created by randomly subsampling 33% of the original merged_BAM file.

merged_BAM_replicate 1( 33% of original merged_BAM) merged_BAM_replicate 2( 33% of original merged_BAM) merged_BAM_replicate 3( 33% of original merged_BAM)

Do you have other suggestions?

chip-seq • 321 views
ADD COMMENTlink written 6 months ago by GK161040

Why do you want to filter out peaks? MACS has internal q-value cutoffs, is there anything wrong with it? By subsampling, you reduce statistical power. Please explain what your goal is.

ADD REPLYlink modified 6 months ago • written 6 months ago by ATpoint7.5k
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: 569 users visited in the last hour