Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-adaptedPipeline-generatedprecheck passjudge pass (gpt-6-luna)
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.

Kriging systems and mean-square optimality

Statement

Use the notation of Simple, ordinary, and universal kriging. For simple kriging, if Σ is positive definite, the unique mean-square minimizing weight vector is λSK=Σ−1c0, and its mean-square prediction error is v0−c0⊤Σ−1c0.

For an unknown trend with design matrix P and target vector p0, λ⊤Z is unbiased for every β∈Rq if and only if P⊤λ=p0. Suppose Σ is positive definite and P has full column rank (that is, Py=0 implies y=0 for y∈Rq). Then the unique unbiased mean-square minimizing weights λ∈Rn and Lagrange multiplier ν∈Rq solve (ΣPP⊤0)(λν)=(c0p0). The augmented matrix is invertible under these hypotheses. The resulting mean-square prediction error is v0−λ⊤c0−ν⊤p0. For ordinary kriging, take P=1n and p0=1, giving (Σ1n1n⊤0)(λν)=(c01).

Facts & Assumptions

Given: A square-integrable target Y and observation vector Z from the definition, with covariance matrix Σ, covariance vector c0, and target variance v0.

[F1]

Simple, ordinary, and universal kriging: simple prediction centers Z at its known mean and the target at its known mean.

[F2]

Simple, ordinary, and universal kriging: in a finite trend model, unbiasedness for every β is constrained by P⊤λ=p0.

[F3]

Normal equations for best affine L2 prediction: the best affine L2 predictor exists as a unique almost-sure class, and its centered coefficients satisfy the covariance normal equations.

[F4]

Moments, variance, and covariance on a probability space and Variance and covariance identities for random variables: variance is the expectation of a centered square, and covariance is symmetric and bilinear on finite linear combinations of square-integrable variables.

Proof

Given: The finite second-moment model, Σ≻0, and, for the unknown-trend case, full column rank of P.

1.1givenF1F4algebra

For a centered error vector ε=Z−EZ and centered target error ε0=Y−EY, every unbiased linear predictor has error ε0−λ⊤ε. Expanding its squared expectation by [F4] gives the common objective Q(λ)=v0−2λ⊤c0+λ⊤Σλ. The centered predictor in the simple case has this same error expression.

1.2F2algebra

In the trend model, E[λ⊤Z−Y]=(P⊤λ−p0)⊤β. This expression is zero for every β∈Rq if P⊤λ=p0. Conversely, if it is zero for every β, taking β successively to be each coordinate vector shows every coordinate of P⊤λ−p0 is zero. Thus the unbiased weights are exactly the affine constraint P⊤λ=p0.

1.3givenalgebra

Let M=(ΣPP⊤0) and suppose M(xy)=0. The lower block gives P⊤x=0. Taking the inner product of x with the upper block gives x⊤Σx+x⊤Py=0; the second term is (P⊤x)⊤y=0, so positive definiteness gives x=0. The upper block now gives Py=0, and full column rank gives y=0. Hence the finite square matrix M is injective and therefore invertible.

2.1givenF3step 1.1algebra

In the simple case, [F3] gives the covariance normal equations Σλ=c0 for the minimizing centered affine predictor. Positive definiteness makes the solution unique, λ=Σ−1c0; substituting it into step 1.1 gives Q(λ)=v0−c0⊤Σ−1c0.

2.2step 1.1step 1.2step 1.3algebra

Invertibility gives a unique solution (λ,ν) of the displayed augmented system, and its lower block makes λ feasible. For any other feasible vector λ+h, one has P⊤h=0. Using the upper block Σλ+Pν=c0 in the quadratic from step 1.1 gives Q(λ+h)−Q(λ)=2h⊤(Σλ−c0)+h⊤Σh=h⊤Σh. This is positive for nonzero h and zero for h=0, so λ is the unique constrained minimizer.

3.1step 1.1step 1.2step 2.2algebra∎

Multiplying Σλ+Pν=c0 on the left by λ⊤ and using P⊤λ=p0 yields λ⊤Σλ+ν⊤p0=λ⊤c0. Substitution in Q(λ) from step 1.1 gives Q(λ)=v0−λ⊤c0−ν⊤p0, as claimed.

Depends on

Used by

Dependency tree · two levels

17 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