Description Usage Arguments Value Author(s) See Also Examples
Given output from neqtl
, finds positions for maximum numbers of
eQTL. Will find next highest maximum outside of window
from
peak, for n.peaks
number of maximums
1 |
neqtl.out |
Output from |
lodcolumn |
Numeric value for which LOD column to use, default is first column |
window |
Numeric value of window size, greater than or equal to |
n.peaks |
Number of peaks to find (e.g. for n.peaks=10, the top 10 positions of maximum values) |
data.frame of chr and pos of peak positions
Karl W. Broman and Aimee Teo Broman
1 2 3 4 5 6 7 8 9 | data(fake.f2expr)
fake.f2expr <- calc.genoprob(fake.f2expr)
scan.f2 <- scanone(fake.f2expr,pheno.col=1:100, method="hk")
maxlod.f2 <- maxlod(scan.f2)
sigpos.f2 <- maxlod.sigpos(maxlod.f2,sig.lod=3)
n.f2 <-neqtl(sigpos.f2,chr=scan.f2[,1],pos=scan.f2[,2],win=5)
findpeaks(n.f2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.