 0; % graph of a single node with one localCPD
CPD         = tabularCpdCreate(mix0');
localCPD    = condGaussCpdCreate(mu0, Sigma0, 'prior', mixGauss.cpd.prior); 
mixGaussDgm = dgmCreate(G, CPD, 'localCPDs', localCPD);
%% fit using dgmTrain
fprintf('DGM\n'); 
mixGaussDgm = dgmTrain(mixGaussDgm, 'localev', localEv, 'verbose', true);
%% compare results
assert(approxeq(mixGaussDgm.localCPDs{1}.mu, mixGauss.cpd.mu)); 
assert(approxeq(mixGaussDgm.localCPDs{1}.Sigma, mixGauss.cpd.Sigma)); 
assert(approxeq(mixGaus