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)