statpars.h

Go to the documentation of this file.

 
#pragma once
 
#if __has_include("defs.h")
#include "defs.h"
#endif
 
#include <deque>
 
namespace mlib {
 
class statpars
{
public:
  statpars (int nmax = 0);
  statpars (std::vector<double> vec);
 
  void add (double val);
  void add (double* vals, int count);
  void add (std::vector<double> v);
  void clear ();
  int count ();
  double average ();
  double stdev ();
  double mad ();
  double variance ();
  double skewness ();
  double kurtosis ();
 
private:
  void calculate ();
  std::deque<double> values;
  int nmax;
  bool calc;
  double sum, adev, var, sdev, skew, kurt;
};
 
inline void statpars::add (double val)
{
  if (nmax && values.size () == nmax)
  {
    sum -= values.front ();
    values.pop_front ();
  }
  values.push_back (val);
  sum += val;
  calc = false;
}
 
inline void statpars::add (double* vals, int count)
{
  for (int i = 0; i < count; i++)
    add (vals[i]);
}
 
inline void statpars::add (std::vector<double> vals)
{
  for (double v : vals)
    add (v);
}
 
inline double statpars::average ()
{
  if (!values.empty ())
    return sum / values.size ();
  else
    return 0.;
}
 
inline double statpars::variance ()
{
  if (!calc)
    calculate ();
  return var;
}
 
inline double statpars::stdev ()
{
  if (!calc)
    calculate ();
  return sdev;
}
 
inline double statpars::mad ()
{
  if (!calc)
    calculate ();
  return adev;
}
 
inline double statpars::skewness ()
{
  if (!calc)
    calculate ();
  return skew;
}
 
inline double statpars::kurtosis ()
{
  if (!calc)
    calculate ();
  return kurt;
}
 
inline int statpars::count ()
{
  return (int)values.size ();
}
 
inline void statpars::clear ()
{
  values.clear ();
  calc = false;
  sum = 0.;
}
 
} // namespace mlib