Fit a mixture of Gaussians to the observed data.

gibbsGMM(y, niter = 1000, nburn = 500, priors = NULL)

Arguments

y

A vector of observed pixel data.

niter

The number of iterations of the algorithm to perform.

nburn

The number of iterations to discard as burn-in.

priors

A list of priors for the parameters of the model.

Value

A matrix containing MCMC samples for the parameters of the mixture model.