{"id":"ebb378ce-a77e-4902-8208-ef2448d3e2d3","arxiv_id":"2607.29619","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For cubic-tensor Gibbs measures, joint pseudolikelihood estimation of the coupling and field parameters is √N-consistent under explicit inhomogeneity conditions and provably ill-conditioned in homogeneous ferromagnetic regimes.","lead":"This paper proves when two parameters of a network model with triple interactions can be estimated together from one observed network using pseudolikelihood: roughly, it works when local statistics differ enough across nodes, and fails in homogeneous ferromagnetic regimes. It gives the first rigorous treatment for cubic tensor Gibbs measures, with concrete results for edge-triangle and edge-three-star ERGMs, arithmetic-progression models, and random hypergraphs.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the central theorems are internally consistent; the restrictive assumptions are verified in all applications, and the unproved impossibility remarks are peripheral.","rationale":"The reader's verdict is CONDITIONAL, citing the unproved impossibility remarks and the restrictive assumptions behind Lemmas 3.4(b)/3.5. I agree these are worth flagging, but I do not agree that they are load-bearing for the central claim. The main theorems explicitly assume the conditions that are pointed out, and every application verifies them; the theorems are not claimed for all tensors. The unproved remarks are peripheral and do not support any theorem. The proof of Theorem 1.3 is carefully executed, including the delicate existence argument via Lemma A.2 and the Hessian lower bound; Theorem 1.4's use of Lemma 4.3 is non-circular and the Q_t large-deviation bound is valid; Theorem 1.13's variational argument, including the use of Lemma 3.8 and Lemma 5.5, checks out. I therefore see no reason to change the reader's CONDITIONAL verdict: it is appropriately cautious because of the unproved remarks, but no further adjustment is needed.","tokens_in":48579,"tokens_out":35360,"duration_ms":300358,"concrete_test":"Independently prove the contiguity-based impossibility claim sketched in Remark 1.14 for the complete tensor A_{ijk}=N^{-2}1_{i≠j≠k}, following the argument of [18, Theorem 1.6]. If this extension succeeds, the remarks become valid; if it fails, relabel Remarks 1.14 and 2.6 as conjectures. This check directly addresses the only genuinely unproved assertions, while the main theorems remain unaffected either way.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After a careful pass through the main argument, I find no internal inconsistency or missing step that threatens the paper's central claim. The reader's identified weakest assumption—the use of strong pseudo-regularity (6) and non-degeneracy (7) in Lemma 3.4(b) and Lemma 3.5—is indeed a limitation: the mean-field condition (4) and spectral-gap condition (10) are verified in applications via these lemmas, so tensors with highly heterogeneous positive pair-codegrees or many zero rows are not covered. However, this is a boundary of the stated theorems, not a flaw in them: the theorems explicitly assume (6)–(7), and all four applications in Section 2 verify these conditions (edge-triangle and three-star with K=1, 3-AP with K=4, inhomogeneous hypergraphs with high probability). The omitted proofs that do exist are in Remarks 1.14 and 2.6, where impossibility for any estimator is asserted by 'straightforward extension' of [18, Theorem 1.6] without proof. These remarks are not used in any theorem's proof, so they do not affect the central consistency/ill-conditioning dichotomy. I also checked several delicate algebraic steps, including the identity 1ᵀ(ΣAᵢ²)1=Tr(R²) in Lemma 3.4(a) and the coefficient 1/2 in equation (58), and they are correct.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies joint estimation of the parameters (β,h) of a high-dimensional Gibbs measure dP_{β,h} ∝ exp(β/3 ⟨A,x^{⊗3}⟩ + h⟨x,1⟩) dµ^{⊗N} from a single sample X, using the maximum pseudolikelihood estimator (MPLE). The central result (Theorem 1.3) asserts that if the empirical variance T_N(X) of the local fields is ω_p(N^{-2/7}), then the MPLE exists with probability tending to one and satisfies max(|β̂−β0|,|ĥ−h0|)=O_p(N^{-1/2}T_N^{-1}); in particular, when T_N=Ω_p(1), the estimator is √N-consistent. Theorem 1.4 gives checkable tensor conditions (Tr(R²)=Ω(N) or nonconstant row sums of R) that force T_N=Ω_p(1) when Λ'(h0)≠0. A mean-field analysis (Theorems 1.10, 1.11, 1.13) gives complementary sufficient conditions for consistency and for asymptotic ill-conditioning (T_N=o_p(1)) in homogeneous ferromagnetic regimes. The general results are applied to edge-triangle and edge-three-star ERGMs, cyclic/integer 3-AP models, and inhomogeneous random hypergraphs, yielding a dichotomy: ill-conditioning in homogeneous/ferromagnetic regimes and √N-consistency in strongly antiferromagnetic or heterogeneous regimes.","tokens_in":48871,"tokens_out":27756,"duration_ms":247464,"significance":"If the main theorems are correct, the paper gives a fairly complete qualitative account of joint pseudolikelihood estimation for cubic tensor Gibbs measures: heterogeneity of local fields, measured by T_N, is the key mechanism, and ill-conditioning in homogeneous ferromagnetic regimes is tied to explicit mean-field, regularity, and spectral-gap conditions. The proofs are detailed and mostly self-contained, and the sufficient conditions in Theorems 1.4 and 1.11 are explicit functions of the model, not fitted to data; the applications in Section 2 verify the structural assumptions rather than assuming them. I found no circularity or post-hoc parameter fitting. The main limitations are genuine boundaries of the stated results: the strong pseudo-regularity condition (6) and non-degeneracy condition (7) are used in Lemmas 3.4(b) and 3.5 to verify the mean-field and spectral-gap conditions, and the impossibility remarks in Remarks 1.14 and 2.6 are asserted without proof. These limitations do not affect the central consistency/ill-conditioning theorems as stated.","major_comments":[{"comment":"The rate claimed in Theorem 1.3 is O_p(N^{-1/2}T_N^{-1}), but the proof derives only min(Y_N,r) ≤ c_N/(2η√N T_N) with probability tending to one for an arbitrary diverging sequence c_N. Choosing c_N=N^{1/5} gives min(Y_N,r)=O_p(N^{-3/10}T_N^{-1}), which is weaker than the stated O_p(N^{-1/2}T_N^{-1}). The sentence 'for any diverging sequence ... thus giving Y_N=O_p(1/(√N T_N))' is not valid as written. The argument can be repaired locally by using the defining property of O_p(√N Y_N) with a fixed constant M_ε and probability 1−ε, but as it stands the proof does not deliver the central √N-consistency conclusion.","section":"Appendix A, proof of Theorem 1.3 (consistency part)"}],"minor_comments":[{"comment":"These remarks assert impossibility results for any estimator via a 'straightforward extension' of [18, Theorem 1.6], without proof. They are peripheral and not used in any theorem, but as written they overstate the contribution. Please label them explicitly as conjectural/open problems or provide proofs.","section":"Remarks 1.14 and 2.6"},{"comment":"The phrasing 'conditional conclusions hold on an event with P^A-probability tending to one for X∼P^A_{β0,h0}' is awkward. Clarify whether the conclusions are conditional on the high-probability event for A and then almost sure/probability statements for X, and define the joint probability space.","section":"Theorem 2.5 statement"},{"comment":"In the proof of Theorem 1.3, the notation O_{β0,h0,γ,µ}(δ^{-7/2}N^{-1}) and the '≍' relations hide constants that are not specified. Defining these constants or replacing them with explicit inequalities would improve readability.","section":"Appendix A, notation"},{"comment":"The strong pseudo-regularity condition (6) and non-degeneracy condition (7) are used to reduce the mean-field condition (4) to Tr(R²)=o(N/logN) and to transfer unweighted spectral gaps. This is a real limitation for tensors with highly heterogeneous positive codegrees or many zero rows. The applications in Section 2 all verify these conditions, so the theorems are internally consistent, but the scope of the mean-field and ill-conditioning results should be stated more prominently as conditional on (6)-(7).","section":"Section 3.2, Lemmas 3.4-3.5"}],"recommendation":"major_revision","confidential_remarks":"The rate gap in the proof of Theorem 1.3 is local and fixable, but it is load-bearing because Theorem 1.4 and the √N-consistency claims in Section 2 rely on the exact O_p(N^{-1/2}T_N^{-1}) rate. I recommend major revision rather than rejection: the central ideas and the rest of the proof appear sound, and the authors should be able to correct the argument with a small change. The unproved impossibility remarks should also be softened or proved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a serious theory paper and the main load-bearing claims hold up on inspection. It extends the quadratic pseudolikelihood program (Chatterjee; Ghosal–Mukherjee; Chen–Sen–Wu) to cubic tensor interactions. The central message—joint estimability of (β,h) is controlled by the variance T_N of local fields; enough heterogeneity gives √N-consistency, near-homogeneity gives ill-conditioning—is stated cleanly and proved in detail. I checked several delicate steps (the identity 1ᵀ(ΣAᵢ²)1=Tr(R²), the coefficient 1/2 in (58), the use of strong pseudo-regularity in Lemmas 3.4–3.5) and they are correct. No circularity, no fitted constants.\n\nWhat's new: the tensor setting itself, the mean-field and Gaussian-width machinery, and the applications. The edge-three-star all-parameter ill-conditioning result and the strong-antiferromagnetic √N-consistency for edge-triangle are genuinely new, and the cyclic-versus-integer 3-AP contrast is a nice qualitative payoff. The inhomogeneous hypergraph example shows the scope.\n\nSoft spots, in proportion:\n\n- Remarks 1.14 and 2.6 assert impossibility for any estimator as 'straightforward extension' of [18, Theorem 1.6] without proof. These are not used in any theorem, so the central results don't depend on them, but an unproved impossibility claim shouldn't stand as a statement of fact. Demote to conjectures or put the extension in an appendix.\n\n- The sufficient conditions rest on strong pseudo-regularity (6) and non-degeneracy (7). That is a real boundary: if positive pair-codegrees are heterogeneous, the mean-field criterion isn't verified by the lemmas. But the applications all check (6)–(7), so this is a stated limitation, not a hidden one.\n\n- The N^{-2/7} threshold in Theorem 1.3 is unusual but derived, not fitted. Not a flaw.\n\n- The proofs are long and there is some notational roughness. An independent line-by-line check of Appendices A and B would be worthwhile before the results are treated as settled, but I did not find a gap that threatens the main theorems. The citation pattern looks fine; self-citations are legitimate technique references.\n\nWho it's for: statisticians and probabilists working on ERGMs, hypergraph models, or high-dimensional Gibbs estimation. It deserves a careful referee. I'd send it to review, and I'd want the two unproved remarks addressed in revision.","headline":"Solid first joint-pseudolikelihood analysis for cubic tensor Gibbs measures; central dichotomy is credible, but two unproved impossibility remarks should be demoted to conjectures or proved.","tokens_in":49379,"tokens_out":1866,"would_cite":true,"duration_ms":18603,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["62F12","60F10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Joint pseudolikelihood estimation in cubic-tensor Gibbs models is governed by the variability of local fields: enough inhomogeneity yields √N-consistent estimators, while homogeneity makes the Hessian degenerate.","keywords":["cubic tensor Gibbs measures","pseudolikelihood estimation","exponential random graph models","mean-field approximation","arithmetic progressions","inhomogeneous hypergraphs","joint parameter estimation","ill-conditioning"],"falsifier":"Simulate the edge–triangle ERGM with $\\mathrm{Bernoulli}(1/2)$ edges at $\\beta_0=-10$, $h_0=0$ and compute $T_N$ on each draw; the paper predicts $T_N$ is bounded below by a positive constant with high probability, so observing $T_N \\to 0$ in probability would overturn Theorem 2.1(b). Conversely, simulating the cyclic 3-AP model at $\\beta_0=1$, $h_0=1$ (where Theorem 1.13 predicts $T_N=o_p(1)$) and seeing $T_N$ bounded away from zero would falsify that half of the dichotomy.","tokens_in":48439,"feed_emoji":"📊","tokens_out":6625,"duration_ms":58090,"temperature":0.7,"texified_at":"2026-08-05T21:56:43.809972+00:00","pith_summary":"This paper studies joint estimation of the two parameters $(\\beta,h)$ of a Gibbs measure with cubic tensor interactions from a single high-dimensional observation, a setting that covers dense exponential random graph models, 3-term arithmetic-progression models, and inhomogeneous random hypergraphs. Its central claim is that joint pseudolikelihood estimation succeeds or fails according to the variability of the local fields: if the empirical variance $T_N$ of the local fields is $\\omega_p(N^{-2/7})$, the maximum pseudolikelihood estimator exists with probability tending to one and has error $O_p(N^{-1/2}T_N^{-1})$, which becomes $\\sqrt{N}$-consistency when $T_N$ stays bounded below. Conversely, under mean-field, regularity, and spectral-gap conditions, in the ferromagnetic regime with nonnegative external field the local fields homogenize and the pseudolikelihood Hessian degenerates. The paper shows that this dichotomy is realized in concrete models: edge–three-star ERGMs are ill-conditioned everywhere, edge–triangle and cyclic 3-AP models are ill-conditioned ferromagnetically but $\\sqrt{N}$-consistent under strong antiferromagnetism, while integer 3-AP and inhomogeneous hypergraphs with nonconstant profiles are always $\\sqrt{N}$-consistent.","texify_model":"deepseek-v4-flash","texify_usage":{"total_tokens":7612,"prompt_tokens":966,"completion_tokens":6646,"prompt_tokens_details":{"cached_tokens":0},"prompt_cache_hit_tokens":0,"prompt_cache_miss_tokens":966,"completion_tokens_details":{"reasoning_tokens":5683}},"feed_headline":"Local-field variance decides cubic-model learnability","feed_subtitle":"Dense ferromagnetic ERGMs stay ill-conditioned; strong antiferromagnetism or boundary effects restore √N-consistency.","key_machinery":"The central object is the empirical variance $T_N$ of the local fields $m_i(x)=\\sum_{j,k} A_{ijk}x_jx_k$; it appears in the determinant of the pseudolikelihood Hessian, $|H|=N^2 \\tilde T_N \\approx N^2 T_N$ up to factors of $\\Lambda''$, so it determines the curvature of the estimating equations. The proof machinery is a mean-field approximation for cubic-tensor Gibbs measures: a variational formula for the free energy with Gaussian-width error, a low-complexity description of the conditional mean vector $b(x)$, and a pair of checkable structural conditions — strong pseudo-regularity (all positive pair-codegrees $R_{ij}$ are comparable) and non-degeneracy (all row sums positive with average bounded away from zero) — whic","core_discovery":"On the paper's own terms, the discovery is a mechanism: the local-field variance $T_N(x)=N^{-1}\\sum_i (m_i(x)-\\bar m(x))^2$, with $m_i(x)=\\sum_{j,k} A_{ijk}x_jx_k$, controls both the existence and the rate of the maximum pseudolikelihood estimator. Theorem 1.3 proves that $T_N=\\omega_p(N^{-2/7})$ implies existence and the error bound $\\max(|\\hat\\beta-\\beta_0|,|\\hat h-h_0|)=O_p(N^{-1/2}T_N^{-1})$; Theorem 1.4 gives two simple tensor conditions — $\\operatorname{Tr}(R^2)=\\Omega(N)$ or $\\sum_i (R_i-\\bar R)^2=\\Omega(N)$ — that force $T_N=\\Omega_p(1)$ and hence $\\sqrt{N}$-consistency. Theorem 1.13 complements this: for mean-field, asymptotically regular, well-connected tensors with nonnegative entries and a stochastically nonnegative reference measure, the ferromagnetic regime $\\beta_0>0$,","pith_inferences":["The non-negativity of A is only essential for Theorem 1.13, as the paper remarks; a testable extension is that signed tensors with the same row-sum statistics can restore local-field heterogeneity and break the ferromagnetic ill-conditioning.","The paper leaves open the limiting distribution of the MPLE; if established, it would turn the consistency rates into confidence sets, a natural next step the paper explicitly flags.","The contrast between integer and cyclic 3-AP tensors suggests a general principle: any source of row-sum inhomogeneity (boundary effects, nonconstant kernels) is a robust route to estimability, while translation-invariant well-connected tensors are the hard case."],"forward_implications":["Any cubic tensor model with bounded row sums whose sample has local-field variance at least a small power of N yields a consistent MPLE with a quantitative rate; no such joint guarantee existed for order-three interactions before.","In dense ERGMs, the edge–triangle model is √N-estimable in a sufficiently strong antiferromagnetic regime, giving a concrete parameter region where both the edge and triangle parameters can be recovered from one network observation.","The edge–three-star ERGM is ill-conditioned for every inverse temperature and external field, so a single network observation cannot separate edge and three-star effects through pseudolikelihood anywhere in parameter space.","For 3-AP models, boundary effects in the integer case produce enough inhomogeneity for √N-consistency, whereas the translation-invariant cyclic case is ferromagnetically ill-conditioned but antiferromagnetically consistent.","The mean-field approximation tools (variational free energy, low-complexity conditional means) are established under explicit tensor conditions and can be applied to other dense high-order interaction models."],"fun_headline_variants":["Local-field variance controls cubic tensor pseudolikelihood rates","Ferromagnetic ERGMs ill-conditioned; antiferromagnetism restores √N-consistency","Cubic tensor estimation: local-field variance is the key","Antiferromagnetic regime enables √N-consistency in ERGMs"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the tensor's positive pair-codegrees are all comparable and every vertex has a positive row sum with average bounded away from zero; all applications verify the mean-field and spectral-gap conditions through this assumption, and without it the dichotomy between consistency and ill-conditioning is not established.","fun_headline_variants_meta":{"raw":{"variants":["Local-field variance controls cubic tensor pseudolikelihood rates","Ferromagnetic ERGMs ill-conditioned; antiferromagnetism restores √N-consistency","Cubic tensor estimation: local-field variance is the key","Antiferromagnetic regime enables √N-consistency in ERGMs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000668,"raw_usage":{"total_tokens":2875,"prompt_tokens":731,"completion_tokens":2144,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":475,"completion_tokens_details":{"reasoning_tokens":2065}},"tokens_in":475,"tokens_out":2144,"duration_ms":16172,"temperature":1.0,"reasoning_tokens":2065,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:23:45.457384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the edge–triangle ERGM with $\\mathrm{Bernoulli}(1/2)$ edges at $\\beta_0=-10$, $h_0=0$ and compute $T_N$ on each draw; the paper predicts $T_N$ is bounded below by a positive constant with high probability, so observing $T_N \\to 0$ in probability would overturn Theorem 2.1(b). Conversely, simulating the cyclic 3-AP model at $\\beta_0=1$, $h_0=1$ (where Theorem 1.13 predicts $T_N=o_p(1)$) and seeing $T_N$ bounded away from zero would falsify that half of the dichotomy.","supporting_citations":[],"review_version":1}