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)