{"id":"f47c8480-1fc8-4c3e-a335-42cc7d0d0851","arxiv_id":"2511.11054","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Develops a joint Catoni-type framework that simultaneously estimates parameters and variance in heavy-tailed models via coupled equations, achieving oracle-matching non-asymptotic bounds under 2β-moment assumptions using a Poincaré-Miranda proof technique.","lead":"This paper creates a tuning-free way to estimate both a model parameter and the unknown noise variance at the same time using two linked Catoni-type equations for data with heavy tails. A smart generalist might read it because heavy-tailed noise is common in real measurements and this approach avoids manual tuning while matching the performance of methods that already know the variance.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Poincaré-Miranda establishes existence of a solution inside the target deviation region but does not rule out additional solutions outside it, leaving the estimator potentially ill-defined.","rationale":"The reader correctly flags the Poincaré-Miranda step as the point where classical convex arguments are replaced. The load-bearing issue is narrower: existence inside the ball is shown, but uniqueness or exclusion of exterior solutions is not addressed, which directly affects whether the stated deviation bounds apply to the estimator as defined. This qualification moves the verdict from UNVERDICTED to CONDITIONAL pending clarification of solution selection.","tokens_in":1720,"tokens_out":424,"duration_ms":76987,"concrete_test":"Locate the application of Poincaré-Miranda (presumably in the proof of the main non-asymptotic theorem). Check whether the argument also contains a step proving that the estimating map cannot be zero outside the rectangle of radius equal to the claimed rate; if that step is missing, the claim holds only for existence of a good solution rather than for the estimator defined by the system.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim defines the estimators as solutions to the coupled non-convex Catoni-type equations and asserts that these solutions jointly satisfy sub-Gaussian deviation bounds under only 2β moments. The proof strategy invokes Poincaré-Miranda to guarantee existence inside a rectangle whose side lengths are set to the desired rates. This theorem supplies a zero inside the rectangle once opposing sign conditions hold on each face, but supplies no information about zeros outside the rectangle. Without a separate argument showing that the estimating functions cannot vanish when the joint deviation exceeds the target radius (e.g., by showing that the expectation term dominates the fluctuation term uniformly outside the ball), nothing precludes spurious distant solutions. Under the stated moment assumptions the fluctuation terms remain only polynomially integrable, so controlling the far-field behavior may require extra truncation or growth conditions not implied by the boundary analysis alone. Consequently the deviation bound is guaranteed only for some solution, not necessarily for every solution of the system.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.3","summary":"The paper develops a tuning-free joint robust estimation framework for parametric models with heavy-tailed noise, simultaneously estimating the target parameter and unknown noise variance via a system of two coupled Catoni-type estimating equations. It instantiates the approach for mean estimation, linear regression, and ℓ₂-penalized regression. The central theoretical claim is the derivation of non-asymptotic sub-Gaussian-type joint deviation bounds under only a finite 2β-th moment assumption (β ∈ (1,2]), with rates matching those of oracle procedures that know the variance in advance. The proofs rely on the Poincaré-Miranda theorem to establish existence of solutions to the non-convex system, bypassing classical convex M-estimation arguments.","tokens_in":1942,"tokens_out":610,"duration_ms":30829,"significance":"If the central claims hold, the work would provide a valuable contribution to robust statistics by delivering optimal, variance-adaptive estimators under weak moment conditions without tuning parameters. The methodological innovation of adapting Poincaré-Miranda for joint non-convex estimation could extend to other problems involving heterogeneous parameters, and the oracle-matching rates under 2β moments represent a strong theoretical achievement.","major_comments":[{"comment":"The proof strategy invokes Poincaré-Miranda to guarantee a zero of the coupled estimating functions inside a rectangle whose dimensions are set to the target deviation rates. However, the theorem only ensures existence within the rectangle once opposing sign conditions hold on the faces; it provides no control over possible additional zeros outside the rectangle. Under the stated polynomial integrability of the fluctuation terms, far-field behavior is not automatically dominated, so the non-asymptotic bounds may apply only to some solutions rather than to every solution of the system. This affects the well-definedness of the estimator and the validity of the joint deviation claim.","section":"Proofs of the main deviation bounds (theorems establishing joint sub-Gaussian rates)"},{"comment":"The estimators are defined as solutions to the coupled non-convex, non-linear equations. Without an additional argument showing that no solutions exist outside the target deviation ball (e.g., via uniform domination of the expectation term by the fluctuation term for large deviations), or a constructive selection rule for the solution inside the rectangle, it remains unclear which root is being bounded and how the procedure is implemented in practice.","section":"Section 2 (definition of the joint estimators) and the subsequent theoretical analysis"}],"minor_comments":[{"comment":"Notation for the Catoni function and the coupled equations could be made more explicit when first introduced to aid readability for readers unfamiliar with the original Catoni estimator.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to fit the scope of a theoretical statistics journal. The citation pattern is appropriate, but the authors should explicitly contrast their Poincaré-Miranda approach with any prior uses of topological fixed-point theorems in robust estimation literature."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the thorough review and valuable comments on our manuscript. We appreciate the recognition of the potential contributions and address the major comments point by point below. We plan to make revisions to clarify the well-definedness of the estimators.","responses":[{"response":"We thank the referee for highlighting this important subtlety. The Poincaré-Miranda theorem is invoked solely to guarantee existence of at least one solution inside the target rectangle. We agree that this does not automatically rule out other zeros outside the rectangle. In the revision we will explicitly define the joint estimator as any solution lying inside the rectangle whose existence is assured by the theorem. The deviation bounds are then stated for this defined estimator. A short remark will be added noting that the selection is by construction inside the region of interest. This resolves the well-definedness issue while remaining faithful to the minimal moment assumptions.","revision_made":"yes","referee_comment":"The proof strategy invokes Poincaré-Miranda to guarantee a zero of the coupled estimating functions inside a rectangle whose dimensions are set to the target deviation rates. However, the theorem only ensures existence within the rectangle once opposing sign conditions hold on the faces; it provides no control over possible additional zeros outside the rectangle. Under the stated polynomial integrability of the fluctuation terms, far-field behavior is not automatically dominated, so the non-asymptotic bounds may apply only to some solutions rather than to every solution of the system. This affects the well-definedness of the estimator and the validity of the joint deviation claim."},{"response":"We agree that the current wording leaves ambiguity about which root is intended. We will revise Section 2 to state that the estimator is defined to be a solution of the coupled system that lies inside the rectangle for which existence is guaranteed by Poincaré-Miranda. For implementation we will add a brief discussion indicating that the low-dimensional (two-equation) system can be solved numerically by standard methods such as Newton iteration or a merit-function minimization, initialized at a point scaled to the target deviation rates. This makes both the theoretical object and the practical procedure unambiguous.","revision_made":"yes","referee_comment":"The estimators are defined as solutions to the coupled non-convex, non-linear equations. Without an additional argument showing that no solutions exist outside the target deviation ball (e.g., via uniform domination of the expectation term by the fluctuation term for large deviations), or a constructive selection rule for the solution inside the rectangle, it remains unclear which root is being bounded and how the procedure is implemented in practice."}],"tokens_in":1498,"tokens_out":524,"duration_ms":59514,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper's main contribution is a joint estimation procedure for a target parameter and the noise variance using two coupled Catoni-type equations that require no tuning. They show these achieve oracle rates under only 2β moments for β in (1,2], and they use the Poincaré-Miranda theorem to establish existence of solutions in the right region for the non-convex system in mean estimation, linear regression, and penalized regression cases.","headline":"The paper gives a tuning-free joint Catoni estimator for parameter and variance in heavy-tailed models via coupled equations and uses Poincaré-Miranda to get oracle-matching rates, but the analysis only guarantees existence inside the target region.","tokens_in":2414,"tokens_out":176,"would_cite":false,"duration_ms":25639,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"echoes","rs_module":"Foundation/AlexanderDuality.lean","rs_theorem":"alexander_duality_circle_linking","paper_passage":"We employ the Poincaré–Miranda Theorem to show that the solutions lie within certain geometric regions, such as cylinders or cones, centered around the true parameter values."},{"relation":"unclear","rs_module":"Cost/FunctionalEquation.lean","rs_theorem":"washburn_uniqueness_aczel","paper_passage":"ψ1 satisfying -log(1-x+|x|²/2) ≤ ψ1(x) ≤ log(1+x+|x|²/2)"}],"headline":"Statistical robust M-estimation via Poincaré-Miranda on coupled Catoni equations; no overlap with RS cost or forcing machinery","alignment":"orthogonal","rationale":"The paper's core is non-asymptotic deviation bounds for joint mean/variance (and regression) estimators defined by two coupled non-convex Catoni-type equations, with existence proved by verifying sign conditions on the boundary of a rectangle (or cylinder) and invoking Poincaré-Miranda. This is classical heavy-tailed robust statistics under 2β-moment assumptions. RS derives J-cost, φ-ladder, 8-tick periodicity and dimension-3 from a single distinction via functional equations and topological forcing (AbsoluteFloorClosure, AlexanderDuality, Cost.FunctionalEquation, DimensionForcing). No shared primitives, no cosh/J-cost identities, no parameter-free constant derivations, and no recognition-ladder structure appear. The topological tool is used only for existence of statistical solutions, not for any RS-style structural forcing.","tokens_in":68891,"confidence":"high","tokens_out":378,"duration_ms":14686,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A system of two coupled Catoni-type equations estimates both a parameter and its unknown variance at sub-Gaussian rates under heavy tails, without tuning.","keywords":["heavy-tailed estimation","Catoni estimator","joint estimation","robust statistics","mean estimation","linear regression","non-asymptotic bounds","Poincaré-Miranda theorem"],"falsifier":"Generate data with exactly 2.1 moments and check whether the observed joint deviation of the estimator from the true parameter and variance exceeds the claimed sub-Gaussian bound by more than a small constant factor.","tokens_in":2632,"feed_emoji":"","tokens_out":651,"duration_ms":28689,"temperature":0.7,"pith_summary":"The paper constructs a joint estimation procedure that solves two linked Catoni-type equations to recover the target parameter and the noise variance at the same time. This is done for mean estimation, linear regression, and penalized regression under the sole assumption that the noise has a finite moment of order 2β for β between 1 and 2. The resulting non-asymptotic deviation bounds for both quantities match the rates that would be available if the variance were known in advance. Because the equations are non-convex, the analysis replaces standard convex M-estimation tools with an application of the Poincaré-Miranda theorem to guarantee the existence of suitable solutions and control their joint error.","feed_headline":"Joint Catoni equations match oracle rates under heavy tails","feed_subtitle":"Two coupled equations estimate parameter and variance together and deliver sub-Gaussian deviation bounds with only a 2β-moment assumption.","key_machinery":"The pair of coupled, non-convex Catoni-type estimating equations for the parameter and the variance, whose joint solutions are controlled via the Poincaré-Miranda theorem.","core_discovery":"The central claim is that the coupled system of two Catoni-type estimating equations admits solutions whose joint deviation from the true parameter and true variance satisfies sub-Gaussian-type bounds under a finite 2β-moment condition with β∈(1,2], with rates that match those of oracle procedures knowing the variance in advance.","pith_inferences":["The same topological control might simplify proofs for joint robust estimation in generalized linear models.","Practitioners facing data with unknown scale could replace separate variance estimation and cross-validation steps with this single procedure.","The moment condition 2β with β close to 1 suggests the method remains useful even when tails are only slightly heavier than Gaussian."],"forward_implications":["The same joint rates hold in mean estimation, linear regression, and ℓ2-penalized regression.","The bounds remain valid without knowledge of the variance or any tuning parameters.","The rates are optimal up to absolute constants in the heavy-tailed regime.","The proof strategy applies to other problems that require simultaneous estimation of parameters of different types."],"fun_headline_variants":["Coupled Catoni estimating equations match oracle rates","Tuning-free joint Catoni robust estimation for heavy tails","Sub-Gaussian joint bounds via tuning-free Catoni equations","Catoni coupled system achieves oracle rates under heavy tails"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The non-convex coupled equations must possess solutions whose joint deviations can be bounded using a topological theorem instead of convexity arguments.","fun_headline_variants_meta":{"raw":{"variants":["Coupled Catoni estimating equations match oracle rates","Tuning-free joint Catoni robust estimation for heavy tails","Sub-Gaussian joint bounds via tuning-free Catoni equations","Catoni coupled system achieves oracle rates under heavy tails"]},"model":"grok-4.3","cost_usd":0.010482,"raw_usage":{"total_tokens":4539,"prompt_tokens":639,"num_sources_used":0,"completion_tokens":62,"cost_in_usd_ticks":104815500,"prompt_tokens_details":{"text_tokens":639,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3838,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":639,"tokens_out":62,"duration_ms":39848,"temperature":1.0,"reasoning_tokens":3838,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-21T20:07:36.945345+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Generate data with exactly 2.1 moments and check whether the observed joint deviation of the estimator from the true parameter and variance exceeds the claimed sub-Gaussian bound by more than a small constant factor.","supporting_citations":[],"review_version":1}