model_comp.h

#pragma once
#include <cmath>
 
 
namespace cgv{
    namespace math {
 
template<typename T>
T AIC(int k, const T& L)
{
    return (const T)2.0*(T)k - (T)2*(T)log((double)L);
}
 
 
template<typename T>
T AIC_ls(int k, int n, const T&RSS)
{
    return 2*k + (T)n*log(RSS/(T)n);
}
 
template<typename T>
T AICc(int k, int n, const T& L)
{
    return AIC<T>(k,L) + (T)(2*k*(k+1))/((T)(n-k-1));
}
 
template<typename T>
T AICc_ls(int k, int n, const T&RSS)
{
    return (T)log((T)RSS/(T)n) + (T)(n+k)/(T)(n-k-2) ;
}
 
template<typename T>
T AICu_ls(int k, int n, const T&RSS)
{
    return (T)log((T)RSS/(n-k)) + (T)(n+k)/(T)(n-k-2) ;
}
 
template<typename T>
T BIC(int k, int n, const T& L)
{
    return -(T)2*(T)log((T)L) + (T)k*log((T)n);
}
 
template<typename T>
T BIC_ls(int k, int n, const T& RSS)
{
    return (T)n*log((T)RSS/(T)n) + (T)k*log((T)n);
}
 
    }
}