How statement and proof provenance work
The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.
- Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
- AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
- AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.
These labels describe origin, not correctness: citations and verification chips remain separate evidence.
MCMC central limit theorems, MCSE, and long-run variance
1 · Prerequisites
- Absolute and Conditional Convergence; Rearrangement; Products
- Absolute Continuity and the Sharp Fundamental Theorem of Calculus
- Algebraic Closure, Embeddings, and Separability
- Algebraic Extensions, Extension Degree, and Finite Fields
- Analyticity of Holomorphic Functions; Liouville and Morera
- Arc Length and Rectifiable Curves
- Areas of Elementary Plane Figures
- Binary Operations, Monoids, Groups and Subgroups
- Bounded Variation and the Riemann–Stieltjes Integral
- Brownian Motion Construction and Continuity
- Brownian Motion, Markov Properties and Hitting Times
- Brownian Path Properties
- Central Limit Theorems
- Characteristic Functions Inversion and Continuity
- Compactness
- Compactness in Metric Spaces
- Complete Metrizability, Čech-Completeness, and Baire Category
- Completeness, Completion, and Uniform Continuity
- Complex Coordinate and Path Transport
- Complex Differentiability and the Cauchy–Riemann Equations
- Complex Lp Spaces and Test-Function Conventions
- Complex Power Series and Analytic Functions
- Composition Series, the Jordan–Hölder Theorem and Solvable Groups
- Concentration
- Conditional Distributions and Regular Conditional Probability
- Conditional Expectation
- Congruences, the Integers Modulo n and the Chinese Remainder Theorem
- Conjugate families and Gaussian conditioning
- Connectedness
- Construction of the Natural Numbers
- Construction of the Real Numbers via Cauchy Sequences
- Construction of the Real Numbers via Dedekind Cuts
- Continuity, IVT, EVT, and Uniform Continuity
- Contour Integration
- Convexity
- Cosets, Index and Lagrange's Theorem
- Countability and Uncountability
- Countability Axioms and Cardinal Functions
- Cyclic Groups and Direct Products
- Darboux, L'Hôpital, and Taylor's Theorem
- Density Separability and Convolution in Lᵖ
- Determinants of Matrices over a Commutative Ring
- Diagonalisation and the Minimal Polynomial
- Differentiation of Monotone Functions and the Vitali Covering Theorem
- Direct Matrix Factorisations: LU, Cholesky and QR
- Discrete Time Martingales
- Divisibility, Euclidean Domains, Principal Ideal Domains and Unique Factorisation
- Divisibility, Greatest Common Divisors and Bézout's Identity
- Dual Spaces, Bilinear and Quadratic Forms, and Sylvester's Law of Inertia
- Eigenvalues, Eigenvectors and the Characteristic Polynomial
- Filters and Ultrafilters
- Finite Counting, Factorials and Binomial Coefficients
- Finite Probability and the Probabilistic Method
- Finite Probability Spaces and Random Variables
- Finite-history linear prediction
- Foundations of the Real Numbers for Analysis
- Fourier Transform Convolution and Approximate Identities
- Fubini and Change of Variables
- Function Space Topologies and the Exponential Law
- Fundamental Trigonometric Identities
- Gaussian Elimination, Elementary Matrices and Reduced Row Echelon Form
- Gaussian Orthogonal Decomposition and Quadratic Forms
- Gaussian-process priors and finite posterior laws
- General-state Markov kernels, invariance, and reversibility
- General-state MCMC ergodicity and mixing
- Goursat's Theorem and Cauchy's Theorem in a Convex Domain
- Group Actions, Orbits, Stabilisers and Cayley's Theorem
- Group Homomorphisms and the Isomorphism Theorems
- Hausdorff via the Diagonal
- Ideals, Quotient Rings and the Isomorphism Theorems for Rings
- Improper and Parameter-Dependent Multiple Integrals
- Improper Integrals
- Independence Borel Cantelli and Zero One Laws
- Infinite Product Measures and Kolmogorov Extension
- Inner Product Spaces, Gram-Schmidt, Projections and Adjoints
- Lebesgue Measure on Euclidean Space
- Lebesgue-Stieltjes Measures and Distribution Functions
- Limits of Real Functions
- limsup, liminf, and Subsequential Limits
- Line Integrals and the Gradient Theorem
- Linear Independence, Bases and Dimension
- Linear Transformations, Rank-Nullity and Quotient Spaces
- Loss, Risk and Decision Theory
- Markov Kernels and Markov Chains
- Martingale Inequalities and Convergence
- Matrices, the Matrix of a Linear Map, and Change of Basis
- Measurable Functions and Simple Approximation
- Measure-Preserving Systems and Mixing Criteria
- Measures and Their Basic Properties
- Metric Spaces
- Mixed effects and longitudinal models
- Mixed Partials, Taylor Formulae, and Extrema
- Modes of Convergence Egorov and Lusin
- Modes of Convergence for Random Variables
- Monotone Functions, Discontinuities, and Continuity Sets
- Monotone Sequences, Bolzano-Weierstrass, and Cauchy Completeness
- Normal Subgroups and Quotient Groups
- Normed and Banach Spaces
- Order, Zorn's Lemma, and the Axiom of Choice
- Outer Measure and the Caratheodory Extension Theorem
- Partitions of Unity and Paracompactness
- pi: the Equivalent Characterizations
- Pointwise Extended-Real Function Products
- Poisson and Cox point processes
- Polyhedral selective inference
- Polynomial Rings, the Division Algorithm and Roots
- Power Series and Real-Analytic Functions
- Primes, Euclid's Lemma and the Fundamental Theorem of Arithmetic
- Probability Spaces Random Variables and Expectation
- Product Measures and the Fubini Tonelli Theorems
- Properties of the Integral and the Working FTC
- Radon Measures and the Riesz Markov Kakutani Theorem
- Recurrence Transience and Hitting Times for Markov Chains
- Reeb Stability and Global Foliation Constructions
- Relations, Functions, and Quotients
- Rings, Subrings, Integral Domains and Fields
- Rⁿ as a Normed Space; Vector-Valued Functions
- Roots, Rational Powers, and Classical Inequalities
- SCM, DAGs, adjustment
- Separation Axioms: the Hierarchy
- Sequences and Limits
- Sequences and Series of Functions; Uniform Convergence
- Series: Convergence and the Nonnegative Tests
- Sigma Algebras and Borel Sets
- Signed and Complex Measures Hahn and Jordan
- Simple Field Extensions and the Construction of the Complex Numbers
- Sine, Cosine, and the Definition of Pi
- Splitting Fields
- Standard-Borel Real Codings and Determining Classes
- State-space model laws
- Stationary Markov Chains and Ergodic Limits
- Statistical Experiments, Models and Identifiability
- Stopping Times and Optional Stopping
- Subspaces, Products, and Quotients
- Suprema and Infima
- Sylow's Theorems, p-Groups and Nilpotent Groups
- Symmetric Groups, Cycle Decomposition and the Sign Homomorphism
- Temporal intensity and likelihood
- The Cantor Set, Baire Category, and Measure Zero in ℝ
- The Complex Exponential and Euler's Formula
- The Derivative and the Mean Value Theorems
- The Determinant of a Linear Operator, Cofactors and Cramer's Rule
- The Ergodic Theorems of von Neumann and Birkhoff
- The Exponential Function
- The Fundamental Theorem of Algebra
- The Fundamental Theorem of Finite Abelian Groups
- The Fundamental Theorems of Calculus
- The Galois Correspondence
- The Identity Theorem, the Maximum Principle and the Open Mapping Theorem
- The Inverse and Implicit Function Theorems
- The Lebesgue and Riemann Integrals Compared
- The Lebesgue Integral and the Convergence Theorems
- The Logarithm and General Powers
- The Lᵖ Spaces Holder Minkowski and Riesz Fischer
- The Maximal Function and Lebesgue Differentiation
- The Radon Nikodym Theorem and Lebesgue Decomposition
- The Riemann Integral in Rᵐ and Jordan Content
- The Riemann Integral: Definition and Integrability
- The Spectral Theorem, Positive Operators and Singular Value Decomposition
- The Topology of Euclidean Space
- The Total Derivative in ℝᵐ → ℝⁿ
- The ZFC Axioms and the Basic Set Constructions
- Topological Spaces and Continuity
- Topology of ℝ
- Triangularisation, Generalised Eigenspaces and Jordan Canonical Form
- Urysohn's Lemma and the Tietze Extension Theorem
- Vector Spaces, Linear Subspaces, Span and Direct Sums
- Weak Convergence Tightness and Representation
- Weak Laws and Series of Independent Random Variables
2 · Summary
Summary
This page defines covariance long-run variance, three stationary mixing coefficients, geometric and polynomial total-variation rates, one-step minorization, and the overlapping-batch-means estimator. It proves the kernel conditional-moment lemma, rate-to-mixing bounds, measurable-factor bounds, the degenerate zero-variance reduction, and the Poisson-equation martingale CLT with its covariance identity.
The six MCMC CLT regimes from Jones (2004) are proved locally with their moment, reversibility and rate-envelope hypotheses. Their suppliers prove the bounded and unbounded strong-mixing limits, the rho and uniform-mixing limits, the reversible geometric correlation bound and the integrable geometric envelope. The Flegal–Jones OBM contract retains every numerical batch-size condition and corrects the source's incompatible simultaneous monotonicity requirements for integer batches. Local Gaussian quadratic-form bounds and a pathwise Brownian-approximation comparison prove the OBM reduction, including zero approximation variance and the exact floor-power batch calculation. The full Markov-chain OBM theorem combines the positive-variance strong approximation with the separate zero-variance strong rate, applies under arbitrary initial laws, and identifies its covariance long-run variance through the stationary mixing CLT. A finite second moment alone is not presented as a CLT criterion, and covariance summability is not used to identify a CLT variance without a theorem.
3 · Logical flowchart
4 · Definitions, theorems and proofs
Flegal–Jones overlapping-batch size conditions
Definition
Let and let be positive integers. Say that the sequence satisfies the Flegal–Jones OBM batch-size conditions for exponent , with the integer monotonicity correction described below, when for all sufficiently large , , , and there is a constant such that while there are and with for all , and
These are all the numerical rate, summability and endpoint hypotheses of Flegal–Jones Theorem 2. The source additionally asks that both and be nondecreasing. Those two requirements are omitted here: for an integer sequence with eventually they are incompatible with . Indeed, ratio monotonicity would imply , while integer monotonicity would force the nonnegative integer to be zero. Thus the literal source conditions would make eventually constant. The corrected interface preserves the intended floor-power batch sequences, and the Brownian quadratic-form proof of OBM consistency needs neither monotonicity assumption.
This is a sufficient-condition interface; failure of a condition does not imply inconsistency of the estimator.
For the floor-power sequence , the intended source range and satisfies this corrected contract. Indeed, eventually . Choose to bound the summability series by . The other displayed quantities are bounded respectively by , , , and , all tending to zero. Also , and eventually. These checks establish the numerical contract without asserting the impossible simultaneous monotonicity.
Geometric and uniform ergodicity for MCMC CLTs
Definition
Let be a probability kernel with invariant probability on a measurable space . Write for its -step transition probability and use for probability measures.
The chain is geometrically ergodic if there are a measurable finite function and such that
It is uniformly ergodic if this bound holds with bounded on . These are convergence-rate properties, not conclusions about a particular function's CLT or long-run variance. If a proof needs , that integrability is an additional stated hypothesis unless a cited theorem establishes it.
Long-run variance for a stationary Markov chain
Definition
Let be a stationary time-homogeneous Markov chain with invariant probability , and let be real Borel with and . For , set . When the Markov-chain CLT holds with finite variance parameter and , the covariance long-run variance is
This definition does not infer a CLT from covariance summability, or infer that a weak-limit variance equals the covariance sum. Whenever is used as the variance parameter in a CLT or MCSE claim, that equality must be established by the theorem being applied. For example, the local Poisson-equation theorem proves the equality under its square-integrable solution hypotheses.
Remarks
Under stationarity, the finite variance identity is Absolute summability makes the right side converge to by dominated convergence. Equality between this variance limit and a variance parameter from weak convergence is a separate assertion; it is not part of this definition.
Strong, uniform, and maximal-correlation mixing coefficients
Definition
Let be a process on a probability space whose finite-dimensional laws are invariant under common nonnegative time shifts; call this strict stationarity. For , set and . For , define the strong-mixing coefficient by
Define the uniform-mixing coefficient by
For square-integrable real random variables with positive variances, write For each , let consist of the triples such that , are real and square-integrable, and are measurable with respect to and respectively, and both have positive variance. Define the maximal-correlation coefficient by where square-integrability is in the sense of The space as the quotient by null functions. Thus when there is no admissible pair, as for a constant process. The process is called strongly mixing, uniformly mixing, or rho-mixing, respectively, when the corresponding coefficient tends to zero. These coefficients are distinct conditions; no implication between them is included in this definition.
Polynomial ergodicity with an integrable rate envelope
Definition
Let be a probability kernel with invariant probability . For , say that has polynomial total-variation rate of order with an integrable envelope if there is a measurable with such that
This is the specific envelope condition used with the polynomial-rate CLT cases on this page. It does not claim that every definition of polynomial ergodicity in the literature includes integrability of .
State-dependent one-step minorization
Definition
Let be a probability kernel on . A state-dependent one-step minorization is supplied measurable data and a probability measure on such that
The function need not have a positive uniform lower bound; this is the exact one-step assumption used in the overlapping-batch-means theorem below.
Kernel conditional expectations for square-integrable Markov observables
Statement
Assume AC. Let be a stationary time-homogeneous Markov chain with probability kernel , invariant law , and natural filtration . Let (The space as the quotient by null functions) and fix a finite real measurable representative, also denoted by . Define and on the -full measurable set where the respective integrals are finite, and set these functions to zero on the exceptional set. Then, for each , Consequently, if , then is a square-integrable martingale-difference sequence and
Facts & Assumptions
Given: AC, a stationary chain with kernel and invariant law , its natural filtration, and a real measurable representative with .
AC is the choice principle used by the library's conditional-expectation class convention. (The Axiom of Choice)
For a bounded measurable and a -chain, ; in particular, at the conditional mean is the kernel integral. (Chapman-Kolmogorov equations)
Increasing nonnegative inputs have increasing conditional-expectation versions converging almost surely to the conditional expectation of their limit. (Conditional monotone convergence)
Invariance is the event identity for every measurable . (Invariant and stationary distribution for a Markov kernel)
Nonnegative increasing functions pass through an integral by monotone convergence. (Monotone convergence for the integral)
Kernel integration of a nonnegative product-measurable function is measurable, with infinity permitted. (Measurability of integration against a kernel)
Conditional expectation is linear and preserves order. (Basic algebra and order properties of conditional expectation)
Conditional Jensen gives for square-integrable . (Conditional jensen inequality)
On a probability space, Hölder gives . (Holder's inequality for integrals, including the endpoint cases)
An integrable -measurable factor can be taken out of a conditional expectation when the products are integrable. (Taking out what is known)
Products of square-integrable random variables are integrable by Hölder's inequality. (Holder's inequality for integrals, including the endpoint cases)
Proof
For bounded measurable , [F1] with gives . Kernel measurability [F5] makes the right side measurable.
Write . By [F8], both parts are integrable. For let . These functions are bounded and increase to . Step 1.1 gives the conditional kernel identity for each ; [F2] makes the conditional expectations increase to , while [F4] makes . Invariance [F3], first for nonnegative simple functions and then by [F4], gives ; hence both limits are finite almost surely. Subtracting the two finite identities using [F6] proves .
Apply the same argument to the nonnegative integrable function : the truncations increase to , and invariance gives . Thus . These identities do not depend on the chosen representative: if two representatives differ on a measurable -null set , invariance gives , hence for -almost every ; their kernel integrals agree there, and their evaluations at each agree almost surely.
By [F7] and step 2.2, almost surely, so stationarity and [F3] imply ; therefore . Step 2.1, the preceding square-integrability bound, conditional linearity and [F9] for the -measurable factor give . The product is integrable by the two bounds and [F10], so [F9] also gives its conditional expectation as . Expanding the square and using step 2.2 now yields . Thus the sequence is a square-integrable martingale-difference sequence.
Source notes
Jones §2 defines the kernel operator and states the conditional one-step kernel convention. The bounded-to-square-integrable extension and the exceptional-set version are proved here; no pointwise identity is asserted outside a -full set.
A kernel has a countably generated stable sub-sigma-field containing a given observable
Statement
Let be a probability kernel on an arbitrary measurable space , and let be measurable. There is a countably generated sub-sigma-field such that is -measurable and , restricted to events in , is a probability kernel on . For every , , and , If is invariant for , then is invariant for .
Facts & Assumptions
Given: The measurable space, probability kernel, and real measurable function in the Statement.
Measure kernel and probability kernel: is a probability and is measurable for every .
Dynkin's pi-lambda theorem: a lambda-system containing a pi-system contains the sigma-field it generates.
Monotone convergence for the integral: increasing nonnegative integrands have increasing integrals with the expected limit.
Iterated transition kernels: iterated transitions are defined by kernel integration, with the identity kernel at time zero.
Invariant and stationary distribution for a Markov kernel: invariance means for every measurable .
A measurable function between measurable spaces: a real function is measurable when inverse images of Borel sets are measurable.
Every nonnegative measurable function admits an explicit increasing sequence of simple approximations: nonnegative measurable functions have increasing nonnegative simple approximations.
Proof
Given: The hypotheses in the Statement.
Let be the algebra generated by . Inductively, let be the algebra generated by and all sets with and . These are subsets in by [F1] and [F6]. Each algebra is countable: enumerate its generators and then all finite Boolean expressions in them. This also explicitly enumerates . Put . Since rational half-lines generate the real Borel sigma-field, is -measurable.
For , all rational superlevel sets of belong to by construction. The class of for which is -measurable contains and . It is closed under complements because and under countable disjoint unions because is a measurable increasing limit. Thus is a lambda-system, and [F2] gives . The restriction therefore is a probability kernel.
For every nonnegative -measurable , integration against agrees with integration against : this holds for indicators by restriction, for nonnegative simple functions by finite additivity of the integral, and for general by [F7] and [F3]. Both integrals are allowed to be infinite. Induction in [F4], using this integration identity and step 2.1, gives for every , including .
If is invariant, the same restriction identity for integrals gives for . This is invariance of by [F5]. All generating enumerations in step 1.1 were explicit, so this construction introduces no choice assumption.
Centered Gaussian quadratic forms have dimension-free even-moment bounds
Statement
Assume AC. For every positive integer there is a finite constant , depending only on , such that for every integer , independent standard-normal variables , and real symmetric matrix , In particular . The matrix may be singular or indefinite.
Facts & Assumptions
Given: AC, a positive integer , independent standard-normal coordinates of , and a real symmetric matrix with .
AC is assumed as in The Axiom of Choice and the normal interfaces below.
A self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis
Orthogonal transformations preserve independence and the standard-normal coordinate laws. Orthogonal transformations preserve the standard normal product law
The standard-normal law has density . Standard normal and normal laws
Expectations of integrable products of functions of independent variables factor. Expectations factor over finite products of independent random variables
Proof
By [F1] choose an orthogonal matrix with . By [F2], has independent standard-normal coordinates. Integration by parts in the density [F3], with at both infinities, gives . Every absolute moment of is finite: for an integer , , and this polynomial is integrable against . Thus , , and .
Expand the even power over ordered tuples of indices. By [F4], a term whose tuple has an index appearing exactly once has expectation zero. Group all other terms by the partition of their positions into equal-index blocks. Such a partition has block sizes with . For each fixed partition, taking absolute values and dropping the restriction that different blocks have distinct indices bounds its contribution by . All these quantities are finite by step 1.1.
Put . If , the centered form is identically zero and the bound holds. Otherwise and , so for every integer , . Consequently each partition contribution in step 2.1 is bounded by . There are finitely many partitions of the fixed set of positions; the sum of these moment products over the partitions without singleton blocks is a finite constant depending only on .
Since is even, the expanded power is the absolute th power. Summing the bounds of step 3.1 therefore gives . Finally by symmetry and orthogonal diagonalization. This proves the bound and the expectation assertion, including , singular matrices and matrices of either sign.
Summable conditional projections give a stationary martingale CLT
Statement
Assume AC. Let be an invertible bimeasurable probability-preserving transformation on such that implies . Let be real and centered, put for , and put . Suppose Then the series converges. Put and . The are stationary square-integrable martingale differences for , and for , With , including . If additionally , then
Facts & Assumptions
Given: AC, the invertible probability-preserving ergodic transformation , centered real , its stationary translates, natural past filtration, and the displayed summability assumption.
AC supplies the inherited choices for the following conditional-expectation and limit suppliers. (The Axiom of Choice)
Conditional expectation is linear, preserves expectation and constants, and contracts the real norm. (Basic algebra and order properties of conditional expectation, Conditional lp contraction)
Conditional expectations satisfy the tower property and allow a known measurable factor to be taken out when the input and the required products are integrable. (Tower property of conditional expectation, Taking out what is known)
Under countable choice is complete. (Riesz-Fischer completeness of for )
Square-integrable products satisfy Cauchy–Schwarz. (Cauchy-Schwarz for random variables)
A stationary integrable process with ergodic canonical shift has averages converging almost surely to its expectation. Ergodicity requires invariant events to have probability zero or one. (Birkhoff limit for a stationary integrable process, Ergodicity relative to an invariant measure)
A square-integrable martingale-difference array with terminal conditional variance clock tending to in probability and vanishing unconditional Lindeberg sums has terminal sum converging to . (Martingale central limit theorem)
Dominated convergence passes limits through integrals dominated by an integrable function. (Dominated convergence)
convergence implies convergence in probability; a perturbation tending to zero in probability preserves a weak limit; continuous scaling preserves weak convergence. ( convergence implies convergence in probability, Converging together lemma, Continuous mapping theorem)
Variance and covariance are centered second and product moments. (Moments, variance, and covariance on a probability space)
Proof
For an integrable and any integer , the variable is measurable because the pullback of by is . Every event in is the pullback of an event in , since is invertible and bimeasurable. Changing variables under the probability-preserving bijection proves that the integrals of this candidate and agree over every event. Uniqueness of conditional expectation therefore gives as almost-everywhere classes.
By summability and [F3], the partial sums defining are Cauchy in and have a limit there. Composition with is an isometry, so step 1.1 yields in , with measurable in and . Representatives may be fixed for the countably many indices; discard a countable union of null sets for their identities.
By [F1], conditional expectation is continuous in . Apply [F2]'s tower property to the partial sums in step 2.1 and pass to the limit to obtain . Thus is square-integrable and measurable, with . Also , so the differences are stationary. Summing their defining identity gives , where , , and . In particular in and hence in probability by [F8].
Put . It is nonnegative, integrable, and has expectation by [F1]. Step 1.1 gives . The map intertwines with the canonical one-sided shift. Hence a shift-invariant event for this process pulls back to a -invariant event modulo null sets, whose probability is zero or one by the assumed ergodicity. The process is stationary, integrable and ergodic, so [F5] gives almost surely. The limiting variance is therefore a constant.
If , let for , and let for ; use the filtration for . Step 3.1 makes each row a square-integrable martingale-difference row. The row sums and clocks are eventually constant, and step 4.1 shows their terminal clocks converge almost surely to . For every , stationarity gives by [F7]. Thus [F6] gives . Scaling and the perturbation in step 3.1 give by [F8].
If , stationarity implies for every , so almost surely. Step 3.1 gives in and hence in probability. The constant sequence has weak limit the point mass ; [F8] transfers that limit to . This establishes the asserted CLT in both variance cases.
Suppose the additional covariance summability holds, and put . Stationarity and expansion of the finite square give , by [F7] on counting measure with summable dominator . For , is measurable, so [F2], [F4], and step 3.1 give . Therefore . The identity and [F4] imply , since is bounded. All are centered, so [F9] identifies the two variance limits and proves the stated covariance formula, including .
Source notes
The absolutely summable projection assumption permits an exact coboundary series, so no maximal martingale approximation theorem is imported. The deterministic variance follows from the explicitly assumed ergodicity and its factor argument. The covariance identity is obtained from second moments, independently of convergence in distribution.
Finite moments control increments over short moving windows
Statement
Fix . Let be identically distributed centered real random variables with . Assume either that they are independent, or that they are one-dependent, meaning that and are independent for every . Put and .
There is an event of probability one on which, for every and , with for ,
Facts & Assumptions
Given: , the centered sequence, its common finite th absolute moment, and one of the two dependence hypotheses in the Statement.
Independent centered real variables bounded by satisfy the two-sided Bernstein bound when the sum of their variances is positive; when the sum is zero almost surely. Bernstein inequality for bounded independent sums
For a nonnegative random variable and , . Markov's inequality for random variables
If the sum of event probabilities is finite, only finitely many of the events occur almost surely; independence is unnecessary. First Borel-Cantelli lemma for events
Independence is independence of the generated sigma-algebras; iid means independence together with equality of marginal laws. Independent random elements, Identical distribution and IID families
For centered square-integrable , . Moments, variance, and covariance on a probability space
Proof
First assume independence and fix positive integers and with . Choose strictly between and , put , and, for dyadic integers , put , , and . The second moment is finite since . If , the nonnegative variable has zero expectation and is zero almost surely: for each integer , [F2] bounds by zero, and the union of these events is . Equality of marginal laws and countable subadditivity make every zero simultaneously almost surely. This proves the assertion in this case. Henceforth assume .
The union bound and [F2] give . Since , these probabilities are summable over dyadic . By [F3], eventually no one of the first variables is truncated at level , almost surely along the dyadic horizons. Define and . For each fixed , these variables are independent and centered, , and . Moreover centering of gives .
For an interval of at most indices, [F1] bounds the tail of its sum by ; if its actual variance is zero the bound still holds because the tail is zero. Set , where . Then and , so the bound is at most . There are at most intervals with integer endpoints in . The probability that any interval of length at most violates the bound is therefore at most , summable over dyadic . By [F3], all such interval sums are bounded by eventually almost surely.
On the eventual no-truncation event from step 2.1, every original interval sum differs from its centered truncated sum by at most . Now because , and because . Also . Thus step 3.1 proves the claimed little-o bound at dyadic horizons. These comparisons follow by dividing by : in each expression a strictly negative power of multiplies a fixed power of , which tends to zero.
For , monotonicity of gives , so every admissible interval for is included among those for . Since , step 4.1 proves the assertion at every horizon . Intersect the probability-one events obtained for the countably many integer pairs . For arbitrary real , choose integers and ; when , makes the real-parameter window no longer than the integer-parameter window. The conclusion therefore holds on one event for every in the Statement.
Under one-dependence, each parity subsequence is independent. Indeed, the sigma-algebra of any finite initial collection of even variables is contained in , which is independent of by the given separation; induction on the number of factors proves the finite rectangle factorization for all even variables. The identical argument applies to the odd variables. Both subsequences retain the common centered marginal and its th moment. Apply step 5.1 to each subsequence and intersect their two probability-one events. An interval of at most original indices splits into one contiguous interval of each parity, each with at most terms and with its subsequence endpoints at most . Each component is by the independent result, uniformly over the original intervals; their sum has the same bound. No independence between the parity subsequences is asserted or needed. All windows are nonempty by ; empty component windows contribute zero. This completes both dependence cases, including zero variance and every permitted window parameter.
Remarks
The estimate also controls short calendar-time reward increments once a genuine cycle decomposition is supplied. Suppose cycle lengths are positive integers, centered complete-cycle rewards satisfy the lemma, and each partial-cycle reward is bounded by an identically distributed nonnegative envelope with . Dyadic Markov bounds as in step 1.2, at levels with , give almost surely. An interval of at most calendar steps crosses at most complete or partial cycles. Its full-cycle contribution is controlled by the lemma and its two endpoint pieces by twice the maximum envelope. Thus its centered reward increment is uniformly up to time . This remark is conditional on the actual decomposition and its moment/dependence hypotheses; it asserts no regeneration theorem or Brownian coupling.
The cited full source's Lemma 5 was inspected. Its proof uses Nagaev's inequality for a probability estimate. The proof here instead uses the already proved local Bernstein supplier and the spare moment ; no unproved Nagaev or strong-approximation theorem is used.
Independent Brownian motion straightens a Poisson clock
Statement
Assume AC. Let be a sigma-algebra on a probability space, let , and let be -measurable at each time. Suppose , its paths are nondecreasing right-continuous integer-valued functions finite on bounded intervals, its increments on disjoint time intervals are independent, and has law for . Let be a standard Brownian motion independent of .
On an extension preserving the original probability space there is a standard Brownian motion , independent of , such that The almost-sure bound has a finite random constant and a finite random starting index. In particular is independent of the entire Poisson clock; and are not asserted to be independent.
Facts & Assumptions
Given: AC, , , , and the independent Brownian motion satisfying the Statement.
AC is assumed. The Axiom of Choice
Brownian motion is continuous on a common probability-one event and is centered Gaussian with covariance ; conversely that covariance with common path continuity gives Brownian motion. Brownian motion, Brownian covariance is equivalent to independent stationary normal increments
Finite jointly Gaussian blocks with zero cross-covariance are independent, including singular blocks. Uncorrelated jointly Gaussian blocks are independent
A standard Brownian motion exists, and the countable product construction adjoins countably many independent copies independently of the original space. Existence of continuous Brownian motion, Assuming countable and dependent choice, countable products of arbitrary probability spaces, Coordinate random elements of a countable product are independent
The Poisson mass at with parameter is ; at parameter zero it is the unit mass at zero. The exponential is its convergent power series and satisfies . The scalar Poisson count law, The real exponential function and the number by a power series, The exponential addition formula
For a standard normal and , ; reflection gives for Brownian , . Two-sided Mills bounds for the standard normal tail, Law of the Brownian maximum
Markov's inequality applies to nonnegative variables; summable event probabilities imply that only finitely many events occur almost surely. Markov's inequality for random variables, First Borel-Cantelli lemma for events
A differentiable map with uniformly bounded derivative on a convex open set is Lipschitz there. On a convex open set, a uniform bound implies
Agreement of finite measures on a generating pi-system extends to its sigma-algebra. Dynkin's pi-lambda theorem
Proof
Work on common full-measure events of Brownian continuity and the assumed clock paths. For each integer , [F4] gives . The events decrease in , so continuity of probability shows that is unbounded almost surely after taking the countable union over . By [A1] and [F3], extend the space by independent standard Brownian motions , independent also of the original variables. They will supply continuous bridges; the original space and all its joint laws are preserved by this product extension.
We prove the necessary uniform clock bounds directly. If , summing [F4]'s nonnegative series gives for every real . For , the exponential series and for imply . Exponential Markov and the choice , applied with either sign, therefore give when ; at the deviation is zero. With this implies for a constant and integer . This derivation uses only the Poisson mass series and Markov, not a strong approximation theorem.
Partition the original time axis into macrointervals . For each macrointerval , extract the root normal . For every dyadic cell inside a macrointerval, of length , extract its midpoint normal . This is a countable iid standard-normal family, independent of . To verify independence, the midpoint contrast on has covariance zero with the endpoint increment on and with all Brownian evaluations outside its interior: for the three covariance terms cancel, and for they are . Its variance is one by direct expansion of . Two dyadic cells are either disjoint or nested; in the nested case, the larger cell's contrast uses only evaluations outside the smaller cell's interior or at its endpoints. Hence distinct contrasts and roots have zero covariance. Every finite family is jointly Gaussian by [F1], so repeated use of [F2] proves mutual independence.
For dyadic horizons , write . Applying step 1.2 with and choosing a sufficiently large deterministic makes the probability at most . There are at most integer intervals in . The union bound and [F6] show that, almost surely eventually on dyadic horizons, every such interval satisfies . Choose a deterministic sufficiently large. For every dyadic cell of length , its normalized half-counts and then lie in . For every integer point in the cell, with and , one has . Indeed the interval estimate gives this bound first with an additional term , which is at most another constant times since . The same estimate implies that every integer interval of length at most has count increase at most . Constants here and below are deterministic and may depend on .
For a cell put , with half-counts and when . Prescribe a increment over the count interval of each macrointerval, starting from and concatenating the macrointervals. Given both endpoint values for , if prescribe its count-midpoint value as . When all three count times agree and the midpoint value is that endpoint. When only one half-count is zero the displayed prescription also equals the appropriate existing endpoint. Recurse down to unit cells. For a unit cell of count length , with left count time and prescribed endpoint increment , fill its count interval continuously by for . Zero-count leaves require no filling. Unboundedness from step 1.1 makes these intervals cover . All operations depend measurably on finitely many integer counts and Gaussian evaluations for each fixed : partition by the first integer with and its finite count history. The resulting countable pieces give measurable evaluations, including endpoints.
The tail bound in [F5] gives for . On the number of root and midpoint coefficients in step 2.1 is at most . A union bound with , followed by [F6] with sufficiently large, shows that their maximum absolute value is almost surely eventually along dyadic horizons. Also [F5] and the union bound give, for Brownian , when . Therefore a bridge over length has supremum tail bounded by for , since its supremum is at most twice that of . These estimates also apply to translated Brownian increment paths by [F1].
Fix a deterministic finite count history. The Gaussian split used in step 3.1 is exact: if an interval has total increment and child lengths , , independent gives child increments and . Their variances are and their covariance is . Thus they are independent by [F2], with zero-length increments identically zero. Recursing proves independent normal increments of the prescribed count lengths over all nonzero unit leaves, also across macrointervals. For the independent bridge , [F1] directly gives mean zero, continuous paths, endpoints zero, and covariance . Filling a leaf as in step 3.1 therefore gives centered Gaussian increments whose covariance at offsets is . Different leaves have independent increment processes because their endpoint increments are independent and all bridge motions are independent. Summing these increments shows that, for every finite evaluation list, has centered Gaussian covariance on this fixed count history.
Fix a large dyadic horizon and a macrointerval . If , stop its binary tree at cells with lengths in ; if , stop at the root. For an integer , consider the piecewise-linear interpolant of on this stopped grid, and the piecewise-linear interpolant of on the corresponding count grid, evaluated at . Repeated midpoint interpolation from step 3.1 expands each interpolant in the same root normal and in the midpoint normals of its ancestor chain. In an active cell of length containing in its left half, the coefficients of its midpoint normal are respectively and , with as in step 2.2. This follows by multiplying the left interpolation fraction by . In the right half the coefficients are obtained by replacing by . At a midpoint or endpoint these same formulas agree with the endpoint prescription.
For any finite list of evaluation times, partition the original space into countably many -events specifying a long enough finite integer count history to cover them. The Gaussian and added bridge variables are independent of ; step 4.1 shows that on every such count configuration the resulting evaluation law is the same Brownian Gaussian law. Consequently, for any and any Borel set of that finite evaluation vector, summing over this partition gives , where is its Brownian probability. No regular conditional probability on the original space is needed. The construction has common continuous paths: count intervals have integer endpoints, finitely many nonzero intervals meet each bounded count-time interval, their bridges are continuous, and endpoint prescriptions agree. Thus [F1] makes standard Brownian. Continuous evaluations are limits of rational evaluations, and [F8] extends the factorization from rational cylinders to . Therefore is independent of all of , not just of .
The function has uniformly bounded first derivatives for and : its displayed denominators are bounded away from zero, and differentiating gives only bounded products of these reciprocals and square roots. One may take a slightly larger open convex box still with positive and use [F7]. At its value is . Hence step 2.2 bounds the midpoint coefficient discrepancy by . The normalized right-half prefix differs from its deterministic value by at most another constant times , so the same bound applies there. The root coefficients for a large macrointerval are and ; the same bounded-derivative argument and step 2.2 bound their discrepancy by . By step 3.2, each shared Gaussian coefficient therefore contributes at most to the interpolant difference. At most coefficients occur on an ancestor chain, so the increment-interpolant error is .
We control the discarded bridge pieces uniformly. Each stopped cell has deterministic length at most and, on the count event from step 2.2, count length at most . The Brownian deviation from linear interpolation on a cell is a continuous bridge with the covariance computed in step 4.1, regardless of the endpoint values. Step 3.2 therefore bounds the probability that its supremum exceeds by . There are at most stopped cells. For the bridges the same estimate holds on every count history having these length bounds, by the conditional Gaussian calculation in steps 4.1, 5.1; summing over the count configurations gives the identical unconditional bound on the good count event. Choose large, add the summable failure probability of that count event, and apply [F6]. Thus all discarded and scaled bridge pieces are almost surely eventually. This argument includes zero count length, whose bridge is identically zero.
A completed large macrointerval contributes only its root error by steps 2.2, 3.2. For a macrointerval of length at most , its entire increment supremum and its scaled increment supremum are by the same reflection and count-history argument as step 6.1, applied without subtracting an interpolant. There are at most macrointervals before time . Add the root or small-macrointerval errors of the completed macrointervals, then the interpolant error from step 5.2 and the two leaf errors from step 6.1 in the last one. This gives almost surely along dyadic horizons.
For real , Brownian reflection and a union bound over the unit intervals make uniformly eventually almost surely. Step 2.2 gives . On every such count history, step 5.1 gives Brownian law for , and reflection bounds its increment supremum on each by except on a probability bounded by ; this is summable on dyadic horizons when is large. Integrating over the count configurations and using [F6] shows that these increments are uniformly . Monotonicity of places in that count interval. Add these two bounds to step 7.1 to obtain the required real-time supremum on dyadic horizons. For any real , the next dyadic horizon is less than , and its logarithm differs from by a bounded additive constant, proving the asserted bound at every horizon.
The finite initial time range is absorbed into the random constant: is bounded there by continuity, is bounded there by hypothesis, and is bounded on the corresponding compact count-time interval by continuity. Repeated or zero count times were treated in steps 3.1, 4.1. The restriction is needed for the displayed scaling and the lower large-cell count bounds; no assertion is made at rate zero. Step 5.1 proves the full claimed independence, while steps 7.1, 8.1 prove the approximation. AC was used precisely through [F1]–[F3] and the product extension in step 1.1; all concentration, interpolation, and countable-event arguments are local.
Source notes
The complete source PDF was retrieved. Lemma 2.4's full statement, its auxiliary-randomness clause, and the opening Haar construction were inspected, together with the introduction and proof roadmap. The source's detailed exponential-tail estimates were not represented as independently verified. The proof above supplies its own weaker logarithmic-square bound, which is sufficient for every fixed positive polynomial approximation scale. It invokes no KMT, embedding, renewal approximation, or unproved strong-approximation result.
Gaussian anchored bounded sums have uniform cumulant and Fourier bounds
Statement
Assume AC. Fix . Let , , , and let be independent centered random vectors in , where , almost surely, and . Suppose has covariance and in quadratic form order. There is an absolute constant such that, with , the logarithmic moment generating function has a single analytic branch on the complex Euclidean ball , and Here uses the complex bilinear dot product and . Every real exponential tilt is defined, and its characteristic function satisfies where for real , also outside the analytic ball. Put . For every fixed and integer there is a finite , uniform in , such that, if , All assertions hold with the same , , , and for , where is any real matrix with and . This includes every coordinate prefix, every orthogonal change of coordinates, and their compositions. In particular it covers the normalized joint sum and difference of two independent children by applying the assertion in their joint dimension, provided that dimension is at most four and the joint input satisfies the stated covariance assumptions.
Facts & Assumptions
Given: The hypotheses of the Statement; all vector norms are Euclidean, and all matrix orders are orders of real quadratic forms.
Finite products of integrable functions of independent variables factor in expectation. Expectations factor over finite products of independent random variables.
Dominated convergence interchanges limits and integrals under an integrable absolute majorant. Dominated convergence.
Absolute Fubini and nonnegative Tonelli apply to the probability and Lebesgue integrals below. Fubini's theorem for L^1 functions on a sigma-finite product, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product.
A positive definite Gaussian law has density and a Gaussian with covariance has characteristic function . A positive-definite multivariate normal law has its Gaussian density, Characteristic function of a multivariate normal law.
Locally uniformly convergent holomorphic series are holomorphic, the one-variable identity theorem extends equality from a real interval, and linear changes of Lebesgue variables have their determinant Jacobian. A locally uniformly convergent series of holomorphic functions may be differentiated term by term, Identity theorem for holomorphic functions, A C^1 diffeomorphism satisfies the change-of-variables formula for nonnegative Lebesgue measurable functions.
Euclidean Cauchy–Schwarz bounds , and its integral version follows by integrating and minimizing in real . Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation.
AC is assumed and inherited by the Gaussian suppliers. The Axiom of Choice.
Proof
Fix one centered bounded summand with covariance , and put . Its exponential series converges uniformly on every bounded complex set because , and the expectation of its absolute series is bounded by ; thus termwise integration defines an entire function and its derivatives are the expectations of the corresponding polynomial factors times the exponential. Fix an absolute with . If , then , so the convergent power series for supplies a single analytic logarithm vanishing at zero.
Completing the square in the density in [F4] shows for real and that its real tilt is . The same integral for complex converges absolutely on bounded complex sets: times the Gaussian density is a constant times a translated Gaussian density, which dominates its polynomial derivatives on smaller bounded sets. Its exponential series or differentiated integral therefore gives an entire function, and the real identity extends to complex arguments. To justify that last extension without a several-variable identity theorem, fix all but the first coordinate real and apply the one-variable identity theorem, then successively free each remaining coordinate. Thus the Gaussian cumulant is and has zero third derivative.
On this ball, writing subscripts for directional derivatives, the bounds , , and give , , , , and . These inequalities also hold when is singular: a direction with zero second moment has almost surely.
Every real tilt exists because all bounded-summand moment generating functions and the Gaussian one are finite and strictly positive. Its density with respect to the original independent joint law factors as the product of the separately normalized tilt densities. Hence those summands remain independent under the tilt. Step 1.2 shows that the Gaussian factor retains covariance , while each other factor has a characteristic function of modulus at most one because . Factorization and [F4] therefore give for every real .
Direct differentiation gives . Since , step 2.1 bounds this by for one absolute finite . Choose so that and , and set . Each bounded summand therefore has the asserted derivative bound on , using its own covariance.
Substitute in each requested tail integral, using the constant Jacobian . For , . For every fixed , the function is bounded on by differentiation, so the extracted exponential is at most . The integral of is finite by the normal density normalization in [F4]; its product with is integrable as well, since is bounded and the remaining Gaussian is integrable. Using gives the two tail bounds, with an additional in the weighted bound. Taking the maximum of these constants over makes them uniform.
By independence [F1], the moment generating function of is the product of the factors in steps 1.1 and 1.2. On , its logarithm is the sum of the analytic logarithms just constructed, with value zero at zero. Since the covariances add and , summing step 3.1 gives . This proves the analytic assertion, including the empty bounded sum.
If , the operator norm of is one: is an isometry and by [F5]. Thus remains bounded by and centered, its sum has covariance , and is Gaussian with covariance by [F4]. Independence is preserved by applying measurable maps separately. All earlier steps therefore apply to with the same . The map is orthogonal because it preserves by expansion; consequently its coordinate prefixes have the required row isometry property. This proves all assertions.
Overlapping batch means for MCMC standard errors
Definition
Inputs: Real outputs , an integer , and an integer batch size with . For an MCMC run, for the chosen real-valued observable .
Outputs: The sample mean , the overlapping-batch-means estimate , and its estimated Monte Carlo standard error.
-
Compute the sample mean
-
For each starting index , compute the overlapping batch mean
-
Return the overlapping-batch-means variance estimate
-
Return the estimated Monte Carlo standard error
All sums are finite, the denominator is positive because , and the variance estimate is nonnegative. The procedure always returns a number for the stated finite input. It does not by itself assert a CLT or consistency; consistency requires the hypotheses of the overlapping-batch-means theorem.
A total-variation rate bounds strong-mixing coefficients
Statement
Assume AC. Let be a stationary time-homogeneous Markov chain with invariant probability . Suppose there are a measurable with and such that Then the strong-mixing coefficient of the stationary process satisfies
Facts & Assumptions
Given: AC, a stationary homogeneous Markov chain with kernel and invariant law , an integrable envelope , and the displayed geometric total-variation bound.
AC is the choice principle for the conditional classes used by the Markov-chain interfaces. (The Axiom of Choice)
The strong-mixing coefficient is the supremum of over past and separated future events. (Strong, uniform, and maximal-correlation mixing coefficients)
For every bounded product-measurable future-path functional , the conditional expectation of its value on is a measurable function of . (Markov property for bounded future path functionals)
For bounded measurable , . (Chapman-Kolmogorov equations)
Total variation is . (Total variation distance for probability laws)
A known measurable factor can be taken out of conditional expectation. (Taking out what is known)
Proof
Fix , , , and a future event . By [F2], for the path indicator defining there is a measurable with . Applying [F3] and the tower property gives . Stationarity gives .
For probability measures and measurable , the dyadic simple functions differ from by at most . Their integrals differ by at most , since each summand is an indicator and [F4] bounds its measure difference. Letting proves .
Apply step 1.2 with and . The given rate then yields for every .
Using [F5], [F2], and [F3], The first inequality uses ; the last follows from the preceding pointwise bound and stationarity of . Since were arbitrary, taking the supremum in [F1] proves the claim.
Source notes
The proof is local: the Markov property converts separated future events to a bounded function of the state, and the total-variation bound controls its expectation. The source's factor convention is matched to the repository definition of total variation as a supremum over events.
A measurable observation does not increase stationary mixing coefficients
Statement
Let be a strictly stationary process on a measurable state space and let be a measurable real-valued function. Put . Then is strictly stationary and, for every , where the coefficients are as in Strong, uniform, and maximal-correlation mixing coefficients.
Facts & Assumptions
Given: A strictly stationary process and a measurable real-valued map .
Strict stationarity means every finite-dimensional distribution is invariant under a common nonnegative time shift. (Strong, uniform, and maximal-correlation mixing coefficients)
The coefficients are suprema over past and separated-future events, or over square-integrable variables measurable with respect to those sigma-algebras. (Strong, uniform, and maximal-correlation mixing coefficients)
Proof
For any finite indices and time shift , strict stationarity [F1] makes equal in distribution to . Applying the measurable map coordinatewise preserves equality in distribution, so has the same law as . Thus is strictly stationary.
For every , measurability gives and . Therefore every event pair used in the suprema for and is also allowed for , with the same probability difference. Taking the suprema in [F2] proves and .
Every square-integrable past variable and future variable measurable with respect to the corresponding -sigma-algebras is also measurable with respect to the corresponding -sigma-algebras. Its correlation is unchanged, so the supremum for is over a subset of the pairs allowed for . Hence .
Stationary Markov-chain total-variation convergence implies strong mixing
Statement
Assume AC. Let have a supplied countable algebra generating , and let be a stationary time-homogeneous Markov chain with invariant probability . If then its stationary output process is strongly mixing:
Facts & Assumptions
Given: AC, the countable generating algebra , a stationary homogeneous -chain with invariant law , and the stated -almost-everywhere total-variation convergence.
AC is the choice principle for the conditional classes used by the Markov-chain interfaces. (The Axiom of Choice)
The strong-mixing coefficient is the supremum of over past and separated future events. (Strong, uniform, and maximal-correlation mixing coefficients)
The bounded-future-path Markov property gives a measurable function of the present state for a bounded future-path functional. (Markov property for bounded future path functionals)
For bounded measurable , . (Chapman-Kolmogorov equations)
Total variation is . (Total variation distance for probability laws)
A known measurable factor can be taken out of conditional expectation. (Taking out what is known)
Increasing measurable sets have increasing measures converging to the measure of their union, and symmetric differences satisfy measure subadditivity. (Continuity from below for measures, Measure of a set difference when the smaller set has finite measure)
Dominated convergence applies to measurable functions bounded by an integrable majorant. (Dominated convergence)
Kernel integration makes measurable for bounded measurable , in particular for measurable . (Measure kernel and probability kernel, Iterated transition kernels, Measurability of integration against a kernel)
Proof
Fix , , , and . By [F2], the indicator of this future-path event has conditional expectation for a measurable . Applying [F3] and the tower property gives . Stationarity gives .
For probability measures and measurable , the dyadic simple approximation differs from by at most . Each indicator integral difference is bounded by [F4], so the integrals of differ by at most . Letting gives .
For fixed and , is a finite measure. Let be the sets approximable in -symmetric-difference measure by members of . It contains and is closed under complements because is an algebra. For and , [F6] gives with . Choose with for . The finite union belongs to and by finite subadditivity. Thus is a sigma-algebra containing the generator , so . Approximating each measurable in this way shows : the discrepancy for differs from that for an approximating by at most . Since is countable and each evaluation is measurable by [F8], is measurable; also .
Using [F5] and step 1.1, . The last inequality follows from step 1.2 applied pointwise to and , and then from step 1.3. By hypothesis for -almost every ; since , [F7] makes the integral tend to zero. This bound is independent of , so taking the supremum in [F1] proves .
Source notes
The proof uses the supplied countable algebra to make the total-variation supremum a measurable function of the starting state. The bound on the strong-mixing coefficient then follows by dominated convergence.
Uniform ergodicity bounds the uniform-mixing coefficient
Statement
Assume AC. Let be a stationary time-homogeneous Markov chain with invariant probability . If for constants and then the uniform-mixing coefficient satisfies
Facts & Assumptions
Given: AC, a stationary homogeneous Markov chain with kernel and invariant probability , and the displayed uniform rate.
AC is the choice principle for the conditional classes used by the Markov-chain interfaces. (The Axiom of Choice)
The uniform-mixing coefficient is the supremum of over past events with positive probability and separated future events . (Strong, uniform, and maximal-correlation mixing coefficients)
For every bounded product-measurable future-path functional , the conditional expectation of its value on is a measurable function of . (Markov property for bounded future path functionals)
For bounded measurable , . (Chapman-Kolmogorov equations)
Total variation is . (Total variation distance for probability laws)
A known measurable factor can be taken out of conditional expectation. (Taking out what is known)
Proof
Fix , , with , and . By [F2], there is measurable with . By [F3] and [F5], . Stationarity gives .
For probability measures and measurable , the dyadic simple functions differ from by at most . Their integrals differ by at most , since each summand is an indicator and [F4] bounds its measure difference. Letting gives . Applying this with and proves for every .
Averaging the pointwise bound of step 1.2 conditionally on in step 1.1 gives . Taking the supremum in [F1] proves the result.
Source notes
Uniform total-variation convergence controls every conditional future event uniformly over the current state. No moment assumption on the observable is used in this mixing-coefficient bound.
Zero covariance long-run variance gives a degenerate CLT under absolute summability
Statement
Let be a centered, strictly stationary, square-integrable real process (The space as the quotient by null functions), and put . If and then in and hence in distribution to .
Facts & Assumptions
Given: A centered strictly stationary real process with finite second moment and the stated covariance hypotheses.
Variance is the centered second moment and covariance is the centered product expectation. (Moments, variance, and covariance on a probability space)
Strict stationarity makes every covariance depend only on the lag. (Strong, uniform, and maximal-correlation mixing coefficients)
Convergence in implies convergence in probability. ( convergence implies convergence in probability)
Convergence in probability to a constant implies convergence in distribution to its point-mass law. (Convergence in probability implies convergence in distribution)
Dominated convergence applies to a series indexed by counting measure when its terms have a summable majorant. (Dominated convergence)
Proof
Expanding the finite square of and using centering and [F2] gives For each fixed , the weight tends to and its absolute value is at most . The summable sequence dominates the summands, so [F5] gives .
Since , step 1.1 gives . This is convergence to zero in ; [F3] gives convergence in probability, and [F4] gives convergence in distribution to the point mass at zero, namely .
Poisson-equation martingale decomposition and MCMC CLT
Statement
Assume AC. Let be a stationary time-homogeneous Markov chain with invariant probability , natural filtration , and ergodic canonical path shift. Let be real Borel with and (The space as the quotient by null functions). Suppose there is a measurable such that -almost everywhere, where . Put and . Then for every , and including the degenerate case . If additionally , then so this CLT's variance equals the covariance long-run variance of Long-run variance for a stationary Markov chain.
Facts & Assumptions
Given: AC, a stationary Markov chain with invariant probability , natural filtration , ergodic canonical shift, centered real , and a solution of -almost everywhere.
AC is the choice principle for the conditional-expectation classes used below. (The Axiom of Choice)
For , is square-integrable, has conditional mean zero, and conditional variance . (Kernel conditional expectations for square-integrable Markov observables)
For a stationary ergodic process , its averages converge almost surely to . (Birkhoff limit for a stationary integrable process)
A square-integrable martingale-difference array with terminal conditional variance clock converging to in probability and vanishing unconditional Lindeberg sums has terminal sum converging in distribution to . (Martingale central limit theorem)
For square-integrable real random variables, variance is the second centered moment and covariance is the centered product expectation. (Moments, variance, and covariance on a probability space)
When the covariance series is absolutely convergent and the CLT is established, the covariance long-run variance is the sum . (Long-run variance for a stationary Markov chain)
A known measurable factor can be taken out of conditional expectation when products are integrable; Hölder bounds products of variables by their norms. (Taking out what is known, Holder's inequality for integrals, including the endpoint cases)
If and in probability, then . (Converging together lemma)
Dominated convergence passes limits through the integrals used below. (Dominated convergence)
Convergence to zero in implies convergence to zero in probability. ( convergence implies convergence in probability)
The canonical path law is shifted by the coordinate left shift; ergodic means the shift preserves the law and every invariant event has probability zero or one. (Canonical Markov chain on path space, Measure-preserving transformations and systems, Ergodicity relative to an invariant measure)
A continuous scaling map sends weak convergence to the correspondingly scaled limiting law. (Continuous mapping theorem)
Proof
By [F1], each is a square-integrable martingale difference, so . Since -almost everywhere and each has law , the identity holds at every simultaneously outside a countable union of null events. Therefore almost surely.
Set on the -full set where these terms are finite, and set elsewhere. By [F1], is the conditional variance of , is nonnegative, and . Thus . The process is a measurable factor of the stationary ergodic path process : the inverse image of any shift-invariant event for is shift-invariant for . Hence it too is stationary and ergodic by [F10], and [F2] gives almost surely.
Stationarity and give . The -norm of the endpoint divided by therefore tends to zero, so [F9] gives convergence to zero in probability.
Summing step 1.1 for cancels the intermediate terms and yields , where .
Suppose . For row , set for , set for , and use for every . These are square-integrable martingale-difference rows by step 1.1. Their terminal variance clock is almost surely by step 1.2. For every , stationarity of the transition pairs gives by [F8], since is integrable. Extending each row by zeros makes it an infinite array with the same terminal sum and clock. Therefore [F3] gives ; scaling by using [F11] yields .
If , stationarity and [F1] imply , so almost surely and hence . By step 2.1, . Step 1.3 shows in probability; [F7] then gives convergence to the point mass .
If , step 2.2 gives . Combining this with step 1.3, the decomposition in step 2.1, and [F7] gives .
Assume additionally that the covariance series is absolutely convergent. Stationarity and expansion of the finite square give , where ; [F8] applies to the series because is summable. On the other hand, step 2.1 writes with . For , is -measurable and , so [F1, F6] give . Stationarity gives . Also and [F6] gives . Dividing by therefore shows that its limit is . Steps 3.1 and 3.2 give the CLT, so [F5] identifies its variance with the displayed covariance sum.
Source notes
The decomposition is proved directly by telescoping. The CLT uses the conditional-variance ergodic theorem and the stated martingale-array CLT; finite variance of alone is not used as a CLT criterion. The covariance formula is proved only in the stated Poisson-equation case.
Geometric small-set renewal convergence gives an integrable exponential hitting moment
Statement
Assume the Axiom of Choice. Let be a probability kernel on any measurable space , with invariant probability . Let satisfy , and suppose where is a probability measure. Suppose some and satisfy For a -chain started at , put . Then some makes where . In particular is finite -almost everywhere; finiteness at every state is not asserted.
For the same hypotheses, the split construction also has the following return and occupation properties. Put and . Define , on , and off . On use the kernel and invariant law . Its atom has mass . If has law and , then some satisfies and for every nonnegative measurable on the split space, In particular almost surely and . These additional assertions use exactly the construction and hypotheses above; they do not assert an exponential entrance-time moment from every state.
Facts & Assumptions
Given: AC, the invariant law, small-set minorization, and scalar geometric-convergence bound in the Statement.
The Axiom of Choice: AC is assumed, as required by the path-space and Markov suppliers below.
Small-set minorization and geometric drift: the displayed minorization permits the probability residual on .
Canonical Markov chain on path space: probability kernels on arbitrary measurable spaces have canonical path laws under AC.
Measurability of integration against a kernel: integration of a nonnegative measurable function against a probability kernel is measurable.
Discrete strong Markov property: future bounded path functionals obey the Markov identity at stopping times.
Markov property for bounded future path functionals: bounded future-path expectations are measurable functions of the starting state and satisfy the deterministic-time Markov identity.
Invariant initial law makes a Markov chain stationary: an invariant initial distribution makes every finite-dimensional law stationary.
Monotone convergence for the integral: nonnegative increasing limits pass through integrals.
Tonelli's theorem for nonnegative measurable functions on a sigma-finite product: nonnegative double sums and probability integrals may be interchanged.
The sum of a complex power series is analytic throughout its open disc of convergence and Every complex analytic function is holomorphic: an absolutely convergent complex power series is holomorphic in its open disk of convergence.
A complex power-series sum has complex derivatives of every order, obtained by repeated termwise differentiation: power series can be differentiated repeatedly term by term inside their convergence disk.
All higher complex derivatives exist and satisfy Cauchy's integral formula on an interior circle: a holomorphic function has derivatives of every order locally; these derivatives are continuous because each has a further complex derivative.
Cauchy's inequalities bound the Taylor coefficients by the circle supremum: if on a radius- circle about contained with its disk in the holomorphic domain, then .
Invariant and stationary distribution for a Markov kernel: invariance is equality of the initial law with its one-step image.
Measure kernel and probability kernel: a probability kernel is measurable in its starting point and countably additive in its target event.
Proof
Given: The hypotheses in the Statement, with AC [A1].
Set and . On define ; define for , and otherwise. The split kernel is Nonnegativity follows from [F1], measurability from [F3], and finite summation over the two flag sections defines a probability on every product-measurable event. Thus [F2] realizes the split path. If the initial flag at has probabilities , summing it gives . The same summation at each subsequent time shows that the -coordinate path has the original -chain law. The set is an atom: all its states have the identical outgoing transition obtained by using .
Put . Summing its incoming transitions first over gives and then assigns the new flag according to , proving by [F13]. Let , and let , an ordinary normalized restriction. Write ; its future after time zero is the same from every atom state. Starting in , the next coordinate law is and subsequent coordinate transitions are , so The hypothesis implies for all after increasing a finite .
Let and for , with . Decompose a visit at time by the first return time and use [F4] and the common atom transition to obtain . The possible event contributes zero. For , absolute convergence and [F8] give Here , so the denominator is legitimate. Since decays geometrically, [F9] gives with holomorphic on . For real , ; hence and [F7] gives .
Let and start the split chain from . By [F6], have mean and . Expanding the square of and summing the geometric covariance bound from step 2.1 gives expectation at most for some finite . On this square equals , so . In particular the stationary split chain hits the atom almost surely; no unproved Harris or skeleton-recurrence theorem is needed.
Near the expression is holomorphic and agrees with on its overlap with . Indeed the denominator has value at and remains nonzero on a sufficiently small disk; the quotient is complex differentiable there by the difference-quotient formula for the derivative of a reciprocal. By [F10], for and , , where . Letting and using [F7] and derivative continuity from [F11] gives . Choose a radius whose closed disk about lies strictly within this holomorphic neighborhood and a finite bound for on its circle. [F12] then bounds these nonnegative derivative values by .
For bounded nonnegative , decompose by the last atom visit at times . A last visit at time contributes , by stationarity and [F5]. The remaining paths have no atom visit through time , and their contribution lies between zero and . Thus Let using step 3.2, then extend from bounded to all nonnegative measurable using [F7] and [F8]. This proves the stationary occupation identity . At the atom state has law ; subsequent transitions are independent of which atom state was chosen, so this argument does not require to be constant on .
For , the finite binomial expansion for each integer , followed by [F8], yields Thus some has . The return is almost surely finite by step 3.1. This argument establishes an exponential moment from positivity and analytic continuation; it does not assume in advance that the original series converges beyond the unit disk.
Fix and let , allowing infinity. The bounded path functions increase to , so [F5] and [F7] give measurability and the Markov identity for this nonnegative function by truncation. On , . For , on , the residual time to is ; the Markov identity and step 4.2 give The last moment is finite because and step 5.1 bounds . For the coordinate chain initialized with the prescribed flag, its first visit is no later than its first atom visit . Consequently is measurable by the same bounded-path approximation and . This proves the conclusion, including the null-set qualification.
Strong-mixing covariance bound in terms of tail quantiles
Statement
Let be sub-sigma-algebras of a probability space and put For any real random variable , define its absolute tail quantile by This is a nonnegative nonincreasing Borel function and, for every nonnegative Borel function , If are square-integrable real random variables measurable respectively with respect to , then In particular, if and almost surely, the right-hand side is at most . For complex with these bounds and the same measurability, For two collections of observations separated by time steps in a process , the coefficient of their sigma-algebras is at most as defined in Strong, uniform, and maximal-correlation mixing coefficients. An integral whose upper endpoint is zero is zero.
Facts & Assumptions
Given: The sigma-algebras, quantiles, moment and measurability hypotheses in the Statement.
The strong-mixing coefficient is the supremum of the event-probability differences for the specified separated sigma-algebras. (Strong, uniform, and maximal-correlation mixing coefficients)
Covariance is for square-integrable real variables. (Moments, variance, and covariance on a probability space)
Cauchy–Schwarz bounds the expectation of the absolute product of two square-integrable real variables. (Cauchy-Schwarz for random variables)
A nonnegative random variable is the integral of its threshold indicators, and its expectation is its integrated tail probability. (Layer-cake formulas for random variables)
Nonnegative product-measurable functions on sigma-finite product spaces admit either order of integration. (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product)
Finite measures agreeing on a generating pi-system and on the whole space are equal. (Finite measures agreeing on a generating pi-system and on the whole space are equal)
Proof
For a finite real , write for . The tail tends to zero at infinity, is nonincreasing and is right-continuous. Thus is finite for , is nonincreasing, and is Borel measurable. More precisely, if and only if : if , then belongs to the set defining the infimum; if , right-continuity supplies an with , so that infimum exceeds . Consequently Lebesgue measure on gives for every . The laws of under and of therefore agree on all closed left rays and both have total mass one. These rays form a pi-system and generate the Borel sigma-algebra: their differences give half-open intervals, and countable unions of those give all open intervals. By [F6] the laws agree, proving the stated integral identity. In particular .
Suppose first and . Their first moments and are finite by [F3]. The threshold-indicator identities and [F5] give and . Both integrands have finite integrals, so their difference represents the covariance and the integral triangle inequality applies. The absolute threshold-indicator covariance is bounded by , by , and by , since the joint probability and the product probability both lie between zero and either marginal probability.
Put and . By [F5] and step 1.1, Indeed the minimum is , and each inner threshold set has length its quantile. The right-hand side is finite by Cauchy–Schwarz on and the squared-quantile identity in step 1.1; this particular integral Cauchy–Schwarz inequality follows by integrating and minimizing the resulting quadratic, with the zero-norm case immediate. Step 1.2 therefore proves the covariance bound with constant one for nonnegative .
For real , expand into four terms. Each nonnegative part has a tail quantile at most the corresponding absolute quantile, because its tail probabilities are smaller. The result of step 2.1 bounds each term by ; the triangle inequality gives the asserted factor four. If and , their absolute quantiles are bounded by , so the bound becomes , including or either zero bound.
For complex , expand into the four real bilinear covariance terms of their real and imaginary parts. Each real part is bounded in absolute value by or , so step 3.1 and the complex triangle inequality give . Finally, for observations separated by , both collection sigma-algebras are sub-sigma-algebras of the past and future fields appearing in [F1]; their event supremum is therefore at most . All assertions follow.
Two-sided stationary extension preserves maximal-correlation bounds
Statement
Assume AC. Let be a strictly stationary real sequence, with one-sided coefficients from Strong, uniform, and maximal-correlation mixing coefficients. There is a shift-invariant probability on the coordinate sigma-algebra of such that its coordinate process restricted to has the law of . The invertible left shift preserves . Put and . For every integer , every , and real square-integrable measurable in and measurable in , If , then is ergodic: every event with has .
Facts & Assumptions
Given: AC and the strictly stationary real sequence in the statement.
AC is the inherited choice assumption. (The Axiom of Choice)
Consistent finite-dimensional probability laws on standard-Borel coordinate spaces have a unique probability extension on the coordinate sigma-algebra, assuming AC. (Assuming the Axiom of Choice, Kolmogorov extension for arbitrary families of standard Borel coordinate spaces)
Strict stationarity and use the finite nonnegative past and the entire future at each split. (Strong, uniform, and maximal-correlation mixing coefficients)
Square-integrable products satisfy . (Cauchy-Schwarz for random variables)
Simple functions are dense in on a probability space. (Simple functions with finite-measure support are dense in for )
Two finite measures agreeing on a generating pi-system agree on its generated sigma-algebra. (Dynkin's pi-lambda theorem)
Ergodicity requires invariant measurable events to have probability zero or one. (Ergodicity relative to an invariant measure)
Proof
For finite , choose an integer with and prescribe the law of to be that of . If is another permissible shift, stationarity compares the smaller shift with the larger one, so the law is independent of the choice. Restrictions are consistent because finitely many coordinate sets can be translated by one common shift. Every coordinate space is the standard Borel space , so [F1] gives . Its nonnegative restriction equals the process law of on cylinders and hence everywhere by [F5]. Translating finite sets does not change their prescribed laws, so by the same uniqueness argument. Both and its inverse are measurable because their coordinate maps are coordinates.
For any subset , let be the algebra of events depending on finitely many coordinates in . Every event in can be approximated in symmetric-difference probability by . Indeed, the approximable events are closed under complements. For a countable union, first replace that union by a finite subunion with arbitrarily small lost probability, using continuity from below, then approximate its finitely many constituent events and apply the union bound. Thus the approximable events form a sigma-algebra containing . By [F4], approximate an variable by simple functions, and approximate the finitely many indicator events of each simple function by events. For an indicator the squared error is its symmetric-difference probability; finite summation gives the required arbitrarily small error for a simple function. Consequently every such variable has finite-coordinate approximants.
First suppose and depend on finitely many coordinates in and , respectively. Translate all coordinates by one integer so that the used past coordinates and the split are nonnegative. Their joint distribution is an original-process distribution by step 1.1. The translated variables are measurable in and , so [F2] proves the asserted covariance inequality when both variances are positive. A variable of variance zero equals its mean almost surely and has zero covariance, so the inequality holds then too. For arbitrary , take step 1.2's finite-coordinate approximants. Their means, second moments, and covariances converge by [F3], so passage to the limit proves the same inequality.
Assume , and let satisfy , with . Given , step 1.2 supplies a finite-coordinate event with ; enlarge its coordinate support to for some . For , is in and is in . Step 2.1 applied to their indicators gives , since every indicator variance is at most . Invariance modulo null sets and measure preservation imply . The intersection changes in probability by at most upon replacing both copies of by copies of , and . Hence . Let and then to obtain . Therefore , proving [F6] and the stronger mod-null formulation claimed.
Moderate deviation density bounds for Gaussian anchored sums
Statement
Assume AC and the hypotheses of Gaussian anchored bounded sums have uniform cumulant and Fourier bounds, with , fixed , , , and logarithmic moment generating function . Let denote the density of under the real exponential tilt with parameter . Every such density exists, is continuously differentiable and strictly positive, and There are constants , , independent of and the summand laws, with the following property. If and , put and . Then For put and let . Whenever , For in the same region, These assertions apply separately to every row isometry image in the supplier, including every required coordinate prefix and normalized joint child sum and difference. Their tilt parameter and covariance are taken in that image's own coordinates.
Facts & Assumptions
Given: The hypotheses of the Statement. The inverse square root of a positive definite matrix is its symmetric positive inverse square root, and directional derivatives at a point always keep the chosen exponential tilt fixed.
The bounded and Gaussian sum has an analytic cumulant with on . Under every real tilt its Gaussian component remains independent with covariance , its characteristic function has modulus at most , and its Fourier tails, with or without one linear frequency factor, decay faster than any fixed power of . The same holds for every row isometry image. Gaussian anchored bounded sums have uniform cumulant and Fourier bounds.
Positive definite Gaussian densities have the displayed quadratic form and total mass one; their characteristic functions are the matching Gaussian exponentials. A positive-definite multivariate normal law has its Gaussian density, Characteristic function of a multivariate normal law.
Fubini, Tonelli and dominated convergence license absolutely integrable interchanges and dominated differentiation. Fubini's theorem for L^1 functions on a sigma-finite product, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Dominated convergence.
Real symmetric matrices have orthonormal eigenbases; linear changes of Lebesgue variables multiply by the absolute determinant. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis, A C^1 diffeomorphism satisfies the change-of-variables formula for nonnegative Lebesgue measurable functions.
Cauchy–Schwarz and integration of derivatives along line segments give the norm and Lipschitz estimates below. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation, If is differentiable with integrable then ; and a bounded derivative makes Lipschitz.
A Euclidean Cauchy sequence converges. and for with the Euclidean metric are complete, componentwise from the Cauchy criterion in .
AC is assumed and inherited through the Gaussian and integration suppliers. The Axiom of Choice.
Proof
Under any real tilt, write , where the two terms are independent and is Gaussian of covariance by [A1, F1]. Tonelli and translation of the Gaussian density show that is a probability density. It is strictly positive because its integrand is strictly positive, and is bounded by . The Gaussian density and each of its first derivatives are bounded and continuous, so dominated convergence gives continuous first derivatives of this convolution. This proves the global assertions without a general vector inversion theorem.
Exponential moment differentiation identifies and with the tilted mean and covariance, respectively. To justify the differentiation, the bounded summands permit every polynomial factor, and the Gaussian square completion in [F1] makes all its polynomial moments integrable even after a slightly stronger real exponential; dominated convergence applies. For real , polarization of the symmetric complex valued form gives : apply after replacing by , and minimize the resulting bound over . Integrating on the real segment from zero to therefore gives . Hence for .
We will also need its Fourier formula. Integrating the density of against the character and using its characteristic function in [F2] gives . The constants follow directly by inserting that density and multiplying by . Insert this formula in step 1.1 and use absolute Fubini, dominated by the integrable Gaussian . Thus , and its derivative in direction is the same integral with factor . This differentiation is dominated because is integrable by [F1]. Centering or making an invertible linear change of coordinates gives the identical formulas for those transformed densities.
Fix this and center its law. Its cumulant is . On the original analytic ball still contains , since . Its third derivative is bounded by . Whiten by , writing under this tilt and . Then , , , and, with , it is analytic on and satisfies . The centered Gaussian component of has covariance at least . Its covariance remains unchanged under every further real tilt; consequently all such tilted characteristic functions have modulus at most .
Polarization as in step 1.2 and integration along , , give , , and for real in the ball. For real the same last estimate along the complex segment gives with : this uses only the real direction in each derivative and the factors from the segment parameter.
Put , . If , the map sends the closed ball into itself by step 3.1, and its Lipschitz constant there is at most . Starting at zero, its iterates are Cauchy since consecutive distances decrease geometrically; their limit exists by [F6], lies in the closed ball, and is a fixed point by continuity. Call it . It satisfies , , and . For take . Thus no unproved global inverse assertion is needed for the saddlepoint.
Tilt further by the fixed parameter . Its mean is and covariance . Put ; step 3.1 gives . Its cumulant centered at has third derivative bounded by a fixed multiple of on a complex ball of radius , because lies in the interior of the original analytic ball by step 4.1. Its characteristic modulus on all real frequencies is bounded by the same Gaussian from step 2.2. In what follows reduce the final constant so that and all these interior margins hold.
Let be that last centered characteristic function and let . For , with a sufficiently small fixed , the integrated third derivative estimate gives , , and . The elementary identity therefore gives there. Its integral, and its integral after multiplication by , are bounded by and , respectively. To check the integrability explicitly, is bounded for each fixed , by differentiation in , and the remaining Gaussian is integrable by [F2]. Constants are uniform over and in the displayed neighborhood of .
Outside the ball of step 6.1, both and the Gaussian comparator have unweighted integral and weighted integral . Indeed their global bounds are ; extract and integrate the remaining Gaussian, also with . The extracted exponential is smaller than any prescribed fixed power of because is bounded. This is the same direct Gaussian tail argument as [F1], and uses no smoothness condition on a lattice law. Combining this with step 6.1 and the Fourier formula of step 2.1 gives, for the density of the further tilted , The Gaussian integral in the derivative bound is exact: transform by and rotate onto the first coordinate; the standard Gaussian first absolute moment is , from .
The eigenvalue bounds in step 5.1 give because for , by integrating . Thus is bounded above and away from zero, uniformly, and . Step 7.1 and imply and We choose small enough that all relative errors used here are at most .
The exact exponential tilt density identity is , where is the original density of . It initially follows almost everywhere from the definition of the tilt; both sides are continuous by step 1.1, so it holds everywhere, since a continuous nonzero difference would be nonzero on a ball of positive Lebesgue measure. Step 3.1 and give Indeed , whose error is because is small. Combining this identity with step 8.1 proves . The constant Jacobian for whitening cancels in the density ratio, and gives the stated moderate deviation bound.
Differentiate the exact identity in step 9.1 in direction at , keeping fixed. It gives . By step 7.1, step 8.1, and , reduce further so that ; this is possible because its limiting coefficient is , the error is , and tends uniformly to as . Also step 4.1 gives after the same reduction. Hence for . Whitening identifies this derivative with the derivative in direction in the Statement.
There are only finitely many smallness requirements above: bounds on , , , and fixed multiples depending only on of . Since , one sufficiently small positive enforces all of them whenever and . All constants are uniform over dimensions one through four. Apply the same argument to a row isometry image using [F1]; its Gaussian component retains the same covariance lower bound, and every step applies in its own coordinates. This proves the remaining assertions.
A geometrically ergodic chain has an all-start standard-Borel observation factor
Statement
Assume AC. Let be a probability kernel on an arbitrary measurable space , with invariant probability , finite measurable and satisfying Let be finite measurable. For each specified initial law , there are a nonempty standard-Borel space , a probability kernel on , an invariant probability , an initial probability , and finite measurable functions and such that The complete observation path of a -started -chain has the same law as of a -started -chain. The corresponding stationary observation paths also have the same law, and for every , with infinity allowed. The factor may depend on .
Moreover, any joint law of that observation path with a random element of a standard-Borel space can be lifted to an extension preserving the entire original chain probability space and process.
Facts & Assumptions
Given: AC and the objects and geometric bound in the Statement.
The Axiom of Choice: AC supplies the inherited coding, factorization, disintegration and extension conventions.
A kernel has a countably generated stable sub-sigma-field containing a given observable: a real observable lies in a countably generated kernel-stable sub-sigma-field, with restricted iterates and invariant law.
Standard borel spaces admit bimeasurable real codings, Finite products of standard Borel spaces are standard Borel, and The real line with its Borel sigma-algebra is standard Borel: real finite tuples and standard-Borel spaces have Borel real codings.
Doob-Dynkin factorization through the sigma-algebra generated by a function: a real or extended-real function measurable for the sigma-field of a real code factors through a Borel function of that code.
Disintegration of a joint law on standard borel spaces: standard-Borel joint laws have measurable conditional probability kernels, unique outside a single null set.
Dynkin's pi-lambda theorem: probability measures agreeing on a generating algebra agree everywhere.
Markov property for bounded future path functionals: the original Markov conditional identities hold, and bounded future-observation expectations are measurable.
Ionescu-Tulcea construction of a Markov chain: supplied measurable probability kernels extend an initial probability space while preserving its marginal.
Monotone convergence for the integral and Measurability of integration against a kernel: nonnegative limits and probability-kernel integrals have their usual measurability and integral identities.
The usual real line is Polish, as proved in The real line with its Borel sigma-algebra is standard Borel. Closed subsets of a metric space are , so Under the Axiom of Countable Choice, countable products and subspaces of Polish spaces are Polish makes the closed Cantor space and the countable product Polish under AC. Their Borel measurable spaces are standard Borel by Standard Borel spaces.
Invariant and stationary distribution for a Markov kernel and Iterated transition kernels: invariance is , and iterates are successive kernel integrations.
Proof
Given: The hypotheses in the Statement and AC [A1].
Define This is a probability by [F8], , and : if , then and for every , which also makes . Thus for every . Apply [F2] to code the pair into one real observable, and apply [F1] to it. The resulting countably generated is kernel-stable and contains both and , since the coding has a Borel inverse.
Enumerate generators of and set , taking values in the closed Cantor subset of . Its coordinate digits recover every , so . The target is standard Borel by [F9]; no Borel-image assertion about is needed. By [F3], write and for Borel extended-real functions on . Put , and .
Push the probability to and disintegrate by [F4], obtaining a Borel kernel . For each in a fixed countable generating algebra of , kernel stability and [F3] factor through . Its value is the conditional mean of given under this reference joint law. Hence for -almost every . Intersect over the countable algebra and apply [F5] to obtain equality of these probability measures simultaneously on every Borel , on one -full set. In particular , and because and .
Iterating that kernel equality gives for -almost every , simultaneously for all . At the induction step, avoids the preceding -null exceptional set for -almost every , since ; integration of a nonnegative Borel function of agrees under the two one-step measures by simple approximation and [F8]. The geometric hypothesis therefore implies for -almost every , simultaneously for .
On the total-variation supremum can be taken over its fixed countable generating algebra. Indeed, for any two probabilities, every Borel set is approximable in symmetric-difference measure by that algebra under their finite sum: the approximable sets form a sigma-field, using complements and first approximating a countable union by a finite union. Its probability difference is approximated too. Thus the supremum over the algebra equals the full supremum and is measurable as a function of the starting point. Let be the Borel set where all the bounds in step 4.1 hold and are finite with . It is -full by step 2.1 and the finiteness of . Define and . Absolute continuity makes every full. If , continuity from above gives . Thus is nonempty, Borel, absorbing and -full; no Polish topology on is asserted.
Keep the factor space , and define for and otherwise. This is a Borel probability kernel. Set and ; both assign mass one to , since they are absolutely continuous with respect to . Define and on , and set both to zero on ; these functions are finite measurable. Since is absorbing and , we have . For , the iterates of equal those of , so the geometric bound holds by step 5.1. For , for every , so it holds with . Thus is the required kernel on the nonempty standard-Borel space .
For an original chain with initial law , every has law . Thus step 3.1 and [F6], followed by conditioning on the coded past, show that is a -chain. Since every coded state belongs to almost surely, this is also a -chain. Countably many time and algebra tests suffice; [F5] and bounded approximation extend them to every bounded Borel test. The code chain remains in almost surely and simultaneously for every . The same argument applies with initial law . This proves the complete initial and stationary observation-path equalities and, by pushforward integration, every stated stationary moment identity.
The real observation-path space is standard Borel: its product topology is Polish by [F9], and its Borel sigma-algebra is the coordinate product sigma-algebra, since the real line has a countable base and finite-coordinate basic open sets form a countable base for this product. Disintegrate any specified joint law of this path and a standard-Borel auxiliary random element conditional on the path by [F4]. Evaluate that kernel at the actual original observation path and extend the original probability space by its kernel integral as in [F7]; explicitly, is countably additive by [F8] and has marginal . Its observation/auxiliary joint law is exactly the specified law. This proves the lifting assertion without a conditional kernel on the original arbitrary state space.
Geometric ergodicity and reversibility give exponential rho-mixing
Statement
Assume AC. Let be a probability kernel on with invariant probability , and consider its stationary output process . Assume it is -irreducible (Phi-irreducibility for a general-state Markov kernel), aperiodic (Aperiodicity for a general-state Markov kernel), and Harris recurrent relative to its nonzero -finite irreducibility measure (Harris recurrence relative to an irreducibility measure). Suppose it is geometrically ergodic in the sense of Geometric and uniform ergodicity for MCMC CLTs and reversible with respect to in the symmetric-flow sense of General-state reversibility by symmetric flow measure. Then there are and such that its maximal-correlation mixing coefficient satisfies
Facts & Assumptions
Given: AC, the stationary -chain and all the hypotheses in the Statement.
AC covers the conditional-expectation and path-law interfaces used below. (The Axiom of Choice)
Geometric ergodicity supplies a measurable finite and a common with for every and . (Geometric and uniform ergodicity for MCMC CLTs)
Reversibility means that the probability flow measure is invariant under coordinate interchange. (General-state reversibility by symmetric flow measure)
The square-integrable kernel conditional-expectation identity holds with zero filling on the marginal-null set where an integral is not finite. (Kernel conditional expectations for square-integrable Markov observables)
For square-integrable real variables, . (Cauchy-Schwarz inequality for )
Dominated convergence passes almost-everywhere limits through integrals under an integrable dominating function. (Dominated convergence)
A bounded measurable future-path functional has conditional expectation, given the chain's past through its initial time, equal to a bounded measurable function of the state at that time. (Markov property for bounded future path functionals)
For bounded measurable , , where . (Chapman-Kolmogorov equations)
Nested conditional expectations satisfy the tower property, and a known factor can be taken out when the products are integrable. (Tower property of conditional expectation, Taking out what is known)
Conditional Jensen gives . (Conditional jensen inequality)
The coefficient is the supremum of absolute correlations over square-integrable past and separated future variables of positive variance. (Strong, uniform, and maximal-correlation mixing coefficients)
Proof
On real put and . By [F3], conditional Jensen and stationarity give and . These identities also show that the operator is well defined on almost-everywhere classes. Thus is a contraction preserving , and . For bounded , [F2] gives by integrating against the flow measure. Bounded truncation converges in by [F5], so contraction and [F4] extend this symmetry to all functions. Iterated kernel integration agrees with ; for bounded functions this is [F7], and contraction extends the identity to classes.
For probability measures and , approximate uniformly by . Each indicator integral difference is bounded by the event-supremum TV distance, so on taking limits. Applying this to for nonzero bounded measurable , and separating , gives .
Put for positive integers . If bounded vanishes off and , symmetry and invariance give for . Step 1.2 therefore gives . Write , so . No integral of over the whole state space has been used.
For and , symmetry and [F4] give . Repeating with yields . For the in step 2.1, if the bound is immediate; otherwise and the last display is at most , tending to . Hence that bound holds for every such .
For choose a finite measurable representative and put . Finiteness of gives almost everywhere and , so [F5] gives in . By [F4], , and therefore in . Applying step 3.1 to these centered functions and taking limits through the contraction gives . Iteration on proves for all . This also covers .
Fix and bounded centered real measurable respectively for and . By [F6] there is bounded measurable such that . Stationarity and [F9] give and . The tower property and [F7] give . Taking out , applying [F4] and stationarity, and then step 4.1, yields .
For general centered square-integrable in the same sigma-algebras, truncate each at levels and subtract its expectation. By [F5] and [F4], these bounded centered variables converge in ; [F4] passes their covariances and norms to the limit in step 5.1. Dividing by the positive standard deviations and taking the suprema in [F10] proves . If no positive-variance variables exist, the coefficient is zero by the nonnegative empty-supremum convention. Thus the stated conclusion holds with .
Source notes
Jones, Theorem 2(3), supplies the stated implication. Roberts and Rosenthal, Geometric ergodicity and hybrid Markov chains, Theorem 2 of the author-hosted preprint (published Theorem 2.1), proves the reversible spectral-gap route; its full proof is on PDF pp. 7–10 at https://www.probability.ca/jeff/ftpdir/hybrid.pdf . The proof here instead uses envelope level sets and repeated Cauchy–Schwarz to obtain the operator bound locally. It needs neither an integrable envelope nor a compact or general spectral theorem. The irreducibility, recurrence and aperiodicity hypotheses are retained from the source statement; the explicit common geometric rate and reversibility suffice for this argument.
A rho-mixing CLT for square-integrable stationary sequences
Statement
Assume AC. Let be a centered strictly stationary real sequence with . If then the covariance series converges absolutely and is nonnegative. If , then for , If , Zero covariance long-run variance gives a degenerate CLT under absolute summability gives the degenerate normal limit.
Facts & Assumptions
Given: AC, the centered strictly stationary real sequence , finite second moment, and .
AC is the inherited assumption of the local extension and martingale CLT suppliers. (The Axiom of Choice)
The coefficient bounds absolute correlations between the finite past and separated future. (Strong, uniform, and maximal-correlation mixing coefficients)
The sequence has a two-sided stationary realization with the same nonnegative-index law and the same maximal-correlation upper bounds for the full past; if , its invertible shift is ergodic even for mod-null invariant events. (Two-sided stationary extension preserves maximal-correlation bounds)
Conditional expectation is the orthogonal projection onto the conditioning subspace and preserves ordinary expectation. (Conditional expectation is the l2 orthogonal projection, Basic algebra and order properties of conditional expectation)
For a centered process generated by an invertible probability-preserving ergodic transformation, summability of gives a normal CLT with variance , including zero; under absolute covariance summability this variance equals the covariance series. (Summable conditional projections give a stationary martingale CLT)
Covariance and variance are centered product and square moments. (Moments, variance, and covariance on a probability space)
Proof
Summability of the nonnegative numbers implies . By [F2], realize as the nonnegative restriction in law of the two-sided stationary sequence with an invertible probability-preserving shift that is ergodic for mod-null invariant events. Its full past satisfies for past and future separated by . All second moments of the coordinates agree with those of .
Put for . By [F3], is centered and square-integrable and . If , the bound is immediate. Otherwise apply step 1.1's covariance inequality to and , divide by , and obtain that same bound. Consequently .
Put . Step 1.1 applied to the two coordinates gives for ; if is constant, both sides vanish. Thus . All hypotheses of [F4] now hold by steps 1.1 and 2.1. It follows that and . This includes and proves absolute covariance convergence. Since and have the same nonnegative-index law, the series is exactly and the finite sums have the original laws. Therefore the asserted nondegenerate and degenerate limits both hold for .
Source notes
The proof uses the stronger summable individual-projection criterion, proved locally by an exact martingale-coboundary series. No rho-mixing CLT is imported as a proof premise.
Strong Brownian approximation implies overlapping-batch consistency
Statement
Assume AC. Let be real random variables and , with . Suppose on a probability-space extension with the original process as marginal there are a standard Brownian motion , constants and , and finite random constants such that almost surely Let be deterministic positive integers with eventually and For the exact estimator of Overlapping batch means for MCMC standard errors, applied to , No stationarity, independence of from , or monotonicity of or is required. The conclusion includes .
Facts & Assumptions
Given: AC, the original process and its Brownian approximation on the stated extension, the constants , and the deterministic batch sequence with all displayed numerical hypotheses.
AC is assumed as in The Axiom of Choice and the Brownian and Gaussian interfaces.
Brownian motion starts at zero and has independent normal increments with variance equal to interval length. Brownian motion
A centered Gaussian quadratic form satisfies the dimension-free even-moment bound for every positive integer . Centered Gaussian quadratic forms have dimension-free even-moment bounds
For a nonnegative random variable and , . Markov's inequality for random variables
Events with summable probabilities occur only finitely often almost surely. The first Borel-Cantelli lemma for measures
The standard-normal Mills upper bound , together with symmetry and the normal density, gives for . Two-sided Mills bounds for the standard normal tail
The Euclidean norm is a norm and satisfies . Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation, The reverse triangle inequality in a normed space
The finite OBM algorithm uses the displayed denominator and overlapping batch means. Overlapping batch means for MCMC standard errors
Proof
Work first on the extension. For each sufficiently large write , , for , and for . By [F1], the are independent standard normals. If indicates coordinates , then where . At most windows contain a fixed pair of coordinates, and none does if their separation is at least . Hence , when , , and .
Choose an integer . By [F2] and [F3], for every , . These probabilities are summable since eventually, and . Applying [F4] for each and intersecting the resulting countable family of full-probability events gives almost surely. Since , almost surely. No independence between different is used.
By [F1] and [F5], for . Thus [F4] gives eventually almost surely. Put and . Applying [F6] to the vectors and gives , using . Consequently almost surely.
Set and . On the full event of the given approximation, enlarge its finite random constant to absorb all finitely many early errors. Since is increasing for integers , this gives a finite random with for every sufficiently large . For every window, where . Therefore simultaneously for all windows, and .
Put . By [F6] and , . Hence almost surely. This argument divides neither by nor by , so it also proves the zero-variance case.
By [F7], . This deterministic multiplier tends to one because , and step 5.1 proves the desired limit on the extension. The event that this sequence of finite measurable functions of converges to depends only on the original process. Its probability under the extension equals its probability under the original marginal, so that probability is one for the original process as well. Arbitrarily defining the estimator at any finitely many indices with inadmissible batches does not affect its limit.
Bounded stationary CLT under summable strong mixing
Statement
Assume AC. Let be a centered strictly stationary real sequence, with almost surely for a finite , and suppose Put . Then its covariance series converges absolutely and Moreover including . Writing , every finite subset satisfies The strong-mixing coefficient is that of Strong, uniform, and maximal-correlation mixing coefficients.
Facts & Assumptions
Given: AC, the centered strictly stationary bounded sequence, and summability of its strong-mixing coefficients.
AC is assumed for the characteristic-function convergence criterion and the normal-law supplier. (The Axiom of Choice)
Strict stationarity preserves finite-dimensional laws under nonnegative shifts, and the mixing coefficient is the supremum over separated past and future events. (Strong, uniform, and maximal-correlation mixing coefficients)
Bounded real variables from fields with event coefficient have covariance at most in absolute value; complex variables obey the analogous bound . (Strong-mixing covariance bound in terms of tail quantiles)
Variance and covariance are centered product expectations. (Moments, variance, and covariance on a probability space)
Dominated convergence applies on an arbitrary measure space, including counting measure with a summable majorant. (Dominated convergence)
Cauchy–Schwarz bounds absolute product expectations of square-integrable real variables. (Cauchy-Schwarz for random variables)
For , one has for all real . (Second-order characteristic-function expansion)
Under AC the characteristic function of is , including . (Characteristic function of a normal law)
Under AC, convergence of characteristic functions to that of a specified probability law is equivalent to weak convergence to that law. (Characteristic function criterion for weak convergence)
Proof
Put . By [F1]–[F2], for , so . Expanding the square of a sum indexed by a finite subset gives its diagonal terms and off-diagonal terms for pairs in . For each fixed positive lag there are at most such pairs. Thus the asserted bound holds, also for . For a consecutive block, the same expansion gives . The summable majorant and [F4] imply convergence to the stated , which is nonnegative as a limit of variances.
Extend the coefficient convention by . It is nonincreasing: increasing the separation restricts the future sigma-algebra in [F1]. For ordered indices , write and . If , centering of and [F2] applied to and give ; if , the trivial bound suffices. Splitting off the centered last factor gives the same bound with . Therefore . There are at most orderings, at most choices each for the first index and the middle gap, and pairs with maximum . Hence To justify the last limit, for each fixed , , whereas as ; the unweighted part divided by also tends to zero. Repeated indices are included by the convention at zero.
Set and . By step 1.2, . For sufficiently large , define , interpreting the second entry as if , and ignore the finitely many earlier . Then , , and . Put , , and . Thus , , , and . Retain consecutive blocks of length , separated by discarded blocks of length , and discard the final remainder. Their discarded index set has cardinality at most . Step 1.1 gives , controlling all gaps together without asserting their independence.
Summability and monotonicity imply , since . For large , , so . Let denote the sum in retained block . At a fixed real frequency , put . Both each factor and every partial product have modulus at most one. The partial product through block is measurable before the gap preceding block , so [F2] bounds the discrepancy in factoring its expectation from by . Telescoping these factorizations gives This compares characteristic functions only and requires no independent copies or coupling construction.
Write . These variables are centered and each has the law of by stationarity. Steps 1.1 and 2.1 imply . Also because and . For every , the Lindeberg sum of these block marginals is at most . Applying [F5] first to and then the finite-sum Cauchy–Schwarz inequality yields . Finally . No uniform-integrability hypothesis or deterministic block bound tending to zero has been assumed.
For fixed , [F6] gives , where and . Step 3.2 gives and . The elementary real Taylor bound for follows by integrating the second derivative twice. Since both and are at most one, the product telescoping inequality now gives Thus the product characteristic function tends to .
The discarded normalized sum in step 2.1 tends to zero in and hence in by [F5]. The elementary bound , obtained by integrating the derivative of along the real interval, makes its effect on each characteristic function tend to zero. Combine steps 3.1 and 4.1: for every real . Under [A1], [F7] identifies this with the characteristic function of , and [F8] gives the asserted weak convergence. This remains valid for ; all finite-moment and discarded-set assertions were already proved in step 1.1.
Geometric strong mixing gives an Lp Poisson solution
Statement
Assume AC. Let be a probability kernel on an arbitrary measurable space with invariant probability . Let be its stationary Markov chain, and suppose its strong-mixing coefficients satisfy for , with and . Let be measurable with and for some finite .
For every , put . The kernel operators are understood on -almost-everywhere classes, with finite measurable versions obtained by setting divergent integrals to zero on their -null exceptional sets. Then and the series converges in , with Moreover, writing , one has The covariance sum is finite and nonnegative, and
Facts & Assumptions
Given: AC, the stationary kernel and invariant probability, the geometric strong-mixing bound, centered real , and .
AC is assumed. The Axiom of Choice
For square-integrable variables in sigma-algebras with event mixing coefficient , , and for nonnegative Borel . Separation by bounds by . Strong-mixing covariance bound in terms of tail quantiles, Strong, uniform, and maximal-correlation mixing coefficients
The bounded Markov identity is ; its square-integrable extension identifies the conditional mean with the kernel integral on a full invariant-measure set. Chapman-Kolmogorov equations, Kernel conditional expectations for square-integrable Markov observables
Conditional expectation is an contraction, and on a probability space when . Conditional lp contraction, Finite-measure includes into for
Hölder and its product form bound products in the corresponding spaces; Minkowski bounds the norm of a finite sum by the sum of norms. Holder's inequality for integrals, including the endpoint cases, Generalized Holder inequality puts products into , Minkowski's inequality for integrals, including
Under AC, is complete, consists of almost-everywhere classes of measurable representatives, and norm-convergent sequences have almost-everywhere convergent representative subsequences. Riesz-Fischer completeness of for , The space as the quotient by null functions
Increasing nonnegative measurable functions pass through integration by monotone convergence. Monotone convergence for the integral
Invariance means ; the identity extends from indicators to nonnegative measurable functions by increasing simple approximations. Invariant and stationary distribution for a Markov kernel
Centered second moments are variances, covariance is the centered product expectation, and covariance is bilinear and symmetric. Moments, variance, and covariance on a probability space, Covariance is symmetric and bilinear in finite linear combinations
Proof
Invariance [F7], extended with [F6], gives for integrable , so the absolute kernel integral is finite off a measurable -null set. A change of representative on a -null set does not change its image class, since implies almost everywhere. For , [F3] gives , and [F2] identifies with . Conditional contraction [F3] and stationarity give . The bounded iteration identity in [F2] extends to these classes by replacing with its bounded clipped truncations: their error tends to zero because the integrable tail tends to zero, as follows from [F6] applied to its complementary truncations; conditional and kernel contractions pass this error to zero. Thus is a well-defined linear contraction on , its iterates agree with the kernel iterates, and . This also proves before any covariance estimate is used.
Fix , put , and let and for integers . This function is bounded, belongs to , and satisfies . The centered is in by [F3]. Step 1.1 and bounded conditional multiplication give ; since , this expectation is the covariance of and . Bounded conditional multiplication here follows directly from the defining conditional expectation event identity, first for simple and then by uniform bounded approximation. Thus [F1] is applied only to square-integrable variables, even though .
Put . The quantile identity in [F1] gives and . Since , apply the product inequality [F4] to these two quantile functions in , and then Hölder with , whose norm is . Therefore [F1] and step 2.1 give . If the quantile integral and the right-hand side are zero. Let using [F6]. If the required bound is immediate; otherwise divide by to obtain . The restriction ensures and makes all Hölder exponents just used valid.
The geometric bound makes , including . By step 3.1, [F4] makes a Cauchy sequence in . Completeness [F5] gives . The contraction in step 1.1 gives , while in by step 3.1. Thus as almost-everywhere classes. Choose the finite measurable representative supplied by [F5], setting it to zero on its measurable null exceptional set if necessary, and use the version convention from step 1.1 for . These choices make the asserted almost-everywhere identity precise on the original sigma-algebra without completing the measure.
Apply [F1] directly to , both square-integrable. Their quantiles have common norm . Two applications of [F4], with reciprocal exponents , yield . The exponent is positive because ; the geometric bound therefore gives . At a zero mixing coefficient the quantile integral is zero, so no division by that coefficient is needed.
Centering, stationarity and [F8] give . To pass to its limit, choose any finite ; the terms converge individually to , while the absolute value of either remaining weighted or limiting tail is at most , which tends to zero as by step 4.2. The limit is the claimed finite covariance sum . Each left-hand side is nonnegative, hence . When almost everywhere all bounds and the choice hold directly; no claim of uniqueness of is made. The strict exponent range and all zero-mixing cases were handled in steps 3.1, 4.1, 4.2. AC is used through the Markov/conditional and completeness suppliers, not as an imported covariance or Poisson theorem.
Source notes
The actual Doukhan–Massart–Rio full PDF was retrieved and its printed p. 66 inspected as a rendered page, including the finite-variance restriction of Lemma 1 and the covariance-sum conclusion of Proposition 1. The proof above consumes the fully proved local quantile lemma. Its test variable is bounded before that lemma is applied; the unbounded dual variable, whose exponent is less than two, is never inserted into a result requiring square integrability. The Poisson-series construction, norm convergence and nonnegative variance limit are proved here, not replaced by a source citation.
Conditional quantile couplings for Gaussian anchored sums
Statement
Assume AC and the hypotheses of Gaussian anchored bounded sums have uniform cumulant and Fourier bounds: , , the centered summands are independent, , is Gaussian with covariance , and , where and . Put , with from that supplier, and let be the CDF of . There are constants and , independent of the summand laws, , with the following properties.
For , let be the conditional CDF of given obtained by dividing the joint prefix density by its marginal density. These densities are strictly positive, so this defines a probability kernel at every , including prefixes of probability zero. Its quantile transport is jointly Borel measurable and takes a fresh variable to that conditional law. If , then, with , and also For , the prefix is empty and its norm is zero. The resulting successive coordinate transport sends to the law of and satisfies
If with , write with . The last successive conditional transports give a jointly measurable taking a fresh to the conditional law of given , and when and . The Gaussian reference bridge is independent of its Gaussian parent.
There is also the following scalar unequal-variance form. Let be independent centered scalar sums of independent variables bounded by and independent Gaussian terms, with variances , and with Gaussian component variances at least , respectively. Put Then , , . There is a jointly measurable transport taking a fresh variable to conditional on with when and . The corresponding root transport from to obeys the analogous bound with . The reconstructed children are and .
All these maps are jointly measurable in an additional measurable parameter whenever, on each fixed finite row length, the bounded input laws are probability kernels, their Gaussian covariance and are measurable functions, and the displayed hypotheses hold parameterwise. This includes countable arrays with different row lengths. No claim here asserts the conditional estimates for raw bounded laws without a positive Gaussian component.
Facts & Assumptions
Given: The centered independent Gaussian anchored sum and its parameters as in the Statement; vector norms are Euclidean, and all constants below depend only on .
AC is assumed and inherited by the Gaussian suppliers. (The Axiom of Choice)
The sum has an analytic cumulant with , and every real tilt retains an independent Gaussian component with covariance at least . Each coordinate prefix has the same properties. (Gaussian anchored bounded sums have uniform cumulant and Fourier bounds)
Every tilted density is positive and continuously differentiable, has supremum at most , and on the small tilt and moderate deviation region has its Gaussian density times . Tilted covariance lies between and . The assertions hold in each prefix's coordinates. (Moderate deviation density bounds for Gaussian anchored sums)
Positive definite Gaussian laws have their positive quadratic-form densities, and the standard normal CDF is continuous, strictly increasing, with limits zero and one. (A positive-definite multivariate normal law has its Gaussian density, The standard normal CDF is continuous and strictly increasing)
Nonnegative product integration and dominated convergence apply to the integrals below. (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Dominated convergence)
Integrating a nonnegative product-measurable function against a finite or uniformly sigma-finite kernel gives a measurable function, and finite products of probability kernels are probability kernels. (Measurability of integration against a kernel, Finite products of probability kernels are probability kernels)
Integration of derivatives along real line segments gives the corresponding norm and Lipschitz bounds; Euclidean Cauchy sequences converge. (If is differentiable with integrable then ; and a bounded derivative makes Lipschitz, and for with the Euclidean metric are complete, componentwise from the Cauchy criterion in )
The rationals are countable and dense in the reals. ( is countably infinite, Both and are dense in , and every nonempty open subset of is uncountable)
Proof
Work with a fixed prefix of length and write its dimension as during the proof. Its law satisfies [F1]–[F2]. Set , , and normalize coordinates by . Its cumulant has mean zero, covariance , and third derivatives bounded by . Polarization and [F6], integrating along a real segment, give and . Fix a prefix , put , , and assume sufficiently small. The map sends into itself and has Lipschitz constant at most . Its iterates from zero have geometrically decreasing consecutive distances, hence converge by [F6] to a solution of with . For the empty prefix use . Put and let be the physical last mean under this prefix tilt. The last mixed covariance row at zero vanishes; integrating its Hessian deviation yields . The physical tilt parameter is , so can be kept below .
Let be the joint and prefix densities. The Gaussian convolution representations in [F2] and [F4] give at every : integrate the Gaussian joint density over its last coordinate and then integrate the bounded-summand law. Both densities are everywhere positive. Thus defines a probability density for each and, by integration, a version of the conditional kernel. Denote tilted densities by . Since , their exponential factors cancel exactly, and the density of conditional on is . The centered tilted covariance satisfies by step 1.1. Apply [F2] to its full density at and prefix density at zero. The ratio of the matching Gaussian densities is a centered normal density with variance ; the determinant identity follows by eliminating the off-diagonal blocks in . Consequently for . Here is the density of , and the small constants also ensure all [F2] moderate bounds for the full centered tilt.
In standardized units set and let be the standard normal density. Step 2.1 implies, on , and hence after reducing and . There is also a global bound . Indeed the standardized prefix-tilted law has prefix density at its mean bounded below by a positive constant by [F2]. Tilt its centered full law further by , where . Its full density is bounded above by a constant by the global part of [F2], while integration of its bounded Hessian gives a centered log mgf at most . The density tilt identity therefore gives at every real . Choosing proves the bound. The combined real tilt remains in a fixed smaller analytic ball, where the Hessian bound from step 1.1 applies; the global density bound in [F2] itself permits every real tilt.
We next convert the density estimate into a CDF displacement without a coupling theorem. Write for . Choose a fixed large and then a sufficiently large fixed . Put and , with fixed sufficiently small. Reduce the allowed bounds on and so that for one has and ; this is possible because for and . Also choose within the moderate interval of step 3.1. For , the maps are increasing and their arguments are comparable with . Since , step 3.1 bounds the logarithmic density error at either perturbed argument by , whereas the Gaussian log-density changes by at least a fixed positive multiple of in the favorable direction. The log Jacobians are at most and at least , respectively. Since , first taking large and then large gives and . These comparisons follow by adding the three logarithmic contributions; the Jacobian contribution divided by is at most , and the density-error contribution is at most .
The inverse equations and monotonicity of imply and in the range below. Substituting and in the inequalities of step 4.1 gives, for , Only a smaller range , with , will be retained, so the lower integral has its endpoints in the stated order. The global tail estimate in step 3.1 and the Gaussian tail give : split the first integral at a constant multiple of and integrate the Gaussian and exponential parts. Enlarge the shifts in the displayed inequalities to . Each extra Gaussian interval has mass at least . Choose so small that . For sufficiently small , uniformly for , that interval mass exceeds the discarded tail, since . Thus the full tails have the corresponding inequalities with shift ; here . Reflecting proves the negative-tail version as well.
At step 5.1 bounds the tail-probability difference by . For , step 3.1 bounds the density-ratio error by , where ; integrating over this fixed interval gives a CDF difference at most . The Gaussian density is bounded below by a fixed positive constant on , so a shift of length in either direction dominates that difference, after reducing the allowed bound on to keep the shift within this compact interval. Since , this proves together with step 5.1, for the standardized conditional CDF , All constant choices are finite and made in the indicated order; reducing the common cutoff enforces every smallness requirement. Returning to physical units, the CDF of the centered conditional coordinate satisfies , where .
Every conditional density is positive and continuous by step 2.1, so its CDF is continuous, strictly increasing and has limits zero and one. Its inverse on is therefore unique. For a fresh variable , [F3] implies that is uniform: for , the event is , of probability . Thus the inverse-CDF map has the asserted conditional law. If and , then . Step 6.1 and strict monotonicity imply whenever is moderate. Step 1.1 gives , so implies the required moderate bounds on and . The inequalities and now give .
The bound in Gaussian rather than target values follows by bracketing, not by substituting for in step 7.1. In standardized units fix a Gaussian coordinate and put , , with fixed large . If and are small enough, remain in the moderate region and the standardized version of is at most . Indeed and ; the extra is absorbed into because is small. Step 6.1 therefore gives , so its inverse at lies in . Adding the shift from step 1.1 and restoring units gives under the stated Gaussian-value cutoff.
The maps just constructed are jointly measurable, also under the additional parameter in the Statement. On each fixed row length the product of the bounded input-law kernels exists by [F5]. The Gaussian density is jointly Borel in its positive definite covariance and its argument: its determinant and inverse are continuous functions of the matrix entries by their polynomial and cofactor formulas. Consequently is jointly measurable by [F5], as are the prefix densities, their positive ratios and their CDF integrals. Lebesgue integration here uses its common exhaustion by bounded intervals. The conditional inverse can be written ; by [F7] this equals the continuous strictly increasing inverse and is jointly measurable, since strict sublevel sets are countable unions of measurable sets indexed by . Compose with the continuous . Different finite row lengths are handled on their countable measurable disjoint union. This constructs the kernel versions explicitly at every prefix.
Starting with independent Gaussian coordinates of variance , apply the conditional maps successively to coordinates . Conditional integration of the densities in step 2.1 shows inductively that the first constructed coordinates have density ; their joint density is the product of the successive density ratios, which telescopes to . The full map therefore has the law of . On , every coordinate and prefix meets the target-value cutoff in step 7.1 after reducing by a fixed factor. Summing at most four coordinate bounds gives . For apply only its last conditional maps with the specified prefix and fresh independent Gaussian coordinates; the same telescoping gives the conditional law of and the same coordinate estimates give . Separate bounds on the two norms imply the joint cutoff by reducing . The Gaussian reference has covariance , so its parent and bridge coordinate blocks are independent by the product Gaussian density in [F3]. Joint measurability follows by composing the maps of step 8.2.
For the scalar extension put , and . Independence gives , and , hence . The stated ranges follow from ; the function is increasing in either positive argument, giving . Thus , and . A bounded left increment contributes to , whose norm is at most ; a bounded right increment contributes with the same bound. The joint Gaussian component is obtained from independent child Gaussians by the same linear map. Since their variances are at least , its covariance, by congruence, is at least . Therefore satisfies the original hypotheses with bound and parameter . Apply step 9.1 conditionally on , feeding , and divide the output by . The uniform bounds on convert its estimate and its cutoffs into the asserted scalar bridge bound after changing by fixed factors. Applying the scalar root construction to and rescaling gives the root bound. Finally rearranges to the two stated child formulas. Every coefficient is a continuous function of the positive variances, so the parameter measurability in step 8.2 is preserved. All conclusions follow.
Harris geometric ergodicity admits an integrable total-variation envelope
Statement
Assume the Axiom of Choice. Let be a probability kernel on a measurable space with invariant probability . Assume it is -irreducible (Phi-irreducibility for a general-state Markov kernel), aperiodic (Aperiodicity for a general-state Markov kernel), and Harris recurrent relative to its nonzero -finite irreducibility measure (Harris recurrence relative to an irreducibility measure). Suppose is geometrically ergodic in the sense of Geometric and uniform ergodicity for MCMC CLTs. Then there exist a measurable with and such that
In fact, the same conclusion holds assuming only AC, invariance of , and the everywhere finite measurable common geometric-rate bound in Geometric and uniform ergodicity for MCMC CLTs, on the original arbitrary measurable space. The construction below derives the needed irreducibility on a restricted sigma-field directly from that bound.
Facts & Assumptions
Given: AC and all hypotheses in the Statement. Denote an everywhere finite measurable envelope in the geometric-ergodicity hypothesis by , with rate .
The Axiom of Choice: AC is assumed for the small-set supplier and the canonical path-space and Markov suppliers below.
Geometric and uniform ergodicity for MCMC CLTs: for every and , with finite and measurable.
A kernel has a countably generated stable sub-sigma-field containing a given observable: a real measurable observable belongs to a countably generated sub-sigma-field on which the kernel restricts, with the same iterated transitions on restricted events and the restricted invariant law.
A countably generated phi-irreducible chain has a positive-measure small set: a countably generated irreducible chain has a positive-irreducibility-measure small set , with a minorizing probability satisfying .
Geometric small-set renewal convergence gives an integrable exponential hitting moment: a one-step small set of positive invariant mass, whose minorizing-law return probabilities converge geometrically, has a stationary-integrable exponential entrance-time moment.
Canonical Markov chain on path space: under AC, the original kernel has canonical chains from every initial law.
Discrete strong Markov property: bounded future-path functionals obey the Markov identity at stopping times.
Phi-irreducibility for a general-state Markov kernel: irreducibility relative to a nonzero measure means that each positive-measure measurable set is accessible from every state at some positive time.
Proof
Given: The hypotheses in the Statement and the envelope in [F1], with AC [A1].
Apply [F2] to to obtain a countably generated , the restricted kernel , and invariant probability . If has , [F1] implies for all sufficiently large , for each , because is finite. Thus is -irreducible by [F7]. Applying [F3] on this countably generated space gives an , a with , a probability on with , and such that on . No small-set assertion on the original sigma-field is being made.
Since is finite everywhere, increases to as the positive integer increases. Choose so has both and . Choose a finite integer and put , . Choose . For this gives on , , and on .
For , the restricted-iteration identity in [F2] and [F1] give For the absolute difference is at most . Therefore [F4] applies with constant and rate . It supplies and a measurable , where , such that . Indeed the skeleton path restricted to under the full canonical -chain has the -chain law: its finite-dimensional cylinder probabilities agree by the iterated-kernel identity in [F2]. Thus the same is the entrance-time moment for the original chain of [F5].
Choose , and set and . Fix and . On , decompose into the finitely many events and use [F6] to replace the conditional terminal indicator by . Since , [F1] bounds its difference from by ; when , use . On use the bound . Subtracting on these two events gives For the final inequality, , and on one has ; on , . This argument is for each original measurable event , so taking its supremum controls the full original total variation and does not require measurability of a total-variation supremum.
The measurable set is -null by step 3.1. Define for and for . It is measurable and finite everywhere, and . Off , step 4.1 gives the desired bound with rate ; on , [F1] gives it because . Taking and proves the stronger conclusion from invariance and the everywhere geometric bound alone: the restricted irreducibility used in step 1.1 was derived there directly. The original additional Harris, irreducibility and aperiodicity hypotheses therefore give the same conclusion as a special case. AC was used precisely through [F3] and the path-space and Markov suppliers [F4], [F5] and [F6].
Strong-mixing CLTs under bounded, moment, and logarithmic-moment conditions
Statement
Assume AC. Let be a centered strictly stationary real sequence with , and define . Each of the following is sufficient:
- There is with almost surely and .
- For some , and .
- For some , and , where for and .
In each case the covariance series converges absolutely and is nonnegative, and This includes the degenerate law when , as in Zero covariance long-run variance gives a degenerate CLT under absolute summability. The coefficient and strict stationarity have the meanings in Strong, uniform, and maximal-correlation mixing coefficients.
Facts & Assumptions
Given: AC, the centered strictly stationary square-integrable real sequence, and one of the three sufficient regimes.
AC is assumed for the bounded CLT, normal-law characteristic function and characteristic-function convergence criterion. (The Axiom of Choice)
A centered bounded stationary real sequence with summable strong-mixing coefficients has an absolutely summable covariance series and a CLT with its nonnegative covariance-sum variance, including variance zero. (Bounded stationary CLT under summable strong mixing)
Absolute tail quantiles have the original absolute-value law under uniform Lebesgue measure; square-integrable variables from fields with event coefficient satisfy . (Strong-mixing covariance bound in terms of tail quantiles)
A measurable scalar observation preserves strict stationarity and does not increase the strong-mixing coefficients. (A measurable observation does not increase stationary mixing coefficients)
Dominated convergence applies on arbitrary measure spaces and therefore also to absolutely dominated covariance series. (Dominated convergence)
Tonelli interchanges nonnegative product integrals, including integration against counting measure. (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product)
Holder's integral inequality holds for conjugate finite exponents on every measure space. (Holder's inequality for integrals, including the endpoint cases)
Cauchy–Schwarz bounds absolute product expectations of square-integrable real variables. (Cauchy-Schwarz for random variables)
Under AC the characteristic function of is , including . (Characteristic function of a normal law)
Under AC, weak convergence to a specified probability law is equivalent to pointwise convergence of characteristic functions to that of the law. (Characteristic function criterion for weak convergence)
Zero covariance-sum variance and absolute covariance summability give convergence of normalized sums in and in distribution to . (Zero covariance long-run variance gives a degenerate CLT under absolute summability)
Variance and covariance are centered product expectations. (Moments, variance, and covariance on a probability space)
The coefficient is the event supremum for the past and the future separated by . (Strong, uniform, and maximal-correlation mixing coefficients)
Proof
Write for the absolute tail quantile in [F2], and put and . We first show in every stated regime that and . In regime 1, and [F5] give . This regime also already follows directly from [F1].
In regime 2 put and . Holder's inequality [F6] with exponents and , followed by the quantile identity [F2], gives Tonelli [F5] and the assumed summability yield . Also by its event definition, so , proving unweighted summability.
In regime 3, choose constants and so outside a finite set of lags. Counting the positive integers satisfying shows for , with finite nonnegative constants ; also . The elementary inequality for follows for by maximizing as a function of , whose maximum is ; is immediate. Set and , then integrate and use [F2]: Together with the assumed second moment this proves . Thus the two properties announced in step 1.1 hold in all three regimes, without imposing a power moment in regime 3.
Let . By [F2], [F12] and [F5], Set . Expanding a finite centered sum over and counting at most pairs at each positive lag gives . For consecutive sums the exact expression is . The summable majorant and [F4] prove convergence to the stated , which is nonnegative as a limit of nonnegative variances.
For an integer put , , and . Then . Both processes are strictly stationary by [F3], and is centered and bounded by . Its coefficients are bounded by , whose sum is finite by steps 1.1, 2.1, 2.2. Thus [F1] gives an absolutely convergent covariance series with nonnegative sum , normalized variance convergence to , and .
Put and let be the absolute quantile of . Since , monotonicity of tail probabilities gives , while for by its defining infimum. The centering in does not change any covariance, and . Applying [F2] to the raw, uncentered , followed by [F5], therefore gives Here the integrand is integrable by steps 1.1, 2.1, 2.2, and because is finite almost surely; hence [F4] makes . The same finite covariance expansion as in step 3.1 now gives . This bound uses the raw remainder quantile, without asserting that the centered remainder has that quantile.
The reverse triangle inequality for norms gives by step 4.1. This norm inequality follows by expanding a square and applying [F7] to its mixed term. Using the variance limits in steps 3.1 and 3.2 yields , where the square roots are nonnegative. Therefore .
For fixed real , , obtained by integrating the unit-modulus derivative along the real interval, and [F7] give . By the CLT in step 3.2 and [F9], the second characteristic function tends to , identified by [F8]. Consequently Let and apply step 5.1. Under [A1], [F8]–[F9] prove . The case is also independently supplied by [F10] using absolute covariance summability from step 3.1. All three regimes and the degenerate case follow.
A uniform-mixing CLT for square-integrable stationary sequences
Statement
Assume AC. Let be a centered strictly stationary real sequence with . If then the covariance series converges absolutely and is nonnegative. If , then for , If , Zero covariance long-run variance gives a degenerate CLT under absolute summability gives the degenerate normal limit.
Facts & Assumptions
Given: AC, the centered strictly stationary real sequence , finite second moment, and .
AC is inherited through the rho-mixing CLT. (The Axiom of Choice)
The coefficient bounds for a positive-probability past event and an event in the separated future; is the supremum of their square-integrable variable correlations. (Strong, uniform, and maximal-correlation mixing coefficients)
Real square-integrable random variables satisfy Cauchy–Schwarz. (Cauchy-Schwarz for random variables)
Simple functions are dense in real on a probability space. (Simple functions with finite-measure support are dense in for )
For a centered stationary real sequence of finite second moment, gives absolute convergence and nonnegativity of the covariance variance, and the normal CLT with that variance, including zero. (A rho-mixing CLT for square-integrable stationary sequences)
Variance and covariance are centered second and product moments. (Moments, variance, and covariance on a probability space)
Proof
Fix a gap and a split , put and , and put . Take finite partitions in and in . Write , , , and . For each fixed the sum to zero. If , their positive sum is for , by [F1]. Hence ; when , every and the same bound holds. For each , gives .
For simple and , [F5] and finite Cauchy–Schwarz with weights give , using step 1.1's row and column bounds. The finite Cauchy–Schwarz inequality here is also the probability-space inequality [F2] on a finite set, after normalizing the total weight; if all weights vanish the conclusion is immediate.
For arbitrary real square-integrable -measurable and -measurable , choose simple approximants in their respective spaces using [F3]. Their means, norms and product expectations converge by [F2]; therefore step 2.1's inequality passes to the limit. Apply it to and to get . Taking the supremum over positive-variance variables and then splits in [F1] proves . No reversal of the conditioning direction of has been used.
Summation of step 3.1 gives . The given stationarity, centering, and second-moment assumptions are exactly those of [F4], so that theorem proves absolute convergence and nonnegativity of the displayed covariance series and for positive or zero . Thus every claimed conclusion holds.
Source notes
The phi-to-rho bound is proved directly by finite partitions and approximation. The normal limit is supplied by the locally proved rho-mixing CLT, rather than an external uniform-mixing limit theorem.
Bounded independent vector rows admit a logarithmic cubed Gaussian coupling
Statement
Assume AC. For every there is a finite constant , independent of all row laws and lengths, such that the following holds. Let , , , and let be independent random vectors in with Their distributions need not be identical. There is a coupling with independent such that The coupling can be realized on an extension preserving an already realized raw row as its original marginal. Independence between the raw and Gaussian rows is not asserted.
Facts & Assumptions
Given: AC, , and the independent row with the displayed moment, covariance and almost-sure bounds. All norms are Euclidean.
AC is assumed as in The Axiom of Choice.
A sum of independent centered bounded vectors of norm at most and an independent Gaussian, with total covariance and Gaussian covariance at least , has the uniform analytic cumulant and tilted Fourier bounds of Gaussian anchored bounded sums have uniform cumulant and Fourier bounds. The bounds persist under row isometries and coordinate prefixes in dimensions at most four.
The smooth Gaussian-anchored conditional quantile supplier gives measurable root and conditional bridge transports with their prescribed marginal laws. For fixed positive Gaussian covariance fraction, total covariance , bounded-summand scale , and , their Gaussian-to-target errors are at most at a root and at a split with target parent and target difference , whenever these target norms are at most . Gaussian bridge coordinates are independent of the Gaussian parent. Conditional quantile couplings for Gaussian anchored sums
Finite-dimensional Euclidean spaces are standard Borel, and a joint probability law on two such spaces has a probability-kernel disintegration given either marginal. Finite-dimensional Euclidean space with its Borel sigma-algebra is standard Borel, Disintegration of a joint law on standard borel spaces
A probability-kernel extension depending only on one part of an independent pair preserves independence of that part together with the newly sampled variable from the other part. Conditional independence under kernel extension and conditioning
Markov's inequality applies to nonnegative random variables, and the probability of a finite union is at most the sum of its probabilities. Markov's inequality for random variables, The finite union bound
Gaussian characteristic functions identify their covariance; jointly Gaussian subvectors of zero cross covariance are independent, by factorization of that characteristic function. Characteristic function of a multivariate normal law
Increasing nonnegative integrands may be passed through the integral. Monotone convergence for the integral
A nondegenerate multivariate normal law has its stated Gaussian density. Completing the square in that density gives the standard-normal mgf . A positive-definite multivariate normal law has its Gaussian density
The Euclidean norm obeys the triangle inequality and is convex. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Proof
It suffices to construct a coupling of row laws: by [F3] its joint law can be disintegrated given the raw row, and the measure extends any original realization with that coupling law. This is a probability measure: for disjoint product-measurable events, countable additivity follows by monotone convergence inside the kernel integral. The original marginal remains because the kernel has total mass one. Pad the row to dyadic length , with if , by independent extra entries of the law of . It suffices to control its first entries. Set , where the constant will depend only on .
Every independent row of centered entries, each either one of the raw laws or , obeys a common Bernstein estimate. For a raw coordinate , for , so for exponential-series expansion gives ; the same inequality holds for a standard normal because its log mgf is . For any block of entries, multiply the mgfs and apply [F5] to their exponential, using and both signs, and then the coordinate union bound. Thus for , The constant is uniform in and in the individual row laws.
Consider a compressed dyadic row of length , whose odd entries are Gaussian and whose even entries have their prescribed raw laws. For every aligned dyadic interval of length , both halves contain half Gaussian slots. Its sum has covariance and Gaussian covariance . Its joint normalized parent and difference has covariance and Gaussian covariance . Each bounded contribution to this joint vector has norm at most . Hence [F1] applies with fixed Gaussian fraction , also to all coordinate prefixes, and [F2] supplies root and bridge maps. Rescaling the normalized pair only changes universal constants. No identical-distribution assumption was used: each bounded contribution retains its own original index law.
Put , with a fixed covering [F2], and choose a fixed large . Let be the least dyadic integer at least . If , choose an independent Gaussian root and independent Gaussian split primitives for all dyadic of length greater than . The root map constructs from . Descending the tree, with present, let be the conditional bridge map of given target parent , and put The new primitive is independent of the exposed parents. Its Gaussian conditional law therefore does not depend on . The target split has exactly the prescribed conditional law given by [F2].
At terminal blocks sample the entries of the mixed independent row by their regular conditional law given their sum , with product kernels over distinct terminal blocks. Recover Gaussian entries given sum by fresh independent standard normals and the formula . For the mixed row, independence of the original disjoint blocks implies that their conditional kernels given their sums factor: check the identity on products of block events by the disintegration formula and extend to all events by uniqueness of probability measures. Induction over the binary splits and then these terminal kernels recovers exactly the mixed row product law. For Gaussians, and have zero cross covariance and covariance ; [F6] and induction give the Gaussian product law at each level and after terminal reconstruction. Thus the final two marginal rows have independent entries, despite their coupling.
Apply step 1.2 with to every dyadic mixed sum and to every prefix within a terminal block in either completed marginal row. There are at most such vectors; their unconditional product marginal laws were proved in step 3.1. The union bound gives a good event of failure at most on which each mixed sum satisfies and each terminal residual is at most . At an active split of size , . Choose large enough that both , using the bounds for their ratios to . Their normalized squares are at most . The same holds at the root. Consequently [F2] bounds the root and every split discrepancy by .
At terminal boundaries write . Its endpoint values are zero and . At the midpoint of any active interval , the split formulas give . Convexity of the Euclidean norm shows that adding one level increases the maximal boundary error by at most half the maximal split discrepancy at that level. There are at most levels. Adding the terminal residuals from step 4.1, and using , proves that one round's maximum prefix discrepancy is at most except on probability . If , use instead an independent coupling of the full mixed and Gaussian rows; step 1.2 and a union bound over all their initial prefixes bound their discrepancy by , with the same failure bound. This also covers and arbitrarily large .
Start with an independent Gaussian row of length . Apply the round just proved to that row, putting its designated raw law in every even slot and retaining Gaussian odd slots. This is possible while preserving the existing row: disintegrate the coupling given its Gaussian marginal using [F3] and extend by that kernel. Freeze the newly raw entries. The output mixed marginal is a product law, so the retained Gaussian subsequence is independent of all frozen variables. Repeat on that subsequence using a kernel depending only on it; [F4] preserves independence from the frozen complement. At each stage every newly raw coordinate has the law originally designated for that particular original index. If the active length is less than , replace all its remaining slots by independent raw entries instead; step 5.1's terminal bound applies to the raw and Gaussian marginal prefix sums. After at most stages the final full row has exactly the prescribed product of the possibly different raw laws.
A prefix of the original row meets each active subsequence in an initial prefix, so triangle inequality bounds its total change by the sum of the stage maximum prefix errors. Each is at most and there are at most stages. The total error is therefore at most . Use the same ambient in every stage; since each active length is at most , the total failure probability is at most . Choose large enough, depending only on , that this is at most , including the finitely small lengths by increasing once. As for , . Restricting to the first coordinates and using step 1.1 proves the claimed estimate and the assertion about preserving an existing raw realization.
Source notes
Zaitsev's cited Sections 5–8 give the dyadic conditional-quantile and successive-replacement architecture. The present proof uses half Gaussian anchors and a stopped tree, allowing one extra logarithmic factor and avoiding third-moment matching. It proves only the displayed polynomial-tail logarithmic cubed bound, with constants uniform over the bounded nonidentical row laws; it does not claim the source's sharper exponential bound.
Geometric Markov chains admit delayed adjacent-iid split tours with spare moments
Statement
Assume AC. Let be a probability kernel on a nonempty standard-Borel space , with invariant probability , finite measurable and such that Let be finite measurable, with and for some . For every initial law, on an extension preserving the given chain there are an integer , an almost surely finite delay , a standard-Borel seed space , iid seeds , and measurable finite-word and length functions and with the following properties.
Put and . Almost surely, simultaneously for all , Write , , and . Then and for measurable functions , and For every nonnegative measurable the full-tour occupation identity is For each specified initial law, the same , iid seed law, and post-delay tour law can be used also for the invariant initial law ; only the finite initial segment and coupling change.
In particular is an adjacent-iid centered vector sequence, while the lengths themselves are iid. Full tour words or rewards are not asserted to be iid. No moment of or of the initial finite prefix is required, and the construction uses a useful multistep split chosen from rather than requiring any supplied one-step minorization to have positive stationary mass.
Facts & Assumptions
Given: AC and the standard-Borel chain, geometric bound, invariant law, and finite centered observable in the Statement.
The Axiom of Choice: AC supplies the inherited small-set, path-law, coding, product and disintegration conventions.
Phi-irreducibility for a general-state Markov kernel and A countably generated phi-irreducible chain has a positive-measure small set: on a countably generated irreducible state space there is a positive-measure -small set, with a minorizing probability giving that set positive mass.
Geometric small-set renewal convergence gives an integrable exponential hitting moment: under its scalar geometric-return-probability hypothesis the explicitly specified split atom has an exponential positive return moment and finite mean.
Disintegration of a joint law on standard borel spaces and Finite products of standard Borel spaces are standard Borel: finite standard-Borel joint laws have measurable conditional probability kernels.
Standard borel spaces admit bimeasurable real codings, Standard Borel spaces, and Under the Axiom of Countable Choice, countable products and subspaces of Polish spaces are Polish: standard-Borel spaces have Borel real codings; countable sequences of such coordinates and finite words are standard Borel through Borel subsets of the Polish cube.
Canonical Markov chain on path space, Assuming countable and dependent choice, countable products of arbitrary probability spaces, and Ionescu-Tulcea construction of a Markov chain: canonical kernel paths, iid products and extensions by specified kernels exist under AC.
Measurability of integration against a kernel, Monotone convergence for the integral, and Tonelli's theorem for nonnegative measurable functions on a sigma-finite product: the kernel expectations and nonnegative sums below are measurable, and increasing limits pass through them.
Holder's inequality for integrals, including the endpoint cases and Minkowski's inequality for integrals, including : the indicated Holder and finite-sum bounds hold for .
Invariant and stationary distribution for a Markov kernel and Iterated transition kernels: and iterated transitions are successive kernel integrations; consequently for every by induction.
Dynkin's pi-lambda theorem: equality of probabilities or kernel identities extends from a countable generating algebra to all measurable events.
Proof
Given: The hypotheses in the Statement and AC [A1].
For every -positive measurable set , the geometric bound implies for all sufficiently large , for each . Thus is -irreducible. By [F4], a standard-Borel space has a countable generating algebra, so [F1] gives , with , a probability with , and such that on . Choose an integer so is positive under both and , and restrict to for a finite integer , then normalize. Decrease the minorization constant to . With , the resulting probability satisfies on , and .
The bridge kernel has a jointly measurable uniform sampler. Code its standard-Borel word space into a Borel subset of by [F4]. Its distribution-function evaluations at rational cuts are measurable. The generalized inverse, defined by an infimum over those cuts, is jointly measurable in the endpoints and a uniform variable and has the coded probability law. Its output lies in the coded image almost surely for each endpoint; use a fixed valid word on the Borel complement and decode. For completeness, the word and sequence spaces here are standard Borel: code into a Borel subset ; is Borel in the Polish cube by coordinate tests, and finite words are the Borel union of length-marked sequences with fixed trailing coordinates. The extra uniform sequences have the same property.
Put . Invariance [F8] gives . The geometric bound gives for ; for use . Thus [F2] applies to this exact split. Its atom has a positive return time with for some . The transition out of draws the next skeleton state from independently of the previous state, and assigns its fresh flag with success probability .
Define the substochastic failure kernel , and let be the probability of never obtaining a successful split coin, starting with a fresh coin at . This is measurable by [F5] and [F6], lies in , and satisfies . Integrate against to get , hence . Therefore at every state and for all . Convergence of to implies for every , since for a measurable -valued function its expectation difference is bounded by total variation, by integrating level-set indicators. Integrating this constant against makes it zero. Thus every state and every initial law has a finite first success; its post-transition calendar time is an almost surely finite .
A skeleton core starts from the endpoint just after a success and includes the successive states through the next successful coin, excluding that successful transition's endpoint. A core of length has probability measure The sum over has mass one by step 3.1. Its length has the law of in step 2.1: starting in the split atom, one transition draws its initial state, and the next atom visit occurs at precisely that core's inclusive coin count. Hence its length has an exponential moment. Taking iid copies of this core measure supplies all post-success cores; the first state of the next copy supplies the independent endpoint of the preceding successful transition. Multiplying these core kernels reconstructs the split skeleton law. The analogous initial core starts from the given initial law and ends at its first success, before the endpoint drawn from the first iid core.
To restore the intermediate states, choose for the given initial law . Then and , by the null-set calculation used directly on this displayed sum. Disintegrate the intermediate word conditional on under initial law , using [F3]; for take the unique empty word. Denote the kernel by . Its bridge identity holds for -almost every starting point, simultaneously as a probability identity by [F9]. Every starting distribution encountered in a split skeleton is -absolutely continuous: propagates this through failures, and handles successes. Both success and failure endpoint-pair measures are absolutely continuous with respect to , because their endpoint kernels are bounded by . Thus this bridge identity holds on every transition actually sampled.
Let consist of the th iid core and an independent infinite sequence of bridge uniforms, using [F5]. For a nonterminal core transition both skeleton endpoints lie in . For its terminal successful transition the starting endpoint lies in and the ending endpoint is the first state in ; use a uniform in for its bridge. Concatenate each transition's starting state and its intermediate states, excluding its endpoint. This gives a measurable finite word of length . Concatenate also the finite initial segment from step 4.1, with its terminal bridge ending at the first state's value in . Success and failure endpoint kernels add to , and step 4.2's bridge kernel then restores the original -step path law. Iterating these identities gives the original full -chain law from . Thus these words are the claimed delayed calendar tours in a chain copy, with rewards and absolute envelopes measurable functions of two adjacent seeds.
The invariant identity is . Iteration gives . The remainder's total mass is the stationary probability of no successful coin in its first opportunities, which tends to zero by step 3.1 and monotone convergence. Therefore the potential measure equals . It counts every core state before its fresh coin is sampled, including the state of the terminal successful coin. In particular and .
For nonnegative measurable , set . The event that a skeleton state is counted in the core is determined before its fresh coin and bridge are sampled. Averaging that coin, its endpoint and bridge makes the conditional expected calendar-block cost , because their endpoint mixture is and their restored block law is the original one. Hence [F6] and step 5.2 give This includes the cost inside the terminal successful bridge. Apply the identity to the positive and negative parts of ; follows from its stated th moment and [F7]. Thus .
By step 4.1, for . Starting just after a success, the reconstructed infinite chain has initial law and its original -chain marginal, by step 5.1. Since , invariance gives at every calendar time. For , Minkowski and Holder [F7], first on finite sums, imply The exponential length moment makes the series finite; [F6] passes through the increasing finite sums. This proves every claimed spare moment, while .
The constructed state path, its finite delay, and its iid seed sequence have a joint law on standard-Borel spaces by the explicit countable coding in step 1.2. Disintegrate that law conditional on the full state path by [F3], and evaluate the auxiliary kernel at the actual given chain path. Extend its probability space by [F5]. The original entire process is preserved, and its joint law with the seeds and delay is the constructed one. Thus the word equalities, iid seed law, occupation identity and moment conclusions hold on this extension. The initial segment is almost surely finite by step 3.1; since is finite, its absolute reward is finite, even without an initial-law moment bound.
Each length depends only on , so the lengths are iid with positive finite mean. Each reward and envelope depends on , so their joint sequence with the centered lengths is adjacent-iid, with and all the stated moments. Disjoint seed pairs are independent, giving one-dependence without claiming iid full tours. The same and core law were selected from , and the single bridge version in step 4.2 uses a dominating both the specified initial law and . Thus the same seed law and word functions work from either initial law; only the initial segment and its lifting kernel change. All properties in the Statement follow. The chosen split came from step 1.1 and never used positivity of any supplied one-step minorization, which can therefore be vacuous.
Heterogeneous bounded atom vector rows admit a logarithmic cubed Gaussian coupling
Statement
Assume AC. Fix and . There is a finite such that the following holds uniformly over , integers , and . Let be independent centered vectors with representations where every is finite, all displayed atom and Gaussian terms are independent, every is centered, and every is a centered Gaussian vector, possibly degenerate. Suppose There is a coupling with independent satisfying The constants do not depend on the atom counts or individual laws. Empty atom sums and zero Gaussian terms are allowed.
The coupling laws can be chosen as a measurable probability kernel in a supplied standard Borel parameter : if the atom laws and Gaussian covariance matrices are measurable kernels and functions of , with the hypotheses holding for each , there is a coupling kernel with the stated product marginals and the displayed bound for each . Thus a random environment with these conditional independent atom laws can first be sampled, and this coupling kernel then applied conditionally. This assertion concerns the conditional laws supplied in the hypotheses and does not assume a selection of regular conditional laws separately at each parameter.
Facts & Assumptions
Given: The atom and Gaussian representations, covariance bounds, parameters , and, for the kernel assertion, the supplied measurable atom laws and Gaussian covariance matrices. Matrix inequalities are quadratic form inequalities and all vector norms are Euclidean.
AC is assumed. The Axiom of Choice
Independent sums of centered bounded vectors and a Gaussian, of total covariance and Gaussian covariance at least , have the uniform cumulant and tilted Fourier bounds of Gaussian anchored bounded sums have uniform cumulant and Fourier bounds.
The Gaussian anchored conditional quantile supplier constructs measurable root transports and, for joint dimension , measurable transports of a fresh independent Gaussian residual conditional on the target parent. With covariance , atom bound , and Gaussian covariance fraction , the target value error is bounded by when and . Conditional quantile couplings for Gaussian anchored sums
Real symmetric matrices have an orthonormal eigenbasis. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis
A Gaussian has its covariance characteristic function; zero cross covariance implies independence for Gaussian subvectors by factorization. A positive definite Gaussian has its Gaussian density. Characteristic function of a multivariate normal law, A positive-definite multivariate normal law has its Gaussian density
Integrating a measurable function against a supplied probability kernel is measurable, and finite products of probability kernels are probability kernels. Measurability of integration against a kernel, Finite products of probability kernels are probability kernels
Nonnegative integrals obey Tonelli and monotone convergence. Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Monotone convergence for the integral
A kernel extension depending only on one part of an independent pair preserves independence from its complement. Conditional independence under kernel extension and conditioning
Exponential Markov bounds, the finite union bound, and Euclidean triangle inequality apply. Markov's inequality for random variables, The finite union bound, Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Expectations of finite products of nonnegative functions of independent variables factor. Expectations factor over finite products of independent random variables
Proof
For a coordinate of one centered atom, for every integer . Absolute exponential-series integration is legitimate because ; [F6] permits summing its nonnegative absolute terms, and its positive and negative parts then give the signed series identity. For , centering cancels the linear term and for , so A centered Gaussian coordinate has logarithmic mgf , by completing the square in its positive-variance eigencoordinates using [F3] and [F4]; zero-variance coordinates vanish. Independence multiplies the finitely many atom and Gaussian mgfs. Since their variances add to each coordinate variance of , at most , any sum of entries, also with arbitrary matching Gaussian replacements, has coordinate log mgf at most . For a coordinate threshold use and exponential Markov for both signs. This gives tail at most . The coordinate union bound therefore gives, uniformly in atom counts, This holds for every contiguous interval in any compressed ordering.
Consider a dyadic mixed row whose odd entries are independent matching Gaussians and whose even entries have their raw laws. At a node of entries let its two child covariance matrices be , put , and, for its target child sums , define Independence gives , , and . The symmetric matrix equals and has inverse ; hence it is positive definite without assuming commutation. Since each child has entries, [F3] gives The final inequalities follow by inverting, adding, and reinverting the quadratic form bounds.
The transformed pair has covariance . A raw atom in its left or right child contributes respectively Their norms are at most , since , , and . Each child contains half Gaussian entries, so its Gaussian covariance is at least , at least times its total covariance with . Congruence of quadratic forms therefore makes the Gaussian covariance of at least . All its bounded contributions are the independently transformed original atoms, and all Gaussian contributions combine into an independent Gaussian. Thus [F1] and [F2] apply directly to , including its coordinate prefixes. The conditioning is on its first coordinates. Undoing the separate parent and residual whitening, whose inverse norms are bounded by , gives a transport of a fresh to conditional on , with error when and the two target norms are at most . The root has the analogous transport from with error . There is no additional factor from the atom counts or from .
Temporarily take a dyadic current length , set , and let be the least dyadic integer at least , with to be fixed. Suppose . Start with an independent Gaussian row having the individual covariances . Its root sum and residuals at the active nodes are mutually independent Gaussian vectors of covariances and : the cross covariances vanish, and independence follows recursively from [F4]. Apply the root transport in step 2.1 and then, descending the tree, set and Stop at length . Each fresh residual is independent of the already exposed parents, so every target split has exactly its prescribed conditional law. For Gaussians the child covariance is and the cross covariance is , with the analogous right covariance .
There is an explicit terminal target kernel, also measurable in any supplied environment. In a terminal mixed block let be the product law of its raw positions, let be their sum, and let be the sum of the matching Gaussian-position covariances. Given the target block sum , sample the raw-position vector from The denominator is finite and strictly positive everywhere by [F4]. Fill the Gaussian positions conditional on their remaining sum : for their covariance matrices and fresh independent , put and . The residual vector has sum zero, zero cross covariance with , and covariance blocks ; thus [F4] identifies its Gaussian conditional law. Multiplying the posterior by the total-sum density and integrating in proves the original raw/Gaussian product law, first on rectangles and then by monotone approximation. This verifies the terminal conditional law at every sum, without choosing regular conditional versions. Independent terminal kernels and the split laws of step 3.1 consequently recover exactly the entire mixed row product law by induction through the disjoint binary blocks. The existing Gaussian row is unchanged, and already has its required product law.
Use step 1.1 and a union bound for every dyadic target sum and every within-terminal prefix in the completed target and Gaussian rows. There are at most such vectors. Except on probability , each target sum of size has norm at most . The target residual is , so its norm satisfies the same bound with a different . At active sizes , choose large enough that these target norms are at most ; their ratios to are bounded by . Their squared norms divided by are at most . Step 2.1 therefore bounds the root discrepancy and every residual discrepancy by . The within-terminal prefix norms are at most , since . No conditional concentration inside a terminal block has been assumed.
Write . Its child equations are and . For a prefix ending at a terminal boundary of , let be the sum of its individual covariance matrices. Its discrepancy has the exact expansion If the boundary lies in the left child, ; if it lies in the right child, . To prove the formula, in the left case substitute the child's leading term and cancel ; in the right case add to , obtaining leading coefficient and residual coefficient . Induction from terminal blocks proves the expansion. This cancellation uses the order and never commutes two matrices. If a partial child has of its entries, its coefficient norm is at most . Hence and . There are at most boundary ancestors. The expansion and step 5.1 give terminal-boundary error ; adding its terminal prefix residual gives the same bound for all prefixes.
If , couple a raw or mixed row independently to a matching Gaussian row. Their unconditional concentration bounds from step 1.1 and a union bound over prefixes give maximum discrepancy at most except on probability . This covers the small sizes, including a row of one entry, and large . Thus in all cases a mixed row can be coupled directly to an already realized matching Gaussian row with stage error at most and bad probability at most .
Pad the initial length to dyadic by independent entries; for take . These added laws have empty atom sums and preserve all covariance bounds. Let be the matching independent Gaussian row, saved as the final reference row, and initially mark every position active. Inductively, the current full row consists of independent frozen raw entries and an independent active Gaussian row, all with the designated covariances. When its active length is at least , use that active Gaussian row and fresh terminal randomness in steps 3.1--6.1 to construct the target mixed row, with designated raw laws in even slots and matching Gaussian laws in odd slots. Replace the entire active row by this target row, including its odd coordinates, to obtain from . Freeze the target evens and use the target Gaussian odds as the next active row. The extension kernel depends only on the previous active row, so [F7] makes the entire target independent of the previously frozen complement; its product marginal from step 4.1 makes its new raw evens independent of its Gaussian odds. This proves the induction. At active length below , replace every remaining active coordinate by independent designated raw entries using step 7.1 and stop. After at most stages, the final row has the intended raw product law. The saved row retains its Gaussian product law, although subsequent active Gaussian rows need not equal its corresponding coordinates.
Write . Frozen entries do not change at stage , and every original prefix meets its active positions in an initial compressed prefix. Thus the stage estimate in steps 6.1--7.1 bounds by outside its stage bad event; this includes all changes to target Gaussian odds. The exact identity and triangle inequality therefore give total error at most . Each stage has the Gaussian input marginal and mixed target marginal used in its estimate by step 8.1, so its unconditional bad probability remains at most , where is its active length. With the same ambient , a union bound over at most stages gives . Choose sufficiently large depending only on to make this at most for all . Taking , and restricting to the first entries proves the numerical assertion.
For the parameter assertion, every entry law is a measurable kernel obtained by applying addition to the finite product of its supplied atom and Gaussian kernels, using [F5]. Covariance entries are measurable integrals. In dimensions one and two positive definite inverse square roots are continuous functions of their matrices: in dimension two, for , , as verified in an eigenbasis by [F3]; the semidefinite zero case is defined by zero and the formula extends continuously. Thus the Gaussian starts and all transforms in steps 1.2, 2.1 depend measurably on the parameter. Mixed transformed joint and prefix densities are Gaussian convolution integrals of their raw product kernels, so are jointly measurable by [F5] and strictly positive. Their conditional CDFs are measurable density-ratio integrals. The coordinate quantiles in [F2] are measurable rational infima of those CDFs; hence its root and residual maps are jointly measurable in the parameter and all spatial inputs. The terminal posterior in step 4.1 is likewise a measurable kernel with positive measurable denominator, followed by the explicit Gaussian bridge kernel. Finite kernel products and successive extensions preserve measurability: for one extension define its law by integrating the next kernel against the preceding law; [F5] gives measurable evaluations and [F6] gives countable additivity by monotone convergence of disjoint indicator sums. Iterating this finite construction gives a coupling kernel , and every previous marginal and probability bound holds at each parameter. No pointwise selection of conditional-law versions or coupling laws is used.
Source notes
The source supplies the dyadic conditional-quantile architecture. The heterogeneous covariance formulas in steps 1.2 and 2.1, the stopped Gaussian-anchor replacement, and the parameterwise terminal posterior are proved locally. The bound is deliberately weaker than the source's optimized strong approximation estimates. Whole grouped entries need not be bounded: the boundedness hypothesis concerns their independent centered atoms.
Sufficient MCMC CLT regimes
Statement
Assume AC. Let be countably generated, let be a -irreducible probability kernel (Phi-irreducibility for a general-state Markov kernel, Measure kernel and probability kernel), aperiodic (Aperiodicity for a general-state Markov kernel), and Harris recurrent relative to (Harris recurrence relative to an irreducibility measure) with invariant probability , and let be a stationary -chain whose canonical left shift preserves its path law and whose shift-invariant measurable events have probability zero or one. We take the initial law to be , which implies stationarity by Invariant initial law makes a Markov chain stationary. Let be Borel with . Define only when the covariance series is absolutely convergent and the CLT is established, as in Long-run variance for a stationary Markov chain.
Each of the following six conditions is sufficient for absolute convergence of the covariance series and
- has polynomial total-variation rate of order with an integrable envelope (Polynomial ergodicity with an integrable rate envelope), and is essentially bounded under (there is with -almost everywhere).
- has polynomial total-variation rate of order with an integrable envelope (Polynomial ergodicity with an integrable rate envelope), for some (The space as the quotient by null functions), and .
- is geometrically ergodic (Geometric and uniform ergodicity for MCMC CLTs) and for some (The space as the quotient by null functions).
- is geometrically ergodic (Geometric and uniform ergodicity for MCMC CLTs) and , with .
- is geometrically ergodic (Geometric and uniform ergodicity for MCMC CLTs), reversible with respect to in the symmetric-flow sense (General-state reversibility by symmetric flow measure), and (The space as the quotient by null functions).
- is uniformly ergodic (Geometric and uniform ergodicity for MCMC CLTs) and (The space as the quotient by null functions).
In case 4, the displayed logarithmic moment implies : on one has , and on one has .
When , the zero-variance lemma Zero covariance long-run variance gives a degenerate CLT under absolute summability gives the degenerate normal limit from absolute covariance summability. Finite variance of alone is not one of the sufficient conditions.
Facts & Assumptions
Given: AC, the stationary Harris-ergodic chain and centered real function in the Statement, and any one of its six sufficient regimes.
The Axiom of Choice is assumed by the local probability, conditional-expectation and sequence-CLT suppliers. The Statement already includes this hypothesis.
A measurable observation does not increase stationary mixing coefficients makes stationary and bounds each of its mixing coefficients by the corresponding chain coefficient.
Polynomial ergodicity with an integrable rate envelope gives the pointwise polynomial total-variation bound with . Geometric and uniform ergodicity for MCMC CLTs gives the geometric bound with finite measurable envelope, or with bounded envelope in the uniform case.
The bounded-future Markov property and iterated-kernel identity are Markov property for bounded future path functionals and Chapman-Kolmogorov equations. The event-supremum convention is Total variation distance for probability laws, and known factors can be taken out using Taking out what is known.
Strong-mixing CLTs under bounded, moment, and logarithmic-moment conditions proves absolute covariance summability and the normal limit for bounded observations with summable alpha, for the stated alpha-summability condition, or for exponential alpha with the logarithmic moment.
Harris geometric ergodicity admits an integrable total-variation envelope supplies an integrable geometric envelope under the stated Harris, irreducibility and aperiodicity hypotheses. A total-variation rate bounds strong-mixing coefficients then bounds alpha by that geometric rate.
Geometric ergodicity and reversibility give exponential rho-mixing gives exponential rho for a stationary reversible geometrically ergodic chain under the stated recurrence hypotheses. A rho-mixing CLT for square-integrable stationary sequences proves the covariance-summable CLT when rho is summable.
Uniform ergodicity bounds the uniform-mixing coefficient gives exponential phi under uniform ergodicity. A uniform-mixing CLT for square-integrable stationary sequences proves the covariance-summable CLT when the square roots of phi are summable.
Zero covariance long-run variance gives a degenerate CLT under absolute summability handles zero covariance long-run variance by a finite covariance expansion. Long-run variance for a stationary Markov chain identifies the covariance sum only after absolute convergence and the CLT have been established.
For , , and for the series diverges proves convergence of for every real .
Proof
Given: The hypotheses and any one of the six regimes in the Statement.
Put . This is a centered stationary real sequence by [F1]. It is square-integrable in every case: boundedness suffices in case 1, the power moment suffices in cases 2 and 3 on the probability space, and cases 5 and 6 assume . In case 4, split at : on the bounded part , and on the other part , so the stated logarithmic moment implies . All sequence suppliers below therefore have their centering, stationarity and second-moment premises.
Suppose a pointwise rate is available with . For any past event and separated future event , [F3] supplies with and . Approximating uniformly by dyadic averages of its superlevel indicators gives : each indicator discrepancy is bounded by the defining event supremum. Consequently Taking the suprema and then applying [F1] gives . No measurability of an uncountable total-variation supremum is needed, because its measurable envelope bounds each event discrepancy pointwise.
In case 1, take in step 1.2. Since , grouping bounds each block of by ; [F9] makes the resulting geometric series finite. Thus . The bounded branch of [F4] gives absolute covariance summability and the asserted normal limit for the normalized sums of . This uses the whole approved range .
In case 2, put . Step 1.2 gives . The exact assumption says , so the same dyadic argument and [F9] make these bounds summable. The power-moment branch of [F4] proves the same covariance and distribution conclusions.
In cases 3 and 4, [F5] replaces the possibly nonintegrable original envelope by a finite integrable geometric envelope with some rate . Thus by [F1] and [F5]. In case 3, every positive power of this bound is summable, so the branch of [F4] applies. In case 4, its exponential/logarithmic-moment branch applies directly. Each yields absolute covariance summability and the normal limit without strengthening the given moment assumption.
In case 5, reversibility and the geometric-rate assumptions allow [F6] to give for some . By [F1], , so its series is summable. The rho-mixing theorem in [F6] proves the required covariance-summable CLT.
In case 6, [F7] gives for some . By [F1], . The uniform-mixing theorem in [F7] therefore gives the same conclusions.
In all six cases the covariance series is absolutely convergent, its sum is finite and nonnegative, and the normalized sum converges to the centered normal with precisely that covariance sum. If the sum is zero, [F8] also proves the degenerate limit directly in , so no division by an asymptotic variance is made. Stationarity makes and equal in law. The established limit is therefore exactly the displayed limit in the Statement, and only now does [F8] identify its variance as .
Zero-variance geometric Markov sums have a uniform strong rate
Statement
Assume AC. Let be a probability kernel on an arbitrary measurable space with invariant probability . Suppose its stationary Markov chain has strong-mixing coefficients for , with and . Let be finite measurable with , for a finite , and suppose the absolutely convergent covariance sum is zero: Then for every , under the stationary law, The empty sum at is zero. No ergodicity assumption is needed beyond the stated stationary mixing bound.
If additionally for every , the same assertion holds under every initial probability law, without a moment assumption on that law. In particular, it holds under every initial law whenever has an everywhere finite measurable common geometric total-variation envelope in the sense of Geometric and uniform ergodicity for MCMC CLTs and has the stated moment, centering and zero covariance sum. No positive one-step minorization, irreducibility or Harris hypothesis is added to this last sufficient condition.
Facts & Assumptions
Given: AC, the stationary geometric strong-mixing hypotheses, finite centered real , the zero covariance sum, and any fixed . For the all-start conclusion assume also the stated pointwise total-variation convergence; for the final sufficient condition assume the everywhere geometric envelope instead.
AC is assumed. The Axiom of Choice
Under the stationary geometric strong-mixing and moment hypotheses, there is a finite measurable with almost everywhere; the covariance series is absolutely convergent and . Geometric strong mixing gives an Lp Poisson solution
For , is square-integrable with , where . Kernel conditional expectations for square-integrable Markov observables
On a probability space for ; Cauchy–Schwarz bounds products of variables. Conditional expectation preserves means, and a known factor can be taken out when the products are integrable. Finite-measure includes into for , Cauchy-Schwarz for random variables, Basic algebra and order properties of conditional expectation, Taking out what is known
Monotone convergence passes increasing finite nonnegative sums through expectation; summable event probabilities imply only finitely many occurrences almost surely. Monotone convergence for the integral, First Borel-Cantelli lemma for events
Under AC the canonical Markov path law exists for every initial law. For a bounded product-measurable path functional , is measurable and . An invariant initial law makes the chain stationary. Canonical Markov chain on path space, Markov property for bounded future path functionals, Invariant initial law makes a Markov chain stationary
Total variation is . Total variation distance for probability laws
An invariant kernel with an everywhere finite measurable common geometric total-variation envelope admits such an envelope integrable under , possibly at a slower geometric rate. That integrable bound implies for the stationary chain. Harris geometric ergodicity admits an integrable total-variation envelope, A total-variation rate bounds strong-mixing coefficients, Geometric and uniform ergodicity for MCMC CLTs
Proof
Apply [F1] to obtain with almost everywhere; choose its finite measurable version as in that supplier. By [F3], . Stationarity makes the Poisson identity hold at every simultaneously off a countable union of null events. By [F2], on this common event . Summing gives , with and . No ergodic theorem or central limit theorem has been used.
Put . Stationarity of transition pairs gives for every . If , then is -measurable, its product with is integrable by [F3], and [F2]–[F3] give . Thus . The inequality and stationarity give . Cauchy–Schwarz [F3] and step 1.1 imply . By [F1], the same normalized second moments tend to , so . A nonnegative variable of zero expectation is zero almost surely; applying this to each and taking the countable intersection shows that all vanish simultaneously. Hence for every on one probability-one event.
Fix . Pointwise, the number of positive integers with is at most . Stationarity and [F4], applied first to finite sums and then to their increasing limit, therefore give . Borel–Cantelli [F4], intersected over positive rational , yields almost surely. On that event, for every the sufficiently late terms obey for , while the maximum of the finitely many earlier terms divided by tends to zero. Thus . Step 2.1 gives , proving the stationary conclusion.
For the all-start part, define the original path event . It is product-measurable, since the condition is a countable combination of inequalities for finite measurable sums. It is invariant under deleting the first coordinate: if , the shifted sum is , and the reverse identity is with the second term the shifted sum of length . Since is finite and and tend to one, either normalized convergence implies the other. Iterating gives invariance under deleting any fixed finite prefix. Also, on every real path, is equivalent to : the reverse direction is immediate, and the forward direction follows by discarding a finite initial segment and bounding all later by . Thus it suffices to transfer the event , rather than any Poisson identity outside a stationary full-measure set.
Let . By [F5], is measurable and . The canonical stationary chain has the same path law as the stationary chain already considered, so step 3.1 gives ; this equality also follows by applying [F5] at time zero and taking expectations. For every state and integer , prefix-deletion invariance and [F5] give . These are identities for each canonical starting state, not simultaneous versions of a conditional expectation over an uncountable family of starts.
For any probability measures and measurable , let . Then and [F6] gives . Letting proves the same bound for . Under the additional pointwise total-variation assumption, apply it to and . Step 5.1 gives for every , hence everywhere. For any initial law , [F5] at time zero gives . The deterministic equivalence in step 4.1 transfers the uniform maximal rate as well. No moment or support condition on is needed.
Finally suppose the everywhere common geometric bound in the Statement holds. Since its envelope is finite at every state, it implies in total variation for every . By [F5] the canonical invariant-initial-law chain is stationary. The stronger arbitrary-measurable-space clause of [F7] supplies a finite integrable geometric envelope, and its mixing implication supplies the stationary geometric strong-mixing bound needed by [F1]. Therefore the proof in steps 1.1, 2.1, 3.1, 4.1, 5.1, 6.1 applies to the stated and zero covariance sum under every initial law. The zero function, a one-state chain, empty partial sums and zero martingale variance are all included; an empty state space admits no invariant probability and does not satisfy the hypotheses. Strict is used by [F1] and the comparison. AC is inherited precisely through [F1], [F2], [F3], [F5] and [F7]. This completes both the stationary result and the asserted all-start sufficient conditions.
Remarks
This is the zero-variance rate only. For a moment with , choose strictly between and . The proved rate is stronger than . It therefore supplies the degenerate error estimate needed in a zero-Brownian-coefficient branch without using regeneration, a martingale CLT, KMT, or a nondegenerate SIP. It establishes no positive-variance Brownian coupling.
The Jones and Doukhan–Massart–Rio full PDFs were retrieved, and the source locations specified above were inspected. Their covariance and mixing methods document provenance; the bounded-dual Poisson supplier, zero-martingale variance comparison, integrated-tail rate and deletion-invariant arbitrary-start transfer are all proved locally here or in the explicit mathematical dependencies.
Finite-moment iid vectors admit a strong Gaussian approximation
Statement
Assume AC. Let be iid centered vectors in with for some . For any fixed real numbers , on an extension preserving the original process there are independent Gaussian vectors , where , such that almost surely The extension can also carry a standard two-dimensional Brownian motion with for every integer . Thus, in particular, the partial-sum Brownian approximation error is almost surely. Singular covariance, including , is permitted. Neither nor is required to be independent of .
Facts & Assumptions
Given: AC, the centered iid vectors , their finite th absolute moment, and fixed real exponents .
AC is assumed as in The Axiom of Choice and the construction interfaces below.
For centered iid vectors of dimension one or two with covariance and norm at most , the bounded-vector Gaussian coupling has a bound on the maximum partial-sum discrepancy, with exceptional probability at most for every fixed . Constants are uniform in and the bounded law. Bounded independent vector rows admit a logarithmic cubed Gaussian coupling
A real symmetric matrix has an orthonormal eigenbasis; in particular positive definite matrices have square roots and inverse square roots obtained in that basis. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis
Events with summable probabilities occur only finitely often almost surely. The first Borel-Cantelli lemma for measures
Finite-dimensional Euclidean spaces are standard Borel, and a joint probability law on two standard-Borel spaces has a conditional probability kernel for one coordinate given the other. Finite-dimensional Euclidean space with its Borel sigma-algebra is standard Borel, Disintegration of a joint law on standard borel spaces
Supplied history-dependent probability kernels, with an arbitrary measurable initial probability space, admit an iterated countable extension preserving their finite-prefix laws. Ionescu-Tulcea construction of a Markov chain
A standard normal variable satisfies for , by the Mills upper bound and symmetry. Two-sided Mills bounds for the standard normal tail
Gaussian vectors admit their specified covariance, including singular covariance, by linear transformation of independent standard-normal coordinates; finite independent Gaussian vectors form a jointly Gaussian stack. Multivariate normal law, including singular covariance, A finite stack of independent Gaussian vectors is jointly Gaussian
Independent Brownian copies exist in each finite dimension, and countable products of probability spaces exist under the choice principles supplied by AC. Existence and scaling of -dimensional Brownian motion, Assuming countable and dependent choice, countable products of arbitrary probability spaces
Zero cross-covariance between jointly Gaussian blocks implies their independence. Standard vector Brownian motion is characterized by continuous paths, zero start and independent normal vector increments of the specified variances. Uncorrelated jointly Gaussian blocks are independent, -dimensional Brownian motion
Proof
Put . For each vector in a deterministic finite orthonormal basis of , , so almost surely. Intersect over that basis and the countably many indices to get one full event on which every . If , take and adjoin an independent standard two-dimensional Brownian motion by [F8]; both approximation errors vanish. Otherwise use a deterministic orthonormal eigenbasis of on and work in dimension or , preserving Euclidean norms. Its covariance, still denoted , is positive definite, with smallest eigenvalue by [F2]. The final construction will be mapped back by the same isometry.
Let , , , and . Set and for . The pointwise bound for gives . Since , this is summable in , and [F3] gives no discarded observations in every sufficiently large block almost surely. Also implies . On a block without discarded observations the maximum partial-block difference between and is at most , since .
The truncated covariance satisfies the exact identity . Hence . For all sufficiently large deterministic , its eigenvalues are therefore bounded below by and above by a fixed constant. The standardized truncated vectors are centered iid within each block, with covariance and norm bounded by , after increasing a fixed , because and .
Apply [F1] with separately on each such block, obtaining iid standard Gaussian vectors , . Put . The operator norms of are uniformly bounded, so the maximum partial-sum discrepancy between the truncated row and this Gaussian row is at most , except with probability at most . Here is independent of the varying truncated law and , by the uniform constants in [F1].
Set . If and , then . Taking Frobenius inner product with gives . Treating separately and otherwise dividing, the eigenvalue bounds of step 3.1 and equivalence of these matrix norms in dimension at most two give . This uses positive definiteness only on , with no inverse on covariance-kernel directions.
To realize these couplings while preserving the original process, use [F4] on each finite coupling law to obtain a kernel giving its Gaussian block conditioned on its standardized truncated block. Evaluate that kernel at the actually given . Apply [F5] with initial coordinate the entire original space and each subsequent coordinate its finite Gaussian block, the kernel depending on through that block alone. This extension keeps the original marginal. Independence of the original disjoint blocks implies, by integrating finite block rectangles, that the unconditional paired-block laws factor into the finite coupling laws; thus all across all blocks are independent standard Gaussians. Assign independent Gaussian blocks to the finitely many early blocks before invertibility. Their total discrepancy is finite almost surely. On the resulting extension, [F3] and the summable bounds of step 4.1 give the stated partial-block coupling bound for all sufficiently large blocks almost surely.
Put . Each coordinate of , , is Gaussian with variance at most . By [F6], a union bound over the coordinates and times bounds the probability that its maximum Euclidean norm exceeds by , summable in for fixed sufficiently large . If the discrepancy is zero. Thus [F3] gives a maximum partial-block covariance-replacement error almost surely. Its power is less than because , so this error is .
Combine steps 2.1, 5.1 and 6.1 by triangle inequality. On one full event the maximum partial-block error between and is at most for every , after absorbing finitely many early errors into a finite random constant. Concatenate the independent Gaussian blocks into . For , the global maximum is at most . Because , this bound is . Mapping from back to preserves the error norm and gives covariance equal to the original .
To produce the Brownian driver, retain the independent standardized Gaussian increments in the effective eigenbasis of . If , adjoin independent standard Gaussian increments in its orthogonal complement by [F7] and [F8]. For each coordinate and each integer interval, adjoin an independent Brownian motion by [F8], independent of all existing variables, and form its bridge , . The bridges are centered Gaussian with covariance and zero endpoints. Interpolate each standardized Gaussian random walk on by its linear increment plus , independently for the two coordinates. By [F7], finite stacks are Gaussian; direct covariance calculation gives covariance within each coordinate and zero cross-coordinate covariance. Paths are continuous, including at integer joins. The covariance calculation makes vector increments over disjoint time intervals uncorrelated; finite joint Gaussianity and [F9] make them independent with covariance equal to interval length times , so [F9] identifies standard two-dimensional Brownian motion. Applying the deterministic orthogonal eigenbasis map makes a standard two-dimensional and yields for every integer . Step 7.1 supplies the asserted Brownian approximation and its weaker rate.
Adjacent iid vector functionals admit a finite-moment strong approximation
Statement
Assume AC. Let , let be standard Borel, and let be iid -valued random elements with law . Let be measurable with for some finite . Put , and define Then is positive semidefinite. For every fixed , on an extension preserving the entire original process there are independent centered Gaussian vectors with covariance such that The extension can also carry a standard -dimensional Brownian motion with for every integer . Singular , including , is allowed. Neither the Gaussian row nor is required to be independent of the original process.
Facts & Assumptions
Given: AC, the standard-Borel innovation space, the iid sequence , the centered measurable with its finite th moment, and . All vector norms are Euclidean and matrix inequalities are quadratic-form inequalities.
AC supplies the inherited Gaussian, disintegration and countable-extension choices. (The Axiom of Choice)
In dimension at most two, independent entries that are finite sums of independent centered atoms bounded by , plus an independent Gaussian, and whose covariance matrices lie between and , admit a uniform Gaussian coupling with failure at most for a row of length at most . The coupling is a measurable kernel in a supplied standard-Borel environment. (Heterogeneous bounded atom vector rows admit a logarithmic cubed Gaussian coupling)
Independent centered vectors of dimension at most two, covariance and norm at most , admit a Gaussian coupling of maximum partial-sum error with failure at most , uniformly over all row laws. An existing raw realization can be preserved. (Bounded independent vector rows admit a logarithmic cubed Gaussian coupling)
Bounded independent centered real sums satisfy Bernstein's inequality with variance sum and atom bound : their two-sided tail at is at most , with zero sum when . (Bernstein inequality for bounded independent sums)
Conditional Jensen and Cauchy–Schwarz bound the moments and products below; Tonelli and Fubini allow the indicated nonnegative and absolutely integrable iterated expectations. (Conditional jensen inequality, Cauchy-Schwarz for random variables, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product)
Real symmetric matrices have orthonormal eigenbases. Gaussian characteristic functions specify their covariances; independent Gaussian vectors form jointly Gaussian stacks, and uncorrelated jointly Gaussian blocks are independent. Standard-normal tails obey for . (Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis, Characteristic function of a multivariate normal law, A finite stack of independent Gaussian vectors is jointly Gaussian, Uncorrelated jointly Gaussian blocks are independent, Two-sided Mills bounds for the standard normal tail)
Standard-Borel spaces and their finite products admit disintegration of joint laws; finite Euclidean spaces are standard Borel. Kernel integrals and finite kernel products are measurable. Successive supplied kernels admit a countable extension preserving the initial probability space. (Standard Borel spaces, Finite products of standard Borel spaces are standard Borel, Finite-dimensional Euclidean space with its Borel sigma-algebra is standard Borel, Disintegration of a joint law on standard borel spaces, Measurability of integration against a kernel, Finite products of probability kernels are probability kernels, Ionescu-Tulcea construction of a Markov chain)
Markov's inequality and the finite union bound control exceptional events; summable exceptional probabilities occur only finitely often almost surely. (Markov's inequality for random variables, The finite union bound, The first Borel-Cantelli lemma for measures)
Independent Brownian copies and countable products exist under AC; continuous centered Gaussian paths with the Brownian covariance and independent coordinate increments give standard vector Brownian motion. (Existence and scaling of -dimensional Brownian motion, Assuming countable and dependent choice, countable products of arbitrary probability spaces, -dimensional Brownian motion)
Variance and covariance are centered square and product moments. (Moments, variance, and covariance on a probability space)
Proof
Write , and , defining them arbitrarily on the null sets where the integrals are not finite. By [F4], they are measurable, have finite th moments, are centered, and for independent of law . Expanding the products using these two zero conditional means gives Indeed the variance of is the sum of the covariances of , and , while the lag-one covariance is . This proves both the displayed identity and nonnegativity without any limit theorem.
We derive the finite-horizon iid bound needed for this changing row directly from [F2]. Suppose are iid centered vectors with , and covariance whose eigenvalues lie in a fixed with . Put and , . The tail union probability is at most , and on its complement the maximum sum correction is at most . Also , so . Uniformly for all sufficiently large the eigenvalues of are in . The whitened have covariance and norm at most , so [F2] yields independent standard with maximum coupling error to at most , outside probability . All constants are uniform in the iid law under these moment and eigenvalue bounds.
For a deterministic orthonormal basis of , step 1.1's sum of nonnegative scalar variances implies that each basis component of is zero almost surely and the same component of is zero almost surely. Hence the kernel projection of is , where is the kernel projection of . It is an exact endpoint coboundary simultaneously for all , after removing a countable union of null sets. Its partial sums are . For every , the union and Markov bounds over give probability at most that their maximum exceeds . This is summable, so [F7], first for rational and then for all between successive powers of two, gives maximum . If take and adjoin a Brownian motion by [F8]. Otherwise work in its deterministic positive-covariance range of dimension one or two and later restore this kernel projection.
For independent anchors and interior variables , all of law , put . Its conditional mean given is zero. Its average conditional covariance is To check this, substitute the decomposition from step 1.1 to write . Every cross product vanishes after integrating an endpoint using the zero conditional means of ; the five individual covariance contributions give the formula. Thus , so is positive definite on the range retained in step 2.1. A fixed invertible linear map there makes ; its inverse will be restored at the end, and all subsequent constants may depend on this fixed map.
Choose fixed , and consider an epoch of consecutive triple indices. Put , and . In this epoch replace the normalized active by . Then , , and its means and degenerate component have uniformly bounded th moments by [F4]. They converge in to their untruncated versions: the squared tail is at most , and conditional averaging is an contraction by coordinatewise [F4]. Consequently the average conditional triple-noise covariance tends to . For all sufficiently large deterministic epochs, for fixed . Finitely many earlier epochs can be coupled arbitrarily and contribute only finite errors.
Let , and in the epoch, relabeling its first triple index as zero. Conditional on its anchors , the interior pairs are independent, and the centered triple noises are independent centered vectors bounded by . Define the symmetric matrices and , where these covariances integrate independent . The left and right endpoint terms are independent given the anchors, so Although need not be positive semidefinite, are positive semidefinite, their norms are at most , and . Moreover, for every interval of triple indices, These are identities of measurable anchor functions, by [F4] and [F6].
Divide the epoch into consecutive chunks of length , leaving a shorter final chunk. The entries of are one-dependent functions of adjacent iid anchors, while are iid functions of individual anchors. For a positive semidefinite matrix of norm at most , . Its entrywise second moments here are therefore at most , because the mean matrices in step 5.1 are uniformly bounded. Within each parity the entries of are independent. Apply [F3] separately to both parities and each matrix entry, and to the independent entries of across . For any full chunk , the probability that an entry of either centered covariance sum exceeds a fixed small multiple of is at most . Since for symmetric matrices in dimension , choose the multiplier small relative to and union over at most chunks and their entries. Choosing sufficiently large gives an anchor event of probability at least , for any prescribed , on which for every full chunk, with both sums bounded on both sides. The same entrywise Bernstein calculation gives an upper bound on both sums in the shorter chunk. Here and below the constants may be enlarged independently of and . Since , eventually; if not, the same upper bound handles the entire epoch as one residual chunk.
On a good anchor vector, divide each full-chunk sum of by . Its conditional covariance lies between and , for a fixed , and it is a sum of independent centered atoms of norm at most . Its whole value need not be bounded. Its atom kernels are the measurable pushforwards of under the formula in step 5.1; independence holds conditional on the supplied anchors. Thus [F1], with the ambient length parameter and a fixed atom bound , supplies a measurable coupling of the chunk sums with independent matching Gaussian chunk sums. Undoing the normalization gives maximum chunk-boundary error at most , with conditional failure at most . On bad anchors use independent product coupling instead, preserving both marginals. No theorem about uniformly bounded grouped values is being invoked.
To redistribute covariance without a matrix random clock, start instead with independent standard Gaussian vectors , independent of the anchors, and put . These have the required conditional independent row law even if individual are singular. On good anchors group them into the same full chunks, whose covariances are positive definite. Write for the desired other row. Pad to a dyadic number of chunks with independent identical-covariance Gaussian entries in both rows. All chunk covariances lie between and , and every contiguous whole-chunk interval satisfies by step 5.1's endpoint identity; padded intervals retain that bound. We construct the row from the already given row using the Gaussian tree in the next two steps, leaving the original unchanged.
Fill a Gaussian chunk of total covariance and coupled sum with increments of individual covariances by taking independent and setting . The residual is independent of the original total Gaussian by [F5]; direct covariance calculation shows this is its exact conditional bridge law and its sum is . Thus the completed Gaussian noise row has exactly its conditional independent Gaussian product law. Conditional on good anchors, every within-chunk raw prefix has atom bound and covariance at most . Coordinatewise [F3] and Gaussian tails in [F5], with a union bound over at most prefixes, bound both raw and Gaussian within-chunk norms by , with failure at most after enlarging . Couple a shorter residual chunk independently and use its covariance upper bound in the same calculation. Consequently there is a conditional noise-row coupling , with conditionally independent , of full-prefix error and total failure . Its extension to the actually given follows by [F6] disintegrating the finite joint law given anchors and ; the original innovations are preserved.
At a node of chunks and for either covariance row , write , , and . Thus and have eigenvalues comparable to , and the Gaussian parent and residual are independent by [F5]. Extract their standardized independent Gaussian primitives from the given row, and use the same primitives to construct the root and residuals. For positive definite matrices , the identity , with , yields : solve it by integrating and bound the integral. The endpoint bounds in step 7.2 and the inverse identity imply , , and . A union bound over the Gaussian parents and primitives gives, outside probability , and primitive norms at most . Reconstruct children by and . Their errors satisfy , , with , and the root error is also at most . Every child and root marginal is its exact Gaussian law, so the two completed chunk rows each have their required conditional product law.
The prefix bound uses exact covariance interpolation, not repeated norm bounds on . In any subtree its leading root coefficient is . At a boundary ancestor its residual coefficient is if the boundary is in the left child, and if it is in the right child. For the left case substitute the child's leading term and cancel ; for the right case add the full left error to the partial right error, giving leading coefficient and the stated residual coefficient. This proves the expansion by induction. Each coefficient has norm at most , using the count and eigenvalue bounds on chunk covariances, and there are at most boundary ancestors. Step 8.2 therefore gives prefix error at chunk boundaries. Fill the Gaussian increments by the bridge in step 8.1; the increments were already present. Their within-chunk prefix norms are at most by Gaussian tails and the covariance upper bounds. Handle a short residual chunk independently. This gives full-prefix Gaussian-row error outside probability . On bad anchors use an independent row while retaining the original independent standard . All operations are measurable matrix functions and finite Gaussian kernels.
The couplings of steps 8.1 and 9.1 have the same joint anchor and -Gaussian-row marginal. Disintegrate the second joint law to give its standard row conditional on that common anchor and Gaussian row, and append it to the first coupling. The disintegration formula in [F6] proves that both joint laws, hence both error estimates, are preserved. In particular the original innovations remain unchanged and the row is still independent standard Gaussian conditional on anchors. The two full-prefix errors add. Adding back the conditional means of the triples now gives an iid replacement row since the pairs are iid. Indeed the original mean sum is . The endpoint means have uniformly bounded th moments by step 4.1. Their maximum exceeds with probability at most by [F7]. Therefore the original truncated triple sums and the partial sums of have maximum error with failure at most .
The row is centered and has uniformly bounded th moment: conditional Gaussian moments and coordinatewise [F4] give , and the mean terms satisfy the same bound. Write , , and . The variance of a triple is . Subtracting the covariance of its conditional mean gives by total variance. Since , adding gives This is a direct finite product calculation, independent of convergence in law or an inference from a coupling error. As , , which is positive definite in the active space; hence these iid covariances have fixed upper and positive lower eigenvalue bounds in every sufficiently large epoch.
Step 8.2's square-root bound gives . Use the same and apply Gaussian tails to all its partial sums to bound the covariance correction by with failure at most . This is because . Consequently step 1.2 and this covariance correction give a finite-row iid coupling with matching covariance , error and failure . Apply this result to the varying laws from step 11.1, using their uniform moment and eigenvalue bounds. This uses no fixed-law asymptotic iid theorem and asserts no triangular-law interface absent from a supplier's Statement.
The original normalized active and agree up to their global mean correction unless one of the epoch's observations exceeds , an event of probability at most . The mean correction is at most . Moreover . For the lag-one product difference, use before centering and then include the mean correction. Pointwise the displayed squared tail is at most a constant times the sum of the squared tails of and above , so its expectation is at most . The variance term has the same bound. Since is positive definite in the active space, couple the iid Gaussian triple increments of covariance from step 12.1 to increments of covariance using the same standard normals. Step 8.2's square-root bound and Gaussian tails make the maximum correction at most with failure , because . Thus each sufficiently large epoch of original triple sums can be coupled to independent Gaussian increments of covariance , with error and failure at most Here , and all constants are uniform over epochs.
Retain triple index as a finite initial prefix, coupled arbitrarily, and perform the finite constructions on epochs of triple indices, , using fresh kernels. Each conditional noise coupling and each Gaussian redistribution has the full conditional product marginal given all anchors. Because the original interior pairs are independent conditional on anchors, their finite-epoch kernel products give a product standard law conditional on all anchors. Hence the pairs , and therefore the replacement variables with the deterministic epoch-dependent function in step 10.1, are independent across epochs as well as within them. Each later iid-row coupling is glued through its finite replacement row by [F6]; because these rows are independent, products of their kernels give independent final Gaussian rows across epochs. These statements follow by integrating finite rectangles under the conditional product laws, so no independence is being inferred from an error bound. All anchor sets and row coordinates involved in an epoch are finite standard-Borel products, making these disintegrations legitimate. Apply [F6]'s countable kernel extension with the entire original probability space as its initial coordinate to preserve all original innovations. The failure bounds in step 13.1 are summable over , since . By [F7], all sufficiently large epochs satisfy their error bounds almost surely; finitely many earlier errors can be included in a finite random constant.
Sum the completed epoch errors and the current partial one. The geometric bound gives a global triple-prefix error . Map back from the fixed normalization in step 3.1 and restore the covariance-kernel coboundary from step 2.1. Split each independent Gaussian triple increment into three increments of covariance : for each take a fresh triple , independent across all and of all triple increments, and put . By [F5] these three increments are jointly Gaussian, each has covariance and distinct ones have zero cross covariance, so the entire resulting row is independent with covariance . At nontriple endpoints the original sum changes by at most two , and the Gaussian sum by at most two . Union and Markov bounds on dyadic ranges, using respectively their finite th moments and arbitrary Gaussian moments, show both maxima are , as in step 2.1. This proves the stated full-prefix approximation.
In a fixed orthonormal eigenbasis of , divide each positive-variance Gaussian coordinate by its fixed standard deviation and add independent standard Gaussian coordinates in the zero-variance directions. This yields independent standard normal vector increments realizing by . For each integer interval adjoin independent Brownian copies by [F8] and use their bridges , independently of the existing increments. Interpolate the standard Gaussian random walk linearly on that interval plus its bridge. The paths are continuous and direct covariance calculation gives ; finite Gaussianity and zero cross covariances give independent normal increments by [F5]. Thus [F8] identifies a standard -dimensional Brownian motion, and the fixed eigenbasis map gives at every integer . This construction also applies when by adjoining an independent Brownian motion. It does not alter the original-process marginal or any error already proved.
Source notes
The conditional triple decomposition follows the finite-memory mechanism in Berkes–Liu–Wu. Every coupling theorem used as a premise is proved locally. The first coupling uses independent bounded atoms in a random environment, the covariance redistribution is a direct Gaussian matrix tree, and the changing iid replacement laws are handled by the explicit finite-horizon argument in steps 1.2 and 12.1. No independent-nonidentical strong approximation, matrix random clock, or uniform triangular-law asymptotic theorem is imported.
Strong invariance principle for geometrically ergodic chains
Statement
Let be a geometrically ergodic Markov chain with invariant probability , let be real Borel, and suppose there are with . Suppose also that a state-dependent one-step minorization State-dependent one-step minorization holds. For the Flegal–Jones strong invariance principle convention, assume its asymptotic standard deviation satisfies . Then on a suitable extension of the chain's probability space there are a standard Brownian motion , a finite random constant , and a finite random index such that for all , almost surely.
Facts & Assumptions
Given: The chain, , and positive finite in the Statement. AC is assumed under this item's axiom-use convention. The initial law is arbitrary.
AC is the convention for the Gaussian, path-law, standard-Borel disintegration and probability-space extension suppliers below. The Axiom of Choice
For each specified initial law, an everywhere geometrically ergodic kernel on an arbitrary measurable space has a standard-Borel observation factor with the same initial and stationary real-observation path laws and moments. A joint observation/standard-Borel auxiliary law lifts to an extension of the entire original probability space. A geometrically ergodic chain has an all-start standard-Borel observation factor
On a standard-Borel geometrically ergodic chain, a finite centered observable with stationary th moment has, after an almost surely finite integer delay, actual adjacent-iid full calendar tours: lengths are iid with exponential moment, rewards are centered, and absolute tour envelopes have every moment . The full-tour occupation identity is , where . The same fixed post-delay tour law can be used for the specified and invariant initial laws. Geometric Markov chains admit delayed adjacent-iid split tours with spare moments
Centered adjacent-iid vectors of dimension at most two with finite th moment admit, for every , an original-process-preserving Gaussian partial-sum approximation with uniform error . Its long-run covariance is , and a standard vector Brownian driver is supplied, including singular covariance. Adjacent iid vector functionals admit a finite-moment strong approximation
Centered iid vectors in dimension two with finite th moment, for , have a preserving-extension Gaussian approximation of maximum error , with a standard vector Brownian driver. Singular covariance is allowed, so scalar variables can be embedded as . Finite-moment iid vectors admit a strong Gaussian approximation
If a rate- Poisson clock lies in a sigma-algebra independent of a standard Brownian motion , there is an extended standard Brownian motion independent of that whole sigma-algebra such that almost surely. Independent Brownian motion straightens a Poisson clock
Centered identically distributed independent or one-dependent variables, for , have uniform increments over all windows of length at most , simultaneously for fixed . Finite moments control increments over short moving windows
A finite-window homogeneous Poisson process has, conditional on its total count, independent uniform unordered event positions; its ordered event density with total points in is . The scalar Poisson mass and the convergent exponential series are available. Finite-window homogeneous Poisson construction and ordered event-time density, The scalar Poisson count law, The real exponential function and the number by a power series, The exponential addition formula
Real sequences are standard Borel through their countable Polish product presentation. Standard-Borel joint laws have conditional probability kernels, and supplied kernels extend an arbitrary original probability space without changing its marginal. The real line with its Borel sigma-algebra is standard Borel, Under the Axiom of Countable Choice, countable products and subspaces of Polish spaces are Polish, Disintegration of a joint law on standard borel spaces, Ionescu-Tulcea construction of a Markov chain
Markov's inequality, countable subadditivity and the first Borel–Cantelli lemma give the usual union and summability bounds. Markov's inequality for random variables, Finite and countable subadditivity of measures, First Borel-Cantelli lemma for events
Brownian motion has continuous centered Gaussian paths with covariance . Uncorrelated jointly Gaussian blocks are independent. Independent Brownian copies exist. Reflection and the normal Mills upper bound give when . Brownian motion, Uncorrelated jointly Gaussian blocks are independent, Existence of continuous Brownian motion, Assuming countable and dependent choice, countable products of arbitrary probability spaces, Law of the Brownian maximum, Two-sided Mills bounds for the standard normal tail
An everywhere finite common geometric total-variation envelope admits an integrable geometric envelope, which implies exponential stationary strong mixing. A measurable observation does not increase those coefficients. Harris geometric ergodicity admits an integrable total-variation envelope, A total-variation rate bounds strong-mixing coefficients, A measurable observation does not increase stationary mixing coefficients
A centered stationary sequence with has an absolutely convergent covariance sum and a CLT with that variance, including zero. A perturbation tending to zero in probability preserves a distributional limit; centered normal characteristic functions identify their variances. Strong-mixing CLTs under bounded, moment, and logarithmic-moment conditions, Converging together lemma, Characteristic function of a normal law, Long-run variance for a stationary Markov chain
Proof
Put , , and . The given moment makes , since , and the finite-sum power inequality gives . Choose and set . To obtain exactly the requested observation indices, apply [F1] to the chain starting at , with its actual initial law, and observable . This gives a standard-Borel factor retaining the full path , the stationary path and its moments, and an everywhere finite common geometric bound. Work on this factor until the original-space lift at the end. Apply [F2] there, retaining its finite prefix and integer delay . Define to be the sum of the first centered observations of the retained tail beginning at calendar position , with the empty sum zero. This real-time tail process has cycle lengths , centered rewards , envelopes , and relative boundaries , . Then , , and . The pair is a centered adjacent-iid vector in ; exponential length moments and verify its two coordinate moments. No iid assertion about the full tours is used.
Apply [F3] jointly to . Write its long-run covariance as . Since the lengths are iid, ; the reward and cross entries still include their actual lag-one terms. The Gaussian approximation gives a continuous centered Gaussian vector process of covariance such that and . For now suppose . Put and . Then is standard Brownian. If , is standard Brownian independent of the entire process : each finite stack is jointly Gaussian and every cross-covariance is , so [F10] gives independence on finite rational-time cylinders and continuity extends it to all times. If , the residual Gaussian process is zero on rational times almost surely and hence identically zero; adjoin an independent unused . Thus, with independent , we have and . This regression uses the full , not just the covariance inside a single cycle.
Define and , so and . Choose any fixed and . Apply [F4] to the centered iid scalar Poisson variables , embedded as , obtaining a joint law with a standard Brownian driver and . Their moments of every order follow directly from the Poisson exponential series in [F7]. By [F8], disintegrate this joint law to a kernel for the entire count sequence given the integer driver sequence, and evaluate it at the already present . This preserves all original variables and reproduces the asserted count/driver joint law. The kernel depends only on and is sampled using new randomness independent of everything earlier, so remains independent of the joint pair . Conditional on , insert independent uniform positions in , independently over unit intervals and of earlier variables given the counts. By [F7] these are the exact finite-window Poisson laws. More explicitly, partitioning a unit interval into lengths gives multinomial conditional counts; summing their mass against yields the product . Hence the resulting locally finite simple right-continuous clock has independent Poisson increments of rate , , and remains independent of jointly with .
Two uniform bounds control the remaining calendar errors. By the common th envelope moment and [F9], , summable in . Intersecting Borel–Cantelli events over positive rational and integer yields without any independence requirement. For any standard Brownian motion and fixed , cover by overlapping intervals of length . Every pair of times at distance at most lies in one such interval, and [F10] bounds the oscillation tail at by . There are at most intervals. Taking large and using [F9] gives the uniform modulus for windows of length ; for fixed window length the same proof gives . All fixed integer can be placed on one event. Finally [F6] applies to the centered one-dependent rewards , so their partial-sum increments over at most consecutive cycle indices, among the first indices, are uniformly . A finite random bound for or is covered by their countably many integer choices.
We record the uniform estimates used for this clock and for Brownian motion. Summing the Poisson mass series gives . The exponential series and for give for . Exponential Markov with parameter and either sign therefore bounds the increment deviation tail by . On dyadic horizons , take and union over the at most integer intervals in . A sufficiently large deterministic makes the probabilities summable, so [F9] proves this bound for all such intervals eventually almost surely, for every positive integer by countable intersection. In particular unit counts are uniformly, intervals of integer length at least have count increase at least times their length eventually, and . Also [F10], a union bound and [F9] make all within-unit Brownian oscillations up to time be . These estimates require no independence between evaluation times and the Brownian motion. Consequently the integer coupling in step 3.1 extends to , and step 2.1 gives .
Let and . Brownian reflection on dyadic horizons gives ; thus step 2.1 implies . Together with , this gives , absorbing their finite initial values. Put . Monotonicity gives , except for an integer endpoint convention which is covered by shifting at most two indices. The unit-count bound in step 4.1 therefore gives uniformly for . On every required range of indices, eventually by its little-o bound. Choose an integer with . The lower-increment bound in step 4.1 makes a change of indices change by more than eventually. The corresponding values bracket , using index zero at the lower boundary. Thus . The countable estimates for all deterministic horizon multipliers cover a finite random multiplier as well. This is a direct monotone inverse comparison, not an imported renewal approximation.
Let be the interarrival times of . Their iid exponential law follows from [F7]: for , sum its ordered density over all total counts and integrate the last ordered points above . The simplex volume is , so the first- density is . Increasing covers the whole ordered simplex. The spacing change of variables has determinant one and gives on positive spacings. Hence , and all moments are finite, by elementary exponential integration and integration by parts for integer moments, bounding other positive powers by the next integer power plus one. Apply [F4] to with the same auxiliary exponent , and glue its Brownian driver to the entire already realized arrival sequence by [F8], using a fresh kernel depending only on that sequence. Then . Since , set and use step 4.1's within-unit Brownian bound to obtain . The clock sigma-algebra generated by and their fresh clock randomness remains independent of : each new kernel depended on or on the clock alone and used randomness independent of earlier variables.
Apply [F5] to , and this entire clock sigma-algebra. It supplies a standard Brownian motion independent of that sigma-algebra, hence independent of , and . Define . It is continuous centered Gaussian with independent stationary increments and variance per unit time . Thus, when , for a standard Brownian and . When , positive semidefiniteness forces , so and adjoin an unused standard . The case needs no residual normalization or nonzero Brownian coefficient.
For , the following exact decomposition has eight terms: , where , , , , , , , and . Adding these expressions cancels every intermediate reward, clock and Gaussian evaluation and leaves exactly . This verifies the signs and deterministic drift, including the bounded last term.
Uniformly for , the first term is by the incomplete-cycle bound in [F2] and step 3.2. The second is by the inverse comparison in step 5.1 and [F6]. The third is by step 2.1 at , the bounded fractional-part drift , and the within-unit Brownian bound. The fourth is by the length/Poisson comparison at and the count increase between that point and . The fifth is by step 6.1, and the sixth is by step 7.1. The strict inverse of the simple Poisson clock satisfies , so its distance from is at most one; step 3.2's fixed-window Brownian bound makes the seventh term , while . The indices are by step 5.1, so every approximation was evaluated inside its verified range. Therefore for .
If , the iid lengths have variance zero, hence for every simultaneously almost surely. In particular ; positive semidefiniteness of implies . The joint approximation [F3] gives uniformly for , with a standard Brownian when , and an unused one when . For , , so step 3.2's fixed-window modulus replaces by at error . The incomplete-tour bound is again . Thus . The process is standard Brownian by its continuous Gaussian covariance and increments, and the coefficient is the same . This proves the tail approximation in every covariance case, including a zero residual or zero full Gaussian covariance, without inverting a singular matrix.
Restore the finite initial segment of [F2]. Let be its finite centered reward and its finite integer length. For , the actual sum is . Keep the same standard Brownian motion and use . The fixed-window Brownian estimates in step 3.2 hold on one event for every positive integer window length. On that event choose an integer at least the finite random ; then , regardless of any dependence of and . Finitely many earlier indices contribute a finite maximum. Steps 9.1 and 10.1 therefore give the complete-path bound in the factor construction. Neither a moment of the initial law nor an assertion that a randomly shifted Brownian path is standard has been used.
Identify without making an unsupported second-moment inference from the almost-sure approximation. Use the same factor and fixed post-delay tour law under invariant initial law, as supplied in [F1]–[F2]. Steps 2.1, 3.1, 3.2, 4.1, 5.1, 6.1, 7.1, 8.1, 9.1, 10.1, 11.1 then give a CLT with variance : the approximation error divided by tends to zero almost surely since , and has exactly the centered normal law of variance for each ; apply [F12]'s converging-together result. Independently, [F11] gives exponential stationary strong mixing of the original observations. Their th moment, with , satisfies [F12]'s moment regime because every positive power of an exponential rate is summable. Hence [F12] gives the stationary CLT with its absolutely convergent covariance-sum variance , the asymptotic variance in the Statement's convention. The stationary observation-path equality in [F1] makes these two limits refer to the same law. Uniqueness of centered normal characteristic functions, , gives . The assumed positivity therefore gives . This argument uses no uniform-integrability assumption on the approximation error.
Finally preserve the actual original chain probability space. The factor construction gives a joint law of its real observation sequence and the rational-time coordinates of the final standard Brownian motion . Both are standard-Borel real sequence spaces by [F8]. Use [F1]'s lifting assertion, or disintegrate this joint law and apply its conditional kernel at the actual sequence by [F8]; this preserves the entire original probability space and reproduces that observation/Gaussian-coordinate joint law. The approximation event uses only these observations and integer driver coordinates, so it still has probability one. The rational coordinates extend to a continuous Brownian path on a common measurable full-measure event: for each positive integers , some integer makes all rational pairs in at distance below have coordinate difference below . This is a Borel condition and has probability one under the original Brownian rational law; completeness of gives a unique consistent continuous extension on compact intervals. Every fixed real-time evaluation is the limit along a deterministic rational sequence, hence measurable; on the null failure event set the path to zero. Its finite-dimensional laws agree with standard Brownian laws by continuity, so [F10] identifies genuine standard Brownian motion on the extension. With , steps 11.1, 12.1 give error , and implies . Choosing the finite eventual bound and starting index in this almost-sure comparison gives precisely the in the Statement. AC was used through the stated local Gaussian, factor, tour and extension suppliers; the full moment range, arbitrary initial law, positive-variance convention and original observation indices are retained.
Source notes
Flegal–Jones Lemma 3 supplies the stated rate. Its cited older scalar regenerative route does not itself provide the unrestricted moment-range proof. The argument above instead uses the fully local joint adjacent-iid approximation, delayed multistep tour construction and independent Brownian/Poisson composition. The latter's source manuscript statement and roadmap, and the complete almost-sure proof of Theorem 1 in Bashtova–Shashkin's heavy-tail paper motivating the calendar-time decomposition, were retrieved and inspected as recorded in the owner research notes; no external theorem replaces a proof premise here. The zero-variance rate is proved separately in Zero-variance geometric Markov sums have a uniform strong rate; the Statement here retains its explicit positive-variance hypothesis.
Strong consistency of overlapping batch means
Statement
Let be a geometrically ergodic Markov chain with invariant probability , and let the initial distribution be arbitrary. Let be real Borel and suppose for some . Assume the state-dependent one-step minorization State-dependent one-step minorization: there are measurable and a probability measure such that for all .
Let be positive integer batch sizes with for all sufficiently large , , and . Suppose there is such that and suppose there are and with for every , and
For and , define The strong-consistency conclusion is where is the covariance long-run variance defined by Long-run variance for a stationary Markov chain. The batch sequence satisfies Flegal–Jones overlapping-batch size conditions for the displayed exponent ; these conditions are written out above. This corrects the source's incompatible requirement that both integer and be nondecreasing, retaining every numerical batch-size hypothesis.
As in the intended source corollary, if and , then satisfies these corrected batch-size hypotheses and has the same strong-consistency conclusion. No monotonicity of is asserted.
Facts & Assumptions
Given: The geometrically ergodic chain, its arbitrary initial law, the real measurable observable , the moment exponents , the supplied one-step minorization, and every displayed numerical batch-size hypothesis.
AC is assumed through the chain, Brownian, mixing and strong-approximation interfaces. The Axiom of Choice
The canonical chain exists for each initial law, and choosing its invariant initial law makes it strictly stationary. Canonical Markov chain on path space, Invariant initial law makes a Markov chain stationary
An invariant kernel with an everywhere finite measurable common geometric total-variation bound admits a geometric envelope integrable under , on the original arbitrary measurable space. Harris geometric ergodicity admits an integrable total-variation envelope
An integrable geometric total-variation envelope gives exponentially decaying stationary strong-mixing coefficients. A measurable observation preserves stationarity and does not increase those coefficients. A total-variation rate bounds strong-mixing coefficients, A measurable observation does not increase stationary mixing coefficients
For a centered stationary sequence with a finite th absolute moment, , summability of gives absolute covariance summability, a finite nonnegative covariance-sum variance, and a CLT with that variance, including zero. Strong-mixing CLTs under bounded, moment, and logarithmic-moment conditions
The covariance long-run variance is the stationary covariance sum when absolute summability and the CLT have been established; equality to the CLT variance is a theorem obligation, not a consequence of the definition alone. Long-run variance for a stationary Markov chain
Under the stated geometric-ergodicity, minorization and th moment hypotheses, with positive finite asymptotic standard deviation, the strong invariance principle supplies a Brownian approximation with error on an extension of the original chain. Strong invariance principle for geometrically ergodic chains
Under the everywhere common geometric total-variation bound, a finite centered with zero covariance sum has, for every , the maximal partial-sum rate under every initial law, without a moment condition on that law. Zero-variance geometric Markov sums have a uniform strong rate
A partial-sum approximation by with error , and , implies consistency of the exact OBM estimator when , , , and . It requires neither stationarity nor batch monotonicity and includes zero variance. Strong Brownian approximation implies overlapping-batch consistency
The corrected Flegal–Jones contract retains all numerical rate, summability and endpoint conditions. For and , the integer sequence satisfies that contract; the definition supplies the floor bounds and each numerical limit. Flegal–Jones overlapping-batch size conditions
A standard scalar Brownian motion exists, and the product of its probability space with the original probability space gives an extension preserving the original marginal. Existence and scaling of -dimensional Brownian motion, For sigma-finite factors, the product measure exists, has the rectangle formula, is sigma-finite, and is unique
Proof
Put , and . The moment hypothesis gives , since , so is finite. The elementary bound gives , and . For variance identification consider the canonical chain with initial law , which is stationary by [F1]. The given geometric bound has an everywhere finite measurable common envelope, so the stronger arbitrary-measurable-space clause of [F2] supplies an integrable such envelope. By [F3], this stationary chain and its observation have a strong-mixing bound with and . This construction is only for identifying the stationary covariance parameter; the observed initial law remains arbitrary.
Take . The geometric bound in step 1.1 makes finite. Apply [F4] to the stationary observed sequence. Its covariance series is absolutely convergent, its sum is finite and nonnegative, and its CLT has precisely variance . Thus [F5] identifies , without using the stationary Harris-only formulation of a separate MCMC regime theorem. Write and .
If , the original , the given , geometric ergodicity and supplied minorization meet [F6] with the finite positive standard deviation from step 2.1. Hence, for the original arbitrary initial law, on an extension there are a standard Brownian motion and finite random constants for which almost surely. The same moment exponents and the same covariance variance are used; no stronger moment assumption is introduced.
If , choose a real with , possible because . Apply the all-start geometric-envelope clause of [F7] to ; its moment and zero covariance hypotheses were established in steps 1.1, 2.1. Use the chain indexed from as its time-zero chain, whose initial law is still arbitrary. Then almost surely. Since , this implies . Adjoin a standard Brownian motion on the product extension supplied by [F10] and set . The same approximation as in step 3.1 now holds, with zero Brownian coefficient. No division by or implication from a degenerate CLT to an almost-sure rate is used.
In either variance branch the hypotheses of [F8] hold with this . Indeed gives , while the summability and assumptions are displayed in the Statement. Moreover . The other stated numerical endpoint and growth conditions are retained and impose no additional requirement in this reduction. Thus [F8] gives convergence to for the exact OBM estimator applied to .
Centering the observations by the constant does not change this finite estimator: and , so their difference equals . The indexing in [F8] is exactly after adding one to the starting index, and its denominator is as here. Therefore its estimator is the displayed , and steps 2.1 and 5.1 prove its almost-sure convergence to under the original arbitrary initial law.
Finally, the source-corollary range is . By [F9], satisfies every corrected numerical batch hypothesis, so steps 1.1, 2.1, 3.1, 4.1, 5.1, 6.1 apply to it as well. This verifies the intended floor-power conclusion without asserting the incompatible ratio monotonicity. The proof includes , finite early inadmissible batch indices, and arbitrary initial laws without additional moment hypotheses, as allowed by [F7]–[F8].
5 · Examples, counterexamples and false statements
None yet.
Sources
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Theorem 2
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §2
- Gareth O. Roberts and Jeffrey S. Rosenthal, General State Space Markov Chains and MCMC Algorithms, §§3–4
- Charles J. Geyer, Markov Chain Monte Carlo Lecture Notes, §§4.8–4.9
- Gareth O. Roberts and Jeffrey S. Rosenthal, General State Space Markov Chains and MCMC Algorithms, §5
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo
- Olle Häggström and Jeffrey S. Rosenthal, On Variance Conditions for Markov Chain CLTs
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §§2–4
- Charles J. Geyer, Markov Chain Monte Carlo Lecture Notes, §4.8.2
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §2 and Theorem 9
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Assumption 3 and Theorem 2
- Charles J. Geyer, Markov Chain Monte Carlo Lecture Notes, §§2.2 and 4.8
- Meyn and Tweedie, Markov Chains and Stochastic Stability, measurable kernels and geometric-ergodicity reductions
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Lemma 10 (Whittle's quadratic-form bound)
- Mikhail Gordin and Magda Peligrad, On the functional central limit theorem via martingale approximation, §1 and full martingale construction in §2; the absolutely summable specialization is derived locally
- Elena Bashtova and Alexey Shashkin, Strong Gaussian approximation for cumulative processes with heavy tails, Lemma 5 and its proof
- Florence Merlevède and Emmanuel Rio, Strong approximation for additive functionals of geometrically ergodic Markov chains, Lemma 2.4
- A. Yu. Zaitsev, Estimates of quantiles of smooth conditional distributions and the multidimensional invariance principle, Siberian Mathematical Journal 37 (1996), §§1–2; Gaussian convolution specialization proved here
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, §3 and Theorem 2
- Charles J. Geyer, Markov Chain Monte Carlo Lecture Notes, §4.9.2
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §3
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §§2 and 4.1
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §4.1
- Charles J. Geyer, Markov Chain Monte Carlo Lecture Notes, §4.8
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §4
- Meyn and Tweedie, Markov Chains and Stochastic Stability, Chapters 13 and 15, regeneration and local geometric convergence
- P. Doukhan, P. Massart and E. Rio, The functional central limit theorem for strongly mixing processes, Lemma 1 and Proposition 1
- Mikhail Gordin and Magda Peligrad, On the functional central limit theorem via martingale approximation, stationary realizations in §1 and rho coefficients in §3.3
- A. Yu. Zaitsev, Estimates of quantiles of smooth conditional distributions and the multidimensional invariance principle, Siberian Mathematical Journal 37 (1996), Lemmas 1.1, 1.2, 2.1, and 4.1; Gaussian convolution specialization proved here
- Meyn and Tweedie, Markov Chains and Stochastic Stability, measurable kernels and geometric ergodicity; observation-factor reduction supplied locally
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorem 2, part 3
- Gareth O. Roberts and Jeffrey S. Rosenthal, Geometric ergodicity and hybrid Markov chains, Theorem 2 (published Theorem 2.1)
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorem 7
- Mikhail Gordin and Magda Peligrad, On the functional central limit theorem via martingale approximation, §§1–2 and §3.3
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Appendix B.3 and Lemma 4
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorem 5(1)
- P. Doukhan, P. Massart and E. Rio, The functional central limit theorem for strongly mixing processes, bounded-variable application
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §§3–4
- A. Yu. Zaitsev, Estimates of quantiles of smooth conditional distributions and the multidimensional invariance principle, Siberian Mathematical Journal 37 (1996), Lemmas 6.1 and 8.1; Gaussian anchor specialization and measurable weighted scalar bridges proved here
- A. Yu. Zaitsev, Multidimensional version of the results of Komlos, Major and Tusnady for vectors with finite exponential moments, Sections 3 and 5
- Galin L. Jones, On the Markov Chain Central Limit Theorem, §2, Remark 1
- Sean P. Meyn and Richard L. Tweedie, Markov Chains and Stochastic Stability, second edition, Theorems 14.3.7, 15.3.3, and 15.4.1
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorems 5 and 6
- P. Doukhan, P. Massart and E. Rio, The functional central limit theorem for strongly mixing processes, Lemma 1, Theorem 1 and applications
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorem 8
- Mikhail Gordin and Magda Peligrad, On the functional central limit theorem via martingale approximation, §§1–2 and §3.3; finite-partition phi-to-rho bound proved locally
- A. Yu. Zaitsev, Multidimensional version of the results of Komlos, Major and Tusnady for vectors with finite exponential moments (1998), Sections 5–8; weaker logarithmic bound by Gaussian anchors
- Meyn and Tweedie, Markov Chains and Stochastic Stability, splitting and regeneration; full multistep observation-word construction supplied locally
- Bednorz and Latuszynski, A few remarks on Fixed-Width Output Analysis for Markov Chain Monte Carlo, tour-moment argument
- A. Yu. Zaitsev, Multidimensional version of the results of Komlos, Major and Tusnady for vectors with finite exponential moments (1998), Section 5; heterogeneous covariance interpolation and Gaussian anchor specialization proved here
- Galin L. Jones, On the Markov Chain Central Limit Theorem, Theorem 9 and Corollaries 1–5
- A. Yu. Zaitsev, Multidimensional version of the results of Komlós, Major and Tusnády for vectors with finite exponential moments; finite-moment dyadic truncation reduction
- Istvan Berkes, Weidong Liu and Wei Biao Wu, Komlos-Major-Tusnady approximation under dependence, Annals of Probability 42 (2014), 794–817, Sections 3.2–3.4; adjacent iid specialization and matrix Gaussian redistribution proved locally
- A. Yu. Zaitsev, Multidimensional version of the results of Komlos, Major and Tusnady for vectors with finite exponential moments, Sections 5–8; finite-horizon truncation from the locally proved bounded coupling
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Lemma 3
- Florence Merlevède and Emmanuel Rio, Strong approximation for additive functionals of geometrically ergodic Markov chains, Lemma 2.4 and proof roadmap
- Elena Bashtova and Alexey Shashkin, Strong Gaussian approximation for cumulative processes with heavy tails, Theorem 1 and calendar-time decomposition
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Theorem 2 and Appendix B.3
- James M. Flegal and Galin L. Jones, Batch means and spectral variance estimators in Markov chain Monte Carlo, Corollary 1