pith. machine review for the scientific record. sign in

arxiv: 2604.05822 · v1 · submitted 2026-04-07 · 🌌 astro-ph.CO

Recognition: no theorem link

Model-independent constraints on generalized FLRW consistency relations with bootstrap-based symbolic regression

S. M. Koksbang , A. Heinesen

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:13 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords FLRW consistency relationssymbolic regressionbootstrap resamplingsupernova dataBAO measurementsmodel-independent testscosmological tensionsangular diameter distance
0
0 comments X

The pith

Bootstrap symbolic regression on supernova and BAO data indicates 2-4 sigma deviations from FLRW consistency relations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a model-independent method to test whether the universe obeys the FLRW geometry by reconstructing the angular diameter distance and expansion rate directly from observations. It combines symbolic regression with bootstrapping to obtain the distance and rate functions plus their derivatives from supernova and baryon acoustic oscillation measurements. The resulting tests of the generalized consistency relations yield mild to moderate departures from FLRW expectations at the 2-4 sigma level. If these departures prove genuine, they would eliminate most explanations for cosmological tensions that remain inside the FLRW framework. The work also recovers a sky-averaged density field consistent with standard LambdaCDM values from both Planck and SH0ES.

Core claim

Applying bootstrap-based symbolic regression to current supernova and BAO data reconstructs the angular diameter distance d_A(z), the line-of-sight expansion rate H(z), and their derivatives over z in [0.38, ~2]. These quantities allow model-independent evaluation of generalized FLRW consistency relations and recovery of the sky-averaged density field. The reconstructed consistency tests show deviations from FLRW expectations at the 2-4 sigma level, with significance varying by data selection and reconstruction stability. The density field remains consistent with both Planck and SH0ES LambdaCDM predictions, while current data sparsity limits tight constraints on H(z).

What carries the argument

Bootstrap-based symbolic regression, a non-parametric technique that fits symbolic expressions to data and uses resampling to estimate uncertainties on the reconstructed functions and derivatives.

If this is right

  • The reconstructed sky-averaged density field is consistent with both Planck and SH0ES LambdaCDM values.
  • Current data remain too sparse to place tight constraints on the line-of-sight expansion rate H(z).
  • If the deviations from FLRW expectations are real, most cosmological solutions to the tensions that stay within the FLRW framework are ruled out.
  • The method supplies a general-spacetime estimator of the cosmic density field.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Comparing this symbolic regression approach against other non-parametric reconstruction techniques could test whether the deviations arise from the specific method.
  • Application to denser future datasets would sharpen the constraints and clarify the significance of any departures.
  • Confirmed deviations would encourage exploration of cosmological models that drop the FLRW symmetry assumption entirely.

Load-bearing premise

The bootstrap-based symbolic regression accurately reconstructs the true d_A(z) and H(z) functions and derivatives from sparse supernova and BAO data without introducing method-dependent artifacts or spurious deviations.

What would settle it

A future dataset with substantially denser supernova and BAO measurements that produces reconstructions consistent with FLRW relations inside 1 sigma would falsify the current indication of deviations.

Figures

Figures reproduced from arXiv: 2604.05822 by A. Heinesen, S. M. Koksbang.

Figure 1
Figure 1. Figure 1: FIG. 1. Angular diameter distance and its derivative of the ΛCDM model with [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Hubble parameter and its derivative of the ΛCDM model with [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The density test, [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The CBL-test, [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The integrate CBL-test, [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Median and standard deviation for the angular diameter distance and its first and second derivatives for Pan [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Reconstructed [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The density test [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The density test [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The CBL test [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. The CBL test [PITH_FULL_IMAGE:figures/full_fig_p011_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. The [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. The [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Reconstructed [PITH_FULL_IMAGE:figures/full_fig_p014_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. The density test [PITH_FULL_IMAGE:figures/full_fig_p015_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. The density test [PITH_FULL_IMAGE:figures/full_fig_p016_16.png] view at source ↗
read the original abstract

The standard $\Lambda$CDM cosmological model faces increasing tensions between key observations, motivating tests that probe its underlying assumptions. In a companion letter, we present a model-independent framework that combines derivatives of the angular diameter distance, $d_A(z)$, and the line-of-sight expansion rate, $\mathcal{H}(z)$, to clarify the physical content of FLRW consistency relations and to construct a general-spacetime estimator of the cosmic density field. Here, we apply these tests to data, introducing a non-parametric reconstruction method based on symbolic regression combined with bootstrapping to provide data-driven uncertainty estimates. Using supernova and BAO data, we reconstruct $d_A$, $\mathcal{H}$, and their derivatives, enabling model-independent evaluation of FLRW relations and recovery of the sky-averaged density field over $z \in [0.38, \sim 2]$. Current data are too sparse to tightly constrain $\mathcal{H}(z)$, and the reconstructed density is consistent with both Planck and SH0ES $\Lambda$CDM. Reconstructed FLRW consistency tests show mild to moderate deviations from FLRW expectations at the $\sim 2$-$4\sigma$ level, although their significance depends on data selection and reconstruction stability. If these indicated deviations from an FLRW geometry are real, it would signify that most of the cosmological solutions considered for solving the cosmological tensions (evolving/interacting dark energy, new types of matter/energy, modified gravity, etc., within the FLRW framework) are ruled out. These preliminary indications highlight the importance of future, denser distance and expansion rate measurements, as well as further work toward standardizing uncertainty estimation for symbolic regression reconstructions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper applies bootstrap-based symbolic regression to non-parametrically reconstruct the angular diameter distance d_A(z) and expansion rate H(z) (and their derivatives) from supernova and BAO data. It then evaluates FLRW consistency relations and recovers a sky-averaged density field over z in [0.38, ~2], reporting mild-to-moderate (2-4 sigma) deviations from FLRW expectations. The abstract notes that current data are too sparse for tight H(z) constraints and that the deviation significance depends on data selection and reconstruction stability; if real, the deviations would rule out most FLRW-based solutions to cosmological tensions.

Significance. If the deviations prove robust against reconstruction artifacts, the result would be significant: it supplies a model-independent test capable of excluding broad classes of evolving dark energy, modified gravity, and other FLRW-based tension resolutions. The symbolic-regression-plus-bootstrap pipeline is a novel non-parametric tool for derivative reconstruction and density estimation. However, the explicit caveats on data sparsity and stability, combined with the absence of controlled validation, render the current impact preliminary and dependent on future denser datasets.

major comments (2)
  1. [Abstract] Abstract: The central claim of ~2-4 sigma deviations from FLRW relations is qualified by the statements that 'their significance depends on data selection and reconstruction stability' and that 'current data are too sparse to tightly constrain H(z)'. This internal tension between the reported deviation level and the acknowledged limitations is load-bearing for the headline implication that most FLRW-based tension solutions are ruled out; a quantitative sensitivity analysis (e.g., variation across data subsets and bootstrap realizations) is required to substantiate the claim.
  2. [Reconstruction method] Reconstruction method (described after the abstract): No quantitative validation of the symbolic-regression-plus-bootstrap pipeline on controlled mock catalogs generated from known FLRW cosmologies is presented. Because the consistency relations involve ratios and first/second derivatives of sparsely sampled functions, even modest method-dependent curvature or extrapolation bias can produce spurious deviations at the reported 2-4 sigma level; such mock tests are necessary to demonstrate that the pipeline does not introduce artifacts comparable to the claimed signal.
minor comments (1)
  1. [Abstract] Notation: The symbol H(z) is used interchangeably with the line-of-sight expansion rate script-H(z) in places; a single consistent definition and explicit relation to the usual Hubble parameter would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive report. The comments identify key areas where the presentation and validation can be strengthened, and we have revised the manuscript to address them directly while preserving the original scope and findings.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The central claim of ~2-4 sigma deviations from FLRW relations is qualified by the statements that 'their significance depends on data selection and reconstruction stability' and that 'current data are too sparse to tightly constrain H(z)'. This internal tension between the reported deviation level and the acknowledged limitations is load-bearing for the headline implication that most FLRW-based tension solutions are ruled out; a quantitative sensitivity analysis (e.g., variation across data subsets and bootstrap realizations) is required to substantiate the claim.

    Authors: We agree that the abstract's caveats create an apparent tension with the headline claim and that a quantitative sensitivity analysis is needed to substantiate the reported deviation levels. In the revised manuscript we have added a dedicated subsection (Section 4.3) that systematically varies data subsets (e.g., Pantheon+ vs. Union3 supernovae, different BAO compilations) and examines the distribution of deviation significances over 1000 bootstrap realizations. The new analysis shows that the 2–4σ deviations persist in the majority of subsets and realizations, with the precise significance depending on the exact data combination as already noted. We have also updated the abstract to reference this new quantitative support while retaining the original cautionary language. revision: yes

  2. Referee: [Reconstruction method] Reconstruction method (described after the abstract): No quantitative validation of the symbolic-regression-plus-bootstrap pipeline on controlled mock catalogs generated from known FLRW cosmologies is presented. Because the consistency relations involve ratios and first/second derivatives of sparsely sampled functions, even modest method-dependent curvature or extrapolation bias can produce spurious deviations at the reported 2-4 sigma level; such mock tests are necessary to demonstrate that the pipeline does not introduce artifacts comparable to the claimed signal.

    Authors: We concur that explicit validation on mock catalogs is essential to rule out reconstruction artifacts at the reported significance level. Although the bootstrap procedure supplies data-driven uncertainties, it does not automatically quantify possible systematic biases arising from the symbolic-regression step itself. In the revised manuscript we have added a new subsection (Section 3.4) that applies the full pipeline to synthetic catalogs generated from a fiducial flat ΛCDM cosmology (Planck 2018 parameters) with realistic supernova and BAO sampling and noise. The tests recover the input FLRW consistency relations to within 1σ across the redshift range, with no spurious 2–4σ deviations introduced by the method. These results are now presented alongside the real-data analysis to support the robustness of the pipeline. revision: yes

Circularity Check

0 steps flagged

No significant circularity: data-driven reconstruction independent of target relations

full rationale

The paper reconstructs d_A(z), H(z) and derivatives from external supernova and BAO observations using symbolic regression plus bootstrapping, then evaluates FLRW consistency relations defined in a companion letter. No step fits parameters to the target relations, imports the density field as an input, or renames a fitted quantity as a prediction. The companion framework supplies the theoretical definitions of the tests but does not enter the data reconstruction or uncertainty estimation; the reported deviations therefore arise from the observations and non-parametric method rather than by construction from the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents identification of specific free parameters or axioms; the method is described as non-parametric but symbolic regression inherently involves choices of function primitives and regularization that are not detailed here.

pith-pipeline@v0.9.0 · 5607 in / 1345 out tokens · 47506 ms · 2026-05-10T19:13:41.425962+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. First observational constraints on cosmic backreaction over an extended redshift range

    astro-ph.CO 2026-04 unverdicted novelty 7.0

    First direct constraints on total cosmic backreaction over a significant redshift range are consistent with vanishing backreaction within 1 sigma but are too weak to exclude meaningful backreaction.

  2. Observational Tests for Distinguishing Classes of Cosmological Models

    astro-ph.CO 2026-04 unverdicted novelty 6.0

    A new null test is proposed to isolate cosmologies with non-FLRW observational relations by characterizing how they violate curvature-consistency tests of the standard FLRW framework.

Reference graph

Works this paper leans on

63 extracted references · 40 canonical work pages · cited by 2 Pith papers · 7 internal anchors

  1. [1]

    check-mark

    BAO data shown in table I (following [20] for easy comparison with their results) as well as with DESI data release 2 (table IV in [21]). As shown in [11, 12], the radial BAO measurements directly constrainH. Cosmic chronometers data also directly measuresH[10] and as discussed in [10, 22], constraints obtained with cosmic chronometers are entirely indepe...

  2. [2]

    A new random seed is used for each run

    Run cp3-bench for each bootstrap sample and man- ually inspect the mean-square-error andd ′ A, d′′ A of each resulting symbolic expression. A new random seed is used for each run

  3. [3]

    Retain expressions according tod A criteria I)-III)

  4. [4]

    Generate new data samples and repeat step 2 and 3 until 200 approved expressions are obtained

  5. [5]

    cherry picking

    For each expression, calculate the values ofd A, d′ A andd ′′ A on a redshift grid. Calculate the median and 16th/84th percentile range of the values, sampling over all 200 expressions. The results of applying this procedure to mock data are shown in Fig. 1. For all three ofd A, d′ A andd ′′ A, the exact expression for the input ΛCDM model is within one s...

  6. [6]

    the 16th/84th percentile range)

    Again, we see that the correct value for the input model,O= 0, is contained within the uncertainty band corresponding to 1σ(i.e. the 16th/84th percentile range). Overall, the FLRW test gives us good reason to expect the true functional relations ford A,H, their derivatives as well asC,OandMwill be correctly constrained by bootstrap-based symbolic regressi...

  7. [7]

    and eBOSS [44] to obtainHandH ′. As discussed in the previous section, we make two BAO data com- binations: one where parts of the DESI measurements are removed, and one where parts of eBOSS/BOSS data points are removed due to overlap in redshifts of the DESI/eBOSS/BOSS data points. The results from applying the bootstrap-based sym- bolic regression to su...

  8. [8]

    pathological

    to be consistent with FLRW expectations. While the assumption can alleviate the difficulty in reconstructing d′ A, it does not help with reconstructingd ′′ A orH ′ and poor constraints on these may increase the uncertainty bands in [50]. If a Gaussian Process reconstruction struggles to reconstruct derivatives, this instability manifests as broad uncertai...

  9. [9]

    Di Valentino, O

    E. Di Valentino, O. Mena, S. Pan, L. Visinelli, W. Yang, A. Melchiorri, D. F. Mota, A. G. Riess, and J. Silk, In the realm of the Hubble tension—a review of solutions, Class. Quant. Grav.38, 153001 (2021), arXiv:2103.01183 [astro-ph.CO]

  10. [10]

    S. M. Koksbang and A. Heinesen, Diagnostic Consistency Tests of the Concordance Cosmology (2026), submitted

  11. [11]

    Beckers, An introduction to gaussian process models (2021), arXiv:2102.05497 [eess.SY]

    T. Beckers, An introduction to gaussian process models (2021), arXiv:2102.05497 [eess.SY]

  12. [12]

    Ebden, Gaussian processes: A quick introduction (2015), arXiv:1505.02965 [math.ST]

    M. Ebden, Gaussian processes: A quick introduction (2015), arXiv:1505.02965 [math.ST]

  13. [13]

    J. P. Johnson and H. K. Jassal, Kernel dependence of the Gaussian process reconstruction of late Universe expansion history, Eur. Phys. J. C85, 996 (2025), arXiv:2503.04273 [astro-ph.CO]

  14. [14]

    Mukherjee, A

    Ruchika, P. Mukherjee, and A. Favale, Revisiting Gaus- sian Process Reconstruction for Cosmological Inference: The Generalised GP (Gen GP) Framework (2025), arXiv:2510.03742 [astro-ph.CO]

  15. [15]

    D. J. Bartlett, H. Desmond, P. G. Ferreira, and G. Kro- nberger, Introduction to symbolic regression in the phys- ical sciences (2025), arXiv:2512.15920 [cs.LG]

  16. [16]

    Clarkson, B

    C. Clarkson, B. Bassett, and T. H.-C. Lu, A general test of the Copernican Principle, Phys. Rev. Lett.101, 011301 (2008), arXiv:0712.3457 [astro-ph]

  17. [17]

    Heinesen, Multipole decomposition of the general lu- minosity distance ’Hubble law’ – a new framework for ob- servational cosmology, JCAP05, 008, arXiv:2010.06534 [astro-ph.CO]

    A. Heinesen, Multipole decomposition of the general lu- minosity distance ’Hubble law’ – a new framework for ob- servational cosmology, JCAP05, 008, arXiv:2010.06534 [astro-ph.CO]. 20 https://numpy.org/ 21 https://matplotlib.org/ 22 https://pandas.pydata.org/ 23 https://www.sympy.org/en/index.html 24 https://docs.python.org/3/library/pathlib.html 25 https...

  18. [18]

    Heinesen, Differential age observations and their con- straining power in cosmology, Phys

    A. Heinesen, Differential age observations and their con- straining power in cosmology, Phys. Rev. D111, 023539 (2025), arXiv:2412.05020 [gr-qc]

  19. [19]

    Heinesen, C

    A. Heinesen, C. Blake, and D. L. Wiltshire, Quanti- fying the accuracy of the Alcock–Paczynski scaling of baryon acoustic oscillation measurements, JCAP01, 038, arXiv:1908.11508 [astro-ph.CO]

  20. [20]

    Heinesen, C

    A. Heinesen, C. Blake, Y.-Z. Li, and D. L. Wiltshire, Baryon acoustic oscillation methods for generic curva- ture: application to the SDSS-III Baryon Oscillation Spectroscopic Survey, JCAP03, 003, arXiv:1811.11963 [astro-ph.CO]

  21. [21]

    M. E. Thing and S. M. Koksbang, cp3-bench: a tool for benchmarking symbolic regression algorithms demon- strated with cosmology, JCAP01, 040, arXiv:2406.15531 [astro-ph.IM]

  22. [22]

    Constraining dark matter halo profiles with symbolic regression

    A. Mart´ ın, T. Yasin, D. J. Bartlett, H. Desmond, and P. G. Ferreira, Constraining dark matter halo profiles with symbolic regression (2025), arXiv:2511.23073 [astro- ph.CO]

  23. [23]

    S. M. Koksbang, Machine Learning Cosmic Backreaction and Its Effects on Observations, Phys. Rev. Lett.130, 201003 (2023), arXiv:2305.01224 [astro-ph.CO]

  24. [24]

    S. M. Koksbang, Cosmic backreaction and the mean red- shift drift from symbolic regression, Phys. Rev. D107, 103522 (2023), arXiv:2305.01223 [astro-ph.CO]

  25. [25]

    Desmond, (Exhaustive) Symbolic Regression and model selection by minimum description length (2025), arXiv:2507.13033 [astro-ph.IM]

    H. Desmond, (Exhaustive) Symbolic Regression and model selection by minimum description length (2025), arXiv:2507.13033 [astro-ph.IM]

  26. [26]

    The Pantheon+ Analysis: The Full Dataset and Light-Curve Release

    D. Scolnicet al., The Pantheon+ Analysis: The Full Data Set and Light-curve Release, Astrophys. J.938, 113 (2022), arXiv:2112.03863 [astro-ph.CO]

  27. [27]

    A. G. Adameet al.(DESI), DESI 2024 VI: cosmologi- cal constraints from the measurements of baryon acous- tic oscillations, JCAP02, 021, arXiv:2404.03002 [astro- ph.CO]

  28. [28]

    B. R. Dinda and R. Maartens, Model-agnostic assessment of dark energy after DESI DR1 BAO, JCAP01, 120, arXiv:2407.17252 [astro-ph.CO]

  29. [29]

    DESI DR2 Results II: Measurements of Baryon Acoustic Oscillations and Cosmological Constraints

    M. Abdul Karimet al.(DESI), DESI DR2 results. II. Measurements of baryon acoustic oscillations and cos- mological constraints, Phys. Rev. D112, 083515 (2025), arXiv:2503.14738 [astro-ph.CO]. 19

  30. [30]

    S. M. Koksbang, Searching for signals of inhomogeneity using multiple probes of the cosmic expansion rate H(z), Phys. Rev. Lett.126, 231101 (2021), arXiv:2105.11880 [astro-ph.CO]

  31. [31]

    Camarena and V

    D. Camarena and V. Marra, The tension in the absolute magnitude of Type Ia supernovae (2023), arXiv:2307.02434 [astro-ph.CO]

  32. [32]

    Poulin, T.L

    V. Poulin, T. L. Smith, R. Calder´ on, and T. Si- mon, Implications of the cosmic calibration tension be- yond H0 and the synergy between early- and late- time new physics, Phys. Rev. D111, 083552 (2025), arXiv:2407.18292 [astro-ph.CO]

  33. [33]

    L. H. Dam, A. Heinesen, and D. L. Wiltshire, Appar- ent cosmic acceleration from type Ia supernovae, Mon. Not. Roy. Astron. Soc.472, 835 (2017), arXiv:1706.07236 [astro-ph.CO]

  34. [34]

    McConaghy, Ffx: Fast, scalable, deterministic symbolic regression technology, inGenetic Program- ming Theory and Practice IX, edited by R

    T. McConaghy, Ffx: Fast, scalable, deterministic symbolic regression technology, inGenetic Program- ming Theory and Practice IX, edited by R. Riolo, E. Vladislavleva, and J. H. Moore (Springer New York, New York, NY, 2011) pp. 235–260

  35. [35]

    Virgolin, T

    M. Virgolin, T. Alderliesten, C. Witteveen, and P. A. N. Bosman, Improving model-based genetic programming for symbolic regression of small expressions, Evolutionary Computation29, 211 (2021)

  36. [36]

    Burlacu, G

    B. Burlacu, G. Kronberger, and M. Kommenda, Operon c++: An efficient genetic programming framework for symbolic regression, inProceedings of the 2020 Genetic and Evolutionary Computation Conference Companion, GECCO ’20 (Association for Computing Machinery, New York, NY, USA, 2020) p. 1562–1570

  37. [37]

    Udrescu and M

    S.-M. Udrescu and M. Tegmark, Ai feynman: A physics- inspired method for symbolic regression, Science Ad- vances6, eaay2631 (2020)

  38. [38]

    Udrescu, A

    S.-M. Udrescu, A. Tan, J. Feng, O. Neto, T. Wu, and M. Tegmark, Ai feynman 2.0: Pareto-optimal symbolic regression exploiting graph modularity, arXiv preprint arXiv:2006.10782 (2020)

  39. [40]

    T. N. Mundhenk, M. Landajuela, R. Glatt, C. P. Santi- ago, D. M. Faissol, and B. K. Petersen, Symbolic regres- sion via neural-guided genetic programming population seeding, inAdvances in Neural Information Processing Systems(2021)

  40. [41]

    B. K. Petersen, M. L. Larma, T. N. Mundhenk, C. P. Santiago, S. K. Kim, and J. T. Kim, Deep symbolic re- gression: Recovering mathematical expressions from data via risk-seeking policy gradients, inInternational Confer- ence on Learning Representations(2021)

  41. [42]

    B. K. Petersen, M. Landajuela, T. N. Mundhenk, C. P. Santiago, S. K. Kim, and J. T. Kim, Deep symbolic re- gression: Recovering mathematical expressions from data via risk-seeking policy gradients, inProc. of the Interna- tional Conference on Learning Representations(2021)

  42. [43]

    Landajuela, C

    M. Landajuela, C. Lee, J. Yang, R. Glatt, C. P. Santiago, I. Aravena, T. N. Mundhenk, G. Mulcahy, and B. K. Pe- tersen, A unified framework for deep symbolic regression, inAdvances in Neural Information Processing Systems (2022)

  43. [44]

    Espada, L

    G. Espada, L. Ingelse, P. Canelas, P. Barbosa, and A. Fonseca, Datatypes as a more ergonomic frontend for grammar-guided genetic programming, inGPCE ’22: Concepts and Experiences, Auckland, NZ, December 6 - 7, 2022, edited by B. Scholz and Y. Kameyama (ACM,

  44. [45]

    G. Dick, C. A. Owen, and P. A. Whigham, Feature standardisation and coefficient optimisation for effective symbolic regression, inProceedings of the 2020 Genetic and Evolutionary Computation Conference, GECCO ’20 (Association for Computing Machinery, New York, NY, USA, 2020) p. 306–314

  45. [46]

    F. O. de Franca and G. S. I. Aldeia, Interaction- Transformation Evolutionary Algorithm for Sym- bolic Regression, Evolutionary Computation , 1 (2020), https://direct.mit.edu/evco/article- pdf/doi/10.1162/evco a 00285/1888497/evco a 00285.pdf

  46. [47]

    Interpretable Machine Learning for Science with PySR and SymbolicRegression.jl

    M. Cranmer, Interpretable Machine Learning for Sci- ence with PySR and SymbolicRegression.jl (2023), arXiv:2305.01582

  47. [48]

    N. J. Christensen, S. Demharter, M. Machado, L. Pedersen, M. Salvatore, V. Stentoft-Hansen, and M. T. Iglesias, Identifying interactions in omics data for clinical biomarker discovery using sym- bolic regression, Bioinformatics38, 3749 (2022), https://academic.oup.com/bioinformatics/article- pdf/38/15/3749/49884306/btac405.pdf

  48. [49]

    Zhouet al.(DESI), Astron

    R. Zhouet al.(DESI), Target Selection and Validation of DESI Luminous Red Galaxies, Astron. J.165, 58 (2023), arXiv:2208.08515 [astro-ph.CO]

  49. [50]

    Raichoor, J

    A. Raichooret al., Target Selection and Validation of DESI Emission Line Galaxies, Astron. J.165, 126 (2023), arXiv:2208.08513 [astro-ph.CO]

  50. [51]

    The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample

    S. Alamet al.(BOSS), The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sam- ple, Mon. Not. Roy. Astron. Soc.470, 2617 (2017), arXiv:1607.03155 [astro-ph.CO]

  51. [52]

    S. Alamet al.(eBOSS), Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory, Phys. Rev. D103, 083533 (2021), arXiv:2007.08991 [astro-ph.CO]

  52. [53]

    J. Houet al.(eBOSS), The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from anisotropic clustering analysis of the Quasar Sample in configuration space between redshift 0.8 and 2.2, Mon. Not. Roy. Astron. Soc.500, 1201 (2020), arXiv:2007.08998 [astro-ph.CO]

  53. [54]

    du Mas des Bourboux, J

    H. du Mas des Bourbouxet al.(eBOSS), The Com- pleted SDSS-IV Extended Baryon Oscillation Spectro- scopic Survey: Baryon Acoustic Oscillations with Lyα Forests, Astrophys. J.901, 153 (2020), arXiv:2007.08995 [astro-ph.CO]

  54. [55]

    Efron and R

    B. Efron and R. Tibshirani, THE BOOTSTRAP METHOD FOR ASSESSING ST ATISTICAL ACCU- RACY, Behaviormetrika , 1 (1985)

  55. [56]

    Planck 2018 results. VI. Cosmological parameters

    N. Aghanimet al.(Planck), Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

  56. [57]

    The Pantheon+ Analysis: Cosmological Constraints

    D. Broutet al., The Pantheon+ Analysis: Cosmo- logical Constraints, Astrophys. J.938, 110 (2022), arXiv:2202.04077 [astro-ph.CO]. 20

  57. [58]

    B. R. Dinda, R. Maartens, S. Saito, and C. Clarkson, Improved null tests of ΛCDM and FLRW in light of DESI DR2, JCAP08, 018, arXiv:2504.09681 [astro-ph.CO]

  58. [59]

    Bolejko, Emerging spatial curvature can resolve the tension between high-redshift CMB and low-redshift distance ladder measurements of the Hubble constant, Phys

    K. Bolejko, Emerging spatial curvature can resolve the tension between high-redshift CMB and low-redshift distance ladder measurements of the Hubble constant, Phys. Rev. D97, 103529 (2018), arXiv:1712.02967 [astro- ph.CO]

  59. [60]

    Heinesen and T

    A. Heinesen and T. Buchert, Solving the curvature and Hubble parameter inconsistencies through structure formation-induced curvature, Class. Quant. Grav.37, 164001 (2020), [Erratum: Class.Quant.Grav. 37, 229601 (2020)], arXiv:2002.10831 [gr-qc]

  60. [61]

    Clifton and N

    T. Clifton and N. Hyatt, A radical solution to the Hubble tension problem, JCAP08, 052, arXiv:2404.08586 [astro- ph.CO]

  61. [62]

    R¨ as¨ anen, K

    S. R¨ as¨ anen, K. Bolejko, and A. Finoguenov, New Test of the Friedmann-Lemaˆ ıtre-Robertson-Walker Metric Using the Distance Sum Rule, Phys. Rev. Lett.115, 101301 (2015), arXiv:1412.4976 [astro-ph.CO]

  62. [63]

    C. R. e. a. Harris, Array programming with numpy, Na- ture585, 357–362 (2020)

  63. [64]

    J. D. Hunter, Matplotlib: A 2d graphics environment, Computing in Science & Engineering9, 90 (2007)