bayesint {qtl}R Documentation

Bayesian credible interval

Description

Calculate an approximate Bayesian credible interval for a particular chromosome, using output from scanone.

Usage

bayesint(results, chr, prob=0.95, lodcolumn=1)

Arguments

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).

Details

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.

Value

An object of class scanone, like the input, indicating the position with the maximum LOD, and indicating approximate endpoints for the Bayesian credible interval.

Author(s)

Karl W Broman, kbroman@biostat.wisc.edu

See Also

scanone, lodint

Examples

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)

[Package qtl version 1.06-43 Index]