Hi,

I am working on writing a Gibbs sampler to do hypothesis testing in a Bayesian framework. I have a set of data which give the number of occurrences in a sample of trials. The particular model of interest is a three-layered hierarchical model with the following distributions:

1) The number of occurrences given is drawn from a Binomial distribution, .

2) The distribution of given the hyperparameter is a Beta distribution, .

Here is a known constant between 0 and 1. Note that the shape paramters and are chosen so that the mean is .

3) is drawn from a uniform distribution on where which ensures that the distribution is concave.

If I am understanding the Gibbs sampling procedure correctly, it goes something like this.

Start by generating a value for from the uniform distribution. Then:

1) Given this generate a value for from .

2) Given generate a value for from .

3) Given generate a new value for from .

4) Given this new generate a new value for from .

5) Repeat steps 1-4.

I am having trouble with step 4, since I need the conditional distribution . I should be able to get this from Bayes theorem, since .

I have been unable to calculate the normalization integral that appears in the denominator above:

.

Question 1: Does anyone know of a method to perform this integration of the Beta distribution with respect to analytically? I have been unable to do so or find much on integration of the Beta distribution with respect to the shape parameter.

Question 2: For those readers who are familiar with Gibbs sampling and the Bayesian framework, please feel free to comment on the method of approach I outlined here. I am rather new at this and not entirely confident in the way I am attempting to do the sampling.

Thanks in advance,

basmati