likelihoodStats.cxx

gistfile1.cpp

#include <admodel.h>

 //=================================================================================  
 // Likelihoods for abundance indices or catch
 //=================================================================================  

 dvar_vector nlikelihood(dvector observed, dvar_vector expected, dvector cv, int LikeType)
    {
      RETURN_ARRAYS_INCREMENT();
        dvar_vector logLike;
        switch (LikeType)
          {
            case 0: //not used
            break;
            case 1: //normal
              logLike=elem_div(square(observed-expected), 2*square(elem_prod(cv,observed))); 
            break;
            case 2: //lognormal
              logLike=elem_div(square(log(observed)-log(expected)),2*square(cv)); 
            break;
            default: 
              cout<<"No such likelihood: "<<LikeType<<" /n";
            break;
          }
   
      RETURN_ARRAYS_DECREMENT();
         return(logLike);
    }