I'm wondering how I can calculate the power function of the test statistic
$$T_n=\frac{1}{n_1(n_1-1)}\frac{1}{n_2(n_2-1)}\sum_{k=1}^p \underset{i\not=j}{\sum^{n_1}\sum^{n_1}} \underset{s\not=t}{\sum^{n_2}\sum^{n_2}}\frac{(X_{1ik}-X_{2sk})(X_{1jk}-X_{2tk})}{\hat{\sigma}^2_{1k(i,j)}+\gamma\hat{\sigma}^2_{2k(s,t)}}$$
where $\gamma=n_1/n_2$, $\hat{\sigma}^2_{1k(i,j)}$ is the sample variance of $\{X_{1lk}\}_{l=1}^{n_1}$ excluding $X_{1ik}$ and $X_{1jk}$
which is given in Feng and Sun (2015) [1], page 31, from theorem 1 (same page):
Under Conditions (C1)--(C3), as $p,n \to \infty$, $\frac{T_n-E(T_n)}{\sqrt{\mathrm {var}(T_n)}}\overset{\mathcal{L}}{\longrightarrow} N(0,1).$
Any help will be appreciated.
[1]: L.Feng and F.Sun (2015),
"A note on high-dimensional two-sample test",
Statistics and Probability Letters, 105, 29–36
http://www.sciencedirect.com/science/article/pii/S0167715215001789
[arXiv version arXiv:1502.05455 [stat.ME]
here]