interest in the presence of nuisance parameters is proposed. The approximation is obtained
using a matching prior. The procedure improves the normal first-order approximation and
has several advantages. It does not require the elicitation on the nuisance parameters,
neither numerical integration nor Monte Carlo simulation, and it enables us to perform
accurate Bayesian inference even for small sample sizes. Numerical illustrations are given …