Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-adaptedPipeline-generatedprecheck pass
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.

Self-normalized importance sampling ratio CLT and consistent Monte Carlo standard error

Statement

Let (X,X,μ), p~, Zp, π, q, Q, f, and Wi satisfy the hypotheses and definitions of Self-normalized importance sampling is strongly consistent under finite first moments. Assume also EQ[W12]<∞,EQ[(W1f(X1))2]<∞, and assume AC for the central limit conclusion. For every n≥1, put Sn=∑i=1nWi and define the numerical extension θ^SN,n=∑i=1nWif(Xi)/Sn when Sn>0, and θ^SN,n=0 when Sn=0. This replaces the algorithm's no-estimate status by 0 for all convergence statements below. Set θ=∫f dπ and τ2:=1Zp2EQ ⁣[W12(f(X1)−θ)2]. This is finite under the stated second moments. Then this numerical extension of the estimate in Self-normalized importance sampling satisfies n(θ^SN,n−θ)⇒N(0,τ2). Define its plug-in variance estimate, when Sn>0, by Var⁡^MC(θ^SN,n):=∑i=1nWi2(f(Xi)−θ^SN,n)2Sn2, and set this quantity to 0 if Sn=0. Then, without using AC, nVar⁡^MC(θ^SN,n)⟶τ2almost surely. Define se⁡^MC,n to be the nonnegative square root of this plug-in variance. Then n se⁡^MC,n→τ almost surely. No finite-sample unbiasedness assertion is made.

Facts & Assumptions

Given: The IID sampling and first-moment hypotheses of Self-normalized importance sampling is strongly consistent under finite first moments, the additional second moments, and AC for the CLT branch.

[F1]

Self-normalized importance sampling: The algorithm defines a measurable weight w and, for Sn>0, returns θ^SN,n=∑iWif(Xi)/Sn; when Sn=0, it returns the status no estimate: zero observed target weight.

[F2]

Change of variables for expectation: A nonnegative measurable function of an IID draw with law Q has expectation given by integration against Q, possibly with value +∞; the same identity holds as a finite integral for integrable functions.

[F3]

Multivariate iid central limit theorem: Under AC, IID Rd-valued variables with finite second moment and covariance Σ have sample-mean limit n(Xˉn−m)⇒Nd(0,Σ); singular Σ is allowed.

[F4]

First- and second-order delta method for random vectors: If n(Tn−η)⇒Z and g(η+h)−g(η)=Ah+o(∥h∥), then n(g(Tn)−g(η))⇒AZ.

[F5]

Self-normalized importance sampling is strongly consistent under finite first moments: Under its support and finite first-moment assumptions, the ratio estimate converges almost surely to θ, and its denominator is eventually positive almost surely. For the convergence statement, that theorem interprets the algorithm's no-estimate status as the numeric value 0 on any finite prefix with zero observed target weight.

[F6]

Kolmogorov iid l1 strong law: Every IID integrable real sequence has sample mean converging almost surely to its expectation.

[F7]

Measurable coordinatewise functions preserve independence: Applying one measurable function to each member of an independent family preserves independence.

[F8]

Identical distribution and IID families: A family is IID when it is mutually independent and has a common law.

[F9]

Moments, variance, and covariance on a probability space: Covariance is the expectation of the product of centered square-integrable variables, and variance is the self-covariance.

[F10]

Multivariate normal law, including singular covariance: Under AC, Z∼Nd(0,Σ) exactly when each linear projection u⋅Z is normal with variance uTΣu; singular covariance is allowed.

[F11]

Holder's inequality for random variables: For conjugate exponents p=q=2, square-integrable real random variables R,S satisfy E∣RS∣≤(E∣R∣2)1/2(E∣S∣2)1/2, so RS is integrable.

[F12]

Arithmetic and lattice operations preserve measurability whenever they are defined: Finite products and sums, and absolute values, of measurable real functions remain measurable.

[F13]

A nonnegative integral over a null set vanishes: A nonnegative measurable function has integral zero over a measurable null set.

[A1]

The Axiom of Choice: AC asserts that every family of nonempty sets has a choice function; the CLT branch inherits it through [F3] and [F10], while the strong-law standard-error branch is choice-free.

[F14]

Integrating against a density agrees with integrating the product: for a nonnegative measurable density r and nonnegative measurable h, ∫h d(r dμ)=∫hr dμ, including infinite values. Integrable real and complex functions, and their integrals: an absolutely integrable real function has integral equal to the finite difference of the integrals of its positive and negative parts.

Proof

Given: Let Zi=(Wi,Wif(Xi)), Wˉn=n−1∑iWi, and Vˉn=n−1∑iWif(Xi).

1.1A1F1F2F3F7F8F12F13givenF14

Use [F14] with density r=q for integration against Q and with r=p~/Zp for integration against π. By the support condition, π({q=0})=0. For each i, Xi has law Q, so the given finite first moment applies to Wi∣f(Xi)∣. The integrable case of [F2] gives EQ[Wi∣f(Xi)∣]=∫{q>0}p~(x)∣f(x)∣q(x)q(x) dμ(x)=∫{q>0}p~∣f∣ dμ<∞. By [F13], ∫{q=0}∣f∣ dπ=0, and hence ∫{q=0}p~∣f∣ dμ=Zp∫{q=0}∣f∣ dπ=0. Thus ∫Xp~∣f∣ dμ<∞. For any real measurable h whose absolute integral is finite against r dμ, apply [F14] separately to h+=max⁡(h,0) and h−=max⁡(−h,0): both product integrals are finite because h±≤∣h∣. Subtracting their finite identities yields ∫h d(r dμ)=∫hr dμ, since (hr)±=h±r for r≥0 (with the density convention 0⋅∞=0). The integrable case of [F2] applied to Wif(Xi) yields EQ[Wif(Xi)]=∫{q>0}p~f dμ=∫Xp~f dμ=Zpθ, because the omitted signed integral is zero by its absolute integrability on the π-null set. The nonnegative case of [F2] also gives EQWi=∫{q>0}p~ dμ=Zp, since ∫{q=0}p~ dμ=Zpπ({q=0})=0. The second-moment assumptions make Zi square-integrable. The map x↦(w(x),w(x)f(x)) is measurable by [F1] and [F12]. Applying [F7] to this map proves that the pairs Zi are independent; their laws agree because each Xi has law Q and the same map is applied, so they are IID by [F8]. Therefore [F3] and [A1] give n((Wˉn,Vˉn)−(Zp,Zpθ))⇒N2(0,Σ), where Σ=Cov⁡Q(Z1).

2.1F5F6F7F8F11F12givenstep 1.1algebra

Put Ui=Wi2(f(Xi)−θ)2, Ci=Wi2(f(Xi)−θ), and Bi=Wi2. Their measurability follows from [F12]. The stated second moments and (Wif(Xi)−θWi)2≤2(Wif(Xi))2+2θ2Wi2 give EUi<∞ and EBi<∞; [F11] gives E∣Ci∣≤(EUi)1/2(EBi)1/2<∞. Each transformed sequence is IID by [F7] and [F8], so [F6] applies to all three. Also (Wi) is IID and integrable, with mean Zp by step 1.1; applying [F6] to this sequence gives Wˉn→Zp>0 almost surely. Separately, [F5] gives θ^SN,n→θ almost surely.

2.2F1F4step 1.1givenalgebra

Define the Borel function g(s,t)=t/s for s>0 and g(s,t)=0 for s≤0; it is Borel because the ratio is continuous on the open set s>0 and the other branch is constant on its Borel complement. At η=(Zp,Zpθ), its derivative is A(h1,h2)=(−θh1+h2)/Zp: for h1 small enough that Zp+h1>0, g(η+h)−g(η)−Ah=−h1(h2−θh1)Zp(Zp+h1)=O(∥h∥2). Since Zp>0, this is a valid first-order expansion. Since Wˉn≥0, g(Wˉn,Vˉn) equals the ratio when Sn>0; when Sn=0, it equals the numerical value 0 prescribed by the Statement. Thus g(Wˉn,Vˉn)=θ^SN,n for every n, including the zero-weight branch where [F1] returns a status. Apply [F4] to step 1.1 to obtain the scalar limit AZ, where Z∼N2(0,Σ).

3.1F2F9F10step 1.1step 2.2algebra

By [F10], AZ is centered normal with variance AΣAT. The centered linear combination of Z1 defining this variance is Zp−1(W1f(X1)−θW1), whose mean is zero by [F2]. Therefore AΣAT=Zp−2EQ[W12(f(X1)−θ)2]=τ2, which proves the stated CLT.

4.1F5F6F7F8F11F12step 2.1givenalgebra∎

On the probability-one event from step 2.1 the denominator is eventually positive, and expansion of the squared residual gives nVar⁡^MC=Uˉn−2(θ^SN,n−θ)Cˉn+(θ^SN,n−θ)2BˉnWˉn2. The numerator tends to EU1 and the denominator tends to Zp2, so the limit is Zp−2EU1=τ2. Continuity of the nonnegative square-root map gives n se⁡^MC,n→τ. This argument uses only the choice-free suppliers [F5], [F6], [F7], [F8], [F11], and [F12], not AC.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

71 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources