Driven by several successful applications such as in stochastic gradient descent or in Bayesian computation, control variates have become a major tool for Monte Carlo integration. However, standard methods do not allow the distribution of the particles to evolve during the algorithm, as is the case in sequential simulation methods. Within the standard adaptive importance sampling framework, a simple weighted least squares approach is proposed to improve the procedure with control variates. The procedure takes the form of a quadrature rule with adapted quadrature weights to reflect the information brought in by the control variates. The quadrature points and weights do not depend on the integrand, a computational advantage in case of multiple integrands. Moreover, the target density needs to be known only up to a multiplicative constant. Our main result is a non-asymptotic bound on the probabilistic error of the procedure. The bound proves that for improving the estimate's accuracy, the benefits from adaptive importance sampling and control variates can be combined. The good behavior of the method is illustrated empirically on synthetic examples and real-world data for Bayesian linear regression.
翻译:在一些成功的应用(例如,在随机梯度梯度下降或巴伊西亚计算中)的驱动下,控制变异已成为蒙特卡洛整合的主要工具。然而,标准方法不允许在算法期间粒子的分布演变,正如相继模拟方法那样。在标准的适应重要性抽样框架内,建议采用简单加权最小方格的方法来改进控制变异的程序。该程序的形式是四方形规则,以经调整的二次体积重量来反映控制变异带来的信息。四方位点和重量并不取决于指数,对于多个数子而言,这是一种计算优势。此外,目标密度仅需要知道到一个倍复制的常数。我们的主要结果是一个不依赖程序概率误差的不考虑方形方法。为了提高估计的准确性,可将适应重要性取样和控制变异变的效益结合起来。该方法的良好行为是用合成实例和巴伊斯直线性真实世界数据进行实验性说明。