export(gibbs_sampling,gibbs2)

