Conditional on parameters ρ t-1 , β t and σ 2 t , the posterior predictive distribution for the highest level response (PBE0) is Gaussian. However, the marginal posterior predictive distribution, with these parameters integrated out with respect to their posterior distributions (2) and (3), is not available in closed-form. However, its mean and variance are available, and take the following form: