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)