| bayesint {qtl} | R Documentation |
Calculate an approximate Bayesian credible 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 1, and then calculate the connected interval with density above some threshold and having coverage matching 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 credible interval.
Karl W Broman, kbroman@jhsph.edu
data(hyper) hyper <- calc.genoprob(hyper, step=0.5) out <- scanone(hyper, method="hk") bayesint(out, chr=1) bayesint(out, chr=4) bayesint(out, chr=4, prob=0.99)