This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Manhattan plot of Fisher's Exact Test results

I've been trying to plot FET data (scaffold, position, and p-value) using a Manhattan plot with scaffolds ordered by decreasing size. There were around 77 million SNPs identified and the script I'm using keeps silently failing and not producing an output file or informative error. This is the first time I've done this and I'm wondering how other people handle this situation? I've tried to find other conversations about this online but haven't found anything yet. I'm using a HPC and R script, but am open to anything at this point. Anything helps. Thank you.

manhattanplot

1 answer

This is a common problem for manhattan plots. An easy way around this issue is to just subset p-values below 0.05 and only plot those. This will give you a simple illustration of the important results you can work with in the meantime. When it comes to making a publication quality plot, you can subset a few thousand random p-values above 0.05 per chromosome and this will fill in the background. The vast majority of the non-significant data points will overlap so they become redundant very quickly.

There is a note about filtering and making nice plots in RGraphGallery's page on Manhattan Plots.

Log in to answer this question.