I need to estimate $$ \mathbb{E}_x[f_i(x)] = \int_{\mathbb{R}^n} f_i(x) p(x) dx $$ for many functions $f_i(x)$, where $p(x)$ is the density of a normal distribution. The evaluation of all the functions $f_i(x)$ is expensive. I was looking into Bayesian Quadrature but I cannot find an (approximately) optimal way of choosing sample points.

  • $\begingroup$ How accurate does this estimate need to be? Highly accurate or somewhat accurate? Also: what properties about $f$ do you know? Is it $L$-Lipschitz or smooth in some other sense, or does $f$ vary widely? $\endgroup$ – cdipaolo Aug 18 at 15:41
  • $\begingroup$ The approximation doesn't need to be higly accurate, $f$ is L-Lipschitz. $\endgroup$ – Rosh Aug 18 at 16:23
  • $\begingroup$ in that case, you can always get a Monte-Carlo estimate (depending on your internal definition of “highly accurate.”) Since $f$ is Lipschitz, $f(x)$ is sub-Gaussian and obeys nice concentration properties. $\endgroup$ – cdipaolo Aug 19 at 5:29
  • $\begingroup$ The problem is that Monte-Carlo converges slowly. I was looking to find a way of defining the sampling points (with the corresponding weights) to converge faster than Monte-Carlo. $\endgroup$ – Rosh Aug 19 at 8:27

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.