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)