How to debug for my Gibbs sampler of Bayesian regression in R?

Question

I am working on a Gibbs sampler of three parameters and we know the full conditional distribution of three parameters beta, pi~Beta(sum(z)+0.5, sum(1-z)+0