I am working on a fast version of mirror manhattan plot that essentially repalces geom_point with geom_scattermore. I only modifed few places in the gmirror.R from hudson package.
gmirror_fast(top=gwas.t, bottom = gwas.b, tline=5e-8, bline=5e-8)
I also added an option of using pseduo-log-transformed y-axis. This is helpful for visualizing some traits (e.g. height) with many small p-values (and therefore much larger on -log10(p) than rest of the SNPs).
gmirror_fast(top=gwas.t, bottom = gwas.b, tline=5e-8, bline=5e-8, log_axis=TRUE)
Annotations and highlight are same as in the hudson package.
rsid <- gwas[gwas$pval<5e-30,]$rsid
gmirror_fast(top=gwas.t, bottom = gwas.b, tline=5e-8, bline=5e-8, log_axis=TRUE)
