theta_init function

Initialize the variables for the Gibbs sampler chain