dinvgamma.cpp

Go to the documentation of this file.

#include "statsLib.h"

//inverse-gamma
dvariable dinvgamma(const dvariable& x, const double a, const double b)
{
  return -a*log(b)+gammln(a)+(a+1)*log(x)+b/x;
}