diff --git a/eo/test/CMakeLists.txt b/eo/test/CMakeLists.txt
index 2c0772de..4f194681 100644
--- a/eo/test/CMakeLists.txt
+++ b/eo/test/CMakeLists.txt
@@ -25,8 +25,6 @@ LINK_DIRECTORIES(${EO_BINARY_DIR}/lib)
 ######################################################################################
 
 SET (TEST_LIST
-  t-eoParetoFitness
-  t-eoPareto
   t-eofitness
   t-eoRandom
   t-eobin
@@ -65,8 +63,8 @@ SET (TEST_LIST
   t-eoSyncEasyPSO
   t-eoOrderXover
   t-eoExtendedVelocity
-  # t-eoFrontSorter
   t-eoLogger
+  t-eoIQRStat
   )
 
 
diff --git a/eo/test/t-eoIQRStat.cpp b/eo/test/t-eoIQRStat.cpp
new file mode 100644
index 00000000..60bcff8f
--- /dev/null
+++ b/eo/test/t-eoIQRStat.cpp
@@ -0,0 +1,69 @@
+
+#include <eo>
+#include <es.h>
+#include <utils/eoStat.h>
+
+#include "real_value.h"
+
+typedef eoReal<eoMinimizingFitness> realVec;
+
+double test( eoPop<realVec>& pop, double target_value )
+{
+    eoEvalFuncPtr<realVec, double, const std::vector<double>&> eval( real_value );
+
+    eoPopLoopEval<realVec> pop_eval(eval);
+
+    pop_eval(pop,pop);
+
+    eoInterquartileRangeStat<realVec> iqr_stat(0.0, "IQR");
+
+    iqr_stat( pop );
+
+    std::cout << iqr_stat.longName() << "=" << iqr_stat.value() << " should be " << target_value << std::endl;
+
+    return iqr_stat.value();
+}
+
+int main()
+{
+    eoPop<realVec> pop;
+
+    // fixed test
+    realVec sol1(2,-1);
+    realVec sol2(2,-1);
+    realVec sol3(2,1);
+    realVec sol4(2,1);
+    pop.push_back( sol1 );
+    pop.push_back( sol2 );
+    pop.push_back( sol3 );
+    pop.push_back( sol4 );
+    // on the sphere function everyone has the same fitness of 1
+    if( test(pop, 0) != 0 ) {
+        exit(1);
+    }
+
+    pop.erase(pop.begin(),pop.end());
+
+    // fixed test
+    sol1 = realVec(2,0);
+    sol2 = realVec(2,0);
+    sol3 = realVec(2,1);
+    sol4 = realVec(2,1);
+    pop.push_back( sol1 );
+    pop.push_back( sol2 );
+    pop.push_back( sol3 );
+    pop.push_back( sol4 );
+    if( test(pop, 1) != 1 ) {
+        exit(1);
+    }
+
+    // test on a random normal distribution
+    eoNormalGenerator<double> normal(1,rng);
+    eoInitFixedLength<realVec> init_N(2, normal);
+    pop = eoPop<realVec>( 1000000, init_N );
+    double iqr = test(pop, 1.09);
+    if( iqr < 1.08 || iqr > 1.11 ) {
+        exit(1);
+    }
+}
+