[R] comment in src/nmath/dgamma.c

[Faheem Mitha]

Hi,

In src/nmath/dgamma.c the comment at the top says

  * DESCRIPTION
  *
  *   Computes the density of the gamma distribution,
  *
  *                   1/s (x/s)^{a-1} exp(-x/s)
  *        p(x;a,s) = -----------------------
  *                            (a-1)!
  *
  *   where `s' is the scale (= 1/lambda in other parametrizations)
  *     and `a' is the shape parameter ( = alpha in other contexts).

Maybe I'm missing something, but shouldn't the denominator be Gamma(a) 
(since a is not necessarily an integer)?

                                                               Faheem.