Changeset 900 for library/bdm/stat

Timestamp:

04/11/10 22:05:11 (14 years ago)

Author:

smidl

Message:

particle bug fixing

Location:

library/bdm/stat

Files:

• emix.cpp (modified) (1 diff)
• exp_family.h (modified) (1 diff)

library/bdm/stat/emix.cpp

@@ -244 +244 @@
         vec tmp ( dim ); //second moment
         for ( int i = 0; i < no_factors(); i++ ) {
-                vec pom = factor ( i )->mean();
+                vec pom = factor ( i )->variance();
                 dls ( i )->pushup ( tmp, pow ( pom, 2 ) );
         }
-        return tmp - pow ( mean(), 2 );
+        return tmp;
 }
 vec eprod_base::sample() const {

library/bdm/stat/exp_family.h

@@ -1879 +1879 @@
                 double evallog (const vec &dt) const {return -inf;}
                 vec mean () const {return point;}
-                vec variance () const {return pow(point,2);}
+                vec variance () const {return zeros(point.length());}
                 void qbounds ( vec &lb, vec &ub, double percentage = 0.95 ) const { lb = point; ub = point;}
                 //! access