ik: -3448.96
2	 loglik: -1884.58
3	 loglik: -1746.68
4	 loglik: -1580.99
5	 loglik: -1490.54
6	 loglik: -1473.37
7	 loglik: -1472.62
8	 loglik: -1472.51
</pre><h2>create the initial dgm<a name="5"></a></h2><pre class="codeinput">G           = 0; <span class="comment">% graph of a single node with one localCPD</span>
CPD         = tabularCpdCreate(mix0');
localCPD    = condGaussCpdCreate(mu0, Sigma0, <span class="string">'prior'</span>, mixGauss.cpd.prior);
mixGaussDgm = dgmCreate(G, CPD, <span class="string