R.DBETA(x,a,b,give_log)
x: observation
a: the first shape parameter of the distribution
b: the second 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 beta distribution.