Cut-off value from qvalue package R for Manhattan plot
Entering edit mode
12 weeks ago
b.g.tamang ▴ 20

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,

GWAS qvalue threshold significant plot manhattan rrBLUP • 190 views

Login before adding your answer.

Traffic: 2024 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6