gg_Manhattan_Traits
gg_Manhattan_Traits.Rd
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.