Documentation

TauCeti.Probability.Distributions.Gaussian.MixedMoments

Fourth mixed moments of Gaussian laws #

The fourth centered moment of four linear observables is the sum of the three products of their pairwise covariances. The formula allows repeated observables and singular covariance. For a multivariate Gaussian this gives the coordinate identity E[Zᵢ Zⱼ Zₖ Zₗ] = Sᵢⱼ Sₖₗ + Sᵢₖ Sⱼₗ + Sᵢₗ Sⱼₖ for centered coordinates.

The linear-observable result is proved for Gaussian measures on real Banach spaces, then specialized to the multivariate law. Finite products of affine linear observables are integrable. Polarizing the existing scalar fourth-moment formula avoids any invertibility requirement.

References #

theorem ProbabilityTheory.IsGaussian.integrable_prod_dual {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) [IsGaussian μ] {ι : Type u_2} (s : Finset ι) (L : ι → StrongDual ℝ E) (c : ι → ℝ) :
MeasureTheory.Integrable (fun (x : E) => ∏ i ∈ s, ((L i) x + c i)) μ

Finite products of affine linear observables of a Gaussian measure are integrable, including the empty product.

The 2nth centered moment of a linear observable of a Gaussian measure is Varⁿ (2n - 1)‼, with the variance expressed by the covariance bilinear form.

theorem ProbabilityTheory.IsGaussian.integral_mul_mul_mul_sub_dual {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [MeasurableSpace E] [BorelSpace E] [CompleteSpace E] [SecondCountableTopology E] (μ : MeasureTheory.Measure E) [IsGaussian μ] (L₁ L₂ L₃ L₄ : StrongDual ℝ E) :
∫ (x : E), L₁ (x - ∫ (x : E), id x ∂μ) * L₂ (x - ∫ (x : E), id x ∂μ) * L₃ (x - ∫ (x : E), id x ∂μ) * L₄ (x - ∫ (x : E), id x ∂μ) ∂μ = ((covarianceBilinDual μ) L₁) L₂ * ((covarianceBilinDual μ) L₃) L₄ + ((covarianceBilinDual μ) L₁) L₃ * ((covarianceBilinDual μ) L₂) L₄ + ((covarianceBilinDual μ) L₁) L₄ * ((covarianceBilinDual μ) L₂) L₃

The fourth centered mixed moment of four linear observables of a Gaussian measure is the sum over their three pairings. No independence or nondegeneracy is assumed.

Every shifted coordinate monomial is integrable under every multivariate Gaussian. Taking c = 0 gives raw monomials and taking c = m gives centered monomials. No validity assumption on the covariance parameter is needed.

@[simp]
theorem TauCeti.Probability.integral_mul_mul_mul_sub_multivariateGaussian {ι : Type u_1} [Fintype ι] [DecidableEq ι] (m : EuclideanSpace ℝ ι) {S : Matrix ι ι ℝ} (hS : S.PosSemidef) (i j k l : ι) :
∫ (x : EuclideanSpace ℝ ι), (x.ofLp i - m.ofLp i) * (x.ofLp j - m.ofLp j) * (x.ofLp k - m.ofLp k) * (x.ofLp l - m.ofLp l) ∂ProbabilityTheory.multivariateGaussian m S = S i j * S k l + S i k * S j l + S i l * S j k

The fourth centered coordinate moment of a multivariate Gaussian, including repeated coordinates and positive-semidefinite singular covariance matrices.

@[simp]
theorem TauCeti.Probability.integral_mul_mul_mul_multivariateGaussian {ι : Type u_1} [Fintype ι] [DecidableEq ι] {S : Matrix ι ι ℝ} (hS : S.PosSemidef) (i j k l : ι) :
∫ (x : EuclideanSpace ℝ ι), x.ofLp i * x.ofLp j * x.ofLp k * x.ofLp l ∂ProbabilityTheory.multivariateGaussian 0 S = S i j * S k l + S i k * S j l + S i l * S j k

Isserlis' fourth-moment identity for centered multivariate Gaussian coordinates.

theorem ProbabilityTheory.HasLaw.integral_mul_mul_mul_sub_multivariateGaussian {ι : Type u_1} {Ω : Type u_2} [Fintype ι] [DecidableEq ι] [MeasurableSpace Ω] {P : MeasureTheory.Measure Ω} {X : Ω → EuclideanSpace ℝ ι} {m : EuclideanSpace ℝ ι} {S : Matrix ι ι ℝ} (hX : HasLaw X (multivariateGaussian m S) P) (hS : S.PosSemidef) (i j k l : ι) :
∫ (ω : Ω), ((X ω).ofLp i - m.ofLp i) * ((X ω).ofLp j - m.ofLp j) * ((X ω).ofLp k - m.ofLp k) * ((X ω).ofLp l - m.ofLp l) ∂P = S i j * S k l + S i k * S j l + S i l * S j k

The fourth centered coordinate moment of a random vector with a multivariate Gaussian law. Repeated coordinates and singular covariance matrices are allowed.