Changeset 1383 for library/bdm/estim/arx.h

Timestamp:

08/18/11 16:02:53 (13 years ago)

Author:

sindj

Message:

Program update, kontrola vypoctu provedena, integruje spravne v nedegenerovanem pripade. JS

Files:

• library/bdm/estim/arx.h (modified) (1 diff)

library/bdm/estim/arx.h

@@ -311 +311 @@
         
         void bayes ( const vec &val, const vec &cond ) {
-                // fill window
-                Y.append_col(val);
-                Cond.append_col(cond);          
-                if (Y.cols()>win_length){
-                        // shift the buffer
-                        Y=Y.get_cols(1,Y.cols()-1);
-                        Cond=Cond.get_cols(1,Cond.cols()-1);
+                
+                if(cond.size()>0)
+                {
+                        // fill window
+                        Y.append_col(val);
+                        Cond.append_col(cond);          
+                        if (Y.cols()>win_length){
+                                // shift the buffer
+                                Y=Y.get_cols(1,Y.cols()-1);
+                                Cond=Cond.get_cols(1,Cond.cols()-1);
+                        }
+                        
+                        est._V()=V0;
+                        est._nu()=nu0;
+                        for ( int t = 0; t < Y.cols(); t++ ) {
+                                ARX::bayes ( Y.get_col ( t ), Cond.get_col ( t ) );
+                        }
+                }
+                else
+                {
+                        Y.append_col(val);
+
+                        if (Y.cols()>win_length){
+                                // shift the buffer
+                                Y=Y.get_cols(1,Y.cols()-1);                             
+                        }
+
+                        est._V()=V0;
+                        est._nu()=nu0;
+
+                        for ( int t = 0; t < Y.cols(); t++ ) {
+                                ARX::bayes (Y.get_col ( t ));
+                        }
                 }
                 
-                est._V()=V0;
-                est._nu()=nu0;
-                for ( int t = 0; t < Y.cols(); t++ ) {
-                        ARX::bayes ( Y.get_col ( t ), Cond.get_col ( t ) );
-                }
-                
-        }
+        }
+
         void from_setting(const Setting &set){
                 ARX::from_setting(set);