Pith. sign in

REVIEW 3 major objections 5 minor 29 references

Filtering and Statistical Properties of Unimodal Maps Perturbed by Heteroscedastic Noises

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read For a class of chaotic unimodal maps with two state-dependent noises, the filtering distribution converges to a unique limit independent of the initial guess, carrying CLT, large-deviation, concentration, Gumbel, and Poisson laws.

desk verdict The filter-stability result is real and well-proved under its Main Assumption, but the flagship financial application is not verified because that assumption is never checked for the calibrated parameters. read the letter →

arxiv 2411.13939 v3 pith:MPH7K2HV submitted 2024-11-21 math.ST math.DSmath.PRstat.TH

Fabrizio Lillo , Stefano Marmi , Matteo Tanzi , Sandro Vaienti This is my paper · ORCID
classification math.STmath.DSmath.PRstat.TH MSC 60G3562M2037A5060F0560F10
keywords filteringstabilityunimodalmapsheteroscedasticnoiseobservationalequivariantmeasurerandomconesHilbertprojectivemetricbankleveragesystemicrisk
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

This paper treats a hidden process that evolves by a chaotic unimodal map plus state-dependent dynamical noise, then passes through a second, independent state-dependent observational noise; both noises are heteroscedastic, meaning their size varies with the state. The main claim is filter stability: the conditional distribution of the hidden state given the whole observed history converges, as observations accumulate, to a unique absolutely continuous distribution that does not depend on the initial prior, so different agents starting from different guesses eventually make the same forecasts. The proof uses a random-cone contraction for the filtering cocycle, and the same spectral machinery yields a central limit theorem, a large-deviation principle, concentration inequalities, a Gumbel extreme-value law, and Poisson statistics for the observed process. The authors connect these results to a bank-leverage model of systemic risk, where the hidden state is the leverage of a representative bank.

What carries the argument

The central object is the filtering cocycle $P((x,z),\nu)=(\sigma(x,z),P_z\nu)$ on fibres over the stationary process $(X_n,Z_n)$, where $P_z$ is the Bayesian likelihood update followed by the Markov-kernel prediction. The paper equips each fibre with the cone $V_0(J_z)$ of integrable nonnegative densities supported on the observational-noise support $J_z$, carrying the Hilbert projective metric $\Theta_0$. Birkhoff's theorem (Theorem 3.3) turns a linear map sending $V_0(J_z)$ into a smaller cone $V_c(J_{z'})$ into a strict contraction; the Main Assumption guarantees this happens with positive frequency, so the intersection of iterated images is a unique direction, the equivariant measure. The limit theorems are obtained by perturbing the quasi-compact Perron-Frobenius operator $L$: the perturbed operators $L_z$ and the Keller-Liverani spectral theory give the leading-eigenvalue expansion that produces the CLT, large deviations, Gumbel law, and Poisson statistics.

What would settle it

Simulate the bank-leverage chain with the parameters of Figure 1 using two different initial priors for the hidden state, and compare the filtering distributions after a long common observation sequence; if they remain at positive distance on a recurrent set of observations, the Main Assumption fails for those parameters. A more direct check is to evaluate $|T(J_{z_0})| < \sigma(x)\delta/2$ for a point $z_0$ in the support of $\mu_Z$ and for all $x \in J_{z_0}$; condition A of Lemma 3.7 failing means the paper's sufficient condition is not met.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that filtering is well posed for a class of chaotic unimodal maps with two heteroscedastic noises. The filtering cocycle constructed from the Bayesian update and the Markov-chain prediction admits a unique absolutely continuous equivariant probability measure, and any initial density converges to it at an exponential rate in the Hilbert projective metric (Proposition 3.6). The key mechanism is Birkhoff contraction on random cones: when the observation $z$ lies in a set $I'$ of positive measure for which the transition density $\zeta(x,\cdot)$ is bounded below by $c>0$ on an interval $F_z$ containing the support $J_z$ of the next observational-noise likelihood, the update maps the cone of nonnegative densities on $J_z$ into a strictly smaller cone $V_{c^4}(J_{z'})$, contracting distances. Because such good observations recur with positive frequency, the composed contractions select one direction independent of the starting point. The same framework gives a central limit theorem and large-deviation principle (Prop 4.1), a concentration inequality for the empirical measure (Prop 5.1), a Gumbel law with extremal index 1 (Prop 7.3), and Poisson distribution and point-process statistics (Props 8.1 and 8.2) under the stated assumptions.

Load-bearing premise

The filtering-stability conclusion rests on the assumption that, for a positive fraction of observations, the dynamical noise is strong enough relative to the observational noise that the hidden state's transition density is bounded below by a positive constant on an interval containing the support of the next observation's noise; if that ratio is reversed, the contraction argument gives no uniqueness.

Editorial extensions

If this is right

  • For any two initial priors on the hidden state, the filtering distributions stay close and converge to the same equivariant filter as the observation window grows, so the filtering problem is well posed for these systems.
  • A central limit theorem holds for additive functionals of the observed process, and large-deviation bounds with a convex rate function follow (Prop 4.1).
  • If the modulation term is constant, the empirical measure of the observed process converges to the push-forward measure $\mu'$ with exponential concentration, and its deviations from the stationary measure $\mu$ are exponentially small (Prop 5.1 and its application).
  • Under the threshold and ball-shrinking assumptions, rare events of the observed process obey a Gumbel law with extremal index 1, and visit counts converge to a standard Poisson distribution and a Poisson point process (Props 7.3, 8.1, and 8.2).
  • In the bank-leverage model of systemic risk, when the observational noise is sufficiently small relative to the dynamical noise, these statistical laws apply to the leverage process (Lemma 3.7).

Reading between the lines

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

  • Inference: If condition A of Lemma 3.7 fails for the calibrated bank-leverage parameters, the paper's machinery gives no uniqueness; a numerical comparison of $|T(J_{z_0})|$ with $\min_x \sigma(x)\delta/2$ would delimit the regime where filter stability can be expected and might reveal a transition in the noise ratio.
  • Inference: The GEV detector sketched in Appendix A, where the location parameter should equal $\log(t \int s(x)^{-1}\,d\mu(x))$, could be tested on simulated time series to estimate the modulation average from data; the paper does not carry out that numerical test.
  • Inference: The paper leaves the non-constant-modulation case open for the extremal index and Poisson parameter; a natural test is to compute the quantities $q_k$ numerically for a slowly varying modulation and see whether the compound Poisson law reduces to the standard Poisson law.
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

3 major / 5 minor

Summary. The paper studies a hidden Markov model in which the latent process is a unimodal map perturbed by heteroscedastic Markov-chain noise and the observation is further corrupted by heteroscedastic observational noise. The central result is a filtering theorem: under the Main Assumption of Section 3.2.3, the filtering cocycle admits a unique absolutely continuous equivariant probability measure, so the optimal filter forgets its initial prior (Proposition 3.6). The proof uses Birkhoff contraction in Hilbert metrics on random cones. The paper also claims a CLT and large deviation principle for the observed process (Proposition 4.1), concentration inequalities under a constant modulation assumption (Proposition 5.1), Gumbel extreme-value laws (Propositions 7.1 and 7.3), and Poisson statistics for rare events (Propositions 8.1 and 8.2). The motivating application is a financial leverage model from the authors' previous work.

Significance. If the main theorem is correct, it is a valuable contribution: it gives filter stability for a genuinely nonlinear, non-uniformly hyperbolic hidden Markov model with heteroscedastic noises, and the random-cone method is an elegant way to obtain an equivariant measure. The paper is honest in stating the Main Assumption as a sufficient condition, and the proof of Proposition 3.6 is structurally sound. The advertised financial application is a potentially important selling point, but it is not verified for the calibrated parameters. The EVT and Poisson sections are more conditional and contain admitted gaps, especially when the modulation term is not constant. Overall the paper contains a substantial core result together with several less-developed peripheral claims that need to be either proved or explicitly delimited.

major comments (3)
  1. [§3.2.3, Lemma 3.7, Remark 3.8, Fig. 1] The advertised application to the financial model is not supported by the verified hypotheses. Lemma 3.7 establishes the Main Assumption only when conditions A and B hold, and Remark 3.8 gives the sufficient global bound ε max_x s(x) ≤ min{|I|/10, min_x σ(x)δ/10}. The paper never checks this bound for the parameters used in Fig. 1 (γ0=15.969, α=1.64, Σε=2.7×10^-5), and it does not state the corresponding σ(x), δ, or ε. It is therefore unknown whether the observational noise is small enough relative to the dynamical noise for filter stability to hold in the financial model. The abstract and introduction should either state this small-noise restriction explicitly or provide the missing numerical verification.
  2. [§7.1 and §7.2] Proposition 7.1 is conditional on β>0, but the positivity of the extremal index β is not established in Section 7.1. It is proved only later, in Section 7.2, under the additional Assumption EI that s is constant; the text explicitly says 'When s is not a constant, we do not have a proof.' Thus the Gumbel law for the general heteroscedastic observational-noise model is not demonstrated. In addition, the proof relies on inequality (33), which requires the stationary density h to be strictly positive on the set {x : x+s(x)ε ∈ B_t}; this strict positivity is asserted without proof in Section 7.1.
  3. [§8.1 and §8.2] Propositions 8.1 and 8.2 are stated under Assumption EI, i.e., s is constant, so they do not cover the general heteroscedastic case advertised in the title. Moreover, the proof that l(λ)=1 is only sketched with the phrase 'By the same argument used in the preceding section,' and Proposition 8.2 relies on the unpublished reference [29]. The compound-Poisson intermediate statement does not by itself prove standard Poisson statistics. A complete proof or a published reference is needed before these results are claimed.
minor comments (5)
  1. [§5.1, Eq. (25)] The displayed inequality in Eq. (25) and the subsequent display contain mismatched parentheses and a duplicated event; the notation P×Pε should be used consistently and the parentheses should be balanced.
  2. [§3.1, Eq. (18)] The convention P(Z0=0)=1 appears to conflict with the later use of a stationary process (Xn,Zn) with a nontrivial Z0 distribution; the authors should clarify how these two viewpoints are reconciled.
  3. [§7.1 and §8.1] The symbol W_t is used both for the extreme-value probability in (29)-(30) and for the limiting random variable in Section 8.1; please use distinct notation.
  4. [§8.1] The statement 'Supposing that l(λ) is continuous at 0' should be justified, since identifying the characteristic-function limit requires control of l(λ) on a neighborhood of 0, not merely continuity at a single point.
  5. [References] The reference [18] is likely misspelled as 'Brockner' rather than 'Brocker'; please check.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the filter-stability theorem is proved from an explicit Main Assumption, and citations to [1,2] provide independent spectral groundwork rather than assuming the conclusion.

full rationale

The central claim, Proposition 3.6, is derived rather than assumed: under the Main Assumption in Section 3.2.3, the normalized likelihood operator maps V0(J_{z0}) into V_{c^4}(J_{z1}) and contracts the Hilbert metric on the positive-frequency event E, yielding the unique equivariant measure. The Main Assumption is an explicit hypothesis (lower bound c>0, support inclusions J_{z'} subset F_z) whose verification in Lemma 3.7 is separate from the theorem; none of its conditions is equivalent to existence of the equivariant measure. The paper imports from the authors' prior work [1,2] the compactness/quasi-compactness of the transfer operator L, the uniqueness of the stationary measure of the unobserved Markov chain, and exponential decay of correlations (equations (5) and (9), Remark 1.1). These are parameter-free facts about the chain under Assumption TM and do not include the filtering conclusion or the observational noise; they are therefore independent support, not a circular premise. The limit theorems (Propositions 4.1, 5.1, 7.1, 7.3, 8.1, 8.2) are obtained by standard Keller-Liverani/Kato spectral perturbation around L and do not rename fitted quantities as predictions. The only serious caveat is a scope gap, not a circularity: the sufficient global bound in Remark 3.8 ensuring conditions A and B of Lemma 3.7 is not checked for the parameters used in Figure 1 (gamma0=15.969, alpha=1.64, Sigma_epsilon=2.7e-5), so the abstract's unqualified claim that the results apply to the financial model exceeds the verified hypotheses. No step satisfying the quote-and-reduction standard for circularity is exhibited.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The paper introduces no fitted numerical parameters and no new postulated entities. Theorems depend on the prior stationary-measure machinery from [1,2] (Assumption TM) plus several extra assumptions specific to this paper (Main Assumption, conditions A/B, CI, EI, S/S'). All numerical parameters in the financial example are inherited from earlier calibrations and do not enter the proofs.

assumptions (6)
  • domain assumption Assumption TM: T is a continuous unimodal map of I with unique maximum, transitive, with acip, extended to I~; noise σ differentiable, density g as in (7) with a satisfying (6).
    This class of systems is taken from the authors' prior papers [1,2] and underlies every theorem; the filtering and limit results are proven only for maps and noises satisfying TM.
  • ad hoc to paper Main Assumption (Section 3.2.3): there exists c>0 and a set I' of positive µZ-measure such that for a.e. z∈I' there is an interval Fz with ζ(x,y)>c for x∈Jz, y∈Fz, Jz'⊂Fz and g(z',y)>c for z'∈F'_z, and ζ,g<c^{-1}.
    This abstract condition is introduced solely to make the Birkhoff-cone contraction argument work; it is not implied by Assumption TM.
  • ad hoc to paper Lemma 3.7 conditions A and B: |T(Jz0)|<σ(x)δ/2 for all x∈Jz0, and there exists an interval F' around T(z0) with Jz'⊂O_p(T(Jz0)) for z'∈F'.
    These are the concrete small-observation-noise hypotheses under which the financial model satisfies the Main Assumption; they are not verified for the empirical parameter values cited in the paper.
  • ad hoc to paper Assumption CI (Section 5): the modulation function s is constant, s≡sM.
    The concentration inequality Prop 5.1 holds only after imposing this restriction, which removes the heteroscedasticity of the observational noise.
  • domain assumption Assumptions S and S' (Section 7.1): ∫ψ((Bt-x)/s(x))dx→0 and t∫ψ((Bt-x)/s(x))dµ(x)→τ as t→∞.
    These are standard shrinking-target and scaling conditions needed to apply Keller-Liverani perturbation theory for extremes.
  • ad hoc to paper Assumption EI (Section 7.2): s is constant (taken to be 1) and ψ is non-atomic.
    The proof that the extremal index equals 1 and that the Poisson law holds uses this restriction; the paper's title advertises heteroscedastic noise, but these theorems do not cover it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Filtering and Statistical Properties of Unimodal Maps Perturbed by Heteroscedastic Noises." pith.science (2026). https://pith.science/paper/MPH7K2HV

@misc{pith2026241113939,
  author       = {Pith},
  title        = {Pith review of: Filtering and Statistical Properties of Unimodal Maps Perturbed by Heteroscedastic Noises},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MPH7K2HV}},
  note         = {Machine review of arXiv:2411.13939}
}
read the original abstract

We propose a theory of unimodal maps perturbed by an heteroscedastic Markov chain noise and experiencing another heteroscedastic noise due to uncertain observation. We address and treat the filtering problem showing that by collecting more and more observations, one would predict the same distribution for the state of the underlying Markov chain no matter one's initial guess. Moreover we give other limit theorems, emphasizing in particular concentration inequalities and extreme value and Poisson distributions. Our results apply to a family of maps arising from a model of systemic risk in finance.

Figures

Figures reproduced from arXiv: 2411.13939 by the authors.

Figure 1
Figure 1. Taken from [2]: plot of the deterministic component T(ϕ), γ0 = 15.969, α = 1.64, Σϵ = 2.7 × 10−5 . The value for γ0 is taken from the empirical analysis in [1], Section 7.2, (where it is denoted simply by γ) . The value α = 1.64 corresponds to a VaR constraint of 5% in case of a Gaussian distribution for the returns. The values Σϵ = 2.7 × 10−5 is taken from [25], [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. The image above depicts the typical situation one has under the Main Assumption. The horizontal axis represents a subset of I, while the vertical axis gives the values taken by the transition densities ζ(x, ·). Almost every point x ∈ Jz, where z ∈ I ′ , has transition density larger than c > 0 on Fz, and every point of F ′ z ⊂ Fz has the corresponding set Jz′ contained in Fz. The above entails that when an observati… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references · 28 canonical work pages

  1. [1]

    Lillo, G

    F. Lillo, G. Livieri, S. Marmi, A. Solomko, S. Vaienti, Analysis of bank leverage via dynamical systems and deep neural networks , SIAM Journal on Financial Mathematics, 14, (2), pp. 598-643, (2023)

    work page 2023

  2. [2]

    Lillo, G

    F. Lillo, G. Livieri, S. Marmi, A. Solomko, S. Vaienti, Unimodal maps perturbed by heteroscedastic noise: an application to a financial systems , Journal of Statistical Physics, 190, 156 (2023)

    work page 2023

  3. [3]

    Thunberg, Periodicity versus Chaos in One-Dimensional Dynamics , SIAM Review 43.1 (2001), pp

    H. Thunberg, Periodicity versus Chaos in One-Dimensional Dynamics , SIAM Review 43.1 (2001), pp. 3–30

    work page 2001

  4. [4]

    Kantz, T

    H. Kantz, T. Schreiber, Nonlinear Time Series Analysis, second ed., Cambridge University Press, 2004

    work page 2004

  5. [5]

    de Melo, S

    W. de Melo, S. van Strien, One-Dimensional Dynamics, Vol. 25. Springer Science and Business Media, 2012

    work page 2012

  6. [6]

    Keller, Rare events, exponential hitting times and extremal indices via spectral perturbation , Dynamical Systems 27.1 (2012), pp

    G. Keller, Rare events, exponential hitting times and extremal indices via spectral perturbation , Dynamical Systems 27.1 (2012), pp. 11–27. 25

    work page 2012

  7. [7]

    Keller, C

    G. Keller, C. Liverani, Rare events, escape rates and quasistationarity: some exact formulae , Journal of Statistical Physics 135 (2009), pp. 519–534

    work page 2009

  8. [8]

    Lalley, Removing the noise from chaos plus noise , in: Nonlinear Dynamics and Statistics (Cambridge, 1998), Birkh¨ auser, Boston, MA, 2001, pp

    S.P. Lalley, Removing the noise from chaos plus noise , in: Nonlinear Dynamics and Statistics (Cambridge, 1998), Birkh¨ auser, Boston, MA, 2001, pp. 233–244

    work page 1998

Show all 29 references
  1. [9]

    Lalley, Beneath the noise, chaos , Ann

    S.P. Lalley, Beneath the noise, chaos , Ann. Statist. 27 (2) (1999) 461–479

    1999

  2. [10]

    Schreiber, Extremely simple nonlinear noise-reduction method , Phys

    T. Schreiber, Extremely simple nonlinear noise-reduction method , Phys. Rev. E 47 (1993) 2401–2404

    1993

  3. [11]

    Schreiber, P

    T. Schreiber, P. Grassberger, A simple noise-reduction method for real data, Phys. Lett. A (1991) 411–418

    1991

  4. [12]

    Lalley, A.B

    S.P. Lalley, A.B. Nobel, Denoising deterministic time series , Dyn. Partial Differ. Equ. 3 (4) (2006) 259–279

    2006

  5. [13]

    Judd, Chaotic-time-series reconstruction by the Bayesian paradigm: right resultsby wrong meth- ods, Phys

    K. Judd, Chaotic-time-series reconstruction by the Bayesian paradigm: right resultsby wrong meth- ods, Phys. Rev. E 67 (2003) 026212

    2003

  6. [14]

    McGoff, S

    K. McGoff, S. Mukherjee, N.S. Pillai, Statistical inference for dynamical systems: a review , 2012. arXiv:1204.6265v1 [math.ST

    2012 arXiv

  7. [15]

    Atnip, G

    J. Atnip, G. Froyland, C. Gonzalez-Tokman, S. Vaienti, Compound Poisson Statistics for dynam- ical systems via spectral perturbation , https://arxiv.org/pdf/2308.10798.pdf

    arXiv

  8. [16]

    Chigansky, R

    P. Chigansky, R. Liptser and R. Van Handel, Intrinsic methods in filter stability, in The Oxford handbook of nonlinear filtering, 319–351, Oxford Univ. Press, Oxford

  9. [17]

    Capp´ e, E

    O. Capp´ e, E. Moulines, T. Ryd´ en, Inference in hidden Markov models, Springer Series in Statistics. Springer, New York, (2005)

    2005

  10. [18]

    Brockner, G

    J. Brockner, G. Del Magno, Asymptotic stability of the optimal filter for random chaotic maps , Nonlinearity, (2017) 30 1809

    2017

  11. [19]

    Chigansky, R

    P. Chigansky, R. Liptser and R. Van Handel, Intrinsic methods in filter stability , in The Oxford handbook of nonlinear filtering, 319–351, Oxford Univ. Press, Oxford

  12. [20]

    Th. Caby, D. Faranda, S. Vaienti, P. Yiou, On the computation of the extremal index for time se- ries, Journal of Statistical Physics, 179, 1666-1697, (2020), https://arxiv.org/pdf/1904.04936.pdf

    2020 arXiv

  13. [21]

    Maldonado, Fluctuation Bounds for Chaos Plus Noise in Dynamical Systems, JSP, (2012)

    C. Maldonado, Fluctuation Bounds for Chaos Plus Noise in Dynamical Systems, JSP, (2012)

    2012

  14. [22]

    Paulin, Concentration inequalities for Markov chains by Marton couplings and spectral methods, Electronic Journal of Probability, (2015)

    D. Paulin, Concentration inequalities for Markov chains by Marton couplings and spectral methods, Electronic Journal of Probability, (2015)

    2015

  15. [23]

    Aimino, M

    R. Aimino, M. Nicol, S. Vaienti, Annealed and quenched limit theorems for random expanding dynamical systems, Probability Theory and Related Fields, June 2015, 162, Issue 1, pp 233-274

    2015

  16. [24]

    Lucarini V, Faranda D, Freitas A M, Freitas J M, Holland M, Kuna T, Nicol M, Todd M and Vaienti S, Extremes and Recurrence in Dynamical Systems (New York: Wiley), 2016

    2016

  17. [25]

    and Marmi, S., (2019), When panic makes you blind: A chaotic route to systemic risk , Journal of Economic Dynamics and Control, 100:176–199

    Mazzarisi, P., Lillo, F. and Marmi, S., (2019), When panic makes you blind: A chaotic route to systemic risk , Journal of Economic Dynamics and Control, 100:176–199

    2019

  18. [26]

    Faranda, S

    D. Faranda, S. Vaienti, Extreme Value laws for dynamical systems under observational noise , Physica D, 280-281, (2014), 86-94

    2014

  19. [27]

    Liverani, Decay of correlations, Ann

    C. Liverani, Decay of correlations, Ann. Math., 142 (2), 239-301, (1995)

    1995

  20. [28]

    Kallenberg Random measures, Akademia-Verlag Berlin, (1986) 26

    O. Kallenberg Random measures, Akademia-Verlag Berlin, (1986) 26

    1986

  21. [30]

    Atnip, G

    J. Atnip, G. Froyland, C. Gonzalez-Tokman, S. Vaienti, Thermodynamic Formalism and Per- turbation Formulae for Quenched Random Open Dynamical Systems , to appear in Dissertationes Mathematicae, https://arxiv.org/pdf/2307.00774.pdf. 27

    arXiv

Pith tools

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