gibbs_jitter function

Sample node dates conditional on root time using Gibbs sampling and date prior