LikelihoodFunction.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2014 by Federico Marulli and Alfonso Veropalumbo  *
 *  federico.marulli3@unibo.it                                      *
 *                                                                  *
 *  This program 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 2 of  *
 *  the License, or (at your option) any later version.             *
 *                                                                  *
 *  This program 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 this program; if not, write to the Free      *
 *  Software Foundation, Inc.,                                      *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.       *
 ********************************************************************/
 
#ifndef __LIKELFUN__
#define __LIKELFUN__
 
#include "Model2D.h"
 
 
// ===================================================================================================
 
 
namespace cbl {
 
  namespace statistics {
    
    enum class LikelihoodType {
 
      _NotSet_,
 
      _Gaussian_Error_,
      
      _Gaussian_Covariance_,
 
      _Poissonian_,
      
      _UserDefined_
 
    };
 
    inline std::vector<std::string> LikelihoodTypeNames () {return {"NotSet", "Gaussian_Error", "Gaussian_ColikelihoodTypeiance", "Poissonian", "UserDefined"}; }
 
    inline LikelihoodType LikelihoodTypeCast (const int likelihoodTypeIndex) {return castFromValue<LikelihoodType>(likelihoodTypeIndex);}
 
    inline LikelihoodType LikelihoodTypeCast (const std::string likelihoodTypeName) {return castFromName<LikelihoodType>(likelihoodTypeName, LikelihoodTypeNames());}
 
    inline std::vector<LikelihoodType> LikelihoodTypeCast (const std::vector<int> likelihoodTypeIndeces) {return castFromValues<LikelihoodType>(likelihoodTypeIndeces);} 
 
    inline std::vector<LikelihoodType> LikelihoodTypeCast (const std::vector<std::string> likelihoodTypeNames) {return castFromNames<LikelihoodType>(likelihoodTypeNames, LikelihoodTypeNames());}
 
    struct STR_likelihood_inputs
    {
      std::shared_ptr<data::Data> data;
 
      std::shared_ptr<Model> model;
 
      std::vector<double> xx;
 
      std::vector<double> yy;
 
      std::vector<double> weights1D;
 
      std::vector<std::vector<double>> weights2D;
 
      std::shared_ptr<glob::FuncGrid> interp_function1D;
 
      std::shared_ptr<glob::FuncGrid2D> interp_function2D;
 
      STR_likelihood_inputs (const std::shared_ptr<data::Data> input_data, const std::shared_ptr<Model> input_model, const std::vector<size_t> input_x_index={0, 2}, const int input_w_index=-1); 
      
    };
 
    double LogLikelihood_1D_interpolated (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_2D_interpolated (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_Gaussian_1D_error (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_Gaussian_1D_covariance (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_Gaussian_2D_error (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_Poissonian_1D_ (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    double LogLikelihood_Poissonian_2D_ (std::vector<double> &likelihood_parameter, const std::shared_ptr<void> input);
 
    typedef std::function<double (std::vector<double> &, const std::shared_ptr<void>)> LogLikelihood_function;
 
    typedef std::function<double (std::vector<double> &, const std::shared_ptr<void>)> Likelihood_function;
  }
}
 
#endif