made a test function for the gibbs sampler - bug at line 127 in calculating odds ratio for X(n,k) sample