arXiv:2408.05963v4

Non-asymptotic Estimates for Markov Transition Matrices via Spectral Gap Methods

De Huang, Xiangyuan Li

math.STmath.PR60J1037A2562M05

Abstract

We establish non-asymptotic error bounds for the classical Maximal Likelihood Estimation of the transition matrix of a given Markov chain. Meanwhile, in the reversible case, we propose a new reversibility-preserving online Symmetric Counting Estimation of the transition matrix with non-asymptotic deviation bounds. Our analysis is based on a convergence study of certain Markov chains on the length-2 path spaces induced by the original Markov chain.

AI-generated audit

Audit summary

Audited against arXiv v4

Not a correctness certificate. A “Correct” result may include yellow typos or minor formal corrections that do not affect substantive soundness. It means this audit found no unresolved substantive error under the stated criteria; it does not replace expert scrutiny or formal verification.

Current report

Detailed mathematical audit

Generated August 18, 2026
01Statements4 reported findingsCorrect

The almost-sure convergence, entrywise tail bound, and dimension-free mean-square bound for the joint-count MLE matrix are correct. The reversible symmetric-counting construction and its matrix tail and mean-square bounds are also correct. The spectral statements for the two path-space chains support exactly the gap substitutions used in those estimates. The standing state-space hypothesis should exclude the trivial one-state case so that the displayed spectral-gap definitions are nonempty.

Theorem 2.1Correct

The MLE joint-count estimates hold

Pages 4–5 and 17–21 · Theorem 2.1 and its proof · arXiv:2408.05963v4

Writing RnR_n as the empirical mean of edge indicators on the length-two path chain gives stationary mean DμPD_\mu P. The ergodic theorem yields almost-sure convergence. The cited scalar Bernstein inequality applies with variance at most μ(u)p(u,v)\mu(u)p(u,v), and Theorem 3.3 supplies ηp(P2)ηp(P)/(1+ηp(P))\eta_p(P_2)\geq\eta_p(P)/(1+\eta_p(P)). For the mean-square estimate, stationary covariances reduce to a polynomial in P2Π2P_2-\Pi_2; the factorization P2=STP_2=ST, P=TSP=TS reduces its norm to a polynomial in PΠP-\Pi, and telescoping gives the claimed O((nηp(P))1)O((n\eta_p(P))^{-1}) bound. The alternative absolute-gap and reversible bounds follow from the same polynomial estimate.

Theorem 2.2Correct

The reversible symmetric-counting estimates hold

Page 5 and pages 18 and 22–23 · Theorem 2.2 and its proof · arXiv:2408.05963v4

The orientation coin and one oracle call generate the stated Markov chain on unordered length-two paths. Detailed balance for PP makes its invariant edge law reversible, and the symmetrized matrix observable has expectation DμPD_\mu P. The path-chain absolute gap is η(P)/2\eta(P)/2, so the cited Markov-dependent matrix Bernstein bound gives the displayed operator-norm tail estimate after the constants are rescaled. Nonnegativity of the path-chain spectrum and the stationary covariance calculation give the dimension-free bound EHnDμPF2(4η(P))/(nη(P))\mathbb E\lVert H_n-D_\mu P\rVert_F^2\leq(4-\eta(P))/(n\eta(P)), with the stated initial-density factor in the nonstationary case.

Theorems 3.3 and 3.4Correct

The path-space spectral identities are correct

Pages 11–16 · Theorems 3.3–3.4 · arXiv:2408.05963v4

The rectangular factorizations give P2=STP_2=ST and P=TSP=TS, so the nonzero eigenvalues agree. After weighting by the invariant measures, SS has orthonormal columns and TT has orthonormal rows; consequently P2kP_2^k and Pk1P^{k-1} have the same nonzero singular values. The positive and negative parts of the additive symmetrization identify the reversible spectrum with those of (P+I)/2(P+I)/2 and (PI)/2(P-I)/2, while the conditional-expectation argument yields the claimed IP-gap lower bound. For the unordered path chain, reversing before extension gives T~S~=(P+I)/2\widetilde T\widetilde S=(P+I)/2, proving both spectral-gap identities and uniqueness of the invariant law.

Standing state-space hypothesisMinor formal correction

The spectral-gap definitions require at least two states

Pages 1 and 4 · standing hypotheses and Subsection 2.1.2 · arXiv:2408.05963v4

The paper permits an arbitrary finite state space but defines λ(P)\lambda(P) and ηp(P)\eta_p(P) by a supremum or infimum over nonzero vectors in L2,μ0L^0_{2,\mu}. If Ω=1|\Omega|=1, that set is empty and the displayed gaps and subsequent fractions are not finite quantities under the stated definitions. Add the standing assumption Ω2|\Omega|\geq2. This is a local range correction: the excluded one-state estimation problem is identically zero-error, and every proof and application in the paper otherwise remains unchanged.

02Proofs4 reported findingsCorrect

The central proofs are correct and complete, apart from the harmless one-state scope correction. The weighted rectangular-factorization arguments, stationary covariance sums, polynomial norm bounds, and applications of the two external concentration theorems all satisfy their required hypotheses and yield the displayed constants.

Theorem 3.3Correct and complete

The singular-value and gap comparison proof is complete

Pages 12–16 · proof of Theorem 3.3 · arXiv:2408.05963v4

The identities S1TS1=IS_1^\mathsf TS_1=I and T1T1T=IT_1T_1^\mathsf T=I follow directly from stationarity and the edge measure. They preserve the nonzero singular values in S1(PΠ)k1T1S_1(P-\Pi)^{k-1}T_1, proving the absolute- and pseudo-gap claims. The additive symmetrization is squeezed between its positive and negative Gram components; in the reversible case its square has the positive component as its unique positive square root, giving the asserted spectral decomposition. Finally, conditioning (IP2)h(I-P_2)h on the first coordinate and applying the IP inequality for PP gives the stated lower bound with no missing case.

Theorems 4.1–4.2 and the main tail boundsCorrect and complete

The cited concentration inequalities are applied correctly

Pages 17–18 · Theorems 4.1–4.2 and proofs of Theorems 2.1(2), 2.2(2) · arXiv:2408.05963v4

For the scalar edge indicator, centering, the unit bound, the variance proxy, the IP gap, and the initial-density ratio all match the hypotheses of Huang–Li's Bernstein theorem. For symmetric counting, F~DμP\widetilde F-D_\mu P is Hermitian, has norm at most 22, and has matrix second moment at most 2I2I; the unordered path chain has the required positive absolute gap. Neeman–Shi–Ward's Bernstein theorem and its nonstationary corollary then give the paper's deliberately weakened constants, including the two-sided factor for the operator norm.

Huang–Li, Theorem 2.4
Theorem 2.2(2)Correct and complete

The matrix concentration input supports the displayed dependence

Page 18 · proof of Theorem 2.2(2) · arXiv:2408.05963v4

The source theorem uses the contraction parameter 1ηa1-\eta_a through (1+λ)/(1λ)(1+\lambda)/(1-\lambda) and a linear Bernstein term proportional to (1λ)1(1-\lambda)^{-1}. Substituting ηa(P~2)=η(P)/2\eta_a(\widetilde P_2)=\eta(P)/2, the variance proxy 22, the norm bound 22, and tt for the empirical average yields a bound at least as strong as the one printed in Theorem 2.2(2). The source's nonstationary extension supplies exactly the factor ν/μ\lVert\nu/\mu\rVert_\infty.

Neeman–Shi–Ward, matrix Bernstein theorem
Proofs of Theorems 2.1(3) and 2.2(3)Correct and complete

The covariance-polynomial estimates are complete

Pages 19–23 · Subsection 4.2 · arXiv:2408.05963v4

Stationarity turns each mean-square Frobenius error into a sum of lagged scalar covariances. The sum of squared centered-indicator norms is at most one, so an operator norm controls the entire expression without a dimension factor. For MLE, the path-space factorization and the invertibility of IPI-P on the mean-zero finite-dimensional subspace reduce the polynomial to a telescoping sum bounded by 4/n4/n. For SCE, all nonconstant eigenvalues of P~2\widetilde P_2 lie in [0,1)[0,1), and evaluating the nonnegative covariance polynomial at the largest one gives the stated constant.

03Novelty0 reported findingsNo non-novelty findings

No non-novelty findings.

Detailed audit reportFull reasoning, manuscript locations, and references.
Open report PDF ↗

Author response

Challenge an audit finding

Local workflow preview

A listed author may submit formal evidence that an audit is inaccurate. The response would be considered in a fresh AI re-evaluation; it would not edit the audit automatically.

Paper
arXiv:2408.05963v4
Authors listed
De Huang, Xiangyuan Li
Audit date
August 18, 2026
  1. 01Establish identityMatch an authenticated scholarly identity to this paper.
  2. 02Submit evidenceIdentify the finding and give a formal mathematical response.
  3. 03Re-evaluateA separate agent checks the response and records a disposition.
Recommended production method

Authenticate with ORCID, then require an exact arXiv match

MathAudit should accept the identity only when ORCID OAuth authenticates the claimant's iD and this exact arXiv paper appears in arXiv's public authority feed for that iD. A matching name alone is not sufficient.

ORCID OAuth and arXiv authority-record lookup are not connected in this local prototype.

Email fallback for papers without a linked ORCID

A production fallback could send a one-time link only when the submitted address matches an independently maintained author-contact allowlist for this paper. MathAudit must return the same message for every address so the form cannot reveal which contacts are on that list.

This demonstration does not send, store, or compare the address.

Structured response preview

This form remains unavailable until production identity verification succeeds. Nothing entered here is submitted.

This panel never establishes authorship in the local prototype. A production result should be described narrowly as an authenticated ORCID match or control of a separately allowlisted author-contact mailbox.