knitr::opts_chunk$set(fig.width=8, fig.height=6)

For this demonstration of the qqplotr package we will work with GWAS summary statistics from the PGC's cross disorder analysis.

# devtools::install_github("aaronwolen/pgcxd")
library(pgcxd)
scz

First, let's create a qqplot for only the schizophrenia results:

library(gqq)
qq_plot(scz$pval, highlight = 6)

Next, let's add points for the bipolar results:

bip
<<qqplot-scz>>
qq_plot(bip$pval, qpoints = TRUE, col = "red", highlight = 6)


aaronwolen/gqq documentation built on Feb. 17, 2024, 11:27 a.m.