Learn R Programming

mixAK (version 2.2)

NMixMCMCwrapper: Wrapper to the NMixMCMC main simulation.

Description

This is wrapper to the NMixMCMC main simulation which allows vectorized evaluation and possibly parallel computation.

THIS FUNCTION IS NOT TO BE CALLED BY ORDINARY USERS.

Usage

NMixMCMCwrapper(chain=1,
                scale, prior, inits, Cpar, RJMCMC, CRJMCMC,
                actionAll, nMCMC, keep.chains, PED,
                dens.zero)

Arguments

chain
identification of the chain sampled in a particular call of this function, usually number like 1, 2, ...
Cpar
a list with the following components [object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object]
scale
a list specifying how to scale the data before running MCMC. See argument scale in NMixMCMC
prior
a list specifying prior hyperparameters. See argument prior in NMixMCMC.
inits
a list of length at least chain. Its chain-th component is used. Each component of the list should have the structure of init argument of function NMixMCMC
RJMCMC
a list specifying parameters for RJ-MCMC. See argument RJMCMC in NMixMCMC
CRJMCMC
a numeric vector with parameters for RJ-MCMC.
actionAll
argument for underlying C++ function.
nMCMC
vector giving the length of MCMC etc.
keep.chains
logical. If FALSE, only summary statistics are returned in the resulting object. This might be useful in the model searching step to save some memory.
PED
a logical value which indicates whether the penalized expected deviance (see Plummer, 2008 for more details) will be computed (which requires two parallel chains). Even if keep.chains is FALSE, it is necessary to
dens.zero
small number (1e-300) to determine whether the contribution to the deviance ($-\log$ density) is equal to infinity. Such values are trimmed when computing expected deviance.

Value

  • A list having almost the same components as object returned by NMixMCMC function.

See Also

NMixMCMC.