pith. machine review for the scientific record. sign in

arxiv: 2604.07336 · v1 · submitted 2026-04-08 · 🌌 astro-ph.CO · astro-ph.IM· physics.data-an· stat.AP

Recognition: unknown

The Non-Gaussian Weak-Lensing Likelihood: A Multivariate Copula Construction and Impact on Cosmological Constraints

Authors on Pith no claims yet

Pith reviewed 2026-05-10 17:37 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.IMphysics.data-anstat.AP
keywords weak lensingcopula likelihoodnon-Gaussian likelihoodtwo-point correlation functionsS8 parametercosmological constraintsstage-IV surveys
0
0 comments X

The pith

A copula likelihood built from exact marginals and multivariate dependence matches weak-lensing correlation sampling distributions better than Gaussian approximations, shifting S8 by roughly one standard deviation in 1000 deg² surveys.

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

The paper constructs a non-Gaussian likelihood for two-point correlation functions measured in weak-lensing surveys by combining accurate one-dimensional marginal distributions with a dependence structure taken from the exact multivariate case via a copula. This targets the strong non-Gaussianity that appears on large scales, which future surveys will measure with high precision. When the resulting likelihood is used for full cosmological parameter inference, it produces shifts in the parameter S8 of order one standard deviation for survey areas of 1000 square degrees, while the shifts shrink to negligible levels for areas of 10 000 square degrees. The construction is validated by direct comparison to simulated sampling distributions of the data vector and by sampling the full parameter space of current weak-lensing analyses. The authors conclude that standard Gaussian likelihoods remain adequate for stage-IV survey sizes.

Core claim

A multivariate copula likelihood is constructed by inserting exact one-dimensional marginals and a dependence structure derived from the exact multivariate likelihood; this form reproduces the simulated sampling distribution of weak-lensing correlation functions more accurately than a Gaussian likelihood, especially on large scales, and induces S8 shifts of order one standard deviation for 1000 deg² surveys but negligible shifts for 10 000 deg² surveys.

What carries the argument

The copula likelihood, which combines exact one-dimensional marginal distributions with a dependence structure taken from the exact multivariate likelihood to capture non-Gaussianity in the correlation-function data vector.

If this is right

  • S8 parameter shifts reach about one standard deviation when the non-Gaussian copula likelihood is used on 1000 deg² data.
  • The same shifts become negligible once the survey area reaches 10 000 deg².
  • The copula likelihood agrees more closely with simulated sampling distributions than Gaussian forms, especially at large angular scales.
  • Gaussian likelihoods remain sufficient for stage-IV weak-lensing analyses under typical mask geometries.
  • The size of the shifts depends on the detailed mask geometry and the length of the data vector.

Where Pith is reading between the lines

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

  • Smaller or intermediate-area weak-lensing analyses may need non-Gaussian likelihoods to avoid systematic offsets in S8.
  • The method could be tested on other two-point statistics such as galaxy clustering to check whether similar scale-dependent non-Gaussianity appears.
  • The explicit dependence on mask geometry implies that real-data applications require case-by-case validation of the copula construction.
  • Large-scale modes appear to carry enough non-Gaussian information to affect parameter inference even when smaller scales dominate the signal-to-noise.

Load-bearing premise

The dependence structure extracted from the exact multivariate likelihood can be incorporated into the copula without introducing sampling artifacts and remains valid across the mask geometries and data-vector lengths used in the tests.

What would settle it

Generate many independent realizations of the correlation-function data vector for a 1000 deg² survey mask, compute the posterior for S8 under both the copula and Gaussian likelihoods, and check whether the two posteriors differ by approximately one standard deviation in their means.

Figures

Figures reproduced from arXiv: 2604.07336 by Tilman Tr\"oster, Veronika Oehl.

Figure 1
Figure 1. Figure 1: Comparison of two- and one-dimensional marginals of a set of 106 simulations against the analytically predicted copula likelihood. The angular separations correspond to 𝜃¯ 3 = [2.25◦ , 5.00◦ ] and 𝜃¯ 4 = [5.00◦ , 11.12◦ ] and the redshift bins 𝑆3 and 𝑆5 are the same as in the KiDS-1000 analysis. Relative differences are shown in the upper triangle with respect to the maximum value of the copula PDF. using … view at source ↗
Figure 2
Figure 2. Figure 2: Posteriors comparing the use of a fully Gaussian (dashed line) to the copula likelihood (solid line) using a 10 000 deg2 mask and a mock KiDS- 𝜉 + -data vector excluding scales below 15 arcmin. In each case all redshift-bin combinations are included, while different colors indicate the use of different angular bins. The orange lines show the posteriors when all angular bins are included. Single angular bin… view at source ↗
Figure 3
Figure 3. Figure 3: Two-dimensional posteriors in the Ωm-𝑆8 plane using a 1000 deg2 mask (left panel) and a 10 000 deg2 mask (right panel). For both posterios, a full simulated KiDS-𝜉 + -data vector excluding small scales < 30 arcmin has been used and the 68% and 95% confidence levels are shown as contours. Black dashed lines mark the fiducial cosmology, while dots or solid lines indicate the modes. 0.1 0.2 0.3 0.4 0.5 Ωm 0.5… view at source ↗
Figure 4
Figure 4. Figure 4: Two-dimensional posteriors in the Ωm-𝑆8 plane using a 1000 deg2 mask. A simulated KiDS-𝜉 + -data vector is used for both panels, including angular bins down to 5 arcmin (left panel) and down to 30 arcmin (right panel). Contours show the 68% and 95% confidence levels and the fiducial cosmology is indicated by the black dashed lines. Note that the right panel corresponds to the left panel in [PITH_FULL_IMAG… view at source ↗
Figure 5
Figure 5. Figure 5: Two-dimensional posteriors obtained from sampling the Gaussian (pink contours) and copula likelihood (purple contours) in a full weak lensing parameter space. We show marginals in the Ωm-𝑆8 plane. A simulated KiDS- 𝜉 + -data vector including angular bins down to 30 arcmin and a 1000 deg2 mask were used. Contours show the 68% and 95% confidence levels and the fiducial cosmology is indicated by the black das… view at source ↗
read the original abstract

We present a framework to compute non-Gaussian likelihoods for two-point correlation functions. The non-Gaussianity is most pronounced on large scales that will be well-measured by stage-IV weak-lensing surveys. We show how such a multivariate likelihood can be constructed and efficiently evaluated using a copula approach by incorporating exact one-dimensional marginals and a dependence structure derived from the exact multivariate likelihood. The copula likelihood is found to be in better agreement with simulated sampling distributions of correlation functions than Gaussian likelihoods, particularly on large scales. We furthermore investigate the effect of the non-Gaussian copula likelihood on posterior inference, including sampling the full parameter space of contemporary weak-lensing analyses. We find parameter shifts in $S_8$ on the order of one standard deviation for $1 \ 000 \ \mathrm{deg}^2$ surveys but negligible shifts for areas of $10 \ 000 \ \mathrm{deg}^2$, suggesting Gaussian likelihoods are sufficient for stage-IV surveys, though results depend on the detailed mask geometry and data-vector structure.

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 / 2 minor

Summary. The manuscript presents a copula-based framework for constructing non-Gaussian multivariate likelihoods for weak-lensing two-point correlation functions. It combines exact one-dimensional marginal distributions with a dependence structure extracted from the exact multivariate likelihood, yielding a likelihood that is shown to agree better with simulated sampling distributions than the standard Gaussian approximation, particularly on large scales. Cosmological parameter inference using this likelihood produces S8 shifts of order one standard deviation for 1000 deg² surveys but negligible shifts for 10 000 deg² surveys, leading to the conclusion that Gaussian likelihoods remain adequate for stage-IV analyses, although results depend on mask geometry and data-vector structure.

Significance. If the copula construction proves stable and free of artifacts, the work would offer a practical route to incorporating non-Gaussian effects into weak-lensing likelihoods without requiring full simulation-based methods. The reported negligible S8 impact for large survey areas would support continued reliance on Gaussian approximations in stage-IV pipelines (LSST, Euclid), simplifying analyses while flagging the role of mask and binning choices. The framework itself represents a concrete methodological advance, though the moderate quantitative detail provided limits immediate assessment of the magnitude of improvement.

major comments (2)
  1. The central claim that the copula likelihood matches simulated sampling distributions better than Gaussian and produces only negligible S8 shifts for 10 000 deg² surveys rests on the assumption that the dependence structure extracted from the exact multivariate likelihood can be stably embedded without introducing artifacts. The abstract explicitly notes dependence on mask geometry and data-vector structure, yet provides no quantitative tests of stability under changes in angular binning, mask fragmentation, or scaling of survey area and dimensionality. This is load-bearing for the conclusion that Gaussian likelihoods suffice for stage-IV surveys.
  2. The reported S8 shifts (order 1σ for 1000 deg², negligible for 10 000 deg²) and the claim of improved agreement lack accompanying error budgets, details on the number of mocks used to extract the dependence structure or to build the sampling distributions, or explicit validation that the structure generalizes beyond the specific masks and data-vector lengths employed in the inference tests. Without these, the moderate soundness of the central result cannot be fully evaluated.
minor comments (2)
  1. The abstract would benefit from inclusion of at least one concrete quantitative metric (e.g., a Kolmogorov-Smirnov statistic, log-likelihood ratio, or effective degrees of freedom) quantifying the reported better agreement with simulated distributions.
  2. Clarify the dimensionality of the correlation-function data vectors used in the posterior-sampling tests, as high-dimensional cases are precisely where finite-sample bias in rank-based copula parameters is most acute.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the importance of demonstrating robustness in the copula construction. The comments are constructive and have prompted us to strengthen the presentation of our results. We address each major comment point by point below, with revisions incorporated where the manuscript required additional quantitative support.

read point-by-point responses
  1. Referee: The central claim that the copula likelihood matches simulated sampling distributions better than Gaussian and produces only negligible S8 shifts for 10 000 deg² surveys rests on the assumption that the dependence structure extracted from the exact multivariate likelihood can be stably embedded without introducing artifacts. The abstract explicitly notes dependence on mask geometry and data-vector structure, yet provides no quantitative tests of stability under changes in angular binning, mask fragmentation, or scaling of survey area and dimensionality. This is load-bearing for the conclusion that Gaussian likelihoods suffice for stage-IV surveys.

    Authors: We agree that explicit stability tests are necessary to support the central claim. The manuscript already explores two survey areas (1000 and 10 000 deg²) that differ in mask geometry and data-vector dimensionality, providing an initial indication of scaling behavior. To directly address the referee's concern, we have added a dedicated subsection (now Section 4.4) with quantitative tests: we vary angular binning from 8 to 24 bins and introduce controlled mask fragmentation by randomly masking 10–30% of the survey area. The extracted copula dependence parameters (pairwise rank correlations) vary by less than 7% across these configurations, and the improvement in likelihood agreement relative to the Gaussian case remains consistent. We have also clarified the scaling with survey area by showing how the copula structure is tied to the covariance properties that scale predictably with area. These additions make the robustness explicit. revision: yes

  2. Referee: The reported S8 shifts (order 1σ for 1000 deg², negligible for 10 000 deg²) and the claim of improved agreement lack accompanying error budgets, details on the number of mocks used to extract the dependence structure or to build the sampling distributions, or explicit validation that the structure generalizes beyond the specific masks and data-vector lengths employed in the inference tests. Without these, the moderate soundness of the central result cannot be fully evaluated.

    Authors: We acknowledge that the original manuscript provided insufficient detail on these aspects. In the revised version we now state explicitly that 1500 mocks were used to extract the dependence structure and 8000 mocks to construct the empirical sampling distributions for comparison. We have added error budgets to the reported S8 shifts by including the standard deviation obtained from 200 bootstrap resamples of the mock ensemble. In addition, we include a cross-validation test in which the copula is constructed from one mask geometry and then applied to posterior inference on a second, independent mask with different fragmentation and binning; the resulting S8 shifts remain consistent to within 0.4σ. These changes provide the requested quantitative support and validation of generalization. revision: yes

Circularity Check

0 steps flagged

No circularity: copula uses exact marginals and dependence structure as independent inputs

full rationale

The paper's central construction takes one-dimensional marginals and a dependence structure directly from the exact multivariate likelihood as given inputs, then builds and evaluates the copula approximation. Validation against simulated sampling distributions and posterior shifts on S8 are performed as separate tests on mock data. No equation or claim reduces by construction to a fitted parameter renamed as a prediction, nor to a self-citation chain. The result that Gaussian likelihoods suffice for stage-IV surveys follows from applying the constructed likelihood rather than from tautological re-use of the same quantities.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on standard copula theory and the availability of an exact multivariate likelihood from which the dependence structure is taken; no new free parameters or invented physical entities are introduced in the abstract.

axioms (1)
  • standard math Copula functions allow the joint distribution to be constructed from arbitrary marginal distributions and a separate dependence structure.
    Invoked to justify separating marginals from dependence in the likelihood construction.

pith-pipeline@v0.9.0 · 5496 in / 1320 out tokens · 59260 ms · 2026-05-10T17:37:43.146235+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

44 extracted references · 43 canonical work pages · 2 internal anchors

  1. [1]

    2018, PASJ, 70, S4, doi: 10.1093/pasj/psx066

    Aihara, H., Arimoto, N., Armstrong, R., et al. 2018, PASJ, 70, S4, doi:10.1093/pasj/psx066

  2. [2]

    1996, Monthly Notices of the Royal Astronomical Society, 283, 983, doi:10.1093/mnras/283.3.983

    Amendola, L. 1996, Monthly Notices of the Royal Astronomical Society, 283, 983, doi:10.1093/mnras/283.3.983

  3. [3]

    and Gruen, D

    Amon, A., Gruen, D., Troxel, M., et al. 2022, Physical Review D, 105, doi:10.1103/physrevd.105.023514

  4. [4]

    , keywords =

    Asgari, M., Lin, C.-A., Joachimi, B., et al. 2021, A&A, 645, A104, doi:10.1051/0004-6361/202039070

  5. [5]

    H., Read , M

    Benabed, K., Cardoso, J.-F., Prunet, S., & Hivon, E. 2009, Monthly Notices of the Royal Astronomical Society, 400, 219, doi:10.1111/j.1365-2966.2009.15202.x

  6. [6]

    The Astrophysical Journal Supplement Series , author =

    Chisari, N. E., Alonso, D., Krause, E., et al. 2019, The Astrophysical Journal Supplement Series, 242, 2, doi:10.3847/1538-4365/ab1658

  7. [7]

    A., & Simons, F

    Dahlen, F. A., & Simons, F. J. 2008, Geophysical Journal International, 174, 774, doi:10.1111/j.1365-246x.2008.03854.x

  8. [8]

    doi:10.1093/mnras/stw641 , adsurl =

    Dalal, R., Li, X., Nicola, A., et al. 2023, Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Power Spectra. https://arxiv.org/abs/2304.00701 Dark Energy Survey and Kilo-Degree Survey Collaboration, Abbott, T. M. C., Aguena, M., et al. 2023, DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys. https://arxiv.org/abs/2305.171...

  9. [9]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, Publications of the Astronomical Society of the Pacific, 125, 306, doi:10.1086/670067

  10. [10]

    2009, The European Journal of Finance, 15, 609, doi:10.1080/13518470802604457

    Genest, C., Gendron, M., & Bourdeau-Brien, M. 2009, The European Journal of Finance, 15, 609, doi:10.1080/13518470802604457

  11. [11]

    , keywords =

    Giblin, B., Heymans, C., Asgari, M., et al. 2021, A&A, 645, A105, doi:10.1051/0004-6361/202038850 G´orski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759, doi:10.1086/427976 6 https://github.com/ntessore/wigner Non-Gaussian weak lensing likelihood11 101 102 103 Threshold for Gaussian marginals (arcmin) 1326.4 1326.5 1326.6 1326.7 1326.8 Log-L...

  12. [12]

    2022, Phys

    Hall, A., & Taylor, A. 2022, Phys. Rev. D, 105, 123527, doi:10.1103/PhysRevD.105.123527

  13. [13]

    Likelihood Analysis of CMB Temperature and Polarization Power Spectra

    Hamimeche, S., & Lewis, A. 2008, Phys. Rev. D, 77, 103013, doi:10.1103/PhysRevD.77.103013

  14. [14]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi:10.1038/s41586-020-2649-2

  15. [15]

    and Hartlap, J

    Hilbert, S., Hartlap, J., & Schneider, P. 2011, A&A, 536, A85, doi:10.1051/0004-6361/201117294

  16. [16]

    , keywords =

    Hildebrandt, H., van den Busch, J. L., Wright, A. H., et al. 2021, A&A, 647, A124, doi:10.1051/0004-6361/202039018

  17. [17]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi:10.1109/MCSE.2007.55

  18. [18]

    , keywords =

    Joachimi, B., Lin, C.-A., Asgari, M., et al. 2021, A&A, 646, A129, doi:10.1051/0004-6361/202038831

  19. [19]

    2014, Dependence Modeling with Copulas (New York: CRC Press), doi:10.1201/b17116

    Joe, H. 2014, Dependence Modeling with Copulas (New York: CRC Press), doi:10.1201/b17116

  20. [20]

    2011, A&A, 534, A76, doi:10.1051/0004-6361/201117284

    Keitel, D., & Schneider, P. 2011, A&A, 534, A76, doi:10.1051/0004-6361/201117284

  21. [21]

    2019, A&A, 625, A2, doi:10.1051/0004-6361/201834918

    Kuijken, K., Heymans, C., Dvornik, A., et al. 2019, A&A, 625, A2, doi:10.1051/0004-6361/201834918

  22. [22]

    Euclid Definition Study Report

    Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, Euclid Definition Study Report.https://arxiv.org/abs/1110.3193

  23. [23]

    and Kannawadi, Arun and Reischke, Robert and Shan, HuanYuan , year = 2023, month = nov, journal =

    Li, S.-S., Hoekstra, H., Kuijken, K., et al. 2023a, Astronomy & Astrophysics, 679, A133, doi:10.1051/0004-6361/202347236

  24. [24]

    2023b, Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Two-point Correlation Functions

    Li, X., Zhang, T., Sugiyama, S., et al. 2023b, Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Two-point Correlation Functions. https://arxiv.org/abs/2304.00702

  25. [25]

    2020, Monthly Notices of the Royal Astronomical Society, 499, 2977, doi:10.1093/mnras/staa2948 LSST Science Collaboration, Abell, P

    Lin, C.-H., Harnois-D´eraps, J., Eifler, T., et al. 2020, Monthly Notices of the Royal Astronomical Society, 499, 2977, doi:10.1093/mnras/staa2948 LSST Science Collaboration, Abell, P. A., Allison, J., et al. 2009, LSST Science Book, Version 2.0.https://arxiv.org/abs/0912.0201

  26. [26]

    2010, Journal of Financial Econometrics, 8, 511, doi:10.1093/jfinec/nbp031

    Min, A., & Czado, C. 2010, Journal of Financial Econometrics, 8, 511, doi:10.1093/jfinec/nbp031

  27. [27]

    Nelsen, R. B. 1999, Lecture Notes in Statistics, Vol. 139, An introduction to copulas, 1st edn. (New York, NY: Springer)

  28. [28]

    1989, Journal of the American Statistical Association, 84, 487, doi:10.1080/01621459.1989.10478795

    Oakes, D. 1989, Journal of the American Statistical Association, 84, 487, doi:10.1080/01621459.1989.10478795

  29. [29]

    2025, The Open Journal of Astrophysics, 8, doi:10.33232/001c.144235

    Oehl, V., & Tr¨oster, T. 2025, The Open Journal of Astrophysics, 8, doi:10.33232/001c.144235

  30. [30]

    Sato, M., Ichiki, K., & Takeuchi, T. T. 2010, Phys. Rev. Lett., 105, 251301, doi:10.1103/PhysRevLett.105.251301 —. 2011, Physical Review D, 83, doi:10.1103/physrevd.83.023501

  31. [31]

    Schulten, K., & Gordon, R. G. 1975, Journal of Mathematical Physics, 16, 1961, doi:10.1063/1.522426

  32. [32]

    , keywords =

    Secco, L., Samuroff, S., Krause, E., et al. 2022, Physical Review D, 105, doi:10.1103/physrevd.105.023515

  33. [33]

    Sellentin, E., & Heavens, A. F. 2017, Monthly Notices of the Royal Astronomical Society, 473, 2355, doi:10.1093/mnras/stx2491

  34. [34]

    2018, MNRAS, 477, 4879, doi:10.1093/mnras/sty988

    Sellentin, E., Heymans, C., & Harnois-D´eraps, J. 2018, MNRAS, 477, 4879, doi:10.1093/mnras/sty988

  35. [35]

    Smith, M. S. 2013, in Bayesian Theory and Applications (Oxford University Press), doi:10.1093/acprof:oso/9780199695607.003.0017

  36. [36]

    Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report

    Spergel, D., Gehrels, N., Baltay, C., et al. 2015, arXiv e-prints, arXiv:1503.03757, doi:10.48550/arXiv.1503.03757 St¨olzner, B., Wright, A. H., Asgari, M., et al. 2025, A&A, 702, A169, doi:10.1051/0004-6361/202554893

  37. [37]

    Tessore, A

    Jeffrey, N. 2023, The Open Journal of Astrophysics, 6, doi:10.21105/astro.2302.01942

  38. [38]

    K., & Zimmer, D

    Trivedl, P. K., & Zimmer, D. M. 2007, Foundations and Trends in Econometrics, 1, 1, doi:10.1561/0800000005

  39. [39]

    Uhlemann, O

    Uhlemann, C., Friedrich, O., Boyle, A., et al. 2023, The Open Journal of Astrophysics, 6, doi:10.21105/astro.2210.07819

  40. [40]

    E., Whittaker, L., & Brown, M

    Upham, R. E., Whittaker, L., & Brown, M. L. 2020, MNRAS, 491, 3165, doi:10.1093/mnras/stz3225

  41. [41]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi:10.1038/s41592-019-0686-2

  42. [42]

    2013, A&A, 556, A70, doi:10.1051/0004-6361/201321718

    Wilking, P., & Schneider, P. 2013, A&A, 556, A70, doi:10.1051/0004-6361/201321718

  43. [43]

    doi:10.1051/0004-6361/202554908 , archiveprefix =

    Wright, A. H., St¨olzner, B., Asgari, M., et al. 2025, A&A, 703, A158, doi:10.1051/0004-6361/202554908

  44. [44]

    Zonca, A., Singer, L., Lenz, D., et al. 2019, Journal of Open Source Software, 4, 1298, doi:10.21105/joss.01298 APPENDIX A.CHOICE OF THE THRESHOLD FOR GAUSSIAN MARGINALS The copula approach allows to employ Gaussian marginals next to the computationally expensive exact marginals wherever appropriate. The non-Gaussianity of the likelihood is most prominent...