Quasi-Monte Carlo with One Categorical Variable
研究了当积分变量之一取有限个值时,随机化拟蒙特卡洛估计的误差特性,发现对最小混合成分过采样可提升收敛速度,并给出了最优分配策略。
We study randomized quasi-Monte Carlo (RQMC) estimation of a multivariate integral where one of the variables takes only a finite number of values. This problem arises when the variable of integration is drawn from a mixture distribution as is common in importance sampling and also arises in some recent work on transport maps. We find that when integration error decreases at an RQMC rate that it is then important to oversample the smallest mixture components instead of using a proportional allocation; this can even improve the rate of convergence. The optimal allocations depend on the possibly unknown convergence rate. Designing the sample with an incorrect assumption on the rate still attains that convergence rate, with an inferior implied constant. The penalty for using a pessimistic rate is typically higher than for using an optimistic one. We also find that for the most accurate RQMC sampling methods, it is advantageous to arrange that our n=2m randomized Sobol’ points split into subsample sizes that are also powers of 2.