added a brief description of how we can compute the posterior for h_j by simply numerically marginalizing