images/gamma_dist.gnuplot

   1 _ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x
   2 pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)
   3 dgamma(x, shape, rate) = (x<0)? 0 : (x==0)? ((shape<1)? 1/0 : (shape==1)? rate : 0) : (rate==0)? 0 : exp(_ln_dgamma(x, shape, rate))
   4 gamma_offset = 0.0
   5 gamma_k = 15.7880378786862
   6 gamma_theta = 0.980078674752871
   7 gamma_fit(x) = 100.0 * dgamma(x - gamma_offset, gamma_k, gamma_theta)
   8 #gamma_dist(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma(gamma_k))
   9 #gamma_fit(x) = 100.0 * gamma_dist(x - gamma_offset)
  10