Method sliceStatistics

src/topp/MapStatistics.cpp:106–145

  }

  vector<double> sliceStatistics(const FeatureMap& map, Size begin, Size end) const
  {
    // If we are asked to produce stats for an empty set, return an empty vector.
    if (end <= begin || end > map.size())
    {
      return vector<double>(43);
    }

    Size size = end - begin;
    vector<double> intensities(size);
    vector<double> peak_widths(size);
    vector<double> mz(size);
    vector<double> overall_qualities(size);
    vector<double> mz_qualities(size);
    vector<double> rt_qualities(size);
    double tic = 0.0;

    for (Size i = begin; i < end; ++i)
    {
      intensities[i - begin]  = map[i].getIntensity();
      mz[i - begin]                       = map[i].getMZ();
      peak_widths[i - begin]  = map[i].getWidth();
      rt_qualities[i - begin] = map[i].getQuality(Feature::RT);
      mz_qualities[i - begin] = map[i].getQuality(Feature::MZ);
      overall_qualities[i - begin] = map[i].getOverallQuality();
      tic += map[i].getIntensity();
    }

    vector<double> results;
    results.reserve(43); // 6 7-number stats + tic
    results.push_back(tic);
    results << Math::SummaryStatistics< vector<double> >(intensities);
    results << Math::SummaryStatistics< vector<double> >(mz);
    results << Math::SummaryStatistics< vector<double> >(peak_widths);
    results << Math::SummaryStatistics< vector<double> >(overall_qualities);
    results << Math::SummaryStatistics< vector<double> >(rt_qualities);
    results << Math::SummaryStatistics< vector<double> >(mz_qualities);

    return results;
  }

protected: