README

Please see the Example.R for a basic script to run BEM and MCMC.

The input parameters are as follows:

X		Data(N,P)
K		Maximum Number of Components
workdir 	Location of your executables. 
algorithm	Desired Algorithm, must be either "bem" or "mcmc"
relabel		Do you want MCMC to relabel based on a classification vector?
m0		Prior mean on component means. mu_j~N(m0,gamma*Sigma_j)
nu		Prior degrees of freedom for IW. Sigma_j~IW(nu,nu*eta_j*Phi0)
phi0		Diagonal Elements of Phi0
gamma		See m0
e,f		alpha ~ gamma(e,f)
iter		Number of iterations for MCMC and maximum number of iterations for BEM
seed		See for random numbers
burnin		Number of burnin iterations for MCMC
printout	How often the program prints it's status. -1 Suppresses all output.
percenterr	Percent Error Tolerance in likelihodd for BEM
mu		Starting Values for Component Means (KxP)
Sigma		Starting Values for Component Covariances (PxPxK Array)
pi		Starting Values for Component Weights (K Vector)
Z		Reference classification vector for relabeling
gpu		If true, a GPU binrary will be called
gpuChunk	Chunksize for GPU
gpuNumDev	Number of GPUs in multi GPU MCMC
gpuDevice	Which GPU to use (or start with in multi GPU MCMC)