In addition to recent developments in computing speed and memory, methodological advances have contributed to significant gains in the performance of stochastic simulation. In this paper, we focus on variance reduction for matrix computations via matrix factorization. We provide insights into existing variance reduction methods for estimating the entries of large matrices. Popular methods do not exploit the reduction in variance that is possible when the matrix is factorized. We show how computing the square root factorization of the matrix can achieve in some important cases arbitrarily better stochastic performance. In addition, we propose a factorized estimator for the trace of a product of matrices and numerically demonstrate that the estimator can be up to 1,000 times more efficient on certain problems of estimating the log-likelihood of a Gaussian process. Additionally, we provide a new estimator of the log-determinant of a positive semi-definite matrix where the log-determinant is treated as a normalizing constant of a probability density.
翻译:除了最近在计算速度和记忆方面的发展外,方法上的进步还有助于在进行随机模拟方面取得重大进展。在本文件中,我们侧重于通过矩阵因子化来减少矩阵计算的差异;我们深入了解用于估计大型矩阵条目的现有差异减少方法;通用方法没有利用在矩阵因子化时可能减少的差异;我们展示了在某些重要情况下计算矩阵的平方根因子化如何达到任意性强的随机性能;此外,我们提议了用于追踪矩阵产品和数字性地显示估计器在估计高斯进程日志相似性的某些问题上的效率可高达1 000倍。此外,我们为正半确定性矩阵的日志定值提供了一个新的估计器,将日志-半确定性矩阵作为概率密度的正常常数处理。