Learn R Programming

bayesm (version 3.0-2)

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)

Arguments

y

data array - rows are obs

Bbar

prior mean for mean vector of each norm comp

A

prior precision parameter

nu

prior d.f. parm

V

prior location matrix for covariance priro

a

Dirichlet prior parms

p

prior prob of each mixture component

z

component identities for each observation -- "indicators"

Value

a list containing:

p

draw mixture probabilities

z

draw of indicators of each component

comps

new 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://www.perossi.org/home/bsm-1

See Also

rnmixGibbs