REVIEW 2 major objections 4 minor 44 references
A correlated-observable selection bias can shift the measured Hubble constant by about 1 km/s/Mpc, and the data show weak evidence it is there.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 04:45 UTC pith:EYWP6ENC
load-bearing objection A careful, transparent application of a standard selection-correction model to the Cepheid distance ladder; the H0 shift is real in the model but rests on an acknowledged approximation. the 2 major comments →
Selection effects in correlated observations with application to distance-ladder observations
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central result is the selection-corrected likelihood f(y|x<xlim) = phi(z)/(sigma_y Phi(alpha)) * Phi((alpha - rho z)/sqrt(1-rho^2)), where the sample is truncated in x and the measured variable y is correlated with x. In the weak-correlation regime this reduces to a mean shift proportional to rho sigma_y lambda(alpha), with only higher-order skewness. Because each Cepheid has a different photometric uncertainty, the predicted correction scales with the per-star error, which lets the data constrain the selection amplitude at the same time as the distance scale. Fitting this to the R22 Cepheid data with a single selection parameter lowers H0 by -0.7 to -1.1 km/s/Mpc depending on the distan
What carries the argument
Equation (4), the probit-truncated correlated-Normal likelihood. It is the Bayes-theorem product of the marginal Normal for y and the conditional probability of passing the x cut, expressed through the standard Normal density phi and cumulative distribution Phi. The mean-shift approximation (Eq. 6), E[y|x<xlim] = mu_y - rho sigma_y lambda(alpha), is what is implemented in the distance-ladder fits by adding a term s_i sigma_j to each Cepheid magnitude.
Load-bearing premise
Residual Cepheid selection is adequately described by the mean-shift approximation, adding a term s_i sigma_j to the magnitude model — that is, the selection correction scales linearly with each star's photometric error and other terms in the full likelihood can be neglected.
What would settle it
Take the exactly known bivariate-Normal model of Section 2.1, impose a known cut in x with a known rho, then fit the data with the mean-shift correction alone; if the recovered s_i and H0 shift deviate from the analytic predictions by more than the statistical error, the mean-shift approximation or the assumption of constant rho per host is false.
If this is right
- If residual selection exists in this form, each Cepheid's correction scales with its measured photometric error, so heteroscedastic errors partially break the degeneracy between selection and distance.
- Fitting a single shared selection parameter reduces H0 by -0.7 km/s/Mpc (uniform distance-modulus prior) to -1.1 km/s/Mpc (uniform-r^3 prior).
- The combination of the uniform-r^3 distance prior and the shared selection correction yields H0 = 69.8 +/- 1.2 km/s/Mpc, reducing the Hubble tension.
- The 38-parameter per-host selection model overfits: chi^2 drops by 23.2 for 38 extra parameters, so the paper adopts the single-parameter model as the conservative baseline.
- Measures with smaller photometric errors, such as JWST photometry of extragalactic Cepheids, would shrink the predicted size of this class of selection effect.
Where Pith is reading between the lines
- The same likelihood applies to any survey that cuts on one observable and analyzes another correlated one — e.g., SN Ia host selection, weak-lensing shape measurements, or cluster cosmology — so the machinery is a general tool, not only a H0 correction.
- Because the detection significance is 1.2-1.9σ, the true size of the shift, if any, is largely unknown; the paper's H0 shifts should be read as the sensitivity of the current data to this class of effect, not as a measured systematic value.
- A decisive test would be to forward-model the Cepheid selection in the R22 pipeline and check whether the resulting correction scales with sigma_j as the model assumes; if it does not, the fitted s_i would be an artifact.
- The 2-step correction procedure is not equivalent to joint Bayesian inference; the appendix shows it removes the prior-selection coupling, so future analyses that treat selection corrections as fixed should propagate their uncertainty.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives the likelihood for a measurement y when selection is applied to a correlated observable x, f(y|x<x_lim) = φ(z)/(σ_y Φ(α)) Φ((α−ρz)/√(1−ρ²)) (Eq. 4), and notes that in the weak-correlation, mild-truncation regime this reduces to a nearly Gaussian likelihood with a mean shift proportional to the measurement uncertainty (Eq. 6). It then applies this model to the SH0ES/R22 Cepheid distance ladder by adding per-host or shared selection parameters s_i σ_j to the magnitude model. Fitting to the R22 data, the paper reports weak evidence for nonzero residual selection (1.2–1.9σ depending on the distance prior) and a corresponding reduction in H0 of −0.7 to −1.1 km/s/Mpc for a single shared selection parameter. Larger shifts are found for per-host selection parameters, particularly with a uniform-r³ distance prior, but the paper argues these overfit. The paper also examines the impact of alternative distance priors on the Milky Way Cepheid calibration and the extragalactic distance ladder.
Significance. If the empirical result holds, the paper makes a useful contribution to the Hubble-tension literature by identifying a plausible ~1 km/s/Mpc systematic from correlated-observable selection that is not currently included in distance-ladder error budgets. The theoretical derivation in Section 2 is clean, standard (Heckman-type), and clearly presented, and the analytic MAP solution for the augmented linear model (Eq. 16) is elegant, fast, and validated against MCMC. The paper also contains good empirical practice: posterior predictive checks, explicit discussion of prior sensitivity, and a null test showing that the selection-augmented model does not bias H0 when selection is absent. The reported detection significance is honestly low, and the authors are transparent about the model being an effective description rather than a physical pipeline model. The main weaknesses are that the data analysis uses only the mean-shift approximation to Eq. 4, and the checks performed cannot validate that approximation against the full selection likelihood.
major comments (2)
- [§5.1, Eq. (15)] The data fit replaces the full selection likelihood of Eq. (4) with a simple mean shift s_i σ_j added to a Gaussian likelihood, i.e. Eq. (6) with the variance and higher-order terms omitted. This is valid only in the regime ρ≪1, α≫1, where the corrections in Eq. (7) are negligible. However, the paper never demonstrates that the actual Cepheid selection is in this regime; the fitted s_i = ρλ(α) is a degenerate product, so ρ and α are not separately constrained. The posterior predictive tests in §5.4–§5.8 are generated from the same approximate model, so they cannot establish the adequacy of the mean-shift approximation. Since the reported ΔH0 values are conditional on this approximation, this is a load-bearing assumption. I would like to see either a direct fit using the full likelihood of Eq. (4), or a simulation-based sensitivity study showing that the mean-shift model recovers the same
- [§5.7–§5.8] The null tests are described as posterior predictive tests, but the synthetic data are drawn from the same model that is being fit — namely, the approximate mean-shift model. These tests therefore verify internal consistency of the approximate model, not the external validity of the approximation relative to the true selection process. The paper acknowledges this in §6 ('we cannot be certain that the simple model is appropriate'), but the limitation should be stated directly next to the PPT results in §5.7 and §5.8. As written, a reader could easily overinterpret the PPT agreement as evidence that the selection model is correct.
minor comments (4)
- [General] There are several typos: 'paramters' in the Table 2 caption, 'corelation' in §6, 'difficutly' in §6, and an unclosed parenthesis in the introduction ('Section 4.').
- [§4] The text refers to 'H¨og˚as and M¨ortsell (2025)' and 'H¨og¨as and M¨ortsell (2026)' with inconsistent accent rendering; the reference list should be checked for consistency.
- [Table 2] The notation 'individual (42)' versus 'individual' is not defined until later in the text; consider adding a footnote in the table caption clarifying that (42) includes M31, SMC, and LMC.
- [§5.3] The sentence 'The signal-to-noise is then Z=1ᵀCₛ⁻¹ŝ/√(1ᵀCₛ⁻¹1)' uses notation that is not fully defined (e.g., ŝ, Cₛ). It is clear from context, but a brief definition would improve readability.
Circularity Check
No significant circularity: the selection likelihood is an external (Heckman) result, and the selection parameters are fitted nuisance parameters; residual self-citations are auxiliary and not load-bearing.
full rationale
The paper's derivation chain is self-contained on the statistical side: Eq. (4) is re-derived within the paper (and explicitly credited to Heckman 1979), and Eqs. (6)-(7) follow by Taylor expansion, not by assuming the conclusion. The distance-ladder application (§5.1) adds a correction s_i σ_j that corresponds to the combination ρλ(α) from Eq. (6); s_i is a free parameter estimated from the same data and marginalized over, which is standard joint inference rather than a predicted output. The reported ΔH0 values are posterior shifts from fits, not out-of-sample predictions. The posterior predictive tests (§5.4-5.8) are explicitly described as necessary-consistency checks, and the paper states "Posterior predictive agreement is a necessary condition for a good model but not evidence that the model is uniquely correct." The paper also acknowledges the model is "an effective statistical description, not a physical model of the SH0ES pipeline" and that "we cannot be certain that the simple model is appropriate" — these are model-validity limitations, not circular reasoning. Self-citations (Paradiso et al. 2024, Nguyen et al. 2025, Crespi et al. 2025, Krolewski et al. 2025) appear only as auxiliary references for methods or context and do not carry the argument. No step reduces by construction to its input.
Axiom & Free-Parameter Ledger
free parameters (4)
- shared selection coefficient s (or per-host s_i) =
not tabulated directly; implied mean magnitude correction +0.03 to +0.05 mag
- per-host selection coefficients s_i (38 or 42) =
individual values plotted in Fig. 2; mean +0.03 to +0.05 mag
- hierarchical hyperparameters (μ_s, σ_s) =
not reported
- Cepheid intrinsic scatter σ_intr =
0.06 mag (adopted, not fitted)
axioms (7)
- standard math Heckman selection model: bivariate Normal (x,y) with probit threshold on x (Eqs. 2–5)
- ad hoc to paper Mean-shift approximation Eq. 6 replaces full selection likelihood Eq. 4 in the data fit
- domain assumption Constant selection parameters per host (same α, ρ; only s=ρλ(α) identified)
- domain assumption Selection in visible correlates with NIR photometric uncertainty via crowding
- domain assumption R22 data vector and covariance matrix are correct and complete
- domain assumption Uniform-r3 distance prior for extragalactic hosts is an admissible modeling choice
- domain assumption Adopted intrinsic scatter σ_intr = 0.06 mag for Cepheids
read the original abstract
For over a century, following the work of Eddington, Kapteyn, Malmquist and others, astronomers have wrestled with selection biases when making inferences from samples of objects. Typically, selection is performed on the same observations used to make the measurements of interest. However, selecting objects using one observable while analyzing another can also lead to a selection bias when the observables are correlated. Within a Bayesian framework, unmodelled selection effects correspond to a misspecified generative model. We derive the likelihood for truncated selection in correlated observables and demonstrate its usefulness by searching for residual selection effects of this form in Cepheid variable star brightness measurements used in recent distance-ladder measurements of the Hubble constant H0. We specifically look for the form of bias where the selection correction depends on the photometric uncertainties, which the Cepheid data can constrain while simultaneously measuring H0. We find only weak evidence for non-zero residual selection corrections, at a significance of (1.2sigma) to (1.9sigma), depending on the distance prior adopted. Including a single extra parameter to model the unknown cut-off lowers the recovered H0 by -0.7km/s/Mpc to -1.1km/s/Mpc, again depending on the distance prior. Allowing for a different selection for each host galaxy can decrease H0 further, although this becomes very sensitive to the distance prior applied. While introducing a new selection correction cannot by itself explain the Hubble tension, it may be a component of a multi-faceted solution that includes the choice of priors on distances and other effects.
Figures
Reference graph
Works this paper leans on
-
[1]
Optical Identification of Cepheids in 19 Host Galaxies of Type Ia Supernovae and NGC 4258 with the Hubble Space Telescope. , keywords =. doi:10.3847/0004-637X/830/1/10 , archivePrefix =. 1607.08658 , primaryClass =
-
[2]
Bailer-Jones, C. A. L. and Rybizki, J. and Fouesneau, M. and Demleitner, M. and Andrae, R. , title =. The Astronomical Journal , year =
-
[3]
2025 , eprint =
Desmond, Harry and Stiskalek, Richard and Najera, Jose Antonio and Banik, Indranil , title =. 2025 , eprint =
2025
-
[4]
and Casertano, Stefano and Yuan, Wenlong and Bowers, J
Riess, Adam G. and Casertano, Stefano and Yuan, Wenlong and Bowers, J. Bradley and Macri, Lucas and Zinn, Joel C. and Scolnic, Dan , title =. Astrophysical Journal Letters , volume =. 2021 , doi =. 2012.08534 , archivePrefix =
Pith/arXiv arXiv 2021
-
[5]
Gorski, K. M. and Hivon, E. and Banday, A. J. and Wandelt, B. D. and Hansen, F. K. and Reinecke, M. and Bartelmann, M. , title =. The Astrophysical Journal , year =
-
[6]
Fabricius, C. and Luri, X. and Arenou, F. and Babusiaux, C. and Helmi, A. and Muraveva, T. and Reyl. Gaia Early Data Release 3: Catalogue validation , journal =. 2021 , doi =. 2012.06242 , archivePrefix =
Pith/arXiv arXiv 2021
-
[7]
and Lang, Dustin and Goodman, Jonathan , title =
Foreman-Mackey, Daniel and Hogg, David W. and Lang, Dustin and Goodman, Jonathan , title =. Publications of the Astronomical Society of the Pacific , year =. doi:10.1086/670067 , eprint =
-
[8]
Reassessing the Cepheid-based distance ladder: implications for the Hubble constant , journal =
H. Reassessing the Cepheid-based distance ladder: implications for the Hubble constant , journal =. 2025 , doi =. 2412.07840 , archivePrefix =
Pith/arXiv arXiv 2025
-
[9]
Heckman, James J. , title =. Econometrica , year =. doi:10.2307/1912352 , jstor =
-
[10]
Aghanim, N. and others. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020. doi:10.1051/0004-6361/201833910. arXiv:1807.06209
Pith/arXiv arXiv 2018
-
[11]
and others , title =
Riess, Adam G. and others , title =. The Astrophysical Journal Letters , volume =
-
[12]
and others , title =
Wong, Kenneth C. and others , title =. Monthly Notices of the Royal Astronomical Society , volume =
-
[13]
, title =
Verde, Licia and Treu, Tommaso and Riess, Adam G. , title =. Nature Astronomy , volume =
-
[14]
arXiv e-prints , year = 2025, month = nov, eid =
Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth. arXiv e-prints , year = 2025, month = nov, eid =. doi:10.48550/arXiv.2511.23459 , archivePrefix =. 2511.23459 , primaryClass =
-
[15]
arXiv e-prints , year = 2025, month = nov, eid =
A measurement of H_0 from DESI DR1 using energy densities. arXiv e-prints , year = 2025, month = nov, eid =. doi:10.48550/arXiv.2511.23432 , archivePrefix =. 2511.23432 , primaryClass =
-
[16]
and others , title =
Freedman, Wendy L. and others , title =. ApJ , volume =. 2001 , doi =
2001
-
[17]
and others , title =
Riess, Adam G. and others , title =. ApJ , volume =. 2016 , doi =
2016
-
[18]
, title =
Freedman, Wendy L. , title =. ApJ , volume =. 2021 , doi =
2021
-
[19]
and Pickering, Edward C
Leavitt, Henrietta S. and Pickering, Edward C. , title =. Harvard College Observatory Circular , volume =
-
[20]
Gaia Early Data Release 3: Summary of the contents and survey properties , journal =
-
[21]
A distance to the Large Magellanic Cloud that is precise to one per cent , journal =
Pietrzy. A distance to the Large Magellanic Cloud that is precise to one per cent , journal =
-
[22]
and others , title =
Reid, Mark J. and others , title =. ApJ , volume =
-
[23]
New period-luminosity and period-color relations of classical Cepheids. III. Cepheids in SMC. , keywords =. doi:10.1051/0004-6361:200810550 , archivePrefix =. 0810.1780 , primaryClass =
-
[24]
and Madore, Barry F
Freedman, Wendy L. and Madore, Barry F. , title =. ARA&A , volume =
-
[25]
Phillips, M. M. , title =. ApJL , volume =
-
[26]
A&A , volume =
Tripp, Robert , title =. A&A , volume =
-
[27]
2026 , volume =
The Local Distance Network: A community consensus report on the measurement of the Hubble constant at 1\ journal =. 2026 , volume =
2026
-
[28]
Crowded No More: The Accuracy of the Hubble Constant Tested with High-resolution Observations of Cepheids by JWST. , keywords =. doi:10.3847/2041-8213/acf769 , archivePrefix =. 2307.15806 , primaryClass =
Pith/arXiv arXiv 2041
-
[29]
A Redetermination of the Hubble Constant with the Hubble Space Telescope from a Differential Distance Ladder. , keywords =. doi:10.1088/0004-637X/699/1/539 , archivePrefix =. 0905.0695 , primaryClass =
-
[30]
Extracting distribution parameters from multiple uncertain observations with selection biases. , keywords =. doi:10.1093/mnras/stz896 , archivePrefix =. 1809.02063 , primaryClass =
-
[31]
DESI DR2 results. II. Measurements of baryon acoustic oscillations and cosmological constraints , journal =. doi:10.1103/PhysRevD.112.083515 , archivePrefix=. 2503.14738 , primaryClass =
-
[32]
Marcus H. Physically-motivated priors in the local distance ladder significantly reduce the Hubble tension , year =. 2601.22215 , archivePrefix=
-
[33]
LOO-PIT: A sensitive posterior test. , keywords =. doi:10.1088/1475-7516/2025/01/008 , archivePrefix =. 2410.03507 , primaryClass =
Pith/arXiv arXiv 2025
-
[34]
Annual Review of Astronomy and Astrophysics , year =
Teerikorpi, Pekka , title =. Annual Review of Astronomy and Astrophysics , year =
-
[35]
Statistical Analysis and Data Mining: The ASA Data Science Journal , keywords =
Measurement errors and scaling relations in astrophysics: a review. Statistical Analysis and Data Mining: The ASA Data Science Journal , keywords =. doi:10.1002/sam.11173 , archivePrefix =. 1210.6232 , primaryClass =
-
[36]
Coping with selection effects: a Primer on regression with truncated data. , keywords =. doi:10.1093/mnras/stz320 , archivePrefix =. 1901.10522 , primaryClass =
Pith/arXiv arXiv 1901
-
[37]
Evaluating extensions to LCDM: an application of Bayesian model averaging and selection. , keywords =. doi:10.1088/1475-7516/2024/10/021 , archivePrefix =. 2403.02120 , primaryClass =
Pith/arXiv arXiv 2024
-
[38]
Concerning the Classical Cepheid VI _ C Wesenheit Function's Strong Metallicity Dependence. , keywords =. doi:10.1088/2041-8205/741/2/L36 , archivePrefix =. 1110.1629 , primaryClass =
Pith/arXiv arXiv 2041
-
[39]
Measurements of the Hubble Constant with a Two-rung Distance Ladder: Two Out of Three Ain't Bad. , keywords =. doi:10.3847/1538-4357/ac80bd , archivePrefix =. 2204.10866 , primaryClass =
-
[40]
Absolute Calibration of Cepheid Period-Luminosity Relations in NGC 4258. , keywords =. doi:10.3847/1538-4357/ac51db , archivePrefix =. 2203.06681 , primaryClass =
-
[41]
The Perfect Host: JWST Cepheid Observations in a Background-free Type Ia Supernova Host Confirm No Bias in Hubble-constant Measurements. , keywords =. doi:10.3847/2041-8213/ae0ad6 , archivePrefix =. 2509.01667 , primaryClass =
Pith/arXiv arXiv 2041
-
[42]
JWST Validates HST Distance Measurements: Selection of Supernova Subsample Explains Differences in JWST Estimates of Local H _ 0. , keywords =. doi:10.3847/1538-4357/ad8c21 , archivePrefix =. 2408.11770 , primaryClass =
-
[43]
A Lockdown Perspective on the Hubble Tension (with comments from the SH0ES team). arXiv e-prints , keywords =. doi:10.48550/arXiv.2007.10716 , archivePrefix =. 2007.10716 , primaryClass =
-
[44]
H _ 0 revisited. , keywords =. doi:10.1093/mnras/stu278 , archivePrefix =. 1311.3461 , primaryClass =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.