An approach for finding fully Bayesian optimal designs using normal-based approximations to loss functions
The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo