Inverse second moment of the chi-square distribution: for
ProvedGaussianMatrix.inv_sq_chi_square_momentchi-squaregaussianmomentsprobability
Let and let have independent standard normal coordinates, so that has the chi-square distribution with degrees of freedom. Then is integrable and
where is the standard Gaussian measure on .
This is the second inverse moment of the chi-square law. Applied with degrees of freedom it yields for an standard Gaussian matrix , the diagonal ingredient of the second moments of the inverse Wishart matrix.
Formalization Note. The integrand is written ((∑ j, x j ^ 2)⁻¹) ^ 2; with Lean's convention it is defined everywhere (the event is null). The hypothesis is sharp: for the integral diverges.
Preamble
import Definitions.Def_GaussianMatrix_basic open MeasureTheory ProbabilityTheory open scoped Matrix
Formal statement
namespace GaussianMatrix
theorem inv_sq_chi_square_moment {d : ℕ} (hd : 5 ≤ d) :
Integrable (fun x : Fin d → ℝ => ((∑ j, x j ^ 2)⁻¹) ^ 2)
(Measure.pi fun _ : Fin d => gaussianReal 0 1) ∧
∫ x, ((∑ j, x j ^ 2)⁻¹) ^ 2 ∂(Measure.pi fun _ : Fin d => gaussianReal 0 1)
= 1 / (((d : ℝ) - 2) * ((d : ℝ) - 4)) := by sorry
end GaussianMatrixSource
standard fact: for and , ; see N. L. Johnson, S. Kotz, N. Balakrishnan, Continuous Univariate Distributions, Vol. 1, 2nd ed., Wiley 1994, Ch. 18 (moments of the chi-square distribution) [chapter cited from memory]; the case is Halko–Martinsson–Tropp, SIAM Review 53(2), 2011, Appendix A.