Skip to contents

Creates a manhattan plot.

Usage

gg_Manhattan_Traits(
  folder = "GWAS_Results/",
  traits = list_Traits(folder)[1:2],
  title = NULL,
  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 = F,
  facet = F,
  addQQ = T,
  pmax = NULL,
  models = "MLM",
  highlight.sig = F,
  sig.col = "darkred",
  trait.colors = c("darkgreen", "darkorange3", "steelblue", "darkred", "darkorchid4",
    "burlywood4", "darkseagreen4"),
  chrom.unit = "100 Mbp",
  legend.rows = 1,
  plotHBPvalues = F,
  skyline = NULL
)

Arguments

folder

Folder containing GWAS results.

traits

The traits 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.

addQQ

Logical, whether or not to add a QQ plot.

pmax

A max value for the y-axis.

models

Model to read.

highlight.sig

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

trait.colors

Colors for each trait.

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.

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.