logo

Approximation of Normal Distribution Variance Stabilization from a Binomial Distribution 📂Probability Distribution

Approximation of Normal Distribution Variance Stabilization from a Binomial Distribution

Example 1

If Y=YnY = Y_{n} follows a binomial distribution Bin(n,p)\text{Bin} (n,p), arcsinYnDN(arcsinp,n/4) \arcsin \sqrt{ {{ Y } \over { n }} } \overset{D}{\to} N \left( \arcsin \sqrt{p} , n/4 \right)


Explanation

The binomial distribution Bin(n,p)\text{Bin} (n, p ) converges to the normal distribution N(np,np(1p))N \left( np, np(1-p) \right) when nn \to \infty, so there’s nothing particularly miraculous about the normal distribution itself. However, through such a transformation, it’s possible to obtain a limit distribution with a variance constant, regardless of the parameter pp. This is more of a mathematical trick than a practical application, interesting because it employs the differentiation technique of arcsine, which seemed like it wouldn’t be seen after freshman year in college.

Let’s assume the variance of u(Y/n)u ( Y/n ) with uu applied to Y/nY/n is independent of pp. For sufficiently large nn, since Y/npY/n \approx p, then by expanding Taylor’s series, u(Yn)u(p)+(Ynp)u(p) u \left( {{ Y } \over { n }} \right) \approx u (p) + \left( {{ Y } \over { n }} - p \right) u ' (p)

Taking the expectation on both sides yields the mean of u(Y/n)u (Y/n) as u(p)u (p), and the variance, by its property, is Var(aX+b)=a2Var(X)\Var (aX + b) = a^{2} \Var (X), [u(p)]2p(1p)n \left[ u ' (p) \right]^{2} {{ p ( 1 - p ) } \over { n }}

For this variance to become independent of pp, square of u(p)u ' (p) in the numerator must be canceled out by p(1p)p ( 1 - p). Hence, setting such that u(p)=du(p)dp=cp(1p) u ‘(p) = {{ du(p) } \over { dp }} = {{ c } \over { \sqrt{ p ( 1 - p )} }} erases pp from the variance of u(p)u(p). The solution to this differential equation can be directly found using the differentiation of the arcsine function.

Differentiation of inverse trigonometric functions: (arcsinx)=11x2 \left( \arcsin x \right)' = {{ 1 } \over { \sqrt{1-x^{2}} }}

The solution to the differential equation looks a bit different in the denominator, but if you check, it exactly matches due to the differentiation of p\sqrt{p}. u(p)=2carcsinp u (p) = 2c \arcsin \sqrt{p} Since we can set and solve the differential equation in the same manner regardless of what cc is, it doesn’t matter what it is, but for aesthetics, setting it as c=1/2c = 1/2 gives the appearance introduced in the example.


  1. Hogg et al. (2013). Introduction to Mathematical Statistcs(7th Edition): p318. ↩︎