source: trunk/src/Statistics.cc @ 117

// $Id: Statistics.cc 117 2004-07-19 15:09:14Z peter $

// Thep C++ Tools
#include "Statistics.h"

// System includes 
#include <vector>
#include <algorithm>

#include <cstdlib>
#include <iostream>


namespace theplu {
namespace cpptools {  

  Statistics::Statistics(void) 
      
  {
  }

  double Statistics::median(std::vector<double>& vec) 
  {
    return percentile(vec, 50.0);
  }

  double Statistics::median(std::vector<size_t>& vec) 
  {
    return percentile(vec, 50.0);
  }

  double Statistics::percentile(std::vector<double>& vec, double percentile) 
  {
    sort(vec.begin(), vec.end());
    double j = percentile/100 * vec.size();
    if (percentile==100)
      return vec[vec.size()-1];   
    
    double k = j + 0.5 - static_cast<int> (j+0.5);
    if (k<.500000000001 && k>.499999999999){
      int i = static_cast<int>(j+0.5);
      return (vec[i]+vec[i-1])/2;
    }
    else{
      int i = static_cast<int>(j);
      return  vec[i];
    }
    
  }

  double Statistics::percentile(std::vector<size_t>& vec, double percentile) 
  {
    sort(vec.begin(), vec.end());
    double j = percentile/100 * vec.size();
    if (percentile==100)
      return vec[vec.size()-1];   
    
    double k = j + 0.5 - static_cast<int> (j+0.5);
    if (k<.500000000001 && k>.499999999999){
      int i = static_cast<int>(j+0.5);
      double r = static_cast<double>(vec[i]+vec[i-1])/2;
      return r;
      
    }
    else{
      int i = static_cast<int>(j);
      return  static_cast<double>(vec[i]);
    }
    
  }


}} // of namespace cpptools and namespace theplu