sample.bicreg.qtl.models {BayesQTLBIC}R Documentation

Sampling models for BIC analyses of multiple chromosomes

Description

Models are sampled from each of a set of runs with different sets of independent variables (typically corresponding to multiple chromosomes), according to their posterior probabilities.

Usage

sample.bicreg.qtl.models(chrom.fits,nsim,maxtries=10)

Arguments

chrom.fits list of bicreg.qtl objects from separate fits by chromosome or genomic region
nsim number of models to sample
maxtries maximum number of retries to get nsim unique models

Details

Each of nsim combined models is obtained by randomly sampling one model from each chromosome according to its posterior probability, and combining the x-variables from the sampled models.

Value

A list of models represented as a matrix (similar to the which matrix returned by bicreg.qtl, whose (i,j) element is TRUE if the ith sampled model contains the jth variable

Author(s)

R.D. Ball, (rod.ball@AT@scionresearch.com)

References

Ball, R. D. 2001: Bayesian methods for QTL mapping based on model selection: approximate analysis using the Bayesian Information Criterion. Genetics 159: 1351–1364.

See Also

bicreg.qtl,bicreg.models

Examples

  ## Not run: mWhich200 <- sample.bicreg.qtl.models(chrom.fits,nsim=200)

[Package BayesQTLBIC version 1.0-0 Index]