Changeset 741 for library/bdm

Show
Ignore:
Timestamp:
11/26/09 22:46:01 (15 years ago)
Author:
smidl
Message:

Stress tests are passing now. Missing validate calls are filled...

Location:
library/bdm
Files:
5 modified

Legend:

Unmodified
Added
Removed

  • library/bdm/base/bdmbase.h

    r738 r741  
    437437        //! Array<vec> version of evallogcond 
    438438        virtual vec evallogcond_mat ( const Array<vec> &Yt, const vec &cond ) { 
    439                 bdm_error ( "Not implemented" ); 
    440                 return vec(); 
     439                vec v ( Yt.length() ); 
     440                for ( int i = 0; i < Yt.length(); i++ ) { 
     441                        v ( i ) = evallogcond ( Yt( i ), cond ); 
     442                } 
     443                return v; 
    441444        } 
    442445 

  • library/bdm/estim/arx.cpp

    r738 r741  
    1212 
    1313        dyad.set_subvector ( 0, yt ); 
    14         dyad.set_subvector ( dimy, cond ); 
     14        if (cond.length()>0) 
     15                dyad.set_subvector ( dimy, cond ); 
    1516        // possible "1" is there from the beginning 
    1617 

  • library/bdm/estim/mixtures.h

    r737 r741  
    7777                        Coms ( i ) = ( BMEF* ) Coms0 ( i )->_copy_(); 
    7878                } 
     79                weights.set_parameters(alpha0); 
     80                weights.validate(); 
    7981                build_est(); 
    8082        } 

  • library/bdm/stat/exp_family.h

    r739 r741  
    518518        void set_parameters ( const vec &beta0 ) { 
    519519                est.set_parameters ( beta0 ); 
     520                est.validate(); 
    520521                if ( evalll ) { 
    521522                        last_lognc = est.lognc(); 

  • library/bdm/stat/merger.cpp

    r739 r741  
    204204        ARX A0; 
    205205        A0.set_statistics ( dim, V0 ); //initial guess of Mix: 
     206        A0.validate(); 
    206207 
    207208        Mix.init ( &A0, Smp_ex, Ncoms );