Description Usage Arguments Details Value See Also Examples
View source: R/supportinterval.R
Calculates support interval for QTL based on Wald profile and QTL position
1 | supportinterval(x, chr, lodsupport = 1)
|
x |
Object of class |
chr |
Selected chromosomes |
lodsupport |
Size of support interval; default is 1 LOD |
Computes the x-LOD support interval as the region surrounding a QTL peak in which the Wald profile exceeds the equivalent of x LOD less than the peak value.
A list with two components: a matrix containing lower and upper bounds for the support intervals for each QTL, and the positions of QTL on each chromosome.
1 2 3 4 5 6 | map <- sim.map(len=100, n.mar=11, eq.spacing=TRUE, include.x=FALSE)
sim.ped <- sim.mpped(4, 1, 500, 6, 1)
sim.dat <- sim.mpcross(map=map, pedigree=sim.ped, qtl=matrix(data=c(1, 10, .4, 0, 0, 0, 1, 70, 0, .35, 0, 0), nrow=2, ncol=6, byrow=TRUE), seed=1)
mpp.dat <- mpprob(sim.dat, program="qtl", step=2)
mpq.dat <- mpIM(object=mpp.dat, ncov=0, responsename="pheno")
si <- supportinterval(mpq.dat)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.