fitGammaBin function

Fitting the Gamma Binomial distribution when binomial random variable, frequency and shape parameters are given