p-value for peak-calling using MACS
Hi I am trying to do peak-calling using MACS for some of my H3K4me3 ChIP-seq data. Using the default p-value of 1e-5 I only obtained 500+ peaks. But if I change the p-value to 1e-4 I got 1500+ peaks and by increasing it more to 1e-3 I got 7000+ peaks! My question is how do I know which is the best p-value to use? I would like to get more peaks but I also don't want any false positive peaks. Thanks for your help!
-
-
I am also very interested in knowing this because I am using MACS to do the peak calling too -
Depends how many false positives you're willing to handle. The default is 1E-5 and if your enrichment is good enough you should have many thousand peaks at that level.Comment
-
There is no best number and it completely depends on your data set. Load the data onto a genome visualizer and look at the peaks it is calling. Use your best judgement and validate by qPCR. Fold enrichment is another important and easy to understand parameter to consider.--------------
EthanComment
-
How do you mean? Are you talking about the BAM file from Bowtie or the BED file from MACS? I find that the output from MACS is just a bar. On the other hand the BAM files are quite big and time out on me when I try to visualize them.Comment
-
dgrayson, if you use the -w option when running MACS you will generate a wig file which can be uploaded and visualized on a genome browser. I like to use the -w option with the -S option so that I only generate a single wig file instead of one for each chromosome.Comment
-
Thanks- I've tried that. I'm running MACS in a local Galaxy server which does not produce wig output files (MACS 1.4 in Galaxy). In addition this version of Galaxy does not allow us to view our output in UCSC browser.Comment
-
Maybe you can try using the web-based Galaxy? It allows you to generate a wig file, which is, in my opinion, really important for visualization of your chip datasets.Comment
-
igv bam bai bed files
you can use IGV viewer. view bam ( samtools index to get bai file) of experiment control and the bed file to see an overall view.. bed will be bars of the peaks but bam with index will show tags and enrichment.Comment
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...Today, 05:41 AM -
CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).
Despite this, “CRISPR helped turn genome editing from a specialized technique into...07-31-2026, 11:01 AM
Comment