English

Improved Maximum Entropy Analysis with an Extended Search Space

Computational Physics 2013-01-08 v3 High Energy Physics - Lattice Nuclear Theory Data Analysis, Statistics and Probability

Abstract

The standard implementation of the Maximum Entropy Method (MEM) follows Bryan and deploys a Singular Value Decomposition (SVD) to limit the dimensionality of the underlying solution space apriori. Here we present arguments based on the shape of the SVD basis functions and numerical evidence from a mock data analysis, which show that the correct Bayesian solution is not in general recovered with this approach. As a remedy we propose to extend the search basis systematically, which will eventually recover the full solution space and the correct solution. In order to adequately approach problems where an exponentially damped kernel is used, we provide an open-source implementation, using the C/C++ language that utilizes high precision arithmetic adjustable at run-time. The LBFGS algorithm is included in the code in order to attack problems without the need to resort to a particular search space restriction.

Keywords

Cite

@article{arxiv.1110.6285,
  title  = {Improved Maximum Entropy Analysis with an Extended Search Space},
  author = {Alexander Rothkopf},
  journal= {arXiv preprint arXiv:1110.6285},
  year   = {2013}
}

Comments

18 pages, 6 figures, v3 includes several changes in text and figures, t.b.p. in Journal of Computational Physics, source code at http://www.scicode.org/ExtMEM

R2 v1 2026-06-21T19:27:25.309Z