Esempi: Lezione 9
Contents
Esempi: Lezione 9¶
9.0 Log-likelihood¶
main_00.cpp
/*
c++ -o main_00 main_00.cpp
*/
#include <cmath>
#include <cstdlib>
#include <iostream>
#include <vector>
double esponenziale (double x, double tau)
{
if (tau == 0.) return 1. ;
return exp (-1. * x / tau) / tau ;
}
double loglikelihood (const std::vector<double> & data, double param)
{
double result = 0. ;
for (int i = 0 ; i < data.size () ; ++i) result += log (esponenziale (data.at (i), param)) ;
return result ;
}
int main (int argc, char ** argv)
{
// Riempio un vettore con un set di dati
double dataArr[13] = {0.4,0.2,0.1,1.,1.2,0.8,0.9,0.9,1.5,0.1,0.2,0.1,0.2};
std::vector<double> data(&dataArr[0],&dataArr[0]+13);
// Stampo il valore di log(L)
std::cout << "t=1: log(L)=" << loglikelihood(data,1) << std::endl;
return 0 ;
}