DirichletProcessGaussian function

Create a Dirichlet Mixture of Gaussians

Create a Dirichlet Mixture of Gaussians

This is the constructor function to produce a dirichletprocess object with a Gaussian mixture kernel with unknown mean and variance. The base measure is a Normal Inverse Gamma distribution that is conjugate to the posterior distribution.

DirichletProcessGaussian(y, g0Priors = c(0, 1, 1, 1), alphaPriors = c(2, 4))

Arguments

  • y: Data
  • g0Priors: Base Distribution Priors γ=(μ0,k0,α0,β0)\gamma = (\mu _0, k_0 , \alpha _0 , \beta _0)
  • alphaPriors: Alpha prior parameters. See UpdateAlpha.

Returns

Dirichlet process object

Details

G0(θγ)=N(μμ0,σ2k0)InvGamma(σ2α0,β0)G_0(\theta | \gamma) = N \left(\mu | \mu_0, \frac{\sigma^2}{k_0} \right) \mathrm{Inv-Gamma} \left(\sigma^2 | \alpha_0, \beta_0 \right)

We recommend scaling your data to zero mean and unit variance for quicker convergence.