|
We consider the testing problem, one-sided or two-sided, about the correlation coefficient ρ of a bivariate normal distribution. This is not a new inferential problem and good tests exist in the literature, for example, Lehmann (1994) gives the corresponding UMPU and UMPI tests which are based on the sample correlation coefficient R. However, for the alternative H1:ρ>ρ0 with ρ0≠0, the tests above are not easy to implement since their critical points need to be calculated numerically, due to the complicated density function of R. Significantly establish H1:ρ>ρ0 with ρ0 >>0 may be valuable for decision making in some applications, yet one can not get the conclusion of the classical UMPI immediately ─ this motivates our study here.
Hence, in this work, we aim at constructing better level α tests which are easier to implement for any givenρ0, especially for large |ρ0 |. First, we make use of a new confidence interval forρprovided by Berger and Sun (2008) to construct new tests through the duality relationship between confidence interval and tests. Next, we propose new UMPU tests by adding the equal variance assumption.
By the results of simulations and numerical calculations, we notice first that the power functions of the tests by inverting confidence intervals provided by Berger and Sun (2008) are almost identical to the classical UMPI tests; secondly, when equal variances assumption holds, our test is theoretically proved to have larger power function than the classical UMPI test and numerical results indicate that the improvement rate can be more than 100%, quite substantial, in some practically important cases. When variances are not equal, we provide a modified version of our test. Numerical results indicates that, for moderate to large sample sizes, our modified test may have a slightly larger level than the nominal level, yet, in such cases, it has larger power functions than the UMPI and the improvement can be substantial -- acceptable to many practioners.
Moreover, we study the robustness of these tests under a bivariate skew-normal distribution. By the results of simulations, we notice that the power functions of these tests all fail to retain the claimed level. The discrepancies are large as the skew parameters become larger or as the sample sizes become larger. Hence, these tests are not robust against bivariate skew-normality.
|