Methodological refinement of the submillimeter galaxy cross-correlation function measurements and their uncertainty estimation
Pith reviewed 2026-05-24 09:20 UTC · model grok-4.3
The pith
New full-field pair counting and oversampled bootstrap covariance method produces robust cross-correlation measurements compatible between spectroscopic and photometric lens samples.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors present a new methodological framework that counts pairs using the full field area and estimates the covariance matrix internally via an oversampled bootstrap method on k-means defined patches. This yields cross-correlation function measurements that are compatible between spectroscopic and photometric lens samples. Analysis of the three GAMA fields shows that the G15 field has a stronger signal produced by the rare combination of large-scale structure excesses in both foreground and background samples. The results indicate that differences from the previously used mini-tile approach arise from physical properties of the samples.
What carries the argument
Full-field pair counting combined with internal oversampled bootstrap covariance estimation on automatically k-means clustered patches.
If this is right
- Cross-correlation measurements from spectroscopic and photometric lens samples become compatible.
- The G15 field shows a stronger signal attributable to combined large-scale structure excesses in foreground and background samples.
- Differences from earlier mini-tile methods are due to sample properties rather than the estimation technique.
- Using the full field area reduces statistical uncertainty by incorporating all available data.
Where Pith is reading between the lines
- The method could be tested on additional survey fields to check whether similar anomalies appear and affect parameter constraints.
- Compatibility across lens types suggests future analyses can mix sample types without introducing systematic offsets from the measurement process.
- Further study of the G15 anomaly may reveal details about how rare large-scale structure alignments influence magnification bias signals.
- The internal covariance approach might extend to auto-correlation measurements in comparable multi-field datasets.
Load-bearing premise
The automatic k-means division into at least five patches and the oversampled bootstrap procedure correctly capture the true covariance without bias from the choice of patch number or resampling details.
What would settle it
Repeating the analysis with a different number of k-means patches or a standard jackknife covariance estimator produces results that are incompatible between the spectroscopic and photometric lens samples.
Figures
read the original abstract
In this study, we aim to develop a new methodology to estimate the cross-correlation function and uncertainties and apply it to the analysis of magnification bias in galaxy surveys. We adopt a new methodological framework that uses a statistically rigorous approach to obtain more robust measurements for constraining cosmological parameters. This strategy involves using the full field area to count the number of different pairs for each field and combine them into a single estimation, reducing statistical uncertainty and accounting for the full information available in the data. The covariance matrix was estimated internally using an oversampled bootstrap method. We divided each field into at least five patches, that were defined automatically using a k-mean clustering algorithm. We investigate the robustness of the new methodology by comparing the results from a spectroscopic lens sample with those from a photometric lens sample, finding them to be compatible. We also analyse the cross-correlation function and auto-correlation function for individual fields in the three GAMA fields, comparing both samples. The G15 field was found to have a stronger signal compared to the other fields, suggesting that the stronger cross-correlation is produced by the rare combination of two excesses of large-scale structure in both the foreground and background samples. Our results demonstrate the robustness of the new methodology and suggest that the differences with respect to the mini-tile previously used approached may be due to physical properties of the samples themselves. The identified G15 anomalous signal warrants further investigation into its impact on cosmological parameters.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a refined methodology for measuring the submillimeter galaxy cross-correlation function by counting pairs over the full field area and estimating the covariance matrix via an oversampled bootstrap with automatic k-means patch division (at least five patches per field). Applied to magnification bias analysis in GAMA fields using spectroscopic and photometric lens samples, it reports compatible results between the two samples, attributes the stronger G15 signal to physical large-scale structure excesses rather than systematics, and claims the approach is more robust than prior mini-tile methods.
Significance. If the internal covariance estimation proves unbiased and the sample compatibility reflects genuine robustness rather than shared procedural artifacts, the method could reduce statistical uncertainties in cross-correlation measurements and improve cosmological constraints from galaxy surveys. The data-driven nature of the bootstrap is a potential strength, but the absence of quantitative validation or external benchmarks limits the assessed impact on the field.
major comments (2)
- [Abstract] The central robustness claim rests on compatibility between spectroscopic and photometric lens samples, yet the abstract (and available text) provides no quantitative measures such as chi-squared, p-values, or explicit error budgets to support this; without these, it is impossible to evaluate whether differences from prior methods arise from sample properties or unaccounted covariance biases.
- [Methodology (bootstrap and patch division)] The oversampled bootstrap covariance estimation with k-means patches is load-bearing for the uncertainty quantification and compatibility claims, but the description does not include tests varying patch number, comparisons to jackknife or simulated covariances, or checks against field geometry (e.g., equal-area or contiguous patches); k-means clustering risks correlating patches with large-scale structure, potentially underestimating cosmic variance as hinted by the G15 excess.
minor comments (2)
- [Abstract] The abstract would benefit from including at least one numerical result (e.g., measured correlation amplitude or uncertainty reduction factor) to ground the claims of reduced statistical uncertainty.
- [Introduction/Methods] Notation for the cross-correlation function and covariance matrix should be defined explicitly with equations early in the text to aid readability.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive comments. We address each major comment point by point below, indicating where revisions have been made to the manuscript.
read point-by-point responses
-
Referee: [Abstract] The central robustness claim rests on compatibility between spectroscopic and photometric lens samples, yet the abstract (and available text) provides no quantitative measures such as chi-squared, p-values, or explicit error budgets to support this; without these, it is impossible to evaluate whether differences from prior methods arise from sample properties or unaccounted covariance biases.
Authors: We agree that quantitative measures strengthen the presentation of the compatibility claim. The revised manuscript updates the abstract to report the chi-squared value and p-value for the spectroscopic versus photometric sample comparison, and the results section now includes an explicit discussion of the error budget arising from the covariance estimation. These additions allow direct evaluation of whether observed differences are consistent with sample properties. revision: yes
-
Referee: [Methodology (bootstrap and patch division)] The oversampled bootstrap covariance estimation with k-means patches is load-bearing for the uncertainty quantification and compatibility claims, but the description does not include tests varying patch number, comparisons to jackknife or simulated covariances, or checks against field geometry (e.g., equal-area or contiguous patches); k-means clustering risks correlating patches with large-scale structure, potentially underestimating cosmic variance as hinted by the G15 excess.
Authors: We have added text clarifying that k-means was selected to produce patches of comparable area without manual intervention and that the minimum of five patches per field was chosen to ensure sufficient resampling while respecting field boundaries. A short sensitivity discussion on patch number has been included. However, systematic variation of patch number, jackknife comparisons, and external simulated covariance benchmarks are not performed here, as the paper focuses on the internal data-driven approach and the consistency check provided by the two independent lens samples. The G15 excess is supported by the auto-correlation measurements in both samples, which we interpret as physical large-scale structure rather than a covariance artifact. revision: partial
Circularity Check
No significant circularity; methodology is data-driven and self-contained
full rationale
The paper derives cross-correlation measurements and covariance estimates directly from the survey data via full-field pair counting and oversampled bootstrap on the same fields. No equations reduce the reported results to fitted parameters by construction, no self-citation chains justify the central claims, and the compatibility between lens samples is presented as an empirical outcome rather than a definitional or fitted equivalence. The derivation remains independent of its inputs.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We adopt a new methodology that uses ... full field area to count ... covariance matrix was estimated internally using an oversampled bootstrap method. We divided each field into at least five patches ... k-mean clustering algorithm.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The IC can be estimated ... using the formula IC = sum RfRb(θi) wideal(θi) / sum RfRb(θi)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
2003, Journal of Computer and System Sciences, 66, 671, special Issue on PODS 2001
Achlioptas, D. 2003, Journal of Computer and System Sciences, 66, 671, special Issue on PODS 2001
work page 2003
-
[2]
Adelberger, K. L., Steidel, C. C., Pettini, M., et al. 2005, ApJ, 619, 697
work page 2005
-
[3]
2019, arXiv e-prints, arXiv:1912.02905
Ahumada, R., Allende Prieto, C., Almeida, A., et al. 2019, arXiv e-prints, arXiv:1912.02905
-
[4]
2019, Monthly Notices of the Royal Astronomical Society, 483, 4649
Amvrosiadis, A., Valiante, E., Gonzalez-Nuevo, J., et al. 2019, Monthly Notices of the Royal Astronomical Society, 483, 4649
work page 2019
-
[5]
Bakx, T. J. L. C., Eales, S., & Amvrosiadis, A. 2020, MNRAS, 493, 4276
work page 2020
-
[6]
Bakx, T. J. L. C., Eales, S. A., Negrello, M., et al. 2018, MNRAS, 473, 1751
work page 2018
-
[7]
Baldry, I. K., Alpaslan, M., Bauer, A. E., et al. 2014, MNRAS, 441, 2440
work page 2014
-
[8]
Baldry, I. K., Robotham, A. S. G., Hill, D. T., et al. 2010, MNRAS, 404, 86
work page 2010
- [9]
- [10]
-
[11]
Blain, A. W. 1996, MNRAS, 283, 1340
work page 1996
-
[12]
Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28
work page 2017
-
[13]
Bonavera, L., Cueli, M. M., & Gonzalez-Nuevo, J. 2022, Proceedings of the MG16 Meeting on General Relativity, R. Ru ffini & G. Vereshchagin eds., World Scientific., arXiv:2112.02959
-
[14]
Bonavera, L., Cueli, M. M., González-Nuevo, J., et al. 2021, A&A, 656, A99
work page 2021
-
[15]
M., González-Nuevo, J., Casas, J
Bonavera, L., Cueli, M. M., González-Nuevo, J., Casas, J. M., & Crespo, D. 2024, arXiv e-prints, arXiv:2305.13836
- [16]
- [17]
- [18]
- [19]
-
[20]
S., Pérez-Fournon, I., Amber, S., et al
Bussmann, R. S., Pérez-Fournon, I., Amber, S., et al. 2013, ApJ, 779, 25
work page 2013
-
[21]
Calanog, J. A., Fu, H., Cooray, A., et al. 2014, ApJ, 797, 138
work page 2014
-
[22]
Crespo, D., González-Nuevo, J., Bonavera, L., et al. 2022, A&A, 667, A146
work page 2022
-
[23]
M., Bonavera, L., González-Nuevo, J., et al
Cueli, M. M., Bonavera, L., González-Nuevo, J., et al. 2022, A&A, 662, A44
work page 2022
-
[24]
M., Bonavera, L., González-Nuevo, J., & Lapi, A
Cueli, M. M., Bonavera, L., González-Nuevo, J., & Lapi, A. 2021, A&A, 645, A126
work page 2021
-
[25]
M., González-Nuevo, J., Bonavera, L., et al
Cueli, M. M., González-Nuevo, J., Bonavera, L., et al. 2024, arXiv e-prints, arXiv:2305.13835
- [26]
- [27]
- [28]
- [29]
-
[30]
Fernandez, L., Cueli, M. M., González-Nuevo, J., et al. 2022, A&A, 658, A19
work page 2022
-
[31]
Friedrich, O., Seitz, S., Eifler, T. F., & Gruen, D. 2016, MNRAS, 456, 2662
work page 2016
-
[32]
2012, ApJ, 753, 134 González-Nuevo, J., Cueli, M
Fu, H., Jullo, E., Cooray, A., et al. 2012, ApJ, 753, 134 González-Nuevo, J., Cueli, M. M., Bonavera, L., et al. 2021, A&A, 646, A152 González-Nuevo, J., Lapi, A., Bonavera, L., et al. 2017, J. Cosmology Astropart. Phys., 2017, 024 González-Nuevo, J., Lapi, A., Fleuren, S., et al. 2012, ApJ, 749, 65 González-Nuevo, J., Lapi, A., Negrello, M., et al. 2014,...
work page 2012
-
[33]
J., Abergel, A., Abreu, A., et al
Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3
work page 2010
-
[34]
2001, in Cosmological Physics with Gravitational Lensing, ed
Herranz, D. 2001, in Cosmological Physics with Gravitational Lensing, ed. J. Tran Thanh Van, Y . Mellier, & M. Moniez, 197
work page 2001
-
[35]
Hildebrandt, H., van Waerbeke, L., Scott, D., et al. 2013, MNRAS, 429, 3230
work page 2013
-
[36]
Hunter, J. D. 2007, Computing In Science & Engineering, 9, 90
work page 2007
- [37]
-
[38]
Jarvis, M. 2015, TreeCorr: Two-point correlation functions, Astrophysics Source Code Library, record ascl:1508.007
work page 2015
-
[39]
2001, SciPy: Open source scientific tools for Python
Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open source scientific tools for Python
work page 2001
-
[40]
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64
work page 1993
-
[41]
Liske, J., Baldry, I. K., Driver, S. P., et al. 2015, Monthly Notices of the Royal Astronomical Society, 452, 2087
work page 2015
-
[42]
Maddox, S. J. & Dunne, L. 2020, MNRAS, 493, 2363
work page 2020
-
[43]
Marinoni, C., Bel, J., & Buzzi, A. 2012, J. Cosmology Astropart. Phys., 2012, 036 Ménard, B., Scranton, R., Fukugita, M., & Richards, G. 2010, MNRAS, 405, 1025
work page 2012
-
[44]
2016, ApJ, 823, 17 Article number, page 13 of 14 A&A proofs: manuscript no
Nayyeri, H., Keele, M., Cooray, A., et al. 2016, ApJ, 823, 17 Article number, page 13 of 14 A&A proofs: manuscript no. main
work page 2016
-
[45]
Negrello, M., Amber, S., Amvrosiadis, A., et al. 2017, MNRAS, 465, 3558
work page 2017
-
[46]
Negrello, M., Hopwood, R., De Zotti, G., et al. 2010, Science, 330, 800
work page 2010
-
[47]
Negrello, M., Perrotta, F., González-Nuevo, J., et al. 2007, MNRAS, 377, 1557
work page 2007
-
[48]
M., Gaztañaga, E., & Croton, D
Norberg, P., Baugh, C. M., Gaztañaga, E., & Croton, D. J. 2009, MNRAS, 396, 19
work page 2009
-
[49]
J., Bock, J., Altieri, B., et al
Oliver, S. J., Bock, J., Altieri, B., et al. 2012, MNRAS, 424, 1614
work page 2012
-
[50]
2011, MNRAS, 415, 911 Pérez, F
Pascale, E., Auld, R., Dariush, A., et al. 2011, MNRAS, 415, 911 Pérez, F. & Granger, B. E. 2007, Computing in Science and Engineering, 9, 21
work page 2011
-
[51]
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 Planck Collaboration, Aghanim, N., Akrami, Y ., et al. 2020, A&A, 641, A6
work page 2010
- [52]
- [53]
- [54]
-
[55]
Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses
work page 1992
- [56]
-
[57]
1986, The Annals of Statistics, 14, 1322
Shao, J. 1986, The Annals of Statistics, 14, 1322
work page 1986
-
[58]
2021, A&A, 653, A35 The Dark Energy Survey Collaboration
Sureshkumar, U., Durkalec, A., Pollo, A., et al. 2021, A&A, 653, A35 The Dark Energy Survey Collaboration. 2005, arXiv e-prints, astro
work page 2021
-
[59]
Valiante, E., Smith, M. W. L., Eales, S., et al. 2016, MNRAS, 462, 3146
work page 2016
-
[60]
L., Cooray, A., De Bernardis, F., et al
Wardlow, J. L., Cooray, A., De Bernardis, F., et al. 2013, ApJ, 762, 59
work page 2013
- [61]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.