Prove the normal approximation of beta distribution How should I prove the normal approximation of beta distribution as follows:
Let $\mathrm B_{r_1, r_2}\sim \mathrm{Beta}(r_1, r_2)$, then prove that
$\sqrt{r_1+r_2} (\mathrm B_{r_1, r_2}- \dfrac{r_1}{r_1+r_2}) \to \mathrm N(0, \gamma(1-\gamma))$
where $r_1, r_2 \to \infty$ and $\frac{r_1}{r_1+r_2}\to \gamma$  $(0<\gamma<1)$.
My attempt: By some calculation, I figured out that $\mathbb E(\mathrm B_{r_1, r_2})=\dfrac{r_1}{r_1+r_2}$ and $\operatorname{Var}(\mathrm B_{r_1, r_2})=\dfrac{r_1 r_2}{(r_1 + r_2)^2 (r_1 +r_2 +1)}$.
Therefore, it remains to prove that $\sqrt{r_1 + r_2}\dfrac{\mathrm B_{r_1, r_2}-\mathbb E(\mathrm B_{r_1, r_2})}{\sqrt{\operatorname{Var}(\mathrm B_{r_1, r_2})}}$ converges to $Z \sim \mathrm N(0, 1)$.
Here I think I have to apply CLT, but I don't know how to because the given quantity does not contain sample mean. Does anyone have ideas?
Thanks for your help!

I finally figured out the answer:
Let $\mathrm B_{r_1, r_2} = \dfrac{V_1}{V_1+V_2}$ where $V_1 \sim \chi^2(2r_1 ) = \mathrm{Gamma}(r_1, 2)$ and $V_2 \sim \chi^2(2r_2 ) = \mathrm{Gamma}(r_2, 2)$ and $V_1, V_2$ are independent.
Now we know that $\dfrac{V_1-2r_1}{\sqrt{4r_1}} \to \mathrm N(0, 1)$ as $r_1\to\infty$ and $\dfrac{V_2-2r_2}{\sqrt{4r_2}} \to \mathrm N(0, 1)$ as $r_2\to\infty$. Therefore, $\sqrt{r_1+r_2}(\dfrac{V_1-2r_1}{\sqrt{4r_1(r_1+r_2)}}, \dfrac{V_2-2r_2}{\sqrt{4r_2(r_1+r_2)}})^t \to \mathrm N_2(0, I)$ where $I$ is the identity matrix.
Also, $\sqrt{r_1+r_2}\left(\dfrac{V_1-2r_1}{\sqrt{4r_1(r_1+r_2)}}, \dfrac{V_2-2r_2}{\sqrt{4r_2(r_1+r_2)}}\right)^t\sim \sqrt{n}\left(\dfrac {V_1}{\sqrt{4\gamma}n}-\sqrt{\gamma}, \dfrac {V_2}{\sqrt{4(1-\gamma)}n}-\sqrt{1-\gamma} \right)^t$
Where $r_1+r_2=n, r_1 \sim \gamma n, r_2 \sim (1-\gamma)n$.
Now let $g(x, y)=\dfrac{\sqrt{\gamma}x}{\sqrt{\gamma}x+\sqrt{1-\gamma}y}$, then $g'(x, y)=\left(\dfrac{\sqrt{\gamma (1-\gamma)}y}{\left(\sqrt{\gamma}x+\sqrt{1-\gamma}y\right)^2}, -\dfrac{\sqrt{\gamma (1-\gamma)}x}{\left(\sqrt{\gamma}x+\sqrt{1-\gamma}y\right)^2} \right)^t$.
Therefore, $\sqrt{n}\left( g\left(\dfrac {V_1}{\sqrt{4\gamma}n}, \dfrac {V_2}{\sqrt{4(1-\gamma)}n} \right)-g\left(\sqrt{\gamma}, \sqrt{1-\gamma}\right)  \right)=\sqrt{n}(\mathrm B_{r_1, r_2}-\gamma)\to g'\left(\sqrt{\gamma}, \sqrt{1-\gamma}\right)^t\mathrm N(0, I) = \mathrm N(0, \gamma(1-\gamma))$
Am I right?
