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
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment