hpdchr: Find highest probability density (HPD) region by Chromosome.

In byandell/qtlbim: QTL Bayesian Interval Mapping

Description

Find area under posterior by chromosome for genome-wide HPD region.

Usage

qb.hpdchr(qbObject, level = 0.5, height, hpd, chr, smooth = 3)

Arguments

qbObject	Object of class qb.
level	Value between 0 and 1 of HPD coverage. Ignored if height is supplied.
height	Height of posterior corresponding to level. Determined from level and hpd by default.
hpd	Object of class qb.hpdone. Created internally be default.
chr	Chromosomes to include; default determined by genome-wide HPD region.
smooth	Degree of smoothing.

Details

Determine 100*level percent HPD coverage by chromosome for a given height threshold.

Value

A List with the following elements:

hpd.height	Height of HPD threshold. Name of hpd.height is the level.
chr.posterior	Posterior probability by chromosome as percent.

Author(s)

Brian S. Yandell

References

http://www.qtlbim.org

See Also

qb.scanone

Examples

data(qbExample)

qb.hpdchr(qbExample)

byandell/qtlbim documentation built on Dec. 19, 2021, 12:47 p.m.