jannink_mcmc: Draw samples from a MCMC chain per Jannink & Wu 2003

In fboehm/reducedscan: QTL Scans with Simpler Statistical Models in Multiparental Populations

Description

Draw samples from a MCMC chain per Jannink & Wu 2003

Usage

jannink_mcmc(
  initial_values,
  niter = 10000,
  genoprobs,
  trait,
  allelic_number_prior = "poisson",
  poisson_prior_mean = 2
)

Arguments

initial_values	a named list of length 3, where names are 'configuration', 'effects', and 'residual_variance'.
niter	number of iterations
genoprobs	a genoprobs matrix for QTL
trait	a univariate trait affected by the QTL
allelic_number_prior	Takes either of two values, "poisson" or "uniform".
poisson_prior_mean	lamba, ie, mean, for the untruncated poisson.

Value

a list of length two, where each entry is itself a list of results from each iteration

fboehm/reducedscan documentation built on April 30, 2021, 8:31 p.m.