Changeset 1282 for applications/robust/main.cpp

Timestamp:

02/21/11 09:22:08 (15 years ago)

Author:

sindj

Message:

Prvni pokusy s robustnim AR modelem. JS

Files:

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

applications/robust/main.cpp

@@ -9 +9 @@
 #include "robustlib.h"
 #include <vector>
+#include <iostream>
+#include <fstream>
+
 using namespace itpp;
-
 //using namespace bdm;
 
@@ -18 +20 @@
 int main ( int argc, char* argv[] ) {
         
-
-        emlig* emlig1 = new emlig(emlig_size);
+        
+        
+        
+        vector<vector<vector<string>>> string_lists;
+        string_lists.push_back(vector<vector<string>>());
+        string_lists.push_back(vector<vector<string>>());
+        string_lists.push_back(vector<vector<string>>());
+
+        char* file_strings[3] = {"c:\\ar_normal.txt", "c:\\ar_student.txt", "c:\\ar_cauchy.txt"};
+        
+
+        for(int i = 0;i<3;i++)
+        {       
+                ifstream myfile(file_strings[i]);
+                if (myfile.is_open())
+                {
+                        while ( myfile.good() )
+                        {
+                                string line;
+                                getline(myfile,line);
+                                
+                                vector<string> parsed_line;
+                                while(line.find(',') != string::npos)
+                                {
+                                        int loc = line.find(',');
+                                        parsed_line.push_back(line.substr(0,loc));
+                                        line.erase(0,loc+1);                                    
+                                }                               
+
+                                string_lists[i].push_back(parsed_line);
+                        }
+                        myfile.close();
+                }
+        }
+
+        for(int j = 0;j<string_lists.size();j++)
+        { 
+                
+                for(int i = 0;i<string_lists[j].size();i++)
+                {
+                        vector<vec> conditions;
+                        emlig* emliga = new emlig(2);
+                        for(int k = 1;k<string_lists[j][i].size();k++)
+                        {
+                                vec condition;
+                                //condition.ins(0,1);                           
+                                condition.ins(0,string_lists[j][i][k]);                         
+                                conditions.push_back(condition);
+
+                                //cout << "orig:" << condition << endl;
+
+                                if(conditions.size()>1)
+                                {               
+                                        conditions[k-2].ins(0,string_lists[j][i][k]);
+                                        
+                                }
+
+                                if(conditions.size()>2)
+                                {
+                                        conditions[k-3].ins(0,string_lists[j][i][k]);
+
+                                        //cout << "modi:" << conditions[k-3] << endl;
+
+                                        emliga->add_condition(conditions[k-3]);
+
+                                        /*
+                                        if(k>5)
+                                        {
+                                                cout << "MaxLik coords:" << emliga->minimal_vertex->get_coordinates() << endl;
+                                        }
+                                        */
+                                }
+                        }
+
+                        //emliga->step_me(0);
+                        cout << "MaxLik coords:" << emliga->minimal_vertex->get_coordinates() << endl;
+                        cout << "Step: " << i << endl;
+                }
+
+                cout << "One experiment finished." << endl;
+        }
+
+        
+
+        //emlig* emlig1 = new emlig(emlig_size);
         
         /*
@@ -39 +124 @@
         }*/
         
-        vec condition1 = "1.0 1.0 1.0";
+    /*
+        vec condition1 = "1.0 1.0 1.01";
         emlig1->add_condition(condition1);
 
@@ -45 +131 @@
         emlig1->add_condition(condition2);
 
-        vec condition3 = "0.5 -1.1 1.0";
+        vec condition3 = "0.5 -1.01 1.0";
         emlig1->add_condition(condition3);
 
@@ -51 +137 @@
         emlig1->add_condition(condition4);
 
-        vec condition5 = "-0.3 1.7 1.5";
-        emlig1->add_condition(condition5);
-
+        cout << emlig1->minimal_vertex->get_coordinates();
+
+        //vec condition5 = "-0.3 1.7 1.5";
+        //emlig1->add_condition(condition5);
+
+        // DA SE POUZIT PRO VYPIS DO SOUBORU
         // emlig1->step_me(0);
 
@@ -61 +150 @@
 
         
-
         /*
-        for(int i = 0;i<500;i++)
-        {
-                cout << endl << "Step:" << i << endl;
-                
-                
+        
+        for(int i = 0;i<100;i++)
+        {
+                cout << endl << "Step:" << i << endl;           
 
                 double condition[emlig_size+1];         
@@ -139 +226 @@
 }
 
+
+