acetp_mcmc function

Compute CIs for the ACE(t)-p model