Pith. sign in

REVIEW 4 major objections 4 minor 26 references

Approximation Techniques for the Reconstruction of the Probability Measure and the Coupling Parameters in a Curie-Weiss Model for Large Populations

T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A constant-cost estimator recovers the coupling parameters of a multi-group Curie-Weiss model from sample vote margins, with consistency, asymptotic normality, and large-deviation guarantees.

desk verdict Useful new constant-cost estimator for Curie-Weiss couplings, but Proposition 22 and Theorem 14.4 are wrong as printed; both are fixable, yet the paper needs major revision before acceptance. read the letter →

arxiv 2507.17073 v1 pith:MDYLF7F5 submitted 2025-07-22 math.ST math.PRstat.TH

Miguel Ballesteros , Ivan Naumkin , Gabor Toth This is my paper
classification math.STmath.PRstat.TH MSC 62F1082B2060F0591B12
keywords Curie-Weissmodelcouplingparameterslargepopulationapproximationmaximumlikelihoodestimationdeviationsasymptoticnormalitytwo-tiervotingsystemsoptimalweights
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to estimate the coupling parameters $\beta_\lambda$ of a multi-group Curie-Weiss model---the strengths of within-group interaction in a population of binary voters---without computing the partition function, whose cost grows exponentially with population size. Its solution replaces the exact moment condition of maximum likelihood with sharp large-population approximations, producing an estimator that can be evaluated at constant computational cost for any group size. For parameters away from the critical value $\beta=1$, the estimator converges as the number of samples grows to a bias-corrected value $\tilde{\beta}_N$, which in turn approaches the true parameter as the group grows; it is asymptotically normal and obeys a large deviation principle with exponentially small error probabilities. The authors argue this makes the coupling parameters a practical measure of social cohesion and a usable input for assigning optimal weights in two-tier voting systems.

What carries the argument

The load-bearing object is Proposition 22's asymptotic expansion of the even moments of the group margin: for $\beta<1$, $E S^{2k}/N^k = (1/(1-\beta))^k + O(1/\sqrt{N})$ with an asserted constant $D_{\mathrm{high}}$, and for $\beta>1$, $E S^{2k}/N^{2k} = m(\beta)^{2k} + O((\ln N)^{3/2}/\sqrt{N})$ with an asserted constant $D_{\mathrm{low}}$. These expansions are derived by representing moments as ratios of Laplace integrals and applying saddle-point analysis, and they replace the partition function in the optimality equation $E S^2 = $ sample moment. The estimator is the inverse of the piecewise function $\vartheta^\infty_N$ built from $1/(1-\beta)\cdot(1/N)$ on the high-temperature side and $m(\beta)^2$ on the low-temperature side, with condition (10) guaranteeing a positive gap between the two regimes.

What would settle it

For a fixed $N$, try to instantiate Algorithm 24: condition (10) needs numerical values for $D_{\mathrm{high}}$ and $D_{\mathrm{low}}$, and the paper provides none; demonstrating that no such constants can be extracted from the proof would block the advertised constant-cost procedure. Alternatively, simulate the estimator for a true $\beta$ just inside $I_h$ or $I_l$ with the critical band chosen by any candidate gap and test whether the frequency of `u` responses decays exponentially in $n$ as Proposition 11 predicts; a polynomial decay rate would refute the central claim.

Watch

Extended reading notes

Core claim

The central claim is Theorem 14: for each group with $\beta_\lambda \in I_h \cup I_l$ (the high- and low-temperature intervals separated by an excluded critical band), the estimator $\hat{\beta}^\infty_N$ converges in probability to the deterministic bias-corrected target $\tilde{\beta}_N$ as the number of i.i.d. configurations $n$ tends to infinity; $\tilde{\beta}_N$ converges to the true $\beta_\lambda$ as the group size $N$ tends to infinity; $\sqrt{n}(\hat{\beta}^\infty_N - \tilde{\beta}_N)$ converges in distribution to a centered normal with diagonal covariance; and $\hat{\beta}^\infty_N$ satisfies a large deviation principle whose rate function has its minimum at $\tilde{\beta}_N$, giving exponentially decaying probabilities of large deviation from $\beta$. The estimator is defined by inverting the asymptotic approximations: in the high-temperature regime $T = N/(1-\hat{\beta}^\infty_N)$ and in the low-temperature regime $T = m(\hat{\beta}^\infty_N)^2 N^2$, where $m(\beta)$ is the largest solution of $\tanh(\beta x)=x$; when the sample statistic falls in the middle interval $J_c$ the estimator returns the symbol `u` rather than a number.

Load-bearing premise

The whole procedure presupposes condition (10), which requires positive constants $D_{\mathrm{high}}$ and $D_{\mathrm{low}}$ asserted to exist in Proposition 22 but never quantified, so a user cannot verify the gap between the high- and low-temperature intervals and cannot choose $b_1$ and $b_2$ to run Algorithm 24 or 25.

Editorial extensions

If this is right

  • The coupling parameter of each non-interacting group can be estimated at cost independent of group size, so the method scales to populations of millions once the gap condition (10) holds.
  • Confidence intervals follow from the asymptotic normality: the high-temperature variance tends to $2(1-\beta)^2$ and the low-temperature variance tends to $0$ as $N$ grows.
  • The exponential bounds of Proposition 11 let the estimator double as a classifier: the probability of mistaking a high-temperature group for a low-temperature group decays exponentially in the sample size.
  • Plugging the estimated parameters into Theorem 38 yields an estimator $\hat{w}$ for optimal council weights with the same consistency, normality, and large-deviation properties (Theorem 42).
  • For samples that fall in the critical band, no numeric estimate is reported; the paper's guarantee is only that such samples are exponentially rare for parameters away from $\beta=1$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural next step the paper leaves implicit is a two-stage practical protocol: run $\hat{\beta}^\infty_N$, and when it returns `u`, either enlarge the sample or switch to the exact maximum likelihood estimator, since `u` is informative about proximity to the critical point.
  • Because the method only needs moment asymptotics, the same template should transfer to other mean-field families---e.g., block Ising models with known large-$N$ moment limits---yielding constant-cost estimators for their interaction matrices.
  • The optimal-weight asymptotics imply that in large populations, almost all council weight concentrates in strongly cohesive groups; a testable consequence is that in bodies like a confederal council, estimated optimal weights should order countries by estimated cohesion, not by population alone.
  • Theorem 14's large deviation principle is stated for fixed $N$ as $n$ grows; a double limit in both $n$ and $N$, with the bias $\tilde{\beta}_N-\beta$ controlled, would give a single finite-sample error bound and is not derived in the paper.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper proposes a computationally cheap estimator for the coupling parameters of a multi-group Curie-Weiss model with no inter-group interactions, replacing the exact moment E_{\beta,N}S^2 in the maximum-likelihood equation by large-population asymptotic approximations. The main results are Proposition 11 (exponentially decaying regime-misclassification probabilities), Theorem 14 (consistency to a bias-corrected value, asymptotic normality, and a large deviation principle), and Theorem 38 on asymptotically optimal two-tier voting weights. The central statistical claims are conditional on condition (10), which requires a gap between the high- and low-temperature intervals whose verification depends on constants D_high and D_low asserted in Proposition 22.

Significance. If the technical issues below are fixed, the paper would make a useful contribution: it provides a constant-cost estimator with standard asymptotic guarantees in a model of social/electoral interaction, and it connects those estimates to an explicit voting-weight application. The detailed saddle-point/Laplace proofs in Section 4, the self-contained derivation of the moment asymptotics in Propositions 20 and 22, and the large deviation analysis of S/N in Proposition 30 are genuine strengths. The voting-weight theorem is standard but its proof is helpfully re-derived from the moment asymptotics. However, because the statement of Proposition 22 is incorrect as written and because the rate function in Theorem 14(4) is not the correct contraction, the significance is currently conditional on substantial corrections.

major comments (4)
  1. [Proposition 22] The statement of Proposition 22 is false as written: for beta < 1 it asserts E_{\beta,N} S^{2k} \approx (1/(1-beta))^k N^k, but the proof itself derives the limit (2k-1)!! (1/(1-beta))^k N^k (see the simplified sum over r_1 after Eq. (20)). This is load-bearing, not cosmetic: in the proof of Theorem 14(3)(a), the variance computation V_{\beta,N}(S^2/N) \approx 3(1-\tilde\beta_N)^{-2} - (1-\tilde\beta_N)^{-2} uses the missing factor 3 for k=2, and with the formula as printed the variance would vanish asymptotically. The statement must be corrected and all downstream variance formulas checked accordingly.
  2. [Theorem 14(4), Eq. (44)] The rate function in Eq. (44) is not the contraction of the LDP for T/N^2 through the estimator. By Remark 33, \hat\beta_N^\infty = (\vartheta_N^\infty)^{-1}(T/N^2), so the contraction principle (Theorem 49) gives J(y) = \Lambda^*_{(S/N)^2}((\vartheta_N^\infty)^{-1}(y)) on the two branches of the estimator, with infinity on the gap. The printed expression J(y) = \inf\{\Lambda^*_{(S/N)^2}(x) : x^2 = y\} instead minimizes over x = \pm\sqrt y. For beta in I_l, (\vartheta_N^\infty)^{-1}(y) = m^{-1}(\sqrt y), so the correct rate function has its unique zero at y = m(\tilde\beta_N)^2, whereas the function in (44) has its unique zero at y = m(\tilde\beta_N)^4 because \Lambda^*_{(S/N)^2} vanishes at x = m(\tilde\beta_N)^2. This contradicts the consistency statement in Theorem 14(1) and invalidates the LDP as stated. Please replace (44) with the proper contraction and update Theorem 14(4) and Remark 34 accordingly.
  3. [Definition 7, Remark 8, Algorithms 24-25] The estimator and the definition of the intervals I_h, I_l, J_h, J_l depend on constants D_high and D_low that Proposition 22 only asserts to exist pointwise for each beta; no uniform-in-beta bounds are proved, and no explicit or computable values are supplied. As a result, a user cannot verify condition (10), cannot choose b1 and b2, and cannot certify that the advertised constant-cost algorithm is actually runnable. This is a load-bearing gap in the paper's practical claim. The authors should either prove explicit uniform bounds or reformulate the theorems with a directly checkable condition on N, b1, and b2.
  4. [Theorem 38, proof] In the proof of Theorem 38, the moments of the limiting normal distribution N(0, 1/(1-\beta_\lambda)) are stated as m_k = k!! (1/(1-\beta_\lambda))^{k/2} for even k. The correct formula is (k-1)!! (1/(1-\beta_\lambda))^{k/2}; for k=2 the printed formula gives 2 instead of 1. Since this moment sequence is used to apply the method of moments (Theorem 61), the proof as written is invalid, even though the half-normal limit quoted in the theorem is standard and the result is likely correct. Please correct the moment formula and verify that the method-of-moments hypotheses hold with the corrected sequence.
minor comments (4)
  1. [Definition 10 / Notation 9] The symbol u is introduced as a possible value of the estimator, but the target space ([−∞,∞] ∪ {u})^M is not given a topology, and the LDP in Theorem 14(4) is stated on [−∞,∞]^M; the treatment of the u/gap region in the large deviation statement should be made explicit.
  2. [Proposition 30] The sentence 'I_beta has one minimum at m(beta) = 0 if beta \le 1' is misleading; the minimum is at x = 0, not at the point m(beta), and the phrase should be reworded.
  3. [Figure 1 / Remark 26] Figure 1 is only referenced in Remark 26 and its axes and parameter values are not described; adding a caption with the values of N and the intervals used would improve reproducibility.
  4. [Appendix, Proposition 57] Proposition 57 is stated for a single-group statistic, but it is later used for the multivariate statistic T; the passage from the univariate statement to the coordinate-wise application in Theorem 14 should be spelled out.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the central estimator theorems are derived in-paper from Proposition 22, with only auxiliary self-citations to [2]; the Eq. (44) LDP defect is a correctness error, not a circular reduction.

full rationale

I walked the derivation chain. Proposition 22, the main asymptotic moment approximation, is proved in this paper from Proposition 20 via Hubbard–Stratonovich/Laplace expansions and profile-vector sums, so the consistency and asymptotic-normality claims (Theorem 14.1–14.3) do not reduce to the model's inputs by construction. The estimator is the plug-in inverse of this in-paper approximation; convergence to the bias-corrected value tilde-beta_N and then to beta follows from the WLLN/CLT, Slutsky, the delta method, and Proposition 22, not from a fitted parameter renamed as a prediction. Theorem 14.4 invokes the standard iid-average LDP (Proposition 57) and the contraction principle (Theorem 49). However, the displayed rate function J(y) = inf{Λ*_{(S/N)^2}(x) : x^2 = y} in Eq. (44) is not the contraction through beta-hat-infinity_N = (vartheta-infinity_N)^{-1}(T/N^2); contraction would give J(beta) = Λ*_{(S/N)^2}(vartheta-infinity_N(beta)) on the estimator branches, whose unique zero is at tilde-beta_N. As printed, the infimum over x^2 = y has its unique minimum at y = (E[(S/N)^2])^2 ≈ m(tilde-beta_N)^4, contradicting the estimator's own consistency. This is a substantive mathematical error in the LDP statement, but it is not circularity: the derivation does not make the target result equivalent to its input by construction; it is a wrong transformation. Self-citations to [2] (Proposition 5, Proposition 44, Lemma 48, Lemma 56, Lemma 60) supply supporting technical facts that are not the target results, so they do not make the central claim tautological. The unquantified constants D_high and D_low in condition (10) weaken usability but are not a circular step. Verdict: no significant circularity; score 2 for minor, non-load-bearing self-citation.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central proof adds no fitted constants to data; the ledger contains the user-chosen interval boundaries, the unquantified error constants from Proposition 22, the no-interaction and i.i.d.-sample assumptions, and reliance on companion-paper lemmas.

free parameters (2)
  • b1, b2 = user-specified, 0 < b1 < 1 < b2 satisfying condition (10)
    Define the intervals I_h, I_c, I_l and thereby the estimator in Definition 10. Any values satisfying (10) work in principle, but (10) depends on the unquantified constants D_high and D_low.
  • D_high, D_low
    Existential constants in Proposition 22 bounding the moment approximation errors. Their values are never computed, and they are needed to verify condition (10) and to choose b1, b2 in Algorithms 24 and 25.
assumptions (5)
  • domain assumption The voting population is split into M non-interacting groups; the joint law factorizes over groups (Definition 1).
    All estimation is coordinate-wise. If inter-group interactions exist, the estimator and proofs do not apply.
  • domain assumption Samples consist of n i.i.d. configurations from P_beta,N with fixed group sizes N_lambda.
    The CLT, LLN, and LDP arguments require independent replications of the same population model, an idealization for real election data.
  • domain assumption Cited results from companion paper [2], including Lemma 43, Proposition 44, Lemma 48, Lemma 56, and Lemma 60, are correct.
    These results are invoked without proof in this paper; they are stated as propositions in the authors' related work [2].
  • standard math Standard large deviations theory: Varadhan's lemma, contraction principle, Slutsky, delta method, and continuous mapping (Theorems 49-55).
    Used throughout the proofs of Proposition 11 and Theorem 14; these are classical results treated as background.
  • ad hoc to paper D_high and D_low can be chosen uniformly over the parameter intervals I_h and I_l.
    Definition 7 requires one pair of constants valid for all beta in the respective intervals, but uniformity over beta is not proved explicitly in the paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Approximation Techniques for the Reconstruction of the Probability Measure and the Coupling Parameters in a Curie-Weiss Model for Large Populations." pith.science (2026). https://pith.science/paper/MDYLF7F5

@misc{pith2026250717073,
  author       = {Pith},
  title        = {Pith review of: Approximation Techniques for the Reconstruction of the Probability Measure and the Coupling Parameters in a Curie-Weiss Model for Large Populations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MDYLF7F5}},
  note         = {Machine review of arXiv:2507.17073}
}
read the original abstract

The Curie-Weiss model, originally used to study phase transitions in statistical mechanics, has been adapted to model phenomena in social sciences where many agents interact with each other. Reconstructing the probability measure of a Curie-Weiss model via the maximum likelihood method runs into the problem of computing the partition function which scales exponentially with the population. We study the estimation of the coupling parameters of a multi-group Curie-Weiss model using large population asymptotic approximations for the relevant moments of the probability distribution in the case that there are no interactions between groups. As a result, we obtain an estimator which can be calculated at a low and constant computational cost for any size of the population. The estimator is consistent (under the added assumption that the population is large enough), asymptotically normal, and satisfies large deviation principles. The estimator is potentially useful in political science, sociology, automated voting, and in any application where the degree of social cohesion in a population has to be identified. The Curie-Weiss model's coupling parameters provide a natural measure of social cohesion. We discuss the problem of estimating the optimal weights in two-tier voting systems.

Figures

Figures reproduced from arXiv: 2507.17073 by the authors.

Figure 1
Figure 1. Approximation of Eβ,N S 2 on Ih ∪ Il 33 [PITH_FULL_IMAGE:figures/full_fig_p033_1.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

26 extracted references · 24 canonical work pages

  1. [2]

    Reconstruction of the Probability Measure and the Coupling Parameters in a Curie-Weiss Model

    Miguel Ballesteros, Rams´ es H. Mena, Arno Siri-J´ egousse, and Gabor Toth. Reconstruction of the prob- ability measure and the coupling parameters in a curie-weiss model. 2025. arXiv:2505.21778

    work page Pith review arXiv 2025

  2. [1]

    Detection of an Arbitrary Number of Communities in a Block Spin Ising Model

    Miguel Ballesteros, Rams´ es H. Mena, Jos´ e Luis P´ erez, and Gabor Toth. Detection of an arbitrary number of communities in a block spin ising model. 2023. arXiv:2311.18112

    work page Pith review arXiv 2023

  3. [3]

    Welfarist evaluations of decision rules for boards of representatives

    Claus Beisbart and Luc Bovens. Welfarist evaluations of decision rules for boards of representatives. Soc. Choice Welf., 29(4):581–608, 2007

    work page 2007

  4. [4]

    Exact recovery in the ising blockmodel

    Quentin Berthet, Philippe Rigollet, and Piyush Srivastava. Exact recovery in the ising blockmodel. Ann. Statist., 47(4):1805–1834, 2019

    work page 2019

  5. [5]

    Estimation in spin glasses: A first step

    Sourav Chatterjee. Estimation in spin glasses: A first step. Ann. Stat., 35(5):1931–1946, 2007

    work page 1931

  6. [6]

    Joint parameter estimations for spin glasses

    Wei-Kuo Chen, Arnab Sen, and Qiang Wu. Joint parameter estimations for spin glasses. 2024. arXiv:2406.10760

    arXiv 2024

  7. [7]

    Modelling society with statistical mechanics: An application to cultural contact and immigration

    Pierluigi Contucci and Stefano Ghirlanda. Modelling society with statistical mechanics: An application to cultural contact and immigration. Qual. Quant. , 41:569–578, 2007

    work page 2007

  8. [8]

    Inverse problem robustness for multi-species mean-field spin models

    Micaela Fedele, Cecilia Vernia, and Pierluigi Contucci. Inverse problem robustness for multi-species mean-field spin models. J. Phys. A: Math. Theor. , 46(6), 2013

    work page 2013

Show all 26 references
  1. [9]

    Felsenthal and Mosh´ e Machover

    Dan S. Felsenthal and Mosh´ e Machover. Minimizing the mean majority deficit: The second square-root rule. Math. Soc. Sci. , 37(1):25–37, 1999

    1999

  2. [10]

    Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction

    Sacha Friedli and Yvan Velenik. Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction. Cambridge University Press, 2017

    2017

  3. [11]

    Parameter evaluation of a simple mean-field model of social interaction

    Ignacio Gallo, Adriano Barra, and Pierluigi Contucci. Parameter evaluation of a simple mean-field model of social interaction. Math. Models Methods Appl. Sci. , 19(supp01):1427–1439, 2009

    2009

  4. [12]

    Beitrag zur Theorie des Ferromagnetismus

    Ernst Ising. Beitrag zur Theorie des Ferromagnetismus. Zeitschrift f ¨ur Physik , 31:253–258, 1925

    1925

  5. [13]

    Statistical Problems in Ferromagnetism, Antiferromagnetism and Adsorption

    Pieter Willem Kasteleijn. Statistical Problems in Ferromagnetism, Antiferromagnetism and Adsorption . 1956

    1956

  6. [14]

    On penrose’s square-root law and beyond

    Werner Kirsch. On penrose’s square-root law and beyond. Homo oecon., 24(3/4):357–380, 2007

    2007

  7. [15]

    The Fate of the Square Root Law for Correlated Voting , chapter Voting Power and Procedures, pages 147–158

    Werner Kirsch and Jessica Langner. The Fate of the Square Root Law for Correlated Voting , chapter Voting Power and Procedures, pages 147–158. Springer Cham, 2014. 59

    2014

  8. [16]

    Optimal weights in a two-tier voting system with mean-field voters

    Werner Kirsch and Gabor Toth. Optimal weights in a two-tier voting system with mean-field voters. arXiv:2111.08636, November 2021

    2021

  9. [17]

    Collective bias models in two-tier voting systems and the democracy deficit

    Werner Kirsch and Gabor Toth. Collective bias models in two-tier voting systems and the democracy deficit. Math. Soc. Sci. , 119:118–137, 2022

    2022

  10. [18]

    Optimal apportionment.J

    Yukio Koriyama, Antonin Mac´ e, Rafael Treibich, and Jean-Fran¸ cois Laslier. Optimal apportionment.J. Political Econ., 121(3):584–608, 2013

    2013

  11. [19]

    On the democratic weights of nations

    Sascha Kurz, Nicola Maaser, and Stefan Napel. On the democratic weights of nations. J. Political Econ., 125(5):1599–1634, 2017

    2017

  12. [20]

    Exact recovery in block spin Ising models at the critical line

    Matthias L ¨owe and Kristina Schubert. Exact recovery in block spin Ising models at the critical line. Electron. J. Stat., 14:1796–1815, 2020

    2020

  13. [21]

    A note on the direct democracy deficit in two-tier voting

    Nicola Maaser and Stefan Napel. A note on the direct democracy deficit in two-tier voting. Math. Soc. Sci., 63(2):174–180, 2012

    2012

  14. [22]

    Crystal statistics I

    Lars Onsager. Crystal statistics I. A two-dimensional model with an order-disorder transition. Phys. Rev., 65(3-4):117–149, 1944

    1944

  15. [23]

    A conditional Curie-Weiss model for stylized multi-group binary choice with social interaction

    Alex Akwasi Opoku, Kwame Owusu Edusei, and Richard Kwame Ansah. A conditional Curie-Weiss model for stylized multi-group binary choice with social interaction. J. Stat. Phys. , 2018

    2018

  16. [24]

    The elementary statistics of majority voting

    Lionel Penrose. The elementary statistics of majority voting. Journal of the Royal Statistical Society , 109(1):53–57, 1946

    1946

  17. [25]

    Correlated Voting in Multipopulation Models, Two-Tier Voting Systems, and the Democracy Deficit

    Gabor Toth. Correlated Voting in Multipopulation Models, Two-Tier Voting Systems, and the Democracy Deficit. PhD Thesis, FernUniversit ¨at in Hagen, April 2020. doi:10.18445/20200505-103735-0

    2020 doi

  18. [26]

    Square Root Voting System, Optimal Threshold and π, chapter Voting Power and Procedures, pages 127–146

    Karol Zyczkowski and Wojciech Slomczynski. Square Root Voting System, Optimal Threshold and π, chapter Voting Power and Procedures, pages 127–146. Springer Cham., 2014. 60

    2014

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.