individual_likelihood.h

#ifndef _INDIVIDUAL_LIKELIHOOD_H_
#define _INDIVIDUAL_LIKELIHOOD_H_

#include <cmath>
#include <math.h>
#include <iomanip>      // std::setprecision
#include <cfloat>

#include "lnmultinomial.h"
#include "data_types/allele.h"
#include "newton-method-theta.h"
#include "models.h"

count_t init_params(Locus&, Allele&, const float_t&);                                       
count_t maximize_analytical(Locus&, Allele&, models&, std::vector <float_t>&, const float_t&, const size_t&);   
count_t maximize_grid      (Locus&, Allele&, models&, std::vector <float_t>&, const float_t&, const size_t&);   
count_t maximize_newton    (Locus&, Allele&, models&, std::vector <float_t>&, const float_t&, const size_t&);   
count_t maximize_restricted_newton    (Locus&, Allele&, models&, std::vector <float_t>&, const float_t&, const size_t&);    
float_t loglikelihood(const Locus&, const Allele&, models&, const count_t&);    

#endif