A New Test for Correlation on Bivariate Nonnormal Distributions

Journal of Modern Applied Statistical Methods Volume 5 Issue Article 8 --06 A New Test for Correlation on Bivariate Nonnormal Distributions Ping Wang Great Basin College, ping.wang@gbcnv.edu Ping Sa University of North Florida, psa@unf.edu Follow this and additional works at: http://digitalcommons.wayne.edu/jmasm Part of the Applied Statistics Commons, Social and Behavioral Sciences Commons, and the Statistical Theory Commons Recommended Citation Wang, Ping and Sa, Ping (06) "A New Test for Correlation on Bivariate Nonnormal Distributions," Journal of Modern Applied Statistical Methods: Vol. 5 : Iss., Article 8. DOI: 0.7/jmasm/7800560 Available at: http://digitalcommons.wayne.edu/jmasm/vol5/iss/8 This Regular Article is brought to you for free and open access by the Open Access Journals at DigitalCommons@WayneState. It has been accepted for inclusion in Journal of Modern Applied Statistical Methods by an authorized editor of DigitalCommons@WayneState.

Journal of Modern Applied Statistical Methods November 06, Vol. 5, No., 8-7. doi: 0.7/jmasm/7800560 Copyright 06 JMASM, Inc. ISSN 58 97 A New Test for Correlation on Bivariate Nonnormal Distributions Ping Wang Great Basin College Elko, NV Ping Sa University of North Florida Jacksonville, FL A new method to conduct a right-tailed test for the correlation on bivariate non-normal distribution is proposed. The comparative simulation study shows that the new test controls the type I error rates well for all the distributions considered. An investigation of the power performance is also provided. Keywords: Correlation, Edgeworth expansion, Cornish-Fisher inverse expansion, type I error rate, power performance Introduction Bivariate data is data that has two variables. In the bivariate case, the study of the relationship between the two variables is at least as important as analyzing each variable individually. The most popular measure of the strength of the linear relation between two variables is the correlation coefficient, denoted by ρ. The Pearson product-moment correlation, r, is the most frequently-used estimator for ρ. Another widely-used estimator is the Spearman s rank correlation, denoted by r s. Tests Based on Pearson Product-Moment Correlation Pearson (896) developed the initial mathematical formulas for the sample correlation coefficient. Let (X i, Y i ), i =,, n be a random sample, the statistic r is given by: r s X i X Yi Y X i X Yi Y XY sxsy Ms. Wang is the Director of the Academic Success and Testing Center. Email her at: ping.wang@gbcnv.edu. Dr. Sa is a Professor in the Department of Mathematics and Statistics. Email her at: psa@unf.edu. 8

WANG & SA where s XY is the sample covariance of X and Y, n is the sample size, s X and s Y are the sample standard deviations, and X and Y are the sample means for the variables X and Y, respectively. The Pearson product-moment correlation r is the maximum likelihood estimator of the parameter ρ when the population has a bivariate normal distribution. Although r is a biased estimator, the bias is negligible when the sample size is large. Researchers have done intensive work on the distribution of r when the population is bivariate normal (Fisher, 95; Stuart & Ord, 99). r can be used to test H 0 : ρ = 0 when the population is a bivariate normal distribution. The test statistic t * r r n r follows the Student s t-distribution with n degrees of freedom under H 0. r can also be used to test H 0 : ρ = ρ 0, for - ρ 0. The sampling distribution of r is complicated and unstable even when the population is bivariate normal. Fisher (9) introduced a remarkable transformation of r, which tends to normality much faster. When the sample size n is moderately large, given r r ln, ln r * * the distribution of r * ρ * approaches to normal with an approximate mean and variance n. Note that n > 50 is an adequate sample size for the n above approximation (see David, 98). To test H 0 : ρ = ρ 0, the test statistic is: z F r 0 0 ln ln r n n 0 () z F has approximately a standard normal distribution under H 0. 9

A NEW TEST FOR CORRELATION Test Based on Spearman Rank Correlation Spearman (90) proposed a rank correlation which can be used to measure the relationship between two variables when the distribution is neither bivariate normal nor transformed to a bivariate normal. The Spearman rank correlation, r s, is a nonparametric version of the Pearson product-moment correlation. Let (R i, R i ), i =,, n be the paired rank data of two variables, r s is given by: r s s R i R Ri R R i R Ri R RR srs R where s RR is the sample covariance of the paired ranks and s R and R, s R and R are the sample standard deviation and the sample mean of the ranks of the two variables, respectively. The Spearman rank correlation r s can be used to test: The test statistic is H 0 : there is no association between the rank pairs t * s r s n r s which follows the Student s t-distribution with n degrees of freedom under H 0. Other Tests on Correlation The test based on r can only be used when the population is bivariate normal or the sample size is relatively large. Although the test based on r s is applicable to the distribution-free case, it is less powerful and limited to test for zero correlation. However, in real world situations, most distributions are not bivariate normal and the sample sizes may not be large. Furthermore, a test of non-zero correlation is often required. It is desired to develop methods to meet these needs. Beasley et al. (007) proposed two new approaches to test a non-zero by using the bootstrapping method. Their methods do not require any knowledge of the population. One is the hypothesis-imposed univarite sampling bootstrap (HI) and the other one is the observed-imposed univariate sampling bootstrap (OI). Two tests 0

WANG & SA are conducted on populations with various combinations of normal and skewed variates with ρ 0. and the sample size n 0. Their study demonstrated that although OI is preferable to HI under the significance level of 0.05, the type I error rates are still slightly inflated. Also, the simulated populations are limited to the combinations between normal and skewed populations. The methods are not evaluated under the situations that both variables are non-normal. Another drawback of these two methods is that they are computer-intensive methods. Unfortunately, most practitioners do not have the computer programming skills to implement these methods. Beversdorf and Sa (0) proposed tests of correlation for bivariate nonnormal data with small sample sizes. The tests investigated are Fisher s Z transformation z F and the saddlepoint approximation r L. They found that z F and r L have extremely similar performance which could control the type I error rates well when a left-tailed test was performed under all the bivariate non-normal distributions considered. Both methods essentially failed to control the type I error rates when a right-tailed test is desired. The purpose of this study is to develop a new right-tailed test on bivariate non-normal distributions with non-zero correlation. The new test statistic is derived using the Edgeworth expansion and the Cornish-Fisher inverse expansion. Methodology Edgeworth Expansion The Edgeworth expansion was derived by Edgeworth (905), and uses a series to approximate a probability distribution in terms of its cumulants. Let ˆ be an n ˆ be asymptotically normally estimator of an unknown parameter θ, and distributed with mean zero and variance σ. Hall (98) developed the Edgeworth n ˆ as a power series in n. expansion of the distribution function of n ˆ j P u u n p uu n p j uu ()

A NEW TEST FOR CORRELATION where Φ(u), ϕ(u), and p j (u) denote the standard normal distribution function, its density function, and a polynomial function with coefficients depending on cumulants of ˆ, respectively. The inverse of the Edgeworth expansion, obtained by inverting the formula (), is known as the Cornish-Fisher expansion: n ˆ P z n p z n p z n p j z z j where z is the percentile of the standard normal distribution and the p j are polynomials defined in terms of p j s (Hall, 99). Proposed Test Procedure Assume that a bivariate population has finite cumulants and a correlation coefficient ρ. Let κ 0, κ 0, κ 0, κ 0, κ, up to order six be the product cumulants r n r has a limiting normal distribution for the bivariate population. Then * with mean zero and constant variance σ, where σ is of the form (Nakagawa & Niki, 99): 0 0 0 0 0 0 0 0 0 00 0 0 0 Bhattacharya and Ghosh (978) provided the Edgeworth expansion of R, where as r * n r R

WANG & SA P H H n 6 n 7 H u 6 R u u u u u 6 5 H u O n n () where Φ(u) and ϕ(u) denote the standard normal distribution function and its density function, O n n is the big-oh function of order, H (u), H (u), n n H (u), and H 5 (u) are Hermite polynomials with H u u; H u u ; H u u u; H u u 0u 5 u, 5 5 and ν, ν, ν, and ν are parameters such that,,, n n n n are the approximate cumulants of R. The explicit forms of ν, ν, ν, and ν were provided by Nakagawa and Niki (99). Formulas for calculating ν and ν are listed in Appendix A. Formulas for calculating ν and ν are not needed in this study. Nakagawa and Niki (99) applied the inverted Edgeworth expansion to the distribution of R of order /n: P R z z z 6 n 6 n 8 5 z 6 6 8 z O n n () If only order n is required, then () can be reduced to a simpler form:

A NEW TEST FOR CORRELATION P R z B z n O (5) n 6. To test H 0 : ρ = ρ 0 versus H a : ρ > ρ 0, the intuitive decision rule is: where B z reject H : when R z B,i.e. 0 0 reject H : under the significance level, if 0 0 n r 0 z z n 6 (6) Since negative values of B might increase type I errors, the following adjustment is proposed: Define B max z, z z n 6 (7) The decision rule is adjusted to: reject H if n r 0 0 B (8) All the parameters in (6) and (7) can be written in terms of the product cumulants. These product cumulants κ ij are estimated by their corresponding unbiased estimators k ij. Detailed formulas are provided in Appendix B. For the special case of ρ 0 = 0, κ 0 = κ 0 = 0, κ 0 = κ 0 =, and κ pq = 0 for p + q 5, Nakagawa and Niki (99) gave the simplified forms for parameters σ, ν, and ν as follows: 0 0 (9)

WANG & SA To test H 0 : ρ = 0, (6) and (7) are evaluated with the parameters given in (9). Again, all the parameters are estimated by their corresponding unbiased estimators. Simulation Study The simulation study was implemented to evaluate type I error rates, to investigate the power performance, and to compare with the existing Fisher s Z transformation method on the type I error rates. Simulation Description Fleishman (978) proposed a method to generate univariate non-normal random variables with desired coefficients of skewness β and kurtosis γ. Vale and Maurelli (98) extended Fleishman s method to the bivariate non-normal case with a specified correlation coefficient. Five parameters, including two sets of skewness and kurtosis and one correlation coefficient, are required to generate the bivariate non-normal data using Vale and Maurelli method. Seven levels of the skewness, -.0, -., -0.5, 0.0, 0.5,., and.0, and five levels of the kurtosis, 0.0,.0, 0.0,.0, and 5.0, were considered, and combinations were selected. Moreover, five correlation coefficients, 0.0, 0.5, 0.6, 0.75, and 0.9, three significance levels, 0.0, 0.05, and 0.0, and two sample sizes, 5 and 0, were used in the simulation study. Two new methods and the Fisher s Z transformation method were evaluated. The method using (6) was denoted by Zb, and the one using (8) was denoted by Zc. The Fisher s Z transformation method () was shortened as Zf. Both Zb and Zc methods were evaluated with two critical values, z α and t (α,n ). The Algorithm of the Test on Correlation:. Input the desired ρ X,Y and two sets of skewness and kurtosis, (β, γ ) and (β, γ ).. Generate n bivariate non-normal random variates (X, Y) based on the given parameters.. Calculate z F in (), Zb in (6), and Zc in (8).. Compare the tests with their critical values; count one if the test is rejected. 5. Repeat () () 99,999 times. 5

A NEW TEST FOR CORRELATION 6. Calculate the type I error rate, the proportion of the false rejection (out of 00,000) for each test. In the power study, an extra parameter ρ a is input in step () and used to generate the data as the true population correlation. However, all of the test statistics in step () are evaluated under ρ 0. All fo the simulations were run with Fortran 77 for Windows on an IBM T6 Laptop Computer. Simulation Results Type I Error Rate Comparisons Tables - provide the comparative study of the type I error rates on various bivariate non-normal distributions with significance levels 0.0, 0.05, and 0.0 and sample size 0. Comparisons were made among the tests Zf, Zb, and Zc with two critical values, z α and t (α,n ), while Zf only used the critical value, z α. The correlation coefficients 0.00, 0.50, 0.60, 0.75, and 0.90 were targeted during the simulation study. A total of bivariate non-normal distributions with various population conditions were examined. Table shows the results on testing a zero correlation. It can be observed that the Zc method controls the type I error rates well. On the contrary, the Zb method do not control type I error rates at all. Almost all of the type I error rates obtained by the Zb method are slightly inflated except for a few cases. The Zf method can control the type I error rates as long as the skewness and kurtosis are small. Once theses parameters increase, Zf becomes unstable and fails to control the type I error rates in many cases. More specifically, in testing ρ 0 = 0 on a distribution which is bivariate normal or very close to bivariate normal, Zf controls type I error rates a bit better than Zc. However, for the non-normal distributions, Zc is better than Zf in controlling type I error rates. Tables and give the results for right-tailed tests on the non-zero correlation. It is quite interesting to see that the hypothesized value ρ 0 actually affects the type I error rate performance. When n = 0 and ρ 0 = 0.5, both the Zf and Zb methods basically fail to control the type I error rates with very few exceptions. The type I error rates obtained by the Zc method have better performance. However, the cases with controlled type I error rates are restricted to the distributions with small to moderate skewness and kurtosis. When ρ 0 increases to 0.6, the Zc method successfully controls the type I error rates for nearly all the distributions considered 6

WANG & SA with the t critical point. As ρ 0 increases, the type I error rates get more conservative. This tendency can be observed on both Zb and Zc methods. The Zf method fails completely in the right-tailed test on non-zero correlation with only a few exceptions. This result confirms with the study by Beversdorf and Sa (0). Their study shows that Zf can properly control the type I error rates on the left-tailed test but not on the right-tailed test. Therefore, it is fair to conclude that, for the right-tailed test, the only method that can properly control the type I error rates is the Zc method with the t critical point. Due to the similar results in the study, only the moderate sample size 0 and significance levels of 0.05 and 0.0 are reported in the tables. Power Results The power performance of the proposed test is also evaluated. Tables and 5 provide the power performance to test ρ 0 = 0 when ρ a = 0.0, 0., 0., 0.6, and 0.8 with significance levels 0.05 and 0.0. Table 6 provides a small-scale investigation on the power performance to test ρ 0 = 0.55 and ρ a = 0.6 and 0.7. Both the Zf and Zc methods perform well in testing ρ 0 = 0. In testing on an exactly- or nearly-normal distribution, the power from Zf and Zc converges to quickly. When ρ a = 0.6, both achieve a power of 0.99; when ρ a = 0.8, the power rates are essentially. For the distributions with large skewness and kurtosis, the Zc method, which is the only one with controlled type I error rates, converges to more slowly but still reasonably well. A small-scale power study to test non-zero correlation is presented in Table 6. At significance level 0.0, sample size 0, ρ 0 = 0.55 versus ρ a = 0.6 and 0.7, it is observed that the power of Zc steadily increases when ρ a moves away from ρ 0. Conclusions This study proposed a new right-tailed test for the correlation of bivariate nonnormal distributions. This new test adapts the inverse Edgeworth expansion for the nr standardized correlation R by Nakagawa and Niki (99). This newly proposed test can be conducted without any knowledge of the populations. The simulation study shows that this new right-tailed test has the best performance in controlling the type I error rates. The proposed method, along with the t critical point, can be used to test both ρ 0 = 0 and any value of ρ 0 when ρ 0 > 0.5. The power performance of the new test was also evaluated. Zc is as powerful as Zf when testing ρ 0 = 0. To test non-zero correlations, it is meaningless to 7

A NEW TEST FOR CORRELATION compare the two tests since Zf fails to control type I error rates. The power examination of the Zc method shows that the power steadily goes up when ρ a moves away from ρ 0. The new test does have its own limitations. It cannot control the type I error rates well when the population has a small correlation and it is a right-tailed test. In order to better control the type I error rates, a higher-order Edgeworth expansion may be considered. Unfortunately, this might lead to tedious computations when higher-order terms are introduced in the test. References Beasley, W. H., DeShea, L., Toothaker, L. E., Mendoza, J. L., Bard, D. E., & Rodgers, J. L. (007). Bootstrapping to test for nonzero population correlation coefficients using univariate sampling. Psychological Methods, (), -. doi: 0.07/08-989X... Beversdorf, L. M., & Sa, P. (0). Tests for correlation on bivariate nonnormal distributions. Journal of Modern Applied Statistical Methods, 0(), 699-709. Retrieved from http://digitalcommons.wayne.edu/jmasm/vol0/iss/9/ Bhattacharya, R. N., & Ghosh, J. K. (978). On the validity of the formal Edgeworth expansion. The Annals of Statistics, 6(), -5. Available from http://www.jstor.org/stable/958887 David, F. N. (98). Tables of the correlation coefficient. London: Cambridge University Press. Edgeworth, F. Y. (905). The law of error. Transactions of the Cambridge Philosophical Society, 0, 6-65, -. Fisher, R. A. (95). Frequency distribution of the values of the correlation coefficient in samples from an indefinitely large population. Biometrika, 0(), 507-5. doi: 0.07/88 Fisher, R. A. (9). On the probable error of a coefficient of correlation deduced from a small sample. Metron,, -. Fisher, R. A. (90). Moments and product moments of sampling distributions. Proceedings of the London Mathematical Society, s-0(), 99-8. doi: 0./plms/s-0..99 Fleishman, A. I. (978). A method of simulating non-normal distributions. Psychometrika, (), 5-5. doi: 0.007/BF098 Hall, P. (98). Inverting an Edgeworth expansion. The Annals of Statistics, (), 569-576. doi: 0./aos/7666 8

WANG & SA Hall, P. (99). The bootstrap and Edgeworth expansion. New York, NY: Springer-Verlag. Nakagawa, S., & Niki, N. (99). Distribution of sample correlation coefficient for nonnormal populations. Journal of the Japanese Society of Computational Statistics, 5(), -9. doi: 0.58/jjscs988.5. Pearson, K. (896). Mathematical contributions to the theory of evolution. III. Regression, heredity, and panmixia. Philosophical Transactions of the Royal Society of London. Series A, 87, 5-8. Available from http://www.jstor.org/stable/90707 Spearman, C. (90). General intelligence, Objectively determined and measured. The American Journal of Psychology, 5(), 0-9. doi: 0.07/07 Stuart, A., & Ord, J. K. (99). Kendall s advanced theory of statistics (6th ed.) (Vol. ). New York, NY: Halsted Press. Vale, D. C., & Maurelli, V. A. (98). Simulating multivariate non-normal distributions. Psychometrika, 8(), 65-7. doi: 0.007/BF09687 9

A NEW TEST FOR CORRELATION Appendix A: Formulae Used in Edgeworth Expansion of R Let 0 0. Then 8 0 0 0 0 0 0 0 0 00 8 9 60 0 5 0 0 0 0 0 0 8 8 6 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 8 5 9 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 8 9 0 0 0 0 0 0 0 0 0 5 6 0 0 6 0 0 0 0 0 0 0 0 6 0 0 0 0 5 5 6 0 0 0 0 0 0 0 0 0 8 0 0 7 5 6 0 0 5 0 0 0 0 0 9 0 0 0 06 0 0 0 0 0 8 6 0 0 8 50

WANG & SA Appendix B: k-statistics Fisher (90) introduced k-statistics as the unbiased estimator of the m th cumulant n m κ m, i.e. E(k m ) = κ m. Define the power sum of a univariate data as: sm x i i, the first six k-statistics in terms of the corresponding κ m are ( See Stuart & Ord, 99): k s, k ns s, n n k n s nss s, n k n n s n nss n ns nss 6 s, n k n n s n n s s n n s s n n s s n 5 5 5 0 0 5 5 +0n n ss 60nss s, 5 5 5 k6 n 6n n n s 6 6 6 n 6n n n s5s n 6 n n ss n n ss nss s 5n n n s s 0 n n 5n n s 0 n 9n n s s 0 n n s s s 0 n n n s 0 70 60 0 6 s 6 s s 5 s s 0 s 0 s s 0 s s s n 6 5 5 6 0 s 0 s s 70 s s 60 s s 0s 6 7 8 9 0 Stuart and Ord (99) also provided an approach to derive the multivariate k- n r t statistics. Define srt x i i y i, where (x i, y i ), i =,,, n are the bivariate random observations. The following multivariate k-statistics can be derived: k ns s 0s0 n k n s ns s ns s s s n 0 0 0 0 0 5

A NEW TEST FOR CORRELATION k n s nss0 ns0s0 s 0s0 n k n n s n n s s n n s s n n s s n 0 0 0 0 6ns s 6ns s s 6s s 0 0 0 0 0 0 k n n s n n s s n n s s n ns s n 0 0 0 0 6ns s 6ns s s 6s s 0 0 0 0 0 0 k n n s n n s s n n s s n n s n 0 0 n n s0s0 ns0s0 8nss0s0 ns0s0 6s0s0 k s s s s s s s s s s 6 n 0 s s 0 s s s 60 s s s 0 s s s 5 5 5 0 50 0 0 0 5 0 0 0s 5 0 5 0 0 0 6 0 0 6 0 0 0 s s s 0 s s 60 s s 60 s s s 6 0 0 0 7 0 8 0 8 0 0 0 80 s s s 90 s s s 0 s s 9 0 0 9 0 0 0 0 0 0 s s s 0s s 5 0 0 0 0 0 0 k s s s s s 6 s s 8 s s s s 6 n 6 s s s s s 6 s s s s s 0 0 0 0 0 0 5 0 5 0 0 5 0 0 6ss0s0 6ss0s0 8 6sss0 8 6s0s0s0 s s s 6 s s s s s 6 s s 6 0 0 6 0 0 7 0 7 0 0 s s 7 s s s s s s 7 s s 6 8 0 8 0 0 8 0 0 0 9 0 s s s s s s s 8 s s s s 9 0 0 0 9 0 0 0 9 0 0 0 0 0 9 s s s s s s 0s s 0 0 0 0 0 0 0 0 0 5

WANG & SA k s s s s s 6 s s 8 s s s s 6 n 6 s s s s s 6 s s s s s 0 0 0 0 0 0 5 0 5 0 0 5 0 0 6ss0s0 6ss0s0 8 6sss0 8 6s0s0s0 s s s 6 s s s s s 6 s s 6 0 0 6 0 0 7 0 7 0 0 s s 7 s s s s s s 7 s s 8 0 8 0 0 8 0 0 0 9 0 6 s s s s s s s 8 s s 9 0 0 0 9 0 0 0 9 0 0 9 s s s s s s 0s s 0 0 0 0 0 0 0 0 0 s s 0 0 0 k s s s s s s s s s 9 s s 6 n 9 s s s s 6 s s 8 s s s 6 s s 0 0 0 0 0 0 5 0 5 0 0 5 0 6 s s s 6 s s s 8 s s s 6 s s s 6 0 0 0 6 0 6 0 0 6 0 8 s s s 6 s s s 6 s s s s 6 0 0 6 0 0 0 7 7 0 0 6 s s 5 s s s 5 s s s 6 8 0 0 8 0 0 8 0 0 5 s s s 08 s s s 5 s s s s 9 0 0 9 0 0 9 0 0 0 0 5 s s s 7 s s s 7 s s s 9 0 0 0 0 0 0 0 0 0 0 6 s s s 0 s s 0 0 0 0 0 s s 8 0 0 where n n n n n 5 6 5 n n n n n n n n n n n n n n n n n n n n n n and 5

A NEW TEST FOR CORRELATION 5 5 6 7 8 9 n 6n n n n 6n n n n n 7n n n n 5n n n 9n n n n n n n n n n n 0 n 5

WANG & SA Appendix C: Tables Table. Comparison of type I error rates (ρ0 = 0, n = 0) α = 0.0 α = 0.05 α = 0.0 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 zα 0.0985 0.7 0. 0.09 0.0679 0.06 0.00 0.06 0.09 0.0 0 tα 0. 0.079 0.068 0.0580 0.07 0.05 0.0 zα 0.0987 0.68 0.5 0.050 0.0755 0.060 0.00 0.006 0.07 0.0 tα 0. 0.095 0.0697 0.0579 0.07 0.07 0.0 0 zα 0.0 0.5 0.06 0.059 0.0958 0.076 0.0 0.0588 0.008 0.0 0 tα 0.5 0.0998 0.0900 0.0 0.0570 0.0057 0.5 0 zα 0.0988 0.8 0.05 0.0508 0.067 0.0569 0.00 0.0 0.06 0.5 0 tα 0.06 0.0 0.0578 0.05 0.095 0.00.0 0 zα 0.05 0.0 0.095 0.0559 0.05 0.099 0.0 0.06 0.06.0 0 tα 0.0959 0.090 0.088 0.06 0.08 0.0098 0.5 zα 0.0999 0.5 0.7 0.056 0.0709 0.060 0.0 0.075 0.06 0.5 tα 0.68 0.06 0.069 0.058 0.0 0.08. zα 0.055 0.65 0.05 0.056 0.07 0.095 0.0 0.00 0.00. tα 0.09 0.0958 0.067 0.05 0.00 0.0075 5. 0 zα 0.08 0. 0.0 0.0580 0.097 0.069 0.056 0.059 0.0078. 0 tα 0.80 0.096 0.087 0.0 0.05 0.0056 Note: β: skewness; γ: kurtosis; CP: Critical Point 55

A NEW TEST FOR CORRELATION Table (continued). Comparison of type I error rates (ρ0 = 0, n = 0) α = 0.0 α = 0.05 α = 0.0 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc. 0 zα 0.08 0. 0.0 0.0580 0.097 0.069 0.056 0.059 0.0078. 0 tα 0.80 0.096 0.087 0.0 0.05 0.0056. 5 zα 0.06 0.508 0.09 0.059 0.0 0.06 0.0 0.075 0.006. 5 tα 0.56 0.085 0.0975 0.00 0.07 0.00 -. zα 0.050 0.6 0.09 0.056 0.070 0.08 0.06 0.0 0.00 -. tα 0.08 0.097 0.066 0.09 0.097 0.0076 -. 0 zα 0.05 0.6 0.0 0.056 0.095 0.056 0.059 0.05 0.0075 -. 0 tα 0.86 0.0960 0.085 0.09 0.055 0.005 -. 5 zα 0.0 0.507 0.090 0.0589 0.00 0.058 0.005 0.079 0.007 -. 5 tα 0.50 0.0855 0.0978 0.00 0.070 0.00. zα 0.057 0.75 0.0 0.0579 0.0860 0.090 0.05 0.06 0.0090. 0 tα 0.5 0.098 0.0797 0.09 0.06 0.0067. zα 0.06 0.98 0.0 0.056 0.097 0.09 0.055 0.0599 0.0065. 5 tα 0.7 0.0950 0.0888 0.07 0.058 0.005 -. zα 0.07 0.7 0.0 0.056 0.08 0.075 0.0 0.06 0.0087 -. 0 tα 0. 0.0970 0.078 0.08 0.05 0.006 5 -. zα 0.05 0.77 0.00 0.0575 0.0968 0.07 0.06 0.0590 0.0066 -. 5 tα 0. 0.090 0.0908 0.077 0.0576 0.0050 Note: β: skewness; γ: kurtosis; CP: Critical Point 56

WANG & SA Table (continued). Comparison of type I error rates (ρ0 = 0, n = 0) α = 0.0 α = 0.05 α = 0.0 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc. zα 0.09 0.70 0.0877 0.065 0.0779 0.07 0.006 0.0 0.0055.0 tα 0. 0.087 0.077 0.0 0.006 0.006. zα 0.00 0.75 0.095 0.06 0.086 0.09 0.08 0.088 0.0050.0 5 tα 0.7 0.086 0.080 0.06 0.07 0.00. zα 0.08 0.66 0.096 0.06 0.097 0.07 0.0 0.0656 0.005.0 5 tα 0. 0.0875 0.09 0.08 0.06 0.007 -. zα 0.7 0.98 0.0899 0.060 0.0760 0.075 0.006 0.05 0.005 -.0 tα 0. 0.087 0.07 0.08 0.00 0.00 -. zα 0.08 0.8 0.096 0.0607 0.08 0.080 0.088 0.0506 0.005 -.0 5 tα 0.8 0.0869 0.0789 0.08 0.090 0.00 -. zα 0.0 0.7 0.096 0.06 0.097 0.08 0.00 0.065 0.005 -.0 5 tα 0.9 0.0858 0.09 0.06 0.065 0.005.0 5 zα 0.8 0.0 0.08 0.0695 0.096 0.06 0.076 0.060 0.000.0 5 tα 0.5 0.078 0.089 0.087 0.066 0.00 5 -.0 5 zα 0. 0.7 0.087 0.0686 0.09 0.05 0.076 0.068 0.00 -.0 5 tα 0.8 0.076 0.0875 0.075 0.068 0.008 Note: β: skewness; γ: kurtosis; CP: Critical Point 57

A NEW TEST FOR CORRELATION Table. Comparison of type I error rates (n = 0, α = 0.05) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 z α 0.088 0.05 0.055 0.097 0.08 0.005 0.099 0.008 0.0007 0.096 0.0009 0.0000 0.0 0 t α 0.00 0.00 0.089 0.007 0.0085 0.0006 0.0009 0.0000 0.0 z α 0.057 0.057 0.0 0.055 0.06 0.009 0.059 0.08 0.005 0.05 0.00 0.000 0.0 t α 0.05 0.086 0.07 0.008 0.06 0.00 0.00 0.000 0.0 0 z α 0.097 0.77 0.050 0.05 0.65 0.00 0.80 0.07 0.05 0.6 0.07 0.005 0.0 0 t α 0.76 0.0 0.9 0.077 0.06 0.0 0.078 0.00 0.5 0 z α 0.0555 0.0 0.058 0.0550 0.00 0.0058 0.0567 0.006 0.00 0.0590 0.00 0.000 0.5 0 t α 0.00 0.0 0.005 0.005 0.0 0.000 0.00 0.000.0 0 z α 0.089 0.069 0.0 0.09 0.056 0.007 0.008 0.08 0.000 0.0078 0.007 0.0000.0 0 t α 0.056 0.0096 0.06 0.00 0.05 0.0009 0.007 0.0000 0.5 z α 0.0558 0.058 0.05 0.0558 0.066 0.0087 0.057 0.00 0.008 0.0580 0.009 0.000 0.5 t α 0.05 0.09 0.075 0.0079 0.09 0.006 0.0050 0.000. z α 0.0760 0.0998 0.086 0.0809 0.0776 0.098 0.0875 0.0597 0.008 0.09 0.0 0.006. t α 0.098 0.07 0.0790 0.079 0.06 0.0076 0.0 0.005. 0 z α 0.097 0.505 0.059 0.05 0.6 0.009 0.8 0.06 0.06 0.70 0.090 0.00. 0 t α 0.50 0.070 0.88 0.08 0.057 0.05 0.09 0.00 5. 5 z α 0.55 0.06 0.0770 0.75 0.80 0.05 0.00 0.566 0.07 0.68 0.096 0.07. 5 t α 0.065 0.0697 0.89 0.09 0.595 0.0 0.09 0.0 Note: β: skewness; γ: kurtosis; CP: Critical Point 58

WANG & SA Table (continued). Comparison of type I error rates (n = 0, α = 0.05) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. z α 0.076 0.098 0.076 0.080 0.0777 0.099 0.086 0.0609 0.0087 0.096 0.08 0.008 -. t α 0.0968 0.00 0.0790 0.08 0.06 0.008 0.09 0.007 -. 0 z α 0.0955 0.500 0.055 0.060 0.65 0.00 0. 0.07 0.070 0.6 0.085 0.0050 -. 0 t α 0.50 0.058 0.86 0.090 0.097 0.060 0.088 0.009 -. 5 z α 0.580 0.076 0.079 0.7 0.79 0.059 0.99 0.570 0.09 0.0 0.099 0.0 -. 5 t α 0.068 0.07 0.86 0.090 0.60 0.0 0.097 0.08. z α 0.087 0.67 0.05 0.089 0.8 0.00 0.086 0.0609 0.0087 0. 0.066 0.00. 0 t α 0.70 0.05 0.8 0.077 0.06 0.008 0.0667 0.00. z α 0.0996 0.5 0.0668 0. 0.5 0.050 0. 0.07 0.070 0.8 0.7 0.09. 5 t α 0.69 0.0600 0.0 0.08 0.097 0.060 0.7 0.0 -. z α 0.086 0.7 0.050 0.088 0.0 0.005 0.99 0.570 0.09 0.0 0.068 0.007 -. 0 t α 0.7 0.050 0.9 0.076 0.60 0.0 0.0685 0.00 -. z α 0.099 0. 0.0685 0. 0.8 0.080 0.0998 0.85 0.050 0.9 0.70 0.07 -. 5 t α 0.6 0.068 0.8 0.0 0. 0.00 0.7 0.0. z α 0.07 0.60 0.068 0.097 0. 0.09 0.60 0.7 0.08 0.7 0.089 0.008.0 t α 0.585 0.0560 0. 0.08 0.9 0.067 0.09 0.008 5. z α 0. 0.08 0.077 0.0 0.99 0.05 0.099 0.8 0.05 0.69 0.087 0.06.0 5 t α 0.085 0.0669 0.96 0.075 0.0 0.0 0.087 0.055 Note: β: skewness; γ: kurtosis; CP: Critical Point 59

A NEW TEST FOR CORRELATION Table (continued). Comparison of type I error rates (n = 0, α = 0.05) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc. z α 0.6 0.99 0.0757 0.5 0.79 0.058 0.6 0.096 0.09 0.960 0.0680 0.0.0 5 t α 0.98 0.068 0.79 0.080 0. 0.076 0.068 0.07 -. z α 0.06 0.565 0.06 0.089 0.98 0.005 0.5 0.099 0.0 0.7 0.09 0.008 -.0 t α 0.550 0.0558 0.0 0.069 0.096 0.099 0.096 0.008 -. z α 0.5 0.08 0.07 0.9 0.9 0.05 0.9 0.79 0.0 0.68 0.085 0.055 -.0 5 t α 0.08 0.065 0.9 0.08 0.7 0.05 0.087 0.08 -. z α 0.65 0.98 0.07 0.5 0.7 0.056 0.759 0.9 0.0 0.9 0.067 0.0 -.0 5 t α 0.976 0.0659 0.77 0.08 0.5 0.07 0.0676 0.06.0 5 z α 0.68 0.9 0.08 0.85 0.898 0.0588 0.7 0.09 0.06 0.90 0.0 0.056.0 5 t α 0.8 0.0755 0.9 0.057 0.0958 0.005 0.09 0.050 5 -.0 5 z α 0.65 0.86 0.08 0.80 0.9 0.0587 0.0 0.696 0.06 0.57 0.068 0.059 -.0 5 t α 0.7 0.075 0.98 0.05 0.7 0.008 0.07 0.05 Note: β: skewness; γ: kurtosis; CP: Critical Point 60

WANG & SA Table. Comparison of type I error rates (n = 0, α = 0.0) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 z α 0.007 0.05 0.00 0.007 0.07 0.007 0.00 0.07 0.000 0.00 0.00 0.0000 0.0 0 t α 0.075 0.007 0.097 0.00 0.05 0.000 0.00 0.0000 0.0 z α 0.05 0.067 0.0056 0.08 0.056 0.00 0.08 0.05 0.0008 0.0 0.009 0.000 0.0 t α 0.067 0.005 0.056 0.007 0.00 0.0007 0.0050 0.000 0.0 0 z α 0.050 0.768 0.060 0.07 0.65 0.05 0.05 0.6 0.007 0.06 0.050 0.00 0.0 0 t α 0.858 0.06 0.69 0.008 0.95 0.006 0.050 0.009 0.5 0 z α 0.05 0.059 0.006 0.06 0.085 0.007 0.05 0.058 0.000 0.08 0.007 0.0000 0.5 0 t α 0.086 0.000 0.008 0.00 0.067 0.000 0.007 0.0000.0 0 z α 0.08 0.07 0.005 0.0 0.05 0.00 0.008 0.00 0.000 0.00 0.009 0.0000.0 0 t α 0.09 0.000 0.065 0.000 0.05 0.000 0.0097 0.0000 0.5 z α 0.0 0.06 0.005 0.0 0.05 0.00 0.08 0.00 0.0009 0.05 0.0057 0.000 0.5 t α 0.065 0.00 0.05 0.006 0.05 0.0008 0.0057 0.000. z α 0.06 0. 0.0 0.058 0.0 0.008 0.096 0.075 0.00 0.08 0.06 0.0009. t α 0.75 0.0095 0.060 0.0067 0.0778 0.007 0.06 0.0008. 0 z α 0.078 0.75 0.07 0.0 0.597 0.06 0.058 0.59 0.008 0.06 0.058 0.008. 0 t α 0.8 0.0 0.66 0.008 0.9 0.007 0.050 0.005 5. 5 z α 0.08 0.6 0.09 0.0989 0.68 0.060 0.69 0.86 0.087 0.8 0.00 0.0099. 5 t α 0.0 0.05 0.8 0.08 0.866 0.069 0.0 0.009 Note: β: skewness; γ: kurtosis; CP: Critical Point 6

A NEW TEST FOR CORRELATION Table (continued). Comparison of type I error rates (n = 0, α = 0.0) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. z α 0.07 0.6 0.00 0.06 0.00 0.0075 0.000 0.079 0.009 0.06 0.0 0.0008 -. t α 0.7 0.009 0.06 0.006 0.077 0.00 0.0 0.0006 -. 0 z α 0.08 0.780 0.075 0.09 0.605 0.0 0.05 0. 0.0080 0.068 0.059 0.008 -. 0 t α 0.860 0.06 0.676 0.0 0.70 0.007 0.05 0.005 -. 5 z α 0.08 0.9 0.095 0.0977 0.69 0.058 0.06 0.87 0.096 0.6 0.0997 0.009 -. 5 t α 0.7 0.05 0.0 0.06 0.878 0.078 0.00 0.0085. z α 0.086 0.70 0.065 0.08 0.676 0.00 0.06 0.0 0.0069 0.06 0.0689 0.00. 0 t α 0.89 0.06 0.75 0.00 0. 0.006 0.069 0.000. z α 0.000 0.75 0.07 0.057 0.68 0.00 0.0567 0.76 0.05 0.0995 0.77 0.006. 5 t α 0.89 0.006 0.77 0.07 0.07 0.06 0.78 0.0057 -. z α 0.08 0.70 0.06 0.09 0.67 0.05 0.07 0.7 0.007 0.06 0.0698 0.00 -. 0 t α 0.8 0.0 0.75 0.005 0.8 0.0060 0.0699 0.000 -. z α 0.090 0.699 0.0 0.05 0.660 0.099 0.0575 0.7 0.00 0.00 0.8 0.0067 -. 5 t α 0.8 0.00 0.75 0.07 0.00 0.0 0.85 0.0059. z α 0.0 0.707 0.0 0.06 0.560 0.066 0.08 0.080 0.008 0.088 0.008 0.0059.0 t α 0.777 0.077 0.6 0.0 0.0 0.009 0.0 0.005 5. z α 0.068 0.69 0.075 0.05 0.79 0.0 0.0605 0.95 0.08 0.06 0.0865 0.0089.0 5 t α 0.6 0.0 0.5 0.00 0.96 0.06 0.0866 0.0078 Note: β: skewness; γ: kurtosis; CP: Critical Point 6

WANG & SA Table (continued). Comparison of type I error rates (n = 0, α = 0.0) ρ 0 = 0.50 ρ 0 = 0.60 ρ 0 = 0.75 ρ 0 = 0.90 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc. z α 0.068 0. 0.090 0.078 0.079 0.08 0.090 0.65 0.087 0.06 0.0709 0.0090.0 5 t α 0. 0.05 0.5 0.09 0.657 0.065 0.07 0.008 -. z α 0.05 0.756 0.0 0.09 0.566 0.07 0.079 0.088 0.007 0.0809 0.005 0.006 -.0 t α 0.86 0.09 0.66 0.050 0.08 0.009 0.007 0.0059 -. z α 0.066 0.65 0.069 0.05 0.06 0.08 0.059 0.9 0.07 0.0607 0.0880 0.008 -.0 5 t α 0.5 0.05 0.75 0.009 0.9 0.050 0.088 0.0075 -. z α 0.067 0. 0.08 0.0779 0.0 0.08 0.099 0.6 0.079 0.07 0.0708 0.009 -.0 5 t α 0.9 0.0 0.098 0.08 0.6 0.06 0.07 0.0085.0 5 z α 0.09 0.9 0.05 0.058 0. 0.095 0.9 0.87 0.07 0.0 0.080 0.0098.0 5 t α 0.50 0.000 0.85 0.06 0.9 0.097 0.085 0.0088 5 -.0 5 z α 0.096 0.75 0.05 0.060 0.96 0.090 0. 0.87 0.006 0.06 0.09 0.0098 -.0 5 t α 0.6 0.00 0.6 0.05 0.907 0.088 0.00 0.0090 Note: β: skewness; γ: kurtosis; CP: Critical Point 6

A NEW TEST FOR CORRELATION Table. Power performance for test ρ0 = 0 (n = 0, α = 0.05) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 0.098 0.068 0.065 0.778 0.8 0.00 0.7 0.776 0.756 0.9775 0.9785 0.977.0000.0000 0.999 0.0 0 0.06 0.0576 0.976 0.87 0.79 0.6950 0.979 0.967.0000 0.999 0.0 0.09 0.079 0.066 0.8 0. 0.96 0.77 0.789 0.6980 0.9776 0.9798 0.956.0000.0000 0.980 0.0 0.068 0.0565 0. 0.778 0.7 0.6758 0.9768 0.97 0.9999 0.985 0.0 0 0.055 0.099 0.085 0.00 0.89 0.6 0.76 0.79 0.66 0.9809 0.987 0.87.0000 0.9998 0.857 0.0 0 0.09 0.07 0.65 0.9 0.777 0.68 0.98 0.80 0.9997 0.88 0.5 0 0.050 0.065 0.0585 0.8 0.07 0.97 0.70 0.7 0.7089 0.9766 0.977 0.9680.0000 0.9999 0.9977 0.5 0 0.056 0.058 0.888 0.78 0.77 0.688 0.978 0.96 0.9998 0.9969.0 0 0.055 0.0550 0.05 0.655 0.5 0.77 0.659 0.678 0.597 0.909 0.9 0.97 0.9988 0.998 0.9957.0 0 0.09 0.057 0.50 0.9 0.609 0.576 0.95 0.90 0.9978 0.998 0.5 0.05 0.070 0.060 0.8 0.5 0.907 0.709 0.75 0.695 0.976 0.978 0.997.0000 0.9999 0.9808 0.5 0.0650 0.058 0.06 0.79 0.76 0.67 0.975 0.99 0.9999 0.978. 0.0570 0.07 0.050 0.97 0.8 0.57 0.78 0.77 0.65 0.9755 0.9765 0.895.0000 0.999 0.9. 0.0686 0.05 0.090 0.68 0.786 0.69 0.978 0.8809 0.999 0.96. 0 0.0565 0.099 0.067 0.06 0.76 0.568 0.75 0.78 0.6 0.9798 0.9808 0.80.0000 0.9997 0.85. 0 0.086 0.008 0.58 0. 0.767 0.6059 0.9780 0.895 0.9996 0.8. 5 0.0597 0.0 0.059 0. 0.09 0.9 0.7 0.7969 0.5898 0.976 0.9780 0.766 0.9999 0.9997 0.766. 5 0.0988 0.0 0.8 0.7 0.780 0.559 0.97 0.770 0.9996 0.759 -. 0.0556 0.075 0.09 0.96 0.7 0.56 0.708 0.77 0.65 0.9755 0.976 0.89.0000 0.9995 0.9 -. 0.0666 0.06 0.078 0.6 0.78 0.68 0.979 0.8809 0.999 0.9 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 6

WANG & SA Table (continued). Power performance for test ρ0 = 0 (n = 0, α = 0.05) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. 0 0.0578 0.095 0.057 0.050 0.7 0.57 0.78 0.789 0.65 0.9790 0.980 0.8.0000 0.9995 0.85 -. 0 0.0868 0.099 0.5 0. 0.760 0.609 0.977 0.805 0.999 0.8 -. 5 0.0597 0.0 0.06 0.7 0.07 0.88 0.70 0.797 0.5898 0.97 0.9775 0.7606 0.9999 0.9997 0.7667 -. 5 0.0986 0.007 0.86 0.9 0.780 0.5588 0.979 0.70 0.9996 0.75. 0.0565 0.087 0.07 0.0 0.557 0.65 0.759 0.770 0.65 0.9799 0.980 0.87.0000 0.9996 0.890. 0 0.0768 0.0 0.5 0.96 0.755 0.678 0.977 0.8600 0.9995 0.88. 0.0577 0.0956 0.08 0.5 0.8 0.569 0.765 0.807 0.60 0.986 0.9867 0.8.0000 0.9999 0.878. 5 0.0897 0.078 0.67 0. 0.787 0.65 0.987 0.806 0.9999 0.867 -. 0.056 0.08 0.089 0.957 0.508 0.590 0.76 0.7677 0.66 0.980 0.98 0.877.0000 0.9997 0.8895 -. 0 0.078 0.0 0. 0.7 0.79 0.69 0.978 0.8608 0.9996 0.88 -. 0.057 0.095 0.0 0.9 0.89 0.559 0.767 0.80 0.6 0.9858 0.9868 0.86.0000 0.9999 0.888 -. 5 0.0889 0.070 0.656 0.6 0.7867 0.65 0.986 0.87 0.9999 0.889. 0.06 0.0766 0.07 0.067 0.6 0.7 0.75 0.7 0.560 0.9708 0.97 0.80 0.9999 0.9996 0.88.0 0.0708 0.0 0.08 0.98 0.7 0.59 0.9688 0.7957 0.9995 0.875. 0.068 0.089 0.08 0.8 0.59 0.56 0.779 0.77 0.60 0.980 0.980 0.88.0000 0.9997 0.87.0 5 0.079 0.0 0.97 0. 0.7557 0.5806 0.9775 0.8087 0.9996 0.80. 0.069 0.0976 0.08 0.00 0.859 0.9 0.76 0.789 0.600 0.9769 0.9785 0.7878.0000 0.9996 0.789.0 5 0.098 0.00 0.665 0.5 0.776 0.5698 0.979 0.775 0.9995 0.778 -. 0.068 0.0766 0.069 0.09 0.97 0.5 0.77 0.77 0.56 0.97 0.978 0.8 0.9999 0.9995 0.877 -.0 0.075 0.07 0.07 0.96 0.7 0.55 0.9698 0.7967 0.999 0.866 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 65

A NEW TEST FOR CORRELATION Table (continued). Power performance for test ρ0 = 0 (n = 0, α = 0.05) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. 0.06 0.0855 0.088 0.86 0.60 0.88 0.78 0.7755 0.6 0.9806 0.980 0.85.0000 0.9997 0.85 -.0 5 0.080 0.0 0.0 0.59 0.7568 0.58 0.977 0.807 0.9996 0.86 -. 0.068 0.0967 0.078 0.00 0.855 0.95 0.759 0.7878 0.5980 0.9776 0.9785 0.7889.0000 0.9996 0.79 -.0 5 0.098 0.0 0.66 0.6 0.770 0.568 0.975 0.778 0.9996 0.78.0 5 0.0689 0.095 0.0 0.89 0.66 0.9 0.7 0.798 0.556 0.965 0.969 0.758 0.9999 0.999 0.76.0 5 0.0889 0.079 0.89 0.98 0.7 0.59 0.957 0.77 0.9990 0.7509 5 6 -.0 5 0.068 0.099 0.0 0.9 0.686 0.6 0.757 0.7508 0.5 0.965 0.969 0.75 0.9998 0.999 0.7676 -.0 5 0.0875 0.086 0.50 0.96 0.77 0.5 0.957 0.776 0.999 0.7556 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 66

WANG & SA Table 5. Power performance for test ρ0 = 0.00 (n = 0, α = 0.0) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 0.00 0.07 0.00 0.07 0.508 0.6 0.585 0.88 0.5 0.908 0.8966 0.889 0.9996 0.9995 0.987 0.0 0 0.05 0.06 0.99 0.095 0.05 0.908 0.860 0.80 0.999 0.9797 0.0 0.00 0.09 0.080 0.06 0.57 0.85 0.66 0.5000 0.06 0.9090 0.908 0.798 0.9996 0.9995 0.9 0.0 0.08 0.07 0.8 0.096 0.5 0.9 0.878 0.70 0.9990 0.96 0.0 0 0.07 0.0605 0.0087 0.7 0.06 0.077 0.89 0.5508 0.90 0.98 0.90 0.69 0.9997 0.9988 0.7 0.0 0 0.0586 0.006 0.80 0.055 0.96 0.8 0.89 0.5508 0.998 0.68 0.5 0 0.0 0.06 0.080 0.07 0.0 0.5 0.59 0.7 0.5 0.906 0.895 0.87 0.9996 0.999 0.975 0.5 0 0.09 0.00 0.08 0.00 0.0 0.8 0.8597 0.796 0.998 0.967.0 0 0.09 0.069 0.00 0.05 0.096 0.0850 0.05 0.5 0.0 0.897 0.757 0.679 0.999 0.98 0.998.0 0 0.0 0.00 0.079 0.068 0.8 0.5 0.678 0.606 0.97 0.9098 0.5 0.05 0.086 0.07 0.07 0.8 0. 0.6 0.886 0.00 0.9055 0.8997 0.7877 0.9995 0.9990 0.979 0.5 0.09 0.00 0.56 0.09 0. 0.5 0.8665 0.76 0.9986 0.907. 0.08 0.08 0.00 0.0 0.5 0.08 0.7 0.806 0.7 0.90 0.8896 0.686 0.9996 0.997 0.89. 0.05 0.008 0.0 0.067 0.7 0.66 0.85 0.6 0.996 0.79. 0 0.06 0.05 0.0076 0.0 0.98 0.0676 0.9 0.50 0.85 0.9086 0.9087 0.66 0.9997 0.9980 0.799. 0 0.057 0.005 0.706 0.096 0.80 0.5 0.878 0.57 0.997 0.689. 5 0.05 0.07 0.007 0.59 0.5 0.078 0.5077 0.569 0. 0.8869 0.897 0.505 0.9989 0.9977 0.6. 5 0.07 0.00 0.8 0.06 0.555 0.60 0.8670 0.69 0.9965 0.5767 -. 0.05 0.08 0.00 0.08 0.50 0.088 0.7 0.78 0.07 0.90 0.8890 0.68 0.9996 0.997 0.800 -. 0.00 0.0076 0.00 0.060 0.8 0.69 0.85 0.6 0.9960 0.799 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 67

A NEW TEST FOR CORRELATION Table 5 (continued). Power performance for test ρ0 = 0.00 (n = 0, α = 0.0) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. 0 0.06 0.058 0.0077 0.89 0.90 0.0665 0.9 0.57 0.8 0.908 0.9079 0.609 0.9997 0.9980 0.779 -. 0 0.05 0.005 0.695 0.08 0.787 0.6 0.8785 0.59 0.997 0.6858 -. 5 0.0 0.077 0.008 0.578 0.8 0.079 0.50 0.565 0.5 0.8888 0.898 0.5067 0.9988 0.9977 0.608 -. 5 0.070 0.00 0. 0.06 0.57 0.67 0.8695 0.97 0.996 0.575. 0.056 0.07 0.009 0.9 0.70 0.0759 0.80 0.5 0.08 0.9 0.9057 0.65 0.9998 0.9980 0.77. 0 0.0 0.0065 0.5 0.0568 0.598 0.87 0.876 0.59 0.997 0.77. 0.060 0.0586 0.0065 0.8 0.0 0.06 0.58 0.5580 0.80 0.998 0.90 0.669.0000 0.999 0.6999. 5 0.057 0.00 0.807 0.05 0.50 0.9 0.9057 0.556 0.999 0.668 -. 0.0 0.056 0.0089 0.5 0.78 0.0757 0.85 0.5 0.067 0.9 0.906 0.655 0.9998 0.9978 0.770 -. 0 0.0 0.006 0.5 0.0560 0.588 0.77 0.8759 0.5899 0.9969 0.706 -. 0.057 0.058 0.006 0. 0.09 0.067 0.50 0.5586 0.85 0.908 0.9 0.670.0000 0.999 0.700 -. 5 0.0567 0.00 0.80 0.0 0.500 0. 0.907 0.555 0.999 0.668. 0.007 0.0 0.0050 0.5 0.6 0.050 0.85 0.60 0.9 0.888 0.8708 0.59 0.999 0.996 0.698.0 0.09 0.00 0.9 0.050 0.09 0.699 0.896 0.796 0.996 0.660. 0.086 0.050 0.0056 0. 0.796 0.055 0.5095 0.58 0.6 0.99 0.90 0.5867 0.9999 0.9986 0.6977.0 5 0.08 0.006 0.599 0.08 0.6 0.09 0.87 0.55 0.9978 0.666. 0.06 0.066 0.005 0.507 0.5 0.05 0.5075 0.55 0. 0.8979 0.898 0.59 0.9995 0.9978 0.659.0 5 0.066 0.006 0.90 0.08 0.955 0.89 0.868 0.76 0.9969 0.6 -. 0.09 0.009 0.005 0. 0.6 0.086 0.8 0.66 0. 0.8886 0.87 0.595 0.999 0.9965 0.690 -.0 0.09 0.005 0.8 0.00 0.06 0.699 0.88 0.8 0.997 0.6609 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 68

WANG & SA Table 5 (continued). Power performance for test ρ0 = 0.00 (n = 0, α = 0.0) ra = 0.0 ra = 0. ra = 0. ra = 0.6 ra = 0.8 b g Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc Zf Zb Zc -. 0.09 0.050 0.0055 0. 0.806 0.0550 0.5090 0.58 0.6 0.9 0.906 0.5878 0.9999 0.998 0.696 -.0 5 0.08 0.006 0.600 0.085 0.6 0.06 0.876 0.556 0.9976 0.666 -. 0.06 0.0659 0.005 0.99 0.8 0.055 0.5060 0.58 0.06 0.900 0.905 0.56 0.999 0.9978 0.678 -.0 5 0.068 0.005 0.9 0.058 0.956 0.88 0.87 0.775 0.9968 0.66.0 5 0.077 0.065 0.000 0.70 0.099 0.0 0.5080 0.567 0.08 0.87 0.858 0.780 0.9978 0.995 0.659.0 5 0.06 0.005 0.99 0.09 0.70 0.56 0.89 0. 0.99 0.589 5 6 -.0 5 0.08 0.066 0.00 0.79 0.0 0.0 0.5079 0.576 0.0 0.875 0.8589 0.80 0.9978 0.996 0.6 -.0 5 0.06 0.008 0.98 0.080 0.7 0.5 0.89 0.5 0.99 0.585 Note: β: skewness; γ: kurtosis; the Zf, Zb, and Zc results are calculated using the critical points z α and t α as the first and the second number 69

A NEW TEST FOR CORRELATION Table 6 Power performance for test ρ0 = 0.55 (n = 0, α = 0.0) ρα = 0.55 ρα = 0.60 ρα = 0.70 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc 0.0 0 zα 0.099 0.75 0.9 0.005 0.08 0.86 0.658 0.908 0.779 0.0 0 tα 0.5 0.0978 0.6 0.5 0.0 0.07 0.0 zα 0.07 0.5 0. 0.990 0.06 0.795 0.656 0.850 0.678 0.0 tα 0.9 0.0999 0.79 0.5 0.6 0.0 0.0 0 zα 0.8 0.9 0.58 0.5 0.97 0.9 0.660 0.896 0.6 0.0 0 tα 0.855 0.5 0.96 0.690 0.50 0.90 0.5 0 zα 0.075 0. 0.6 0.058 0.07 0.85 0.696 0.99 0.789 0.5 0 tα 0.5 0.005 0.66 0.75 0.5 0..0 0 zα 0.09 0.0 0.098 0.0 0.50 0.56 0.99 0.0 0.00.0 0 tα 0.098 0.08 0.095 0.90 0.756 0.56 0.5 zα 0.06 0.55 0. 0.06 0.06 0.798 0.657 0.870 0.698 0.5 tα 0.0 0.0997 0.77 0. 0.59 0.08. zα 0.80 0.656 0.59 0.6 0.55 0.8 0.69 0.795 0.556. tα 0.556 0. 0.0 0.57 0. 0.057. 0 zα 0.508 0.968 0.87 0.9 0.58 0.969 0.6560 0.85 0.56. 0 tα 0.889 0.7 0.8 0.77 0.66 0. 5. 5 zα 0. 0.88 0.6 0.977 0.76 0.0 0.655 0.797 0.07. 5 tα 0. 0.5 0.57 0.87 0.57 0. Note: β: skewness; γ: kurtosis; CP: Critical Point 70

WANG & SA Table 6 (continued). Power performance for test ρ0 = 0.55 (n = 0, α = 0.0) α = 0.0 α = 0.05 α = 0.0 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc -. zα 0.95 0.656 0.66 0. 0.6 0.86 0.687 0.78 0.55 -. tα 0.558 0.8 0.998 0.550 0.99 0.06 -. 0 zα 0.57 0.97 0.9 0.68 0.55 0.967 0.655 0.886 0.59 -. 0 tα 0.899 0.8 0. 0.76 0.9 0.76 -. 5 zα 0.6 0.75 0.60 0.95 0.696 0.989 0.6558 0.755 0.79 -. 5 tα 0.8 0.509 0.566 0.80 0.509 0.097. zα 0.59 0.958 0.7 0.57 0.7 0.0 0.665 0.55 0.9. 0 tα 0.87 0.55 0.9 0.8 0.80 0.65. zα 0.600 0.97 0.59 0.75 0.87 0.9 0.79 0.57 0.58. 5 tα 0.8 0.80 0.6 0.55 0.586 0.97 -. zα 0.58 0.95 0.7 0.66 0.70 0.095 0.665 0.5 0.90 -. 0 tα 0.879 0.66 0.86 0.88 0.789 0.58 -. zα 0.570 0.59 0.557 0.76 0.0 0.58 0.76 0.575 0.5 -. 5 tα 0. 0.7 0.9 0.87 0.58 0.955. zα 0.55 0.99 0. 0.8 0.569 0.005 0.680 0.908 0.66.0 tα 0.97 0.0 0.95 0.798 0.55 0.88 5. zα 0.7 0.5 0.66 0.790 0.097 0. 0.6958 0.56 0.50.0 5 tα 0.78 0.55 0.97 0.6 0.5009 0.6 Note: β: skewness; γ: kurtosis; CP: Critical Point 7

A NEW TEST FOR CORRELATION Table 6 (continued). Power performance for test ρ0 = 0.55 (n = 0, α = 0.0) α = 0.0 α = 0.05 α = 0.0 β γ CP Zf Zb Zc Zf Zb Zc Zf Zb Zc. zα 0.9 0. 0.578 0.8 0.659 0.00 0.668 0.698 0.76.0 5 tα 0.78 0.79 0.50 0.806 0.5 0.06 -. zα 0.586 0.05 0.59 0.0 0.566 0.006 0.678 0.99 0.685 -.0 tα 0.990 0.6 0.95 0.80 0.58 0.5 -. zα 0.700 0. 0.599 0.790 0.07 0. 0.6958 0.5 0.995 -.0 5 tα 0.65 0.98 0.90 0.09 0.98 0.6 -. zα 0.955 0.6 0.588 0.865 0.68 0.99 0.660 0.66 0.0 -.0 5 tα 0.70 0.8 0.95 0.805 0.9 0.99.0 5 zα 0. 0.57 0.699 0.995 0.760 0.00 0.60 0.70 0.9.0 5 tα 0.7 0.60 0.60 0.8 0.89 0.009 -.0 5 zα 0.0 0.5 0.685 0.996 0.776 0.0 0.6 0.7 0.8 -.0 5 tα 0.6 0.58 0.6 0.8 0.80 0.00 Note: β: skewness; γ: kurtosis; CP: Critical Point 7