# Stein's method

Stein's method is a general method in probability theory to obtain bounds on the distance between two probability distributions with respect to a probability metric. It was introduced by Charles Stein, who first published it in 1972, to obtain a bound between the distribution of a sum of $m$ -dependent sequence of random variables and a standard normal distribution in the Kolmogorov (uniform) metric and hence to prove not only a central limit theorem, but also bounds on the rates of convergence for the given metric.

## History

At the end of the 1960s, unsatisfied with the by-then known proofs of a specific central limit theorem, Charles Stein developed a new way of proving the theorem for his statistics lecture. The seminal paper was presented in 1970 at the sixth Berkeley Symposium and published in the corresponding proceedings.

Later, his Ph.D. student Louis Chen Hsiao Yun modified the method so as to obtain approximation results for the Poisson distribution, therefore the method is often referred to as Stein-Chen method. Whereas moderate attention was given to the new method in the 70s, it has undergone major development in the 80s, where many important contributions were made and on which today's view of the method are largely based. Probably the most important contributions are the monograph by Stein (1986), where he presents his view of the method and the concept of auxiliary randomisation, in particular using exchangeable pairs, and the articles by Barbour (1988) and Götze (1991), who introduced the so-called generator interpretation, which made it possible to easily adapt the method to many other probability distributions. An important contribution was also an article by Bolthausen (1984) on a long-standing open problem around the so-called combinatorial central limit theorem, which surely helped the method to become widely known. {{ safesubst:#invoke:Unsubst||date=__DATE__ |\$B= {{#invoke:Category handler|main}}{{#invoke:Category handler|main}}[citation needed] }}

In the 1990s the method was adapted to a variety of distributions, such as Gaussian processes by Barbour (1990), the binomial distribution by Ehm (1991), Poisson processes by Barbour and Brown (1992), the Gamma distribution by Luk (1994), and many others.

## The basic approach

### Probability metrics

Stein's method is a way to bound the distance of two probability distributions in a specific probability metric. To be tractable with the method, the metric must be given in the form

$(1.1)\quad d(P,Q)=\sup _{h\in {\mathcal {H}}}\left|\int hdP-\int hdQ\right|=\sup _{h\in {\mathcal {H}}}\left|Eh(W)-Eh(Y)\right|$ Here, $P$ and $Q$ are probability measures on a measurable space ${\mathcal {X}}$ , $W$ and $Y$ are random variables with distribution $P$ and $Q$ respectively, $E$ is the usual expectation operator and ${\mathcal {H}}$ is a set of functions from ${\mathcal {X}}$ to the real numbers. This set has to be large enough, so that the above definition indeed yields a metric. Important examples are the total variation metric, where we let ${\mathcal {H}}$ consist of all the indicator functions of measurable sets, the Kolmogorov (uniform) metric for probability measures on the real numbers, where we consider all the half-line indicator functions, and the Lipschitz (first order Wasserstein; Kantorovich) metric, where the underlying space is itself a metric space and we take the set ${\mathcal {H}}$ to be all Lipschitz-continuous functions with Lipschitz-constant 1. However, note that not every metric can be represented in the form (1.1).

In what follows we think of $P$ as a complicated distribution (e.g. a sum of dependent random variables), which we want to approximate by a much simpler and tractable distribution $Q$ (e.g. the standard normal distribution to obtain a central limit theorem).

### The Stein operator

We assume now that the distribution $Q$ is a fixed distribution; in what follows we shall in particular consider the case where $Q$ is the standard normal distribution, which serves as a classical example of the application of Stein's method.

$(2.1)\quad E({\mathcal {A}}f)(Y)=0{\text{ for all }}f\quad \iff \quad Y{\text{ has distribution }}Q.$ We call such an operator the Stein operator. For the standard normal distribution, Stein's lemma exactly yields such an operator:

$(2.2)\quad E\left(f'(Y)-Yf(Y)\right)=0{\text{ for all }}f\in C_{b}^{1}\quad \iff \quad Y{\text{ has standard normal distribution.}}$ thus we can take

$(2.3)\quad ({\mathcal {A}}f)(x)=f'(x)-xf(x)$ We note that there are in general infinitely many such operators and it still remains an open question, which one to choose. However, it seems that for many distributions there is a particular good one, like (2.3) for the normal distribution.

There are different ways to find Stein operators. But by far the most important one is via generators. This approach was, as already mentioned, introduced by Barbour and Götze. Assume that $Z=(Z_{t})_{t\geq 0}$ is a (homogeneous) continuous time Markov process taking values in ${\mathcal {X}}$ . If $Z$ has the stationary distribution $Q$ it is easy to see that, if ${\mathcal {A}}$ is the generator of $Z$ , we have $E({\mathcal {A}}f)(Y)=0$ for a large set of functions $f$ . Thus, generators are natural candidates for Stein operators and this approach will also help us for later computations.

### Setting up the Stein equation

To make this statement rigorous we could find a function $f$ , such that, for a given function $h$ ,

$(3.1)\quad E({\mathcal {A}}f)(W)=Eh(W)-Eh(Y),$ so that the behavior of the right hand side is reproduced by the operator ${\mathcal {A}}$ and $f$ . However, this equation is too general. We solve instead the more specific equation

$(3.2)\quad ({\mathcal {A}}f)(x)=h(x)-Eh(Y),\qquad {\text{for all }}x,$ which is called Stein equation. Replacing $x$ by $W$ and taking expectation with respect to $W$ , we are back to (3.1), which is what we effectively want. Now all the effort is worth only if the left hand side of (3.1) is easier to bound than the right hand side. This is, surprisingly, often the case.

If $Q$ is the standard normal distribution and we use (2.3), the corresponding Stein equation is

$(3.3)\quad f'(x)-xf(x)=h(x)-Eh(Z),\qquad {\text{for all }}x,$ which is just an ordinary differential equation.

### Solving the Stein equation

Now, in general, we cannot say much about how the equation (3.2) is to be solved. However, there are important cases, where we can.

Analytic methods. We see from (3.3) that equation (3.2) can in particular be a differential equation (if $Q$ is concentrated on the integers, it will often turn out to be a difference equation). As there are many methods available to treat such equations, we can use them to solve the equation. For example, (3.3) can be easily solved explicitly:

$(4.1)\quad f(x)=e^{x^{2}/2}\int _{-\infty }^{x}[h(s)-Eh(Y)]e^{-s^{2}/2}ds.$ Generator method. If ${\mathcal {A}}$ is the generator of a Markov process $(Z_{t})_{t\geq 0}$ as explained before, we can give a general solution to (3.2):

$(4.2)\quad f(x)=-\int _{0}^{\infty }[E^{x}h(Z_{t})-Eh(Y)]dt.$ ### Properties of the solution to the Stein equation

After showing the existence of a solution to (3.2) we can now try to analyze its properties. Usually, one tries to give bounds on $f$ and its derivatives (which has to be carefully defined if ${\mathcal {X}}$ is a more complicated space) or differences in terms of $h$ and its derivatives or differences, that is, inequalities of the form

$(5.1)\quad ||D^{k}f||\leq C_{k,l}||D^{l}h||,$ for some specific $k,l=0,1,2,\dots$ (typically $k\geq l$ or $k\geq l-1$ , respectively, depending on the form of the Stein operator) and where often $||\cdot ||$ is taken to be the supremum norm. Here, $D^{k}$ denotes the differential operator, but in discrete settings it usually refers to a difference operator. The constants $C_{k,l}$ may contain the parameters of the distribution $Q$ . If there are any, they are often referred to as Stein factors or magic factors.

In the case of (4.1) we can prove for the supremum norm that

$(5.2)\quad ||f||_{\infty }\leq \min\{{\sqrt {\pi /2}}||h||_{\infty },2||h'||_{\infty }\},\quad ||f'||_{\infty }\leq \min\{2||h||_{\infty },4||h'||_{\infty }\},\quad ||f''||_{\infty }\leq 2||h'||_{\infty },$ where the last bound is of course only applicable if $h$ is differentiable (or at least Lipschitz-continuous, which, for example, is not the case if we regard the total variation metric or the Kolmogorov metric!). As the standard normal distribution has no extra parameters, in this specific case, the constants are free of additional parameters.

Note that, up to this point, we did not make use of the random variable $W$ . So, the steps up to here in general have to be calculated only once for a specific combination of distribution $Q$ , metric $d$ and Stein operator ${\mathcal {A}}$ . However, if we have bounds in the general form (5.1), we usually are able to treat many probability metrics together. Furthermore as there is often a particular 'good' Stein operator for a distribution (e.g., no other operator than (2.3) has been used for the standard normal distribution up to now), one can often just start with the next step below, if bounds of the form (5.1) are already available (which is the case for many distributions).

### An abstract approximation theorem

We are now in a position to bound the left hand side of (3.1). As this step heavily depends on the form of the Stein operator, we directly regard the case of the standard normal distribution.

Now, at this point we could directly plug in our random variable $W$ which we want to approximate and try to find upper bounds. However, it is often fruitful to formulate a more general theorem using only abstract properties of $W$ . Let us consider here the case of local dependence.

Using basically only Taylor expansion, it is possible to prove that

$(6.1)\quad \left|E(f'(W)-Wf(W))\right|\leq ||f''||_{\infty }\sum _{i=1}^{n}\left({\frac {1}{2}}E|X_{i}X_{A_{i}}^{2}|+E|X_{i}X_{A_{i}}X_{B_{i}\setminus A_{i}}|+E|X_{i}X_{A_{i}}|E|X_{B_{i}}|\right)$ Note that, if we follow this line of argument, we can bound (1.1) only for functions where $||h'||_{\infty }$ is bounded because of the third inequality of (5.2) (and in fact, if $h$ has discontinuities, so will $f''$ ). To obtain a bound similar to (6.1) which contains only the expressions $||f||_{\infty }$ and $||f'||_{\infty }$ , the argument is much more involved and the result is not as simple as (6.1); however, it can be done.

$(6.2)\quad d_{W}({\mathcal {L}}(W),N(0,1))\leq 2\sum _{i=1}^{n}\left({\frac {1}{2}}E|X_{i}X_{A_{i}}^{2}|+E|X_{i}X_{A_{i}}X_{B_{i}\setminus A_{i}}|+E|X_{i}X_{A_{i}}|E|X_{B_{i}}|\right).$ Proof. Recall that the Lipschitz metric is of the form (1.1) where the functions $h$ are Lipschitz-continuous with Lipschitz-constant 1, thus $||h'||\leq 1$ . Combining this with (6.1) and the last bound in (5.2) proves the theorem.

Thus, roughly speaking, we have proved that, to calculate the Lipschitz-distance between a $W$ with local dependence structure and a standard normal distribution, we only need to know the third moments of $X_{i}$ and the size of the neighborhoods $A_{i}$ and $B_{i}$ .

### Application of the theorem

$(7.1)\quad d_{W}({\mathcal {L}}(W),N(0,1))\leq {\frac {5E|X_{1}|^{3}}{n^{1/2}}}.$ ## Connections to other methods

• Lindeberg's method. Lindeberg (1922) introduced in a seminal article a method, where the difference in (1.1) is directly bounded. This method usually also heavily relies on Taylor expansion and thus shows some similarities with Stein's method.