The posterior for a,b,e ð Þ can then be obtained by extending the parameter vector to a,b,e,t 1 ,t 2 , . . . ð Þ [H, i.e. including the unobserved event times as parameters, and using MCMC to explore the augmented parameter space H [34, 35] (for recent applications see e.g. [21, [36] [37] [38] ).