Function get_iqr

src/myFunc.cpp:226–265

Calculate inter quartile range across individuals ---------------------------------------------------------------------------

// Calculate inter quartile range across individuals
// ---------------------------------------------------------------------------
float get_iqr(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 IQR
  float lower_quartile, upper_quartile;

  if( ndatapoints % 2 == 1){
    float fl = (ndatapoints-1)/float(4.0);
    float fu = 3*(ndatapoints-1)/float(4.0);

    if( fmod(fl,float(1.0)) == 0 ){
      lower_quartile = data[ (int) fl ];
      upper_quartile = data[ (int) fu ];
    }
    else{
      lower_quartile = (data[ (int) fl ] + data[ (int) fl + 1 ])/float(2.0);
      upper_quartile = (data[ (int) fu ] + data[ (int) fu + 1 ])/float(2.0);
    }
  }
  else{
    double fl = (ndatapoints)/float(4.0);
    double fu = 3*(ndatapoints)/float(4.0);

    if( fmod(fl,1.0) == 0 ){
      lower_quartile = data[ (int) fl ];
      upper_quartile = data[ (int) fu ];
    }
    else{
      lower_quartile = (data[ (int) fl ] + data[ (int) fl + 1 ])/float(2);
      upper_quartile = (data[ (int) fu ] + data[ (int) fu + 1 ])/float(2);
    }
  }

  return upper_quartile - lower_quartile;
}

void readChrNamesInBed(const std::string &targetBed, std::vector<std::string>&names_bed){
    ifstream file(targetBed.c_str());