View source: R/find.gene.pseudomarker.R
find.gene.pseudomarker | R Documentation |
Pull out the pseudomarker that is closest to the position of each of a series of genes.
find.gene.pseudomarker(cross, pmap, geneloc, where = c("prob", "draws"))
cross |
An object of class |
pmap |
A physical map of the markers in |
geneloc |
A data frame specifying the physical locations of the genes.
There should be two columns, |
where |
Indicates whether to pull pseudomarkers from the genotype
probabilities (produced by |
We first convert positions (by interpolation) from those contained within
cross
to physical coordinates contained in pmap
. We then use
qtl::find.pseudomarker()
to identify the closest pseudomarker to
each gene location.
We also include the positions of the pseudomarkers, and we print a warning message if pseudomarkers are > 2 Mbp from the respective gene.
A data frame with columns chr
(the chromosome) and
pmark
(the name of the pseudomarker). The third column pos
contains the Mbp position of the pseudomarker. The final column is the
signed distance between the gene and the pseudomarker. The rownames
indicate the gene names.
Karl W Broman, broman@wisc.edu
qtl::find.pseudomarker()
,
qtl::find.pseudomarkerpos()
, plotEGclass()
,
disteg()
, calc.locallod()
data(f2cross, expr1, genepos, pmap)
library(qtl)
# calc QTL genotype probabilities
f2cross <- calc.genoprob(f2cross, step=1)
# find nearest pseudomarkers
pmark <- find.gene.pseudomarker(f2cross, pmap, genepos, "prob")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.