update_all_y function

Gibbs sampling for yy in xCx model