plot.eqtl: Plot an eqtl Object

In GeneticTools: Collection of Genetic Data Analysis Tools

Description

The function offers informative plots for an eqtl object.

Usage

 ## S3 method for class 'eqtl'
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.

Author(s)

Daniel Fischer

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)

GeneticTools documentation built on Jan. 15, 2017, 11:21 a.m.