arXiv:2408.05963v4
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
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
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.
The MLE joint-count estimates hold
Pages 4–5 and 17–21 · Theorem 2.1 and its proof · arXiv:2408.05963v4
Writing as the empirical mean of edge indicators on the length-two path chain gives stationary mean . The ergodic theorem yields almost-sure convergence. The cited scalar Bernstein inequality applies with variance at most , and Theorem 3.3 supplies . For the mean-square estimate, stationary covariances reduce to a polynomial in ; the factorization , reduces its norm to a polynomial in , and telescoping gives the claimed bound. The alternative absolute-gap and reversible bounds follow from the same polynomial estimate.
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 makes its invariant edge law reversible, and the symmetrized matrix observable has expectation . The path-chain absolute gap is , 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 , with the stated initial-density factor in the nonstationary case.
The path-space spectral identities are correct
Pages 11–16 · Theorems 3.3–3.4 · arXiv:2408.05963v4
The rectangular factorizations give and , so the nonzero eigenvalues agree. After weighting by the invariant measures, has orthonormal columns and has orthonormal rows; consequently and have the same nonzero singular values. The positive and negative parts of the additive symmetrization identify the reversible spectrum with those of and , while the conditional-expectation argument yields the claimed IP-gap lower bound. For the unordered path chain, reversing before extension gives , proving both spectral-gap identities and uniqueness of the invariant law.
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 and by a supremum or infimum over nonzero vectors in . If , that set is empty and the displayed gaps and subsequent fractions are not finite quantities under the stated definitions. Add the standing assumption . 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.
The singular-value and gap comparison proof is complete
Pages 12–16 · proof of Theorem 3.3 · arXiv:2408.05963v4
The identities and follow directly from stationarity and the edge measure. They preserve the nonzero singular values in , 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 on the first coordinate and applying the IP inequality for gives the stated lower bound with no missing case.
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, is Hermitian, has norm at most , and has matrix second moment at most ; 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 ↗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 through and a linear Bernstein term proportional to . Substituting , the variance proxy , the norm bound , and 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 .
Neeman–Shi–Ward, matrix Bernstein theorem ↗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 on the mean-zero finite-dimensional subspace reduce the polynomial to a telescoping sum bounded by . For SCE, all nonconstant eigenvalues of lie in , and evaluating the nonnegative covariance polynomial at the largest one gives the stated constant.
03Novelty0 reported findingsNo non-novelty findings
No non-novelty findings.