MultiStatistics.h

 
#pragma once
 
#if defined(MultiStatistics_RECURSES)
#error Recursive header files inclusion detected in MultiStatistics.h
#else // defined(MultiStatistics_RECURSES)
#define MultiStatistics_RECURSES
 
#if !defined MultiStatistics_h
#define MultiStatistics_h
 
// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
 
namespace DGtal
{
 
// class MultiStatistics
class MultiStatistics
{
    // ----------------------- Standard services ------------------------------
public:
 
 
    ~MultiStatistics();
 
 
 
    MultiStatistics( const unsigned int size, const bool storeSamples = false );
 
 
 
     static void read( std::istream & in, MultiStatistics & samples,
                       const std::vector<unsigned int> & indices );
 
 
 
    // ----------------------- Accessors ------------------------------
  public:
 
    inline unsigned int nb() const;
 
 
    inline unsigned int samples( const unsigned int k ) const;
 
    inline double mean( const unsigned int k ) const;
 
    inline double variance( const unsigned int k ) const;
 
    inline double unbiasedVariance( const unsigned int k ) const;
 
    inline double max( const unsigned int k ) const;
 
 
    inline unsigned int maxIndice( unsigned int k ) const;
 
 
 
    inline double min( const unsigned int k ) const;
 
 
 
    inline unsigned int minIndice(const unsigned int k ) const;
 
 
 
    inline double value( const unsigned int k, const unsigned int i ) const;
  
 
 
    void addValue( unsigned int k, double v );
 
  
    template <class Iter>
    void addValues(const  unsigned int k, Iter b, Iter e );
 
    void terminate();
 
 
    void init( unsigned int size, bool storeSamples );
 
 
    void clear();
  
 
    void erase();
 
  double covariance( const unsigned int x, const unsigned int y, 
                     const unsigned int s = 0, unsigned int e = 0 ) const;
  
    std::pair<double,double> linearRegression( const unsigned int x, 
                                               const unsigned int y ) const;
 
 
 
    double median( const unsigned int k );
 
 
 
 
 
    // ----------------------- Interface --------------------------------------
public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Protected Datas ------------------------------
private:
 
    // ------------------------- Private Datas --------------------------------
private:
 
    unsigned int myNb;
 
    unsigned int* mySamples;
 
    double* myExp;
 
    double* myExp2;
 
    double* myVar;
 
    double* myUnbiasedVar;
    
    double* myMax;
 
    unsigned int* myIndiceMax;
 
    double* myMin;
 
    unsigned int* myIndiceMin;
 
    bool myStoreSamples;
 
    std::vector<double>* myValues;
 
  bool myIsTerminate;
  
    // ------------------------- Hidden services ------------------------------
protected:
 
    MultiStatistics();
 
private:
 
    MultiStatistics ( const MultiStatistics & other );
 
    MultiStatistics & operator= ( const MultiStatistics & other );
 
    // ------------------------- Internals ------------------------------------
private:
 
}; // end of class MultiStatistics
 
 
std::ostream&
operator<< ( std::ostream & out, const MultiStatistics & object );
 
 
} // namespace DGtal
 
 
// Includes inline functions.
#if !defined(BUILD_INLINE)
#include "DGtal/math/MultiStatistics.ih"
#endif
 
 
//                                                                           //
 
#endif // !defined MultiStatistics_h
 
#undef MultiStatistics_RECURSES
#endif // else defined(MultiStatistics_RECURSES)