Stein's method for functions of multivariate normal random variables

Abstract

By the continuous mapping theorem, if a sequence of d-dimensional random vectors (Wn)n≥1 converges in distribution to a multivariate normal random variable 1/2Z, then the sequence of random variables (g(Wn))n≥1 converges in distribution to g(1/2Z) if g:Rd→R is continuous. In this paper, we develop Stein's method for the problem of deriving explicit bounds on the distance between g(Wn) and g(1/2Z) with respect to smooth probability metrics. We obtain several bounds for the case that the j-component of Wn is given by Wn,j=1nΣi=1nXij, where the Xij are independent. In particular, provided g satisfies certain differentiability and growth rate conditions, we obtain an order n-(p-1)/2 bound, for smooth test functions, if the first p moments of the Xij agree with those of the normal distribution. If p is an even integer and g is an even function, this convergence rate can be improved further to order n-p/2. These convergence rates are shown to be of optimal order. We apply our general bounds to some examples, which include the distributional approximation of asymptotically chi-square distributed statistics; the approximation of expectations of smooth functions of binomial and Poisson random variables; rates of convergence in the delta method; and a quantitative variance-gamma approximation of the D2* statistic for alignment-free sequence comparison in the case of binary sequences.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…