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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.