StatCumProbDistDynCalc.h

Go to the documentation of this file.

#ifndef StatCumProbDistDynCalc_h
#define StatCumProbDistDynCalc_h

#include <vector>

namespace Isis {
  class StatCumProbDistDynCalc{
    //class uses the P^2 Algorithim to calculate equiprobability cell histograms from a stream of data without storing the data
    //  see "The p^2 Algorithim for Dynamic Calculations of Quantiles and Histograms Without Storing Observations"
  public:
    StatCumProbDistDynCalc(unsigned int nodes=20);  //individual qunatile value to be calculated
    ~StatCumProbDistDynCalc() { }; //empty destructor

    void addObs(double obs); //

    double cumProb(double value); //given a value return the cumulative probility
    double value(double cumProb); //given a cumulative probibility return a value
    double max(); //return the largest value so far
    double min(); //return the smallest values so far
    void initialize(unsigned int nodes=20); //resets the class to start a new dynamic calculation

  private:

    unsigned int m_nCells;     //the number of cells in the histogram


    unsigned int m_nQuan;     //the number of quantiles being calculated (m_cells+1)



    std::vector<double> m_quan;        //the target quantile



    std::vector<double> m_nIdeal;      //ideal positions of quantiles



    std::vector<int> m_n;



    std::vector<double> m_q;



    unsigned int m_nObs;  //the number of observations
  };

} //end namespace Isis

#endif