Monte Carlo integration function to find the nonparametric adjustments
int.eprior(sdat, g.hat, d.hat)
result