Changeset 1687 for trunk/yat

Timestamp:

Dec 30, 2008, 11:00:24 PM (12 years ago)

Author:

Peter

Message:

fixes #455

File:

• trunk/yat/statistics/KolmogorovSmirnov.cc (modified) (1 diff)

trunk/yat/statistics/KolmogorovSmirnov.cc

@@ -143 +143 @@
     double f1=0;
     double f2=0;
-    // Peter, we should take care of ties!
     while(iter!=data_.end()){
-      if (iter->second)
-        f1 += iter->first.second;
-      else
-        f2 += iter->first.second;
-      res.push_back(f1/sum_w1_-f2/sum_w2_);
-      ++iter;
+      size_t count=0;
+      double value = iter->first.first;
+      while (iter!=data_.end() && iter->first.first==value) {
+        if (iter->second) {
+          f1 += iter->first.second;
+        }
+        else {
+          f2 += iter->first.second;
+        }
+        ++count;
+        ++iter;
+      }
+      res.resize(res.size()+count, f1/sum_w1_-f2/sum_w2_);
     }
   }