bmapqtl.options {bim} | R Documentation |
This function allows the users to set and examine the options and priors for reversible jump MCMC QTL model search. The options are stored in global hidden variable .bmapqtl.options.
bmapqtl.options(...,reset=FALSE)
reset |
reset options to default settings |
... |
any options can be defined, using `name = value'. The
available options are listed below.
|
Initial mean and variance init
are used directly if the
variance is positive. For values c(a,-b)
use mean 0 and
variance .5s^2.
Prior arguments have coded meaning as follows, with a
and
b
non-negative. For the grand mean, c(a,b)
signifies
Normal(a,b)
, while c(a,-b)
stands for
Normal
(bar{Y},bS^2), with bar{Y} and
S^2 the phenotype mean and variance, respectively. Additive and
dominance priors coded as c(a,b)
have Normal(a,b)
; as
c(a,-b)
have Normal
(0,aS^2); as c(-a,-b)
have
Normal
(0,β S^2) with βsimBeta(2,10)
if a=b=0
or βsimBeta(a,b)
otherwise. Variance
has prior inverse Gamma: IG(a,b)
for c(a,b)
, or
IG(a,\eqn{s^2})
for c(a,-b)
.
A list (in any case) with the previous values of the options changed, or all options when no arguments were given.
Hao Wu, hao@jax.org
bmapqtl.options() bmapqtl.options(niter=100000, by=1000)