Description Usage Arguments Details Value Author(s) Examples
Expression values are grouped according to genotype groups and then visualized with boxplots.
1 | genotypePlot(snp,gene,eqtl,geneAnnot=NULL, ylab=NULL, xlab=NULL, mainlab=FALSE)
|
snp |
|
gene |
|
eqtl |
An |
geneAnnot |
|
ylab |
Optional x-axis label. |
xlab |
optional y-axis label |
mainlab |
|
This function plots the expression values of the genotype groups of a certain SNP that can be given in the
snp
option. The expression values are specified in the gene
option.
The eqtl
object is the output of an eQTL
run and carries the required genotype information.
A Figure.
Daniel Fischer
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | # Make the example data available
data(Xgene)
data(genotData)
data(annotTrack)
# We need to have the gene annotation in bed format (Please notice the change to the
# official convention, this is on high priority of the ToDo list of the package to change
# this.)
## Not run:
annotBed <- gtfToBed(annotTrack)
# Perform a basic cis-eQTL with the minimum required input linear model:
lm.myEQTL <- eQTL(gex=Xgene,geno=genotData, xAnnot=annotBed,method="LM",windowSize=1)
# Plot the genotypes
genotypePlot(snp="SNP377", gene="MYBPC1", eqtl=lm.myEQTL, ylab="Expression values",
xlab="Genotypes")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.