Automatic approximation of the marginal likelihood in non-Gaussian hierarchical models