Changeset 104

Timestamp:

05/12/08 17:35:02 (16 years ago)

Author:

smidl

Message:

new mixture model test

Files:

• tests/testSmp.cpp (modified) (3 diffs)

tests/testSmp.cpp

@@ -1 +1 @@
 #include <itpp/itbase.h>
 #include <stat/libEF.h>
+#include <stat/emix.h>
 
 using namespace itpp;
@@ -51 +52 @@
         eG.set_parameters(a,b);
         
-        Smp = eG.sample(N);
+        cout << eG.eval(a)<<endl;
+        Smp = eG.sampleN(N);
 
         vec g_mu = elem_div(a,b);
@@ -64 +66 @@
         Smp=mG.samplecond(mu0,lik,N);
         disp(mu0,pow(mu0,2.0)/k,Smp);
+        
+        cout << "======= EMix ======== " << endl;
+        emix eMix(rv);
+        Array<epdf*> Coms(2);
+        Coms(0) = &eG;
+        Coms(1) = &eN;
+        
+        eMix.set_parameters(vec_2(0.5,0.5), Coms);
+        vec smp = eMix.sample();
+        Smp = eMix.sampleN(N);
+        disp(eMix.mean(),zeros(2),Smp);
+
+        cout << "======= MMix_triv ======== " << endl;
+        mmix_triv mMix(rv,rv,&eMix);
+        
+        Smp = mMix.samplecond(mu0,lik,N);
+        disp(eMix.mean(),zeros(2),Smp);
 
         //Exit program: