STAT_Object.h

Go to the documentation of this file.

#ifndef STAT_Object_H
#define STAT_Object_H

#include "CORE_Object.h"
#include "CORE_Array.h"

DEFINE_SPTR(STAT_Object);
class STAT_Object : public CORE_Object {
    SP_OBJECT(STAT_Object);
    // ATTRIBUTES
  
private:

  
  
    // ASSOCIATIONS
  
  
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    STAT_Object(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~STAT_Object(void);   
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------
    
    // --------------
    // SAVER & LOADER
    // --------------

   


    template<typename T>
    inline static tReal quantile(const T *d,const int& size,const tReal& q) {
        if (size == 0) return 0;
        if (size == 1) return d[0];
        if (q <= 0) return *std::min_element(d, d + size);
        if (q >= 1) return *std::max_element(d, d + size);
        
        double pos = (size - 1) * q;
        unsigned int ind = (unsigned int) pos;
        double delta = pos - ind;
        vector<tReal> w(size); std::copy(d, d + size, w.begin());
        std::nth_element(w.begin(), w.begin() + ind, w.end());
        tReal i1 = *(w.begin() + ind);
        tReal i2 = *std::min_element(w.begin() + ind + 1, w.end());
        return i1 * (1.0 - delta) + i2 * delta;
    }
    template<typename T>
    inline static tReal quantile(T *d,const int& size,const tReal& q) {
        //cout << "sample order is changed...\n";
        if (size == 0) return 0;
        if (size == 1) return d[0];
        if (q <= 0) return *std::min_element(d, d + size);
        if (q >= 1) return *std::max_element(d, d + size);
        
        double pos = (size - 1) * q;
        unsigned int ind = (unsigned int) pos;
        double delta = pos - ind;
        
        std::nth_element(d, d + ind, d+size);
        tReal i1 = *(d + ind);
        tReal i2 = *std::min_element(d + ind + 1, d+size);
        return i1 * (1.0 - delta) + i2 * delta;
    };
    
    template<typename T>
    inline static tReal quantile(const CORE_Array<T>& a,const tReal& q) {
        return quantile(a.getValues(),a.getSize(),q);
    }
};

#endif