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 Data ------------------------------
private:
 
    // ------------------------- Private Data --------------------------------
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)