We propose a Monte-Carlo-based method for reconstructing sparse signals in the formulation of sparse linear regression in a high-dimensional setting. The basic idea of this algorithm is to explicitly select variables or covariates to represent a given data vector or responses and accept randomly generated updates of that selection if and only if the energy or cost function decreases. This algorithm is called the greedy Monte-Carlo (GMC) search algorithm. Its performance is examined via numerical experiments, which suggests that in the noiseless case, GMC can achieve perfect reconstruction in undersampling situations of a reasonable level: it can outperform the $\ell_1$ relaxation but does not reach the algorithmic limit of MC-based methods theoretically clarified by an earlier analysis. The necessary computational time is also examined and compared with that of an algorithm using simulated annealing. Additionally, experiments on the noisy case are conducted on synthetic datasets and on a real-world dataset, supporting the practicality of GMC.
翻译:我们提出一种基于蒙特-卡洛的重建高维环境中微小线性回归的微弱信号的方法。这种算法的基本想法是明确选择变量或共变来代表特定数据矢量或响应,并在能量或成本功能下降时接受随机生成的选择更新。这种算法被称为贪婪的蒙特-卡洛(GMC)搜索算法。它的性能通过数字实验来审查,它表明在无噪音的情况下,GMC可以在没有噪音的情况下在合理水平的抽样情况中实现完美的重建:它可以超过$\ell_1美元,但不能达到早期分析从理论上澄清的以MC为基础的方法的算法极限。必要的计算时间也经过审查,并与使用模拟Annealing算法的算法进行比较。此外,关于噪音案例的实验是在合成数据集和真实世界数据集上进行的,支持GMC的实用性。