Estimate SIR epidemic parameters via Gibbs sampling