bayesint {qtl}R Documentation

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.

Details

We take 10^LOD, rescale it to have area approximately 1, and the look for an interval whose area is 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 probability interval.

Author(s)

Karl W Broman, kbroman@jhsph.edu

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)

[Package qtl version 1.03-14 Index]