| bayesint {qtl} | R Documentation |
Calculate an approximate Bayesian credible interval for a particular
chromosome, using output from scanone.
bayesint(results, chr, prob=0.95, lodcolumn=1)
results |
Output from scanone. |
chr |
A chromosome ID. |
prob |
Probability coverage of the interval. |
lodcolumn |
An integer, or vector of 3 integers, indicating which of the LOD score columns should be plotted (generally this is 1). |
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@biostat.wisc.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)