DOQTL (version 1.8.0)

plot.doqtl: Plot a QTL

Description

Given a genome scan produced by scanone, create a plot of the LOD score across all chromosomes.

Usage

"plot"(x, stat.name = c("lod", "neg.log10.p"),  sig.thr = NULL, sig.col = "red", ...)

Arguments

x

DOQTL object containing the LOD score and model coefficients.

stat.name

Character string containing the name of the mapping statistic.

sig.thr

Numeric matrix containing significance thresholds. Columns must be labelled "A" and "X", in that order and should contain thresholed produced by get.sig.thr().

sig.col

Colors to use when plotting the significance thresholds above. Must be the same length as the nuber of rows in sig.thr.

...

Additional arguments to pass to plot.

Value

Creates a QTL plot.

Examples

  ## Not run: 
#     qtl = scanone(pheno = pheno, pheno.col = 1, probs = probs, snps = snps)
# 	plot(qtl)
#   ## End(Not run)