Calculate mean and standard deviation from a vector of samples in C++ using Boost

by

I don’t know if Boost has more specific functions, but you can do it with the standard library.

Given std::vector<double> v, this is the naive way:

#include <numeric>

double sum = std::accumulate(v.begin(), v.end(), 0.0);
double mean = sum / v.size();

double sq_sum = std::inner_product(v.begin(), v.end(), v.begin(), 0.0);
double stdev = std::sqrt(sq_sum / v.size() - mean * mean);

This is susceptible to overflow or underflow for huge or tiny values. A slightly better way to calculate the standard deviation is:

double sum = std::accumulate(v.begin(), v.end(), 0.0);
double mean = sum / v.size();

std::vector<double> diff(v.size());
std::transform(v.begin(), v.end(), diff.begin(),
               std::bind2nd(std::minus<double>(), mean));
double sq_sum = std::inner_product(diff.begin(), diff.end(), diff.begin(), 0.0);
double stdev = std::sqrt(sq_sum / v.size());

UPDATE for C++11:

The call to std::transform can be written using a lambda function instead of std::minus and std::bind2nd(now deprecated):

std::transform(v.begin(), v.end(), diff.begin(), [mean](double x) { return x - mean; });

More Related Contents:

  1. Sorting zipped (locked) containers in C++ using boost or the STL
  2. Magic number in boost::hash_combine
  3. Cut set of a graph, Boost Graph Library
  4. How to find number of words in a phrase with spaces removed by checking v. dictionary
  5. Fast exact bigint factorial
  6. Generating combinations in c++
  7. How to link C++ program with Boost using CMake
  8. Finding duplicates in O(n) time and O(1) space
  9. Cross-platform way to get line number of an INI file where given option was found
  10. Fast String Hashing Algorithm with low collision rates with 32 bit integer [closed]
  11. Emulate “double” using 2 “float”s
  12. How can I compose output streams, so output goes multiple places at once?
  13. How to link to dynamic boost libs?
  14. Is there a way to test whether a C++ class has a default constructor (other than compiler-provided type traits)?
  15. CMake is not finding Boost
  16. Searching in a sorted and rotated array
  17. Best way for interprocess communication in C++
  18. `std::list::sort()` – why the sudden switch to top-down strategy?
  19. Rolling median algorithm in C
  20. Easiest way of using min priority queue with key update in C++
  21. C++ Boost ASIO simple periodic timer?
  22. Creating Library with backward compatible ABI that uses Boost
  23. Efficiently reading a very large text file in C++
  24. Edges on polygon outlines not always correct
  25. Creating JSON arrays in Boost using Property Trees
  26. version `CXXABI_1.3.8′ not found (required by …) [duplicate]
  27. C++ range/xrange equivalent in STL or boost?
  28. How to find Longest Common Substring using C++
  29. Generate N random numbers within a range with a constant sum
  30. How do I “normalize” a pathname using boost::filesystem?

Leave a Comment