GeneticTools (version 0.3.1)

plot.eqtl: Plot an eqtl Object

Description

The function offers informative plots for an eqtl object.

Usage

"plot"(x, file = NULL, which = NULL, sig = 0.01, verbose = TRUE, centered = TRUE, log = FALSE, x2 = NULL, annot = NULL, track = NULL, trackAnnot = FALSE, trackOrder = NULL, mc.cores = 1, ...)

Arguments

x

Object of class eqtl.

file

Store set of graphics under that file name.

which

Specifies for which genes should the plot be created.

sig

Chosen significance level.

verbose

Logical, extended feedback of the function.

centered

Logical, plot should be centered around center gene.

log

Logical, y-axis scale is log(base=10)-scaled.

x2

Comparison values of a second eqtl object, see details.

annot

Logical, plot annotation track.

track

Gene annotations in bed format.

trackAnnot

Gene annotations in bed format.

trackOrder

Logical, shall the annotation track be ordered.

mc.cores

Amount of cores for parallel computing.

...

Additional plotting parameters.

Details

This function plots the test results of an eqtl object. Typically is the tested gene in the center and the p-values of associated SNPs are visualized. Monomorphic SNPs and those that were missing are separately plotted. Test results that are smaller than the value given to sig are marked in red. The y.axis can be switched to log10 scale by setting the logical parameter log=TRUE in that case are bars instead of dots plotted. If the y-axis is on log-scale it is also possible to give a second eqtl object to the function and plot the test results for both. The annotation feature is currently under development and only available in limited form.

Examples

# Perform eQTL (single location, one gene):
  data(Xgene) 
  data(genotData)
  data(annotTrack)
## Not run: 
#   annotBed110 <- gtfToBed(annotTrack[1:10,])
# 
#   lm.myEQTL <- eQTL(gex=Xgene,geno=genotData, xAnnot=annotBed110,method="LM")
#   
#   plot(lm.myEQTL)
#   
#   dir.myEQTL <- eQTL(gex=Xgene,geno=genotData, xAnnot=annotBed110,method="directional")
#   plot(lm.myEQTL, x2=dir.myEQTL, log=TRUE, sig=2)
# ## End(Not run)