Added a Gibbs sampler update rule. This should be able to more accurately describe uncertainties. It does *not* currently play well with the way that statistical errors are computed in a Pidrixter (bootstrapping). This will result in double counting the uncertainties if both are used.