rmixGibbs: Gibbs Sampler for Normal Mixtures w/o Error Checking
Description
rmixGibbs
makes one draw using the Gibbs Sampler for a mixture of multivariate normals.Usage
rmixGibbs(y, Bbar, A, nu, V, a, p, z, comps)
Arguments
y
data array - rows are obs
Bbar
prior mean for mean vector of each norm comp
A
prior precision parameter
V
prior location matrix for covariance priro
p
prior prob of each mixture component
z
component identities for each observation -- "indicators"
comps
list of components for the normal mixture
Value
- a list containing:
- pdraw mixture probilities
- zdraw of indicators of each component
- compsnew draw of normal component parameters
Warning
This routine is a utility routine that does not check the
input arguments for proper dimensions and type.Details
rmixGibbs
is not designed to be called directly. Instead, use rnmixGibbs
wrapper function.References
For further discussion, see Bayesian Statistics and Marketing
by Allenby, McCulloch, and Rossi, Chapter 5.
http://gsbwww.uchicago.edu/fac/peter.rossi/research/bsm.html