Hi all, I am trying to plot manhattan plot from my GWAS object from rrBLUP. rrBLUP does display Manhattan plots with significant threshold line using qvalue package from Bioconductor. I am instead trying to make my own custom Manhattan plot (using qqman R package) since it is a lot customizable with significant threshold line. For that I am using the same qvalue package as follows:

qobj<- qvalue (p, fdr.level = 0.05)

The qobj has following components:

names (qobj)

[1] "call" "pi0" "qvalues" "pvalues" "lfdr" "fdr.level"
[7] "significant" "pi0.lambda" "lambda" "pi0.smooth"

I can use "significant" component for downstream analysis. But how can I get that single significant threshold value to plot in my custom Manhattan plot?

Your help is much appreciated.

Thank you,