Polynomial Stein operators: a noncommutative algebra perspective

TL;DR

This paper explores the algebraic structure of polynomial Stein operators using noncommutative algebra, specifically the Weyl algebra, providing a complete description for the Gaussian case and insights into their characterising properties.

Contribution

It introduces a novel algebraic perspective on polynomial Stein operators by linking them to the Weyl algebra, offering a complete Gaussian case characterization and analyzing their distributional properties.

Findings

Complete description of Gaussian polynomial Stein operators as a vector space.

Identification of the Gaussian Stein operators as a principal right ideal in the Weyl algebra.

Insight into the non-characterising nature of many polynomial Stein operators for Gaussian distributions.

Abstract

In this paper, we make a novel connection between Stein's method and noncommutative algebra by viewing polynomial Stein operators (Stein operators with polynomial coefficients) as elements of the first Weyl algebra. Through this connection we study the algebraic structure of classes of polynomial Stein operators. In the case of the standard Gaussian distribution, we provide a complete description of the corresponding class of polynomial Stein operators by (i) identifying it as a vector space over $\mathbb{R}$ with an explicit given basis and (ii) by showing that this class is a principal right ideal of the first Weyl algebra generated by the classical Gaussian Stein operator $\partial -x$, with $\partial$ denoting the usual differential operator. We also study the characterising property of polynomial Stein operators for the standard Gaussian distribution, and give examples of general classes of polynomial Stein operators that are characterising, as well as classes that are not characterising unless additional distributional assumptions are made. By appealing to a standard property of Weyl algebras, we shown that the non-characterising property possessed by a wide class of polynomial Stein operators for the standard Gaussian distribution is a consequence of a general result that is perhaps surprising from a probabilistic perspective: the intersection between the class of polynomial Stein operators for any two target distributions with holonomic densities or holonomic characteristic functions is non-trivial.