aic.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
/*****************************************************************************/
#include "tpcmodels.h"
/*****************************************************************************/
 
/*****************************************************************************/
double aicSS(
  double ss,
  const unsigned int n,
  const unsigned int k
) {
  if(!(ss>=0.0) || n<1 || (n-k)<2) return(nan(""));
 
  int dr=n-k-1; 
  int dv=2*k*(k+1);
  double aic=0.0, bias_adj=0.0, css;
  if(dr>0) bias_adj=((double)dv)/((double)dr); else bias_adj=0.0;
  if(ss<1.0e-50) css=1.0e-50; else css=ss; /* Because log(0) is an error */
  aic= n*log(css/(double)n) + 2.0*(double)k + bias_adj;
  return(aic);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int parFreeNr(
  const unsigned int n,
  double *pLower,
  double *pUpper
) {
  if(n==0 || pLower==NULL || pUpper==NULL) return(0);
  unsigned int nf=0;
  for(unsigned int i=0; i<n; i++) {
    double range=pUpper[i]-pLower[i];
    if(range>1.0E-10) nf++;
  }
  return(nf);
}
/*****************************************************************************/
 
/*****************************************************************************/