Bayesian Chi2 goodness-of-fit test

Hello all, I'm trying to implement a GoF measure I found in the literature but haven't had much success replicating even the simplest situation (simulated binomial data).

I've posted the full details with relevant equations and R code over at Cross Validated.

Would be very grateful if anyone can give it a look.