rshrinkage function

Sample from the model-specific posterior of the shrinkage factor t = g / (1 + g)