Changeset 1336 for applications/robust/main.cpp

Timestamp:

04/20/11 20:11:40 (13 years ago)

Author:

sindj

Message:

Sampling dokoncen, zbyva logaritmovat normalizacni faktor a testovat. JS

Files:

• applications/robust/main.cpp (modified) (3 diffs)

applications/robust/main.cpp

@@ -134 +134 @@
         vector<vector<string>> strings;
 
-        char* file_strings[3] = {"c:\\ar_cauchy_single.txt", "c:\\ar_normal_single.txt","c:\\ar_student_single.txt"};
+        char* file_strings[3] = {"c:\\ar_student_single.txt", "c:\\ar_student_single.txt","c:\\ar_cauchy_single.txt"};
 
         for(int i = 0;i<3;i++)
@@ -162 +162 @@
                 vector<vec> conditions;
                 //emlig* emliga = new emlig(2);
-                RARX* my_rarx = new RARX(2,70);
+                RARX* my_rarx = new RARX(3,20,true);
 
                 for(int k = 1;k<170;k++)
@@ -189 +189 @@
                                         
                                 
-                                if(k>5)
+                                if(k>10)
                                 {
                                         //my_rarx->posterior->step_me(0);
 
-                                        my_rarx->posterior->sample_mat(1);
-
+                                        mat samples = my_rarx->posterior->sample_mat(50);
+                                        
+                                        
                                         cout << "MaxLik coords:" << my_rarx->posterior->minimal_vertex->get_coordinates() << endl;
 
-                                        ofstream myfile;
-                                        char fstring[80];
-                                        strcpy(fstring,file_strings[j]);
-                                        strcat(fstring,"_res.txt");
-
-                                        myfile.open(fstring,ios::app);
-                                        myfile << my_rarx->posterior->minimal_vertex->get_coordinates()[0];
-                                        if(k!=strings[j].size()-1)
+                                        for(int s = 0;s<samples.rows();s++)
                                         {
-                                                myfile << ",";
+                                                
+                                                double avg_parameter = samples.get_row(s)*ones(samples.cols())/samples.cols();                                          
+                                                
+                                                ofstream myfile;
+                                                char fstring[80];
+                                                strcpy(fstring,file_strings[j]);
+
+                                                char es[5];
+                                                strcat(fstring,itoa(s,es,10));
+
+                                                strcat(fstring,"_res.txt");
+                                                
+
+                                                myfile.open(fstring,ios::app);
+                                                
+                                                //myfile << my_rarx->posterior->minimal_vertex->get_coordinates()[0];
+                                                myfile << avg_parameter;
+                                                
+                                                if(k!=strings[j].size()-1)
+                                                {
+                                                        myfile << ",";
+                                                }
+                                                else
+                                                {
+                                                        myfile << endl;
+                                                }
+                                                myfile.close();
                                         }
-                                        else
-                                        {
-                                                myfile << endl;
-                                        }
-                                        myfile.close();
                                 }                                       
                         }