R.DGAMMA(x,shape,scale,give_log)
x: observation
shape: the shape parameter of the distribution
scale: the scale parameter of the distribution
give_log: if true, log of the result will be returned instead
This function returns the probability density function of the gamma distribution.