R coding of fiducial estimation using Gibbs sampler