qtl (version 1.01-9)

bayesint: Bayesian probability interval

Description

Calculate a Bayesian probability interval for a particular chromosome, using output from scanone.

Usage

bayesint(results, chr, prob=0.95)

Arguments

results

Output from scanone.

chr

A chromosome ID.

prob

Probability coverage of the interval.

Value

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

Details

We take $10^{LOD}$, rescale it to have area approximately 1, and the look for an interval whose area is the target probability.

See Also

scanone, lodint

Examples

data(listeria)
listeria <- subset(listeria,chr=c(5,13))
listeria <- calc.genoprob(listeria, step=0.5)
out <- scanone(listeria, model="np")
bayesint(out, 5)
bayesint(out, 13)
bayesint(out, 13, 0.99)