diff --git a/mo/src/continuator/moFitnessVarianceStat.h b/mo/src/continuator/moFitnessVarianceStat.h
new file mode 100644
index 000000000..3769c5983
--- /dev/null
+++ b/mo/src/continuator/moFitnessVarianceStat.h
@@ -0,0 +1,115 @@
+/*
+  <moBestFitnessStat.h>
+  Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
+
+  Sébastien Verel, Arnaud Liefooghe, Jérémie Humeau
+
+  This software is governed by the CeCILL license under French law and
+  abiding by the rules of distribution of free software.  You can  use,
+  modify and/ or redistribute the software under the terms of the CeCILL
+  license as circulated by CEA, CNRS and INRIA at the following URL
+  "http://www.cecill.info".
+
+  As a counterpart to the access to the source code and  rights to copy,
+  modify and redistribute granted by the license, users are provided only
+  with a limited warranty  and the software's author,  the holder of the
+  economic rights,  and the successive licensors  have only  limited liability.
+
+  In this respect, the user's attention is drawn to the risks associated
+  with loading,  using,  modifying and/or developing or reproducing the
+  software by the user in light of its specific status of free software,
+  that may mean  that it is complicated to manipulate,  and  that  also
+  therefore means  that it is reserved for developers  and  experienced
+  professionals having in-depth computer knowledge. Users are therefore
+  encouraged to load and test the software's suitability as regards their
+  requirements in conditions enabling the security of their systems and/or
+  data to be ensured and,  more generally, to use and operate it in the
+  same conditions as regards security.
+  The fact that you are presently reading this means that you have had
+  knowledge of the CeCILL license and that you accept its terms.
+
+  ParadisEO WebSite : http://paradiseo.gforge.inria.fr
+  Contact: paradiseo-help@lists.gforge.inria.fr
+*/
+
+#ifndef moFitnessVarianceStat_h
+#define moFitnessVarianceStat_h
+
+#include <continuator/moStat.h>
+
+/**
+ * Statistic that incrementally computes the variance of the fitness of the solutions during the search
+ */
+template <class EOT>
+class moFitnessVarianceStat : public moStat<EOT, typename EOT::Fitness>
+{
+public :
+    typedef typename EOT::Fitness Fitness;
+    using moStat< EOT , typename EOT::Fitness >::value;
+
+    /**
+     * Default Constructor
+     * @param _reInitSol when true the best so far is reinitialized
+     */
+    moFitnessVarianceStat(bool _reInitSol = true)
+    : moStat<EOT, typename EOT::Fitness>(Fitness(), "var"),
+      reInitSol(_reInitSol), firstTime(true),
+      nbSolutionsEncountered(0), currentAvg(0), currentVar(0)
+    { }
+    
+    /**
+     * Initialization of the best solution on the first one
+     * @param _sol the first solution
+     */
+    virtual void init(EOT & _sol) {
+        if (reInitSol || firstTime)
+        {
+            value() = 0.0;
+            nbSolutionsEncountered = currentAvg = currentVar = 0;
+            firstTime = false;
+        }
+        /*else if (firstTime)
+        {
+          value() = 0.0;
+          firstTime = false;
+        }*/
+    	  operator()(_sol);
+    }
+
+  /**
+   * Update the best solution
+   * @param _sol the current solution
+   */
+  virtual void operator()(EOT & _sol) {
+      ++nbSolutionsEncountered;
+      double x = _sol.fitness();
+      double oldAvg = currentAvg;
+      currentAvg = oldAvg + (x - oldAvg)/nbSolutionsEncountered;
+      if (nbSolutionsEncountered > 1) // <- not really necessary
+      {
+        //value() = (value()/nbSolutionsEncountered + _sol.fitness())/(nbSolutionsEncountered+1);
+        double oldVar = currentVar;
+        currentVar = oldVar + (x - oldAvg) * (x - currentAvg);
+        value() = currentVar/nbSolutionsEncountered;
+      }
+  }
+
+    /**
+     * @return name of the class
+     */
+    virtual std::string className(void) const {
+        return "moFitnessVarianceStat";
+    }
+
+protected:
+    bool reInitSol;
+    bool firstTime;
+    double
+      nbSolutionsEncountered
+    , currentAvg
+    , currentVar // actually the var * n
+    ;
+
+};
+
+#endif // moFitnessVarianceStat_h
diff --git a/mo/test/CMakeLists.txt b/mo/test/CMakeLists.txt
index 6e8846be9..319ffb0c9 100644
--- a/mo/test/CMakeLists.txt
+++ b/mo/test/CMakeLists.txt
@@ -95,7 +95,7 @@ set (TEST_LIST
 		t-moStatistics
 		t-moIndexedVectorTabuList
 #		t-moRndIndexedVectorTabuList
-		t-moStdFitnessNeighborStat
+		t-moFitnessVarianceStat
 		t-moStdDevEstimator
 		t-moDynSpanCoolingSchedule
 		t-moTriki
diff --git a/mo/test/t-moStdFitnessNeighborStat.cpp b/mo/test/t-moFitnessVarianceStat.cpp
similarity index 81%
rename from mo/test/t-moStdFitnessNeighborStat.cpp
rename to mo/test/t-moFitnessVarianceStat.cpp
index d562dffa9..f1082864d 100644
--- a/mo/test/t-moStdFitnessNeighborStat.cpp
+++ b/mo/test/t-moFitnessVarianceStat.cpp
@@ -4,7 +4,7 @@
 
 #include <eo>
 //#include "eoReal.h"
-#include "continuator/moStdFitnessNeighborStat.h"
+#include "continuator/moFitnessVarianceStat.h"
 #include "neighborhood/moRealNeighbor.h"
 #include "neighborhood/moRealNeighborhood.h"
 
@@ -15,8 +15,8 @@ typedef moRealNeighbor< EOT > Neighbor;
 
 int main(int ac, char** av)
 {
-    moNeighborhoodStat<Neighbor> nhStat
-    moStdFitnessNeighborStat<Neighbor> stat(nhStat);
+    //moNeighborhoodStat<Neighbor> nhStat
+    moFitnessVarianceStat<Neighbor> stat();
     EOT solution(2, 5);
     stat(solution);
     //assert(stat.value() == 1);