On algebraic Stein operators for Gaussian polynomials

From MaRDI portal
Publication:6330828

DOI10.3150/22-BEJ1460arXiv1912.04605MaRDI QIDQ6330828FDOQ6330828


Authors: Ehsan Azmoodeh, Dario Gasbarra, Robert E. Gaunt Edit this on Wikidata


Publication date: 10 December 2019

Abstract: The first essential ingredient to build up Stein's method for a continuous target distribution is to identify a so-called extit{Stein operator}, namely a linear differential operator with polynomial coefficients. In this paper, we introduce the notion of extit{algebraic} Stein operators (see Definition ef{def:algebraic-Stein-Operator}), and provide a novel algebraic method to find emph{all} the algebraic Stein operators up to a given order and polynomial degree for a target random variable of the form Y=h(X), where X=(X1,dots,Xd) has i.i.d. standard Gaussian components and hinmathbbK[X] is a polynomial with coefficients in the ring mathbbK. Our approach links the existence of an algebraic Stein operator with extit{null controllability} of a certain linear discrete system. A exttt{MATLAB} code checks the null controllability up to a given finite time T (the order of the differential operator), and provides all extit{null control} sequences (polynomial coefficients of the differential operator) up to a given maximum degree m. This is the first paper that connects Stein's method with computational algebra to find Stein operators for highly complex probability distributions, such as H20(X1), where Hp is the p-th Hermite polynomial. Some examples of Stein operators for Hp(X1), p=3,4,5,6, are gathered in the Appendix and many other examples are given in the Supplementary Information.













This page was built for publication: On algebraic Stein operators for Gaussian polynomials

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6330828)