set_initial_gibbs_values function

Set initial values for the Gibbs sampler