Research Article  Open Access
Chioneso Show Marange, Yongsong Qin, "An Empirical Likelihood RatioBased Omnibus Test for Normality with an Adjustment for Symmetric Alternatives", Journal of Probability and Statistics, vol. 2021, Article ID 6661985, 18 pages, 2021. https://doi.org/10.1155/2021/6661985
An Empirical Likelihood RatioBased Omnibus Test for Normality with an Adjustment for Symmetric Alternatives
Abstract
An omnibus test for normality with an adjustment for symmetric alternatives is developed using the empirical likelihood ratio technique. We first transform the raw data via a jackknife transformation technique by deleting one observation at a time. The probability integral transformation was then applied on the transformed data, and under the null hypothesis, the transformed data have a limiting uniform distribution, reducing testing for normality to testing for uniformity. Employing the empirical likelihood technique, we show that the test statistic has a chisquare limiting distribution. We also demonstrated that, under the established symmetric settings, the CUSUMtype and Shiryaev–Roberts test statistics gave comparable properties and power. The proposed test has good control of type I error. Monte Carlo simulations revealed that the proposed test outperformed studied classical existing tests under symmetric shorttailed alternatives. Findings from a real data study further revealed the robustness and applicability of the proposed test in practice.
1. Introduction
The empirical likelihood (EL) methodology was introduced in [1, 2] and has been widely studied as a nonparametric approximation of the parametric likelihood approach (e.g., [3–6]). Thus, it utilizes the concept of the likelihoods in a distributionfree manner in approximating optimal parametric likelihoodbased techniques. The method provides a versatile approach that may be applied to perform inference for a wide variety of statistical applications. An area with substantial new development in the use of the EL methods is hypothesis testing. Various researchers have proposed goodnessoffit (GoF) tests for continuous distributions based on the EL for a wide range of hypothesis tests, which includes exponentiality [7, 8], logistic [9], uniformity [10], and normality [11, 12].
From the various proposed EL testing procedures as well as in the current statistical practice, it is evident that the problem of testing composite hypotheses of normality is undeniably the most common research focus in GoF testing. The continued growing need for normality tests is attributed to the frequent use and applications of normally distributed data in various areas of pure and applied statistical practices. Although it is difficult to propose a test for normality competing with the highly efficient family of Shapiro–Wilk tests (e.g., [13–16]), the proposed ELbased normality tests have proved to be superior under certain alternative distributions [12]. Of these tests, the momentbased tests seem to have gained more traction due to their flexibility, simplicity, power properties, and convenient use of omnibus tests in accessing the normality of underlying continuous distributions.
To test for normality, Dong and Giles [11] proposed an omnibus test statistic by directly utilizing the EL methodology outlined by Owen [17]. They utilized the first four moment constraints that characterize the normal distribution. After outlining drawbacks of the test proposed by Dong and Giles [11], Shan et al. [12] proposed a cumulative sum (CUSUM) type simple and exact empirical likelihood ratiobased (SEELR) test statistic for normality which unlike that of Dong and Giles [11] has good control of type I error (also see [18]) and can be easily implemented in a wide range of statistical packages. The test by Shan et al. [12] is an omnibus test that makes use of standardized sample observations using the Lin and Mudholkar [19] jackknife transformation. In their study, Shan et al. [12] reported that power of their proposed omnibus test is comparable to wellknown existing tests and oftentimes outperforms these tests under certain alternatives, mostly asymmetric distributions.
Just like some tests for normality, the test proposed by Shan et al. [12] suffers the loss of power under several symmetric alternatives. It is a challenge to propose an omnibus test that has high power than the classical Jarque–Bera’s tests [20, 21] and the D’Agostino–Pearson test [22] in detecting departures from normality in alternatives that exhibit the symmetric nature of the normal distribution. Through the utilization of various mathematical and statistical properties that characterize the normal distribution (for example, see [19, 23, 24]) one can remedy such shortfalls in GoF tests. One such remedy is transformation to uniformity, which has several benefits that include increasing the power of a test under certain alternatives (for example, see [8, 25]). For a datadriven omnibus test for symmetry, Fang et al. [25] utilized a bootstrapping approach coupled with the probability integral transformation, and under the null hypothesis, the transformed data had a limiting uniform distribution. For superior power under symmetric alternatives, their proposed test required only oddordered orthogonal moments of the transformed data in constructing the test statistic.
The use of the probability integral transformation in the development of GoF tests of normality has been widely used especially in empirical distribution function (EDF) based tests. Rosenblatt [26] first introduced the concept. Thus, the EDF tests make use of the probability integral transformation . If is the distribution function of , the random variable is uniformly distributed between 0 and 1. Given observations , the values are computed. The most commonly used EDF tests for normality are the Anderson–Darling test [27, 28] and the Lilliefors test which is well known as the modified Kolmogorov–Smirnov test [29] and the Cramér–von Mises test [30]. In addition to the use of the probability integral transformation, several approaches have been used to construct GoF tests for the composite hypothesis of normality. In this study, we adopted the EL methodology to propose a new omnibus test for normality by exploiting different forms of characterizing the normal distribution. The purpose of this paper is to use a jackknife characterization due to Shan et al. [12], as in Lin and Mudholkar [19], followed by a probability integral transformation (see [8, 25]) for developing a goodnessoffit test for normality. Here we consider the approach to obtaining a GoF test statistic by combining two wellknown characterizations, individually powerful against different classes of alternatives. However, following the works of Fang et al. [25], we restrict attention to symmetric alternatives. Power comparisons are conducted with some of the most widely known EDFbased tests, wellknown and powerful momentbased tests, and the powerful classical SW tests.
2. Test Development
Consider an unknown continuous distribution with nonordered random variables denoted by that are assumed to be independent and identically distributed . The intention is to test whether the observations are consistent with a normal distribution. Thus, we intended to test whether to accept or reject the following null hypothesis:where and are unknown parameters. We then proposed to use standardized random variables of the sample observations. To achieve this, we adopted a jackknife transformation technique by deleting one observation at a time following Lin and Mudholkar [19] works (also see [12]). Thus, we transformed our observations usingwhere , , and . It should be noted according to Shan et al. [12] as gets large, the standardized data points become asymptotically independent while under the null hypothesis they are distributed according to a distribution with df, which as grows approaches the standard normal. In addition to this transformation, we then further adopted the probability integral transformation (see [25] as well as [8]). The probability integral transformation then transformed the standardized random variables into independent uniformly distributed random variables . That is, under the null hypothesis, the transformed data follow the uniform distribution asymptotically. From the proposed transformation, are uniformly distributed on , where the density function of is given bywhere = the lowest value of and = the highest value of . The moment of the uniform distribution is defined by
From the uniformly transformed observations, we then proposed to test for the following null hypothesis:versus the alternative that are from a nonuniform distribution defined on (0, 1). In accordance with the EL method, under , we considered unbiased empirical moment equations utilizing the raw moments of the uniform distribution, which are given by
The composite hypothesis for the ELR test was then given by
The nonparametric empirical likelihood function corresponding to the given hypotheses in equation (7) is expressed aswhere the unknown probability parameters ’ are attained under and . Under , the EL function given in equation (8) is maximized with respect to the ’ subject to two constraints:
Following this, the weights of ’ are identified aswhere , for . If we then use the Lagrangian multipliers technique, it can be shown that the maximum EL function under can be expressed aswhere in equation (11) is a root of
Under the alternative hypothesis, is not required to identify the weights, in order to maximize the EL function but only, . Thus, under , the nonparametric EL function is given by
Now let us consider to be −2 loglikelihood test statistic for the hypotheses . It should be noted that, under , minus two times the log LLR has an asymptotic limiting distribution [1]. Considering the null and alternative hypotheses, the test statistic is given by
We then proposed to reject the null hypothesis using two different test statistics. Firstly, we considered the cumulative sum (CUSUM) type statistic given by
Secondly, we considered the common alternative to the CUSUMtype statistic, which is to utilize the Shiryaev–Roberts (SR) statistic (for example, see [31] among others). In our case, the classical SR statistic was of the formwhere is the test threshold and is percentile of the distribution. The set are integer values representing the moment constraints that will maximize the test statistic. Our proposed test statistic (equation (15)) is developed utilizing approaches introduced by Vexler and Wu [32], and various authors have demonstrated that the SR statistic (see equation (16)) and the CUSUMtype statistic have almost equivalent optimal statistical properties due to their common nullmartingale basis [31]. The choice of is also vital in momentbased test statistics. Fang et al. [25] utilized the probability integral transformation and recommended that, under the test of symmetry, only oddordered moments of the transformed data are required in the construction of the test statistic. They further alluded that the use of oddordered moments has several benefits that include power against most symmetric alternatives and robustness and performs well under small sample sizes among others.
For our proposed test statistics, we decided to then conduct an extensive Monte Carlo simulation exercise in order to empirically evaluate the suitable choice of that will give us optimal power under symmetric alternatives in testing for normality. Following the work of Fang et al. [25] as well as Shan et al. [12], we estimated the powers of the test statistics for different alternatives and definitions of oddordered moments of . Table 1 displays a subset of the Monte Caro simulation results. We also considered additional alternatives based on samples of sizes n = 20, 50, and 100 at . We used sizeadjusted critical values for each test statistic, and power for each test was computed using 5,000 replications. The results were that both the CUSUMtype and the Shiryaev–Roberts proposed test statistics with showed an average power that was greater than that of all other cases under symmetric shorttailed alternatives. In addition, the CUSUMtype test statistics with showed an average power that was greater than that of all other cases under symmetric longtailed alternatives. Since our proposed test is meant to perform superior under symmetric alternatives, we recommend to use since the proposed test statistics with are simple and provide relatively high levels of power under symmetric alternatives. The results of our Monte Carlo simulation experiments are consistent with that of Fang et al. [25]. In this article, we denoted CSELR for the CUSUMtype test statistic and SRELR for the Shiryaev–Roberts proposed test statistic. A schematic algorithm of the testing procedure is shown in Figure 1.
 
Bold represents the most superior CStype statistic and italicized represents the most superior SR statistic. 
3. Monte Carlo Simulation Procedures
We utilized the R statistical package for all the simulation procedures. Firstly, sizeadjusted critical values for the proposed test statistics were determined. In order to achieve this, we used 50,000 replications, and without loss of generality, data were simulated from a standard normal distribution at stipulated sample sizes and levels. Only samples of sizes 20 to 100 were considered (see Table 2). This was entirely motivated by the need to utilize samples that commonly arise in practice.

For power comparisons, twelve selected competitor tests were considered. The choice of these tests was guided by potential competitor tests, thus tests developed using similar characterization techniques as well as wellknown powerful classical normality tests. Three broad categories of these competitor tests were established. These include EDFbased tests, momentbased tests, and the Shapiro–Wilkbased tests. For the EDFbased tests, we opted for the Anderson–Darling (AD) test [27, 28], the Lilliefors (LL) test which is well known as the modified Kolmogorov–Smirnov test [29], the Cramér–von Mises (CVM) test [30], and the test [33]. Momentbased tests included the Jarque–Bera’s (JB) test [20, 21], the robust Jarque–Bera’s (JB) test [20, 21], the kurtosis test [14], the skewness test [14], the D’Agostino–Pearson (DP) test [22], and the simple and exact empirical likelihood ratio (SEELR) test based on moment relations [12]. Lastly, the other categories of competitor tests consisted of the classical wellknown and powerful Shapiro–Wilk (SW) test [13] and the Shapiro–Francia’s (SF) test [15] which is a modification of the SW test. Most of these competitor tests have proved to be powerful against a wide range of alternatives including symmetric ones [33–39].
In terms of the alternative distributions, we considered distributions that cover a wide range of symmetric alternative distributional properties. Following Esteban et al. [40] and Torabi et al. [33], we considered the following alternative distributions, which can be classified into two broad sets of symmetric alternative distributions: (1) symmetric shorttailed distributions and (2) symmetric longtailed distributions. In order to evaluate our proposed tests under asymmetric alternatives, we opted for a third set (3) of asymmetric alternatives:(1)Set 1: symmetric shorttailed distributions:(i)The beta distribution with parameters (3, 3), (2, 2), (1, 1), and (0.5, 0.5)(ii)The uniform distribution, with a = 0 and b = 1(iii)The logitnormal distribution with and (iv)The truncated standard normal distribution at a and b, i.e., (−2, 2) and (−1, 1)(v)Tukey’s lambda distribution with , and 1.25(2)Set 2: symmetric longtailed distributions:(i)Student’s distribution with 2, 4, 7, and 15 degrees of freedom.(ii)The Cauchy distribution with and (iii)The logistic distribution with parameters (location) and (scale)(iv)The double exponential distribution (also known as the Laplace distribution) with parameters (location) and (scale)(v)Tukey’s lambda distribution with , and −0.25(3)Set 3: asymmetric distributions:(i)The gamma distribution with parameters (2, 1)(ii)The Weibull distribution with parameters (2, 1)(iii)The skewed normal distribution with parameters (0, 1, 5)(iv)The skewed Cauchy distribution with parameters (0, 2, 5)(v)The beta distribution with parameters (2, 1) and (3, 1.5)
For power simulation, 10,000 samples each of size , and 100 were obtained under the various alternative distributions. Power was computed by considering the number of times the test rejected the null hypothesis over the total number of replications. A numerical bootstrap study on real data was conducted to assess the robustness and applicability of the proposed tests. However, it was necessary to first assess the type I error control of our proposed tests before the power study.
3.1. Type I Error Control
Here, we provide the values of type I error rates along with the associated standard errors of the proposed tests for , and 0.10. In order to compute these quantities, for each nominal alpha, we generated 500,000 random samples from a standard normal distribution, each corresponding to sample size , and 100. The results presented in Table 3 show that the proposed tests control type I error very well. Figure 2 includes plots of the simulated type error rates only for for all the sample sizes considered. The plots for the empirical cumulative probability function of the simulated values for , , and were omitted since their plots were more or less the same as those for other sample sizes and = 0.05, respectively. It is evident that the plots produced the expected appearance in all the simulated scenarios. That is, the plots show close to the level of simulated type I error rates. The closeness of the estimated probabilities of type I error to the nominal value attests that the GoF test does perform as expected. These results were extended in order to evaluate the type I error control when simulating from a normal distribution with varying parameters of and . We considered various scenarios which include , , , , and for samples of sizes at alpha levels of , and 0.10 (see Table 4). Similarly, as observed in Table 3, the estimated probabilities of type I error were close to the respective nominal values which shows that the GoF test does perform as expected. It is important to note that various alternative methods that can also be used to assess the closeness of the simulated type I error rates to the nominal size alpha are available in the literature. The most popular one is based on the central limit theorem and it was described in detail by Batsidis et al. [41]. Once the type I error rates were examined, we then proceeded to evaluate the powers of the proposed tests to determine how well they would detect departures from normality and to see their power performance compared to those of the selected competing tests.

(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
 
Note. Monte Carlo simulations were conducted using 500,000 replications. 
3.2. Monte Carlo Power Simulation Results
Results for the Monte Carlo power comparisons are presented. Bold numbers in all tables represent the two most superior tests under the respective simulated scenarios. From Tables 5 and 6, we found out that when the alternative distributions are shorttailed and symmetric, our proposed tests performed quite well. Under these symmetric alternative cases, our proposed tests (SRELR and CSELR) significantly outperformed all other studied tests. Tests based on LL, JB, RJB, and have the least power as compared to other tests. In general, the tests based on SRELR, CSELR, DP, SW, and are the most powerful under these symmetric shorttailed alternative distributions.


For symmetric longtailed alternatives (see Tables 7 and 8), the tests based on RJB, , SF, and JB are more superior, whereas the tests based on , LL, and SEELR are the least powerful. Our proposed tests performed slightly lesser than the DP test but were comparable to the SW test. It is important to note that, in all of the cases under these symmetric longtailed alternatives, our proposed tests outperformed all the EDFbased tests.


For the considered asymmetric alternatives (see Table 9), the tests based on SEELR, SW, SF, and AD are more superior, whereas the tests based on , RJB, and SRELR are the least powerful. Our proposed test based on CSELR performed slightly lesser than the JB test but was comparable to the LL test. It is important to note that, in all of the cases under these asymmetric alternatives, our proposed test based on CSELR outperformed the SRELRbased test.

In order to get a clearer visualisation of the performance of the different normality tests, the ranking procedure was used. Tables 10 to 12 contain the ranking of all the tests considered in this study according to the average powers computed from the values in Tables 5–8 and 9, respectively. The rank of power is based on the set of alternative distributions and sample sizes, respectively. Using average powers, we can select the tests that are, on average, most powerful against the alternatives from the given sets of alternatives. It should be noted that, under all the symmetric simulated scenarios, our proposed tests (SRELR and CSELR) were comparable in power.



From Table 10, it can be clearly seen that our proposed tests (SRELR and CSELR) are the most powerful tests for both small and moderate sample sizes under symmetric shorttailed alternatives. This is followed rather closely by the DP test. The results of the total rank based on all sample sizes (i.e., to 100) show that our proposed tests (SRELR and CSELR) are overly the most superior tests for symmetric shorttailed distributions.
For symmetric longtailed alternatives (see Table 11), generally the RJB test was the most powerful in both small and moderate sample sizes. Our proposed tests had comparable power with the AD test under small samples for symmetric longtailed alternatives. However, under moderate sample sizes, our proposed tests were slightly more powerful than the DP and SW tests. Lastly, considering all the sample sizes under symmetric longtailed alternatives, our proposed tests were comparable to the SW test.
Lastly, under asymmetric alternatives (see Table 12), our proposed test based on CSELR performed better than the SRELRbased test. It is also important to note that our related test, the SEELR, outperformed all other tests under these considered asymmetric alternatives. It is also important to note that, unlike some of the competitor tests, our proposed tests were consistent in power under all alternative distributions for all simulated scenarios.
4. Real Data Study
We used the snowfall dataset to examine the applicability of the proposed test on real data. The snowfall dataset consists of 63 snow precipitation values that were recorded from the year 1910 to 1972. The dataset has been extensively used in various statistical applications; see, for example, Thaler [42], Carmichael [43], Tukey [44], and Parzen [45] to illustrate and compare various statistical techniques. The snowfall dataset is presented as follows: 126.4, 82.4, 78.1, 51.1, 90.9, 76.2, 104.5, 87.4, 110.5, 25.0, 69.3, 53.5, 39.8, 63.6, 46.7, 72.9, 79.6, 83.6, 80.7, 60.3, 79.0, 74.4, 49.6, 54.7, 71.8, 49.1, 103.9, 51.6, 82.4, 83.6, 77.8, 79.3, 89.6, 85.5, 58.0, 120.7, 110.5, 65.4, 39.9, 40.1, 88.7, 71.4, 83.0, 55.9, 89.9, 84.8, 105.2, 113.7, 124.7, 114.5, 115.6, 102.4, 101.4, 89.8, 71.5, 70.9, 98.3, 55.5, 66.1, 78.4, 102.5, 97.0, 110.0.
The snowfall data is well known to be consistent with the normal distribution. We plotted a histogram and a QQ plot in order to examine the hypothesis for the normality of the snowfall data (see Figure 3).
(a)
(b)
From the plots, it is clearly visible that the snowfall data are consistent with a normal distribution. Following the ideas introduced by Stigler [46], we conducted a bootstrap type study to empirically examine the proposed test based on the two statistics, CSELR and SRELR. The approach was to use a sample of size 60 by randomly selecting from the snowfall data and then test for normality at 0.05 level of significance. We repeated this strategy 10000 times, and the bootstrap type procedure showed that the proposed CSELR test had a value of 0.7755, while the SRELR had a value of 0.1451. In order to further examine the normality of the snowfall data, we repeated the bootstrap type study using the AD, CVM, JB, and SW tests. The values that were obtained, that is, 0.6862 for the AD test, 0.6921 for the CVM test, 0.5702 for the JB test, and 0.6650 for the SW test were all suggestive for one to conclude that the snowfall data are indeed normally distributed. Thus, the values obtained from the traditional tests as well as our proposed tests show to be reliable in illustrating the normality of the snowfall data. Thus, our proposed test statistics have demonstrated that they are indeed applicable when applied on some reallife data.
5. Conclusion
By utilizing the EL methodology and exploiting the mathematical properties and different forms of transforming the normal distribution, we have developed simple and powerful tests for normality against symmetric alternatives. The proposed tests are consistent and control type I error very well, which is consistent with what has been reported in other studies which looked at ELbased GoF tests (see, for example, [8, 12, 18]). They outperformed other common traditional tests under symmetric shorttailed alternatives. The proposed tests also performed quite well under symmetric longtailed alternatives where they were found to be comparable to the SW test and outperformed all the considered EDF tests. The application of our proposed tests on real data revealed the applicability as well as the robustness of the proposed tests in practice. It would be desirable to develop an ELRbased test for normality that outperforms the classical tests under most alternative distributions that occur in practice. This might be the case after certain modifications and improvements that include further exploring the EL methodology as well as other forms of characterizing the normal distribution. The researchers are currently looking at exploiting the use of EDF in developing an empirical likelihood momentbased EDF test for normality. Thus, combining the characterization of EDFbased tests and EL omnibus tests can potentially improve power under small to moderate sample sizes.
Appendix
#Required packages library(emplik) library(zipfR) #Generate standardized data genedata < −function(n){ #Generate data s < −runif(n, 0, 1) x1 < −s for(k in 1 : n){ s[k] < −(x1[k] − mean(x1)) sqrt(n/(n −1))/(sd(x1[−k])) } return(s) } #Generate transformed data genedata1 < −function(n) { < −genedata(n) x < −pnorm(, 0, 1) return(x)} #Moment function for uniform distribution momentFU < −function(k, a, b){ z < −(b ^ (k + 1) −a ^ (k + 1))/((k + 1) (b − a)) } #Compute test statistic teststatistic < −function(x) { #CUSUMtype statistic k3 = el.test(x 3, m3)$2LLR k5 = el.test(x 5, m5)$2LLR return(max(k3, k5)) #Shiryaev–Roberts statistic #k3 = exp(el.test(x 3, m3)$2LLR #k5 = exp(el.test(x 5, m5)$2LLR #return(sum(k3, k5)) } n < −50 #sample size a < −0 #lower limit for moment function b < −1 #upper limit for moment function #Critical values #Critical values for CUSUM statistic CriticalValue < −0.2685595 #n = 50 ########################################## #Critical values for SR statistic #CriticalValue < −2.459526 #n = 50 MC < −10000 #number of replications power < − c() m3 = momentFU(3, a, b) m5 = momentFU(5, a, b) for(i in 1 : MC) { x < −genedata1(n) power[i] < −teststatistic(x) } # Power for the test under alternative PowerELR = (mean(power > CriticalValue))
Data Availability
The data used to demonstrate the applicability of our proposed tests in practice are presented in this article and can also be obtained from respective authors cited in the “Real Data Study” section. All other data were simulated using R and the source code is available in the Appendix.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would wish to extend their gratitude to Professor Albert Vexler for his insightful comments on ResearchGate. This study was funded through the research incentive funds received by the main author from the Department of Higher Education and Training (DHET) of South Africa.
References
 A. B. Owen, “Empirical likelihood ratio confidence intervals for a single functional,” Biometrika, vol. 75, no. 2, pp. 237–249, 1988. View at: Publisher Site  Google Scholar
 A. Owen, “Empirical likelihood ratio confidence regions,” The Annals of Statistics, vol. 18, no. 1, pp. 90–120, 1990. View at: Publisher Site  Google Scholar
 P. Hall and B. L. Scala, “Methodology and algorithms of empirical likelihood,” International Statistical Review / Revue Internationale de Statistique, vol. 58, no. 2, pp. 109–127, 1990. View at: Publisher Site  Google Scholar
 T. DiCiccio, P. Hall, and J. Romano, “Empirical likelihood is Bartlettcorrectable,” The Annals of Statistics, vol. 19, no. 2, pp. 1053–1061, 1991. View at: Publisher Site  Google Scholar
 J. Qin and J. Lawless, “Empirical likelihood and general estimating equations,” The Annals of Statistics, vol. 22, no. 1, pp. 300–325, 1994. View at: Publisher Site  Google Scholar
 J. Yu, A. Vexler, and L. Tian, “Analyzing incomplete data subject to a threshold using empirical likelihood methods: an application to a pneumonia risk study in an ICU setting,” Biometrics, vol. 66, no. 1, pp. 123–130, 2010. View at: Publisher Site  Google Scholar
 J. H. J. Einmahl and I. W. McKeague, “Empirical likelihood based hypothesis testing,” Bernoulli, vol. 9, no. 2, pp. 267–290, 2003. View at: Publisher Site  Google Scholar
 C. S. Marange and Y. Qin, “A momentbased empirical likelihood ratio test for exponentiality using the probability integral transformation,” Journal of Applied Statistics, vol. 46, no. 15, pp. 2786–2803, 2019. View at: Publisher Site  Google Scholar
 H. Alizadeh Noughabi, “Empirical likelihood ratiobased goodnessoffit test for the logistic distribution,” Journal of Applied Statistics, vol. 42, no. 9, pp. 1973–1983, 2015. View at: Publisher Site  Google Scholar
 A. Vexler and G. Gurevich, “Empirical likelihood ratios applied to goodnessoffit tests based on sample entropy,” Computational Statistics & Data Analysis, vol. 54, no. 2, pp. 531–545, 2010. View at: Publisher Site  Google Scholar
 L. B. Dong and D. E. A. Giles, “An empirical likelihood ratio test for normality,” Communications in Statistics  Simulation and Computation, vol. 36, no. 1, pp. 197–215, 2007. View at: Publisher Site  Google Scholar
 G. Shan, A. Vexler, G. E. Wilding, and A. D. Hutson, “Simple and exact empirical likelihood ratio tests for normality based on moment relations,” Communications in StatisticsSimulation and Computation, vol. 40, no. 1, pp. 129–146, 2010. View at: Publisher Site  Google Scholar
 S. S. Shapiro and M. B. Wilk, “An analysis of variance test for normality (complete samples),” Biometrika, vol. 52, no. 34, pp. 591–611, 1965. View at: Publisher Site  Google Scholar
 S. S. Shapiro, M. B. Wilk, and H. J. Chen, “A comparative study of various tests for normality,” Journal of the American Statistical Association, vol. 63, no. 324, pp. 1343–1372, 1968. View at: Publisher Site  Google Scholar
 S. S. Shapiro and R. S. Francia, “An approximate analysis of variance test for normality,” Journal of the American Statistical Association, vol. 67, no. 337, pp. 215216, 1972. View at: Publisher Site  Google Scholar
 J. P. Royston, “An Extension of Shapiro and Wilk’s W Test for Normality to Large Samples,” Applied Statistics, vol. 31, no. 2, pp. 115–124, 1982. View at: Publisher Site  Google Scholar
 A. B. Owen, Empirical Likelihood, Chapman and Hall, New York, NY, USA, 2001.
 C. S. Marange and Y. Qin, “An empirical likelihood ratio based comparative study on tests for normality of residuals in linear models,” Metodoloski Zvezki, vol. 16, no. 1, pp. 1–16, 2019. View at: Publisher Site  Google Scholar
 C.C. Lin and G. S. Mudholkar, “A simple test for normality against asymmetric alternatives,” Biometrika, vol. 67, no. 2, pp. 455–461, 1980. View at: Publisher Site  Google Scholar
 C. M. Jarque and A. K. Bera, “Efficient tests for normality, homoscedasticity and serial independence of regression residuals,” Economics Letters, vol. 6, no. 3, pp. 255–259, 1980. View at: Publisher Site  Google Scholar
 C. M. Jarque and A. K. Bera, “A test for normality of observations and regression residuals,” International Statistical Review / Revue Internationale de Statistique, vol. 55, no. 2, pp. 163–172, 1987. View at: Publisher Site  Google Scholar
 R. D’Agostino and E. S. Pearson, “Tests for departure from normality. Empirical results for the distributions of and ,” Biometrika, vol. 60, no. 3, pp. 613–622, 1973. View at: Google Scholar
 P. Hall and A. H. Welsh, “A test for normality based on the empirical characteristic function,” Biometrika, vol. 70, no. 2, pp. 485–489, 1983. View at: Publisher Site  Google Scholar
 H. Oja, “New tests for normality,” Biometrika, vol. 70, no. 1, pp. 297–299, 1983. View at: Google Scholar
 Y. Fang, Q. Li, X. Wu, and D. Zhang, “A datadriven smooth test of symmetry,” Journal of Econometrics, vol. 188, no. 2, pp. 490–501, 2015. View at: Publisher Site  Google Scholar
 M. Rosenblatt, “Remarks on a multivariate transformation,” The Annals of Mathematical Statistics, vol. 23, no. 3, pp. 470–472, 1952. View at: Publisher Site  Google Scholar
 T. W. Anderson and D. A. Darling, “Asymptotic theory of certain “goodness of fit” criteria based on stochastic processes,” The Annals of Mathematical Statistics, vol. 23, no. 2, pp. 193–212, 1952. View at: Publisher Site  Google Scholar
 T. W. Anderson and D. A. Darling, “A test of goodness of fit,” Journal of the American Statistical Association, vol. 49, no. 268, pp. 765–769, 1954. View at: Publisher Site  Google Scholar
 H. W. Lilliefors, “On the KolmogorovSmirnov test for normality with mean and variance unknown,” Journal of the American Statistical Association, vol. 62, no. 318, pp. 399–402, 1967. View at: Publisher Site  Google Scholar
 H. Cramér, “On the composition of elementary errors,” Scandinavian Actuarial Journal, vol. 1928, no. 1, pp. 13–74, 1928. View at: Publisher Site  Google Scholar
 A. Vexler, A. Liu, and M. Pollak, “Transformation of changepoint detection methods into a ShiryayevRoberts form,” Department of Biostatistics, The New York State University at Buffalo, Buffalo, NY, USA, 2006, Technical Report. View at: Google Scholar
 A. Vexler and C. Wu, “An optimal retrospective change point detection policy,” Scandinavian Journal of Statistics, vol. 36, no. 3, pp. 542–558, 2009. View at: Publisher Site  Google Scholar
 H. Torabi, N. H. Montazeri, and A. Grané, “A test for normality based on the empirical distribution function,” SORT: Statistics and Operations Research Transactions, vol. 40, no. 1, pp. 0055–0088, 2016. View at: Google Scholar
 H. C. Thode, Testing For Normality, vol. 164, CRC Press, Boca Raton, FL, USA, 2002.
 B. Yazici and S. Yolacan, “A comparison of various tests of normality,” Journal of Statistical Computation and Simulation, vol. 77, no. 2, pp. 175–183, 2007. View at: Publisher Site  Google Scholar
 N. M. Razali and Y. B. Wah, “Power comparisons of shapirowilk, Kolmogorovsmirnov, lilliefors and andersondarling tests,” Journal of Statistical Modeling and Analytics, vol. 2, no. 1, pp. 21–33, 2011. View at: Google Scholar
 B. W. Yap and C. H. Sim, “Comparisons of various types of normality tests,” Journal of Statistical Computation and Simulation, vol. 81, no. 12, pp. 2141–2155, 2011. View at: Publisher Site  Google Scholar
 H. A. Bayoud, “Tests of normality: new test and comparative study,” Communications in Statistics  Simulation and Computation, pp. 1–22, 2019. View at: Publisher Site  Google Scholar
 D. K. Wijekularathna, A. B. W. Manage, and S. M. Scariano, “Power analysis of several normality tests: a Monte Carlo simulation study,” Communications in Statistics  Simulation and Computation, pp. 1–17, 2019. View at: Publisher Site  Google Scholar
 M. Esteban, M. Castellanos, D. Morales, and I. Vajda, “Monte Carlo comparison of four normality tests using different entropy estimates,” Communications in Statistics: Simulation and Computation, vol. 30, no. 4, pp. 761–785, 2001. View at: Publisher Site  Google Scholar
 A. Batsidis, N. Martin, L. Pardo, and K. Zografos, “A necessary power divergencetype family of tests for testing elliptical symmetry,” Journal of Statistical Computation and Simulation, vol. 84, no. 1, pp. 57–83, 2014. View at: Publisher Site  Google Scholar
 H. Thaler, “Nonparametric probability density estimation and the empirical characteristic function,” Statistics Department, State University of New York at Buffalo, Buffalo, NY, USA, 1974, Unpublished Ph.D. thesis. View at: Google Scholar
 J. P. Carmichael, “The autoregressive method,” Statistical Science Division, State University of New York at Buffalo, Buffalo, NY, USA, 1976, Unpublished Ph.D. thesis. View at: Google Scholar
 Turkey, Exploratory Data Analysis, AddisonWesley Publishing Co, Reading, MA, USA, 1977.
 E. Parzen, “Nonparametric Statistical Data Modeling,” Journal of the American Statistical Association, vol. 74, no. 365, pp. 105–121, 1979. View at: Publisher Site  Google Scholar
 S. M. Stigler, “Do robust estimators work with real data?” The Annals of Statistics, vol. 5, no. 6, pp. 1055–1098, 1977. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 Chioneso Show Marange and Yongsong Qin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.