Function get_median

src/myFunc.cpp:123–135

Calculate median across individuals ---------------------------------------------------------------------------

// Calculate median across individuals
// ---------------------------------------------------------------------------
float get_median(const std::vector<float> & myvector) {
  vector<float> data (myvector);
  int ndatapoints = data.size();
  if (ndatapoints==0) {
    cerr << "Error: zero values to calculate medians..\n";
    exit(-1);
  }
  sort(data.begin(),data.end());
  // Get median
  float median_value = ndatapoints % 2 == 1 ? data[(ndatapoints-1)/2] : (data[ndatapoints/2 - 1] + data[ndatapoints/2])/float(2.0);
  data.clear();
  return median_value;
}

float get_medianNotNA(const std::vector<float> & myvector) {
  vector<float> data (myvector);