Sparse model identification enables nonlinear dynamical system discovery from data. However, the control of false discoveries for sparse model identification is challenging, especially in the low-data and high-noise limit. In this paper, we perform a theoretical study on ensemble sparse model discovery, which shows empirical success in terms of accuracy and robustness to noise. In particular, we analyse the bootstrapping-based sequential thresholding least-squares estimator. We show that this bootstrapping-based ensembling technique can perform a provably correct variable selection procedure with an exponential convergence rate of the error rate. In addition, we show that the ensemble sparse model discovery method can perform computationally efficient uncertainty estimation, compared to expensive Bayesian uncertainty quantification methods via MCMC. We demonstrate the convergence properties and connection to uncertainty quantification in various numerical studies on synthetic sparse linear regression and sparse model discovery. The experiments on sparse linear regression support that the bootstrapping-based sequential thresholding least-squares method has better performance for sparse variable selection compared to LASSO, thresholding least-squares, and bootstrapping-based LASSO. In the sparse model discovery experiment, we show that the bootstrapping-based sequential thresholding least-squares method can provide valid uncertainty quantification, converging to a delta measure centered around the true value with increased sample sizes. Finally, we highlight the improved robustness to hyperparameter selection under shifting noise and sparsity levels of the bootstrapping-based sequential thresholding least-squares method compared to other sparse regression methods.
翻译:稀疏模型识别能够从数据中发现非线性动力系统。然而,控制稀疏模型识别中的虚假发现具有挑战性,尤其在低数据量和高噪声条件下。本文对集成稀疏模型发现方法展开理论研究,该方法在准确性和噪声鲁棒性方面已展现出实证优势。具体而言,我们分析了基于自助法的序贯阈值最小二乘估计器,证明了这种自助法集成技术能以指数级误差收敛速率执行可证明正确的变量选择过程。此外,我们表明,与通过MCMC的昂贵贝叶斯不确定性量化方法相比,集成稀疏模型发现方法能以计算高效方式实现不确定性估计。通过合成稀疏线性回归与稀疏模型发现的多项数值研究,我们展示了收敛性质及其与不确定性量化的关联。稀疏线性回归实验表明,相较于LASSO、阈值最小二乘及基于自助法的LASSO方法,基于自助法的序贯阈值最小二乘法在稀疏变量选择中具有更优性能。在稀疏模型发现实验中,我们证明基于自助法的序贯阈值最小二乘法能提供有效的不确定性量化,其估计结果随样本量增加收敛至以真值为中心的狄拉克测度。最后,我们强调与其他稀疏回归方法相比,基于自助法的序贯阈值最小二乘法在噪声水平与稀疏度变化条件下对超参数选择具有更优的鲁棒性。