ABC.hpp

/*
    Copyright 2011 Stéphane De Mita
    
    This file is part of the EggLib library.

    EggLib is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    EggLib is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with EggLib.  If not, see <http://www.gnu.org/licenses/>.
*/


#ifndef EGGLIB_ABC_HPP
#define EGGLIB_ABC_HPP

#include <vector>
#include <string>


namespace egglib {

    class ABC {
        
        public:
        
        
            enum TransformMode { NONE, LOG, TAN };


            ABC();

            
            ~ABC();

          
            void number_of_statistics(unsigned int ns);


            void add_fname(const char* fname, unsigned int number_of_params);

            
            unsigned int number_of_samples() const;


            void get_threshold(double tolerance);


            double sd(unsigned int index) const;


            void obs(unsigned int index, double value);

                        
            unsigned int rejection(const char* outfname, bool exportlabels=false);
            
            unsigned int regression(const char* infname, 
                                    const char* outfname,
                                    TransformMode mode,
                                    const char* header = "");
            
            
           double threshold() const; 
            
        private:
        
            ABC(const ABC& src) {}
            ABC& operator=(const ABC& src) {return *this;}
            
            void init();
            void clear();
            
            double *_sd;
            double *_obs;
            unsigned int _nstats;
            unsigned int _nsam;
            double _threshold;
            std::vector<std::string> _fnames;
            std::vector<unsigned int> _nparams;

    };

}

#endif