| bayesint {qtl} | R Documentation |
Calculate a Bayesian probability interval for a particular chromosome, using output from scanone.
bayesint(results, chr, prob=0.95)
results |
Output from scanone. |
chr |
A chromosome ID. |
prob |
Probability coverage of the interval. |
We take 10^LOD, rescale it to have area approximately 1, and the look for an interval whose area is the target probability.
An object of class scanone, like the input, indicating the
position with the maximum LOD, and indicating approximate endpoints
for the Bayesian probability interval.
Karl W Broman, kbroman@jhsph.edu
data(listeria) listeria <- subset(listeria,chr=c(5,13)) listeria <- calc.genoprob(listeria, step=0.5) out <- scanone(listeria, model="np") bayesint(out, 5) bayesint(out, 13) bayesint(out, 13, 0.99)