Changeset 1014 for applications/bdmtoolbox/tutorial/userguide/mixef_basic.m

Timestamp:

05/27/10 23:07:57 (14 years ago)

Author:

smidl

Message:

Mixtures example

Files:

• applications/bdmtoolbox/tutorial/userguide/mixef_basic.m (modified) (4 diffs)

applications/bdmtoolbox/tutorial/userguide/mixef_basic.m

@@ -5 +5 @@
 com.rgr = RV({},[],[]);
 
-Data = [randn(2,50) randn(2,50)+10*ones(2,50)];
+Data = [randn(2,50) randn(2,50)+[5*ones(1,50); 10*ones(1,50)]];
 
-[Mix0,P0]=mixef_init(Data,com);
+[Mix0,P0]=mixef_init(Data,com,4);
 % show predictor
 
@@ -16 +16 @@
 hold on
 epdf_2dplot(Pred);
+hold on
+for i=1:length(Pred.pdfs)
+    mea=epdf_mean(Pred.pdfs{i});
+    plot(mea(1),mea(2),'+r');
+end
 
 % Do Quasi Bayes
@@ -32 +37 @@
     % normalize weights
     w = exp(log_w_nn-max(log_w_nn));
-    w = w/sum(w)
+    w = w/sum(w);
     
     for c=1:n
@@ -56 +61 @@
 epdf_2dplot(PredQB);
 
+PredQBc = bm_epredictor(MixQBc);
+figure(3);
+hold off
+plot(Data(1,:), Data(2,:),'.');
+hold on
+epdf_2dplot(PredQBc);
+
 %% relations between components
 Batta_dist=zeros(n);