bayesint: Find a Bayesian Credible Interval around a QTL.

View source: R/bayesint.R

bayesintR Documentation

Find a Bayesian Credible Interval around a QTL.

Description

This function normalizes the area under the QTL curve on the given chromosome and finds a region that is 95

Usage

  bayesint(qtl, chr, prob = 0.95, expandtomarkers = TRUE)

Arguments

qtl

data.frame: four columns with SNP ID, Chr, position, and LOD in each column.

chr

character: the chromosome on which the QTL lies.

prob

numeric: must be between 0 and 1

expandtomarkers

boolean: if TRUE, expand the QTL interval to the nearest flanking markers. Default = TRUE.

Value

Data frame with the SNP ID, Chr, position and LOD for the left and right side of the interval and the maximum QTL.

Author(s)

Daniel Gatti

References

Saunak, S. (2001) A Statistical Framework for Quantitative Trait Mapping. Genetics, 159 (1), 371–387.

See Also

scanone, scanone.perm

Examples

  ## Not run: bayesint(qtl, 1)

dmgatti/DOQTL documentation built on April 7, 2024, 10:35 p.m.