Run hierarchical MCMC for batch model.
kAn integer value specifying the number of latent classes.
hyperparamsAn object of class 'Hyperparameters' used to specify the hyperparameters of the model.
thetathe means of each component and batch
sigma2the variances of each component and batch
nu.0the shape parameter for sigma2
sigma2.0the rate parameter for sigma2
pimixture probabilities which are assumed to be the same for all batches
mumeans from batches, averaged across batches
tau2variances from batches, weighted by precisions
datathe data for the simulation.
data.meanthe empirical means of the components
data.precthe empirical precisions
zlatent variables
zfreqtable of latent variables
probzn x k matrix of probabilities
logpriorlog likelihood of prior: log(p(sigma2.0)p(nu.0)p(mu))
logliklog likelihood: ∑ p_k Φ(θ_k, σ_k)
mcmc.chainsan object of class 'McmcChains' to store MCMC samples
batcha vector of the different batch numbers
batchElementsa vector labeling from which batch each observation came from
modesthe values of parameters from the iteration which maximizes log likelihood and log prior
label_switchlength-one logical vector indicating whether label-switching occurs (possibly an overfit model)
mcmc.paramsAn object of class 'McmcParams'
is_mendelianinteger vector equal in length to the number of trios
.internal.constraintConstraint on parameters. For internal use only.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.