update_all_y_Sj function

Gibbs sampling for yy in xUx model