Changeset 964 for library/bdm/stat/emix.cpp

Timestamp:

05/21/10 00:44:04 (14 years ago)

Author:

smidl

Message:

Corrections in ARX and PF

Files:

• library/bdm/stat/emix.cpp (modified) (1 diff)

library/bdm/stat/emix.cpp

@@ -53 +53 @@
         //non-central moment
         vec mom2 = zeros ( dim );
+        vec mom1 = zeros ( dim );
+        
         for ( int i = 0; i < w.length(); i++ ) {
-                mom2 += w ( i ) * ( component( i )->variance() + pow ( component ( i )->mean(), 2 ) );
+                vec vi=component( i )->variance();
+                vec mi=component ( i )->mean();
+                mom2 += w ( i ) * ( vi + pow ( mi, 2 ) );
+                mom1 += w(i) * mi;
         }
         //central moment
-        return mom2 - pow ( mean(), 2 );
+        return mom2 - pow ( mom1, 2 );
 }