Randomized Dimension Reduction for Monte Carlo Simulations
提出一种无偏估计算法,通过随机降维降低蒙特卡洛模拟的方差,在相同计算成本下方差可降低d倍,适用于高维函数和马尔可夫链期望估计,实验显示对不连续函数也有效。
We present a new unbiased algorithm that estimates the expected value of f(U) via Monte Carlo simulation, where U is a vector of d independent random variables, and f is a function of d variables. We assume that f does not equally depend on all its arguments. Under certain conditions, we prove that, for the same computational cost, the variance of our estimator is lower than the variance of the standard Monte Carlo estimator by a factor of order d. Our method can be used to obtain a low-variance unbiased estimator for the expectation of a function of the state of a Markov chain at a given time step. We study applications to volatility forecasting and time-varying queues. Numerical experiments show that our algorithm dramatically improves on the standard Monte Carlo method for large values of d and is highly resilient to discontinuities. This paper was accepted by Baris Ata, stochastic models and simulation.