I've read several sources of information that describe the process of Bayesian Monte Carlo Quadrature but am just not understanding the details enough to be able to implement it.
For instance two sources:
I get some high level ideas like:
- It uses prior information better than other monte carlo algorithms such as Quasi Monte Carlo.
- It can converge faster as a result, for integrands that are smooth.
- It gives you a "confidence interval" so you can get an idea of how much variance / noise there is.
- It seems to somehow use gaussian functions as basis functions.
- every data point somehow fits the sum of gaussian basis functions to the actual function better?
Honestly, other than that it is pretty hazy.
Can anyone explain in a simple way how you could integrate a function over a specific domain, using Bayesian Monte Carlo Quadrature? Or any links to drop dead simple explanations, or implementations would be helpful too.