Skip to contents

Creates a manhattan plot.

Usage

gg_Manhattan(
  folder = "GWAS_Results/",
  trait = list_Traits(folder)[1],
  title = trait,
  threshold = NULL,
  sug.threshold = NULL,
  chrom = NULL,
  markers = NULL,
  labels = markers,
  vlines = markers,
  vline.colors = rep("red", length(vlines)),
  vline.types = rep(1, length(vlines)),
  vline.legend = T,
  facet = F,
  addQQ = T,
  pmax = NULL,
  models = c("MLM", "FarmCPU", "BLINK", "MLMM", "GLM", "CMLM", "SUPER"),
  model.colors = c("darkgreen", "darkorange3", "steelblue", "darkred", "darkorchid4",
    "burlywood4", "darkseagreen4"),
  highlight.sig = F,
  sig.col = "darkred",
  chrom.colors = rep(c("darkgreen", "darkgoldenrod3"), 30),
  chrom.unit = "100 Mbp",
  legend.rows = 1,
  plotHBPvalues = F,
  skyline = NULL
)

Arguments

folder

Folder containing GWAS results.

trait

The trait to read.

title

A title for the plot.

threshold

Significant Threshold.

sug.threshold

Suggested threshold.

chrom

Chromosomes to plot. Use if you want to plot a single chromosome.

markers

Markers to be labelled.

labels

Labels to be used for markers.

vlines

Markers which will be used as a location for a vertical lines.

vline.colors

colors for each vertical line.

vline.types

lty for each vertical line.

vline.legend

Logical, whether or not to add a legend for the vlines.

facet

Logical, whether or not to produce a facetted or multi-model plot. Default is `facet = F`

addQQ

Logical, whether or not to add a QQ plot

pmax

A max value for the y-axis.

models

Models to read.

model.colors

Colors for each model. Used if `facet = F`.

highlight.sig

Logical, whether or not to highlight significant associations with a black circle. Used if `facet = F`.

sig.col

Color for significant assoctiations. Used if `facet = T`.

chrom.colors

Colors for each chromosome. Used if `facet = T`.

chrom.unit

Unit for the x-axis. Can be one of c("kbp","100 kbp","Mbp","100 Mbp","Gbp").

legend.rows

Number of rows for the legend.

plotHBPvalues

Logical, if TRUE, H.B.P.values be uses.

skyline

Which skyline type to use. Can be "NYC" or "Kansas". If left NULL, it will use the highest P.value.

Value

A manhattan plot.