Omnibus tests of multinormality based on skewness and kurtosis

Measures of univariate skewness and kurtosis have long been used as a test of univariate normality, several omnibus test procedures based on a combination of the measures having been proposed, see Pearson, D’Agestino and Bowman (1977) and Mardia (1979). Mardia (1970) proposed measures of multivariate skewness and kurtosis, and constructed a test of multinormality based on these measures. we obtain the correlation between these measures and propose several omnibus tests using the two measures. The performances of these tests are compared by means of a Monte Carlo study.     

Recurrence formula on the joint distribution of sample skewness and kurtosis under normality

Abstract

Two recurrence relations with respect to sample size are given concerning the joint distribution of skewness and kurtosis of random observations from a normal population: one between the probability density functions and the other between the product moments. As a consequence, the latter yields a recurrence formula for the moments of sample kurtosis. The exact moments of Jarque-Bera statistic is also given.     

Measures of multivariate skewness and kurtosis for tests of nonnormality

Measures of multivariate skewness and kurtosis are proposed that are based on the skewness and kurtosis of individual components of standardized sample vectors. Asymptotic properties and small sample critical values of tests for nonnormality based on these measures are provided. It is demonstrated that the tests have favorable power properties. Extensions to time series data are pointed out.     

Intriguing yet simple skewness: kurtosis relation in economic and demographic data distributions,pointing to preferential attachment processes

In this paper, we propose that relations between high-order moments of data distributions, for example, between the skewness (S) and kurtosis (K), allow to point to theoretical models with understandable structural parameters. The illustrative data concern two cases: (i) the distribution of income taxes and (ii) that of inhabitants, after aggregation over each city in each province of Italy in 2011. Moreover, from the rank-size relationship, for either S or K, in both cases, it is shown that one obtains the parameters of the underlying (hypothetical) modeling distribution: in the present cases, the 2-parameter Beta function, itself related to the Yule–Simon distribution function, whence suggesting a growth model based on the preferential attachment process.     

Modelling skewness and kurtosis with the BCPE density in GAMLSS

This paper illustrates the power of modern statistical modelling in understanding processes characterised by data that are skewed and have heavy tails. Our particular substantive problem concerns film box-office revenues. We are able to show that traditional modelling techniques based on the Pareto–Levy–Mandelbrot distribution led to what is actually a poorly supported conclusion that these data have infinite variance. This in turn led to the dominant paradigm of the movie business that ‘nobody knows anything’ and hence that box-office revenues cannot be predicted. Using the Box–Cox power exponential distribution within the generalized additive models for location, scale and shape framework, we are able to model box-office revenues and develop probabilistic statements about revenues.     

Graphical comparison of normality tests for unimodal distribution data

A methodology is proposed to compare the power of normality tests with a wide variety of alternative unimodal distributions. It is based on the representation of a distribution mosaic in which kurtosis varies vertically and skewness horizontally. The mosaic includes distributions such as exponential, Laplace or uniform, with normal occupying the centre. Simulation is used to determine the probability of a sample from each distribution in the mosaic being accepted as normal. We demonstrate our proposal by applying it to the analysis and comparison of some of the most well-known tests.     

Bias-reduced estimates for skewness, kurtosis, L-skewness and L-kurtosis

Estimates based on L-moments are less non-robust than estimates based on ordinary moments because the former are linear combinations of order statistics for all orders, whereas the later take increasing powers of deviations from the mean as the order increases. Estimates based on L-moments can also be more efficient than maximum likelihood estimates. Similarly, L-skewness and L-kurtosis are less non-robust and more informative than the traditional measures of skewness and kurtosis. Here, we give nonparametric bias-reduced estimates of both types of skewness and kurtosis. Their asymptotic computational efficiency is infinitely better than that of corresponding bootstrapped estimates.     

Comparison of non-parametric and semi-parametric tests in detecting long memory

The first two stages in modelling times series are hypothesis testing and estimation. For long memory time series, the second stage was studied in the paper published in [M. Boutahar et al., Estimation methods of the long memory parameter: monte Carlo analysis and application, J. Appl. Statist. 34(3), pp. 261–301.] in which we have presented some estimation methods of the long memory parameter. The present paper is intended for the first stage, and hence completes the former, by exploring some tests for detecting long memory in time series. We consider two kinds of tests: the non-parametric class and the semi-parametric one. We precise the limiting distribution of the non-parametric tests under the null of short memory and we show that they are consistent against the alternative of long memory. We perform also some Monte Carlo simulations to analyse the size distortion and the power of all proposed tests. We conclude that for large sample size, the two classes are equivalent but for small sample size the non-parametric class is better than the semi-parametric one.     

X2 and its components as tests of normality for grouped data

We consider testing for an unobservable normal distribution with unspecified mean and variance. It is only possible to observe the counts in groups with boundaries specified before sighting the data. On the basis of a small power study, we recommend the usual X² test be used as an omnibus test, augmented by informal examination of the first two non-zero components of X². We also recommend use of maximum likelihood and method of moments estimation.     

Misspecification tests for periodic long memory GARCH models

Distributional theory for Quasi-Maximum Likelihood estimators in long memory conditional heteroskedastic models is not formally defined, even asymptotically. Because of that, this paper analyses the real size and power of the likelihood ratio and the Lagrange multiplier misspecification tests when periodic long memory GARCH models are involved. The performance of these tests is studied by means of Monte Carlo simulations with respect to the class of generalized long memory GARCH models. For this class of models, analytical derivatives are developed. An application to the USD/JPY exchange rate is also provided.     

Stringency-based ranking of normality tests

Many parametric statistical inferential procedures in finite samples depend crucially on the underlying normal distribution assumption. Dozens of normality tests are available in the literature to test the hypothesis of normality. Availability of such a large number of normality tests has generated a large number of simulation studies to find a best test but no one arrived at a definite answer as all depends critically on the alternative distributions which cannot be specified. A new framework, based on stringency concept, is devised to evaluate the performance of the existing normality tests. Mixture of t-distributions is used to generate the alternative space. The LR-tests, based on Neyman–Pearson Lemma, have been computed to construct a power envelope for calculating the stringencies of the selected normality tests. While evaluating the stringencies, Anderson–Darling (AD) statistic turns out to be the best normality test.     

Correctings tests for trend in proportions for skewness

We develop the score test for the hypothesis that a parameter of a Markov sequence is constant over time, against the alternatives that it varies over time, i.e., θ_t = θ + U_t; t = 1,2,…, where {U_t; t = 1,2,...} is a sequence of independently and identically distributed random variables with mean zero and variance σ^z _u and θ is a fixed constant. The asymptotic null distribution of the test statistic is proved to be normal. We illustrate our procedure by examples and a real life data analysis.     

A re-appraisal of tests for normality

Tests for normality can be divided into two groups - those based upon a function of the empirical distribution function and those based upon a function of the original observations. The latter group of statistics test spherical symmetry and not necessarily normality. If the distribution is completely specified then the first group can be used to test for ‘spherical’ normality. However, if the distribution is incompletely specified and F‘‘x_i - x’/s’ is used these test statistics also test sphericity rather than normality. A Monte Carlo study was conducted for the completely specified case, to investigate the sensitivity of the distance tests to departures from normality when the alternative distributions are non-normal spherically symmetric laws. A “new” test statistic is proposed for testing a completely specified normal distribution     

Comparisons of various types of normality tests

Normality tests can be classified into tests based on chi-squared, moments, empirical distribution, spacings, regression and correlation and other special tests. This paper studies and compares the power of eight selected normality tests: the Shapiro–Wilk test, the Kolmogorov–Smirnov test, the Lilliefors test, the Cramer–von Mises test, the Anderson–Darling test, the D'Agostino–Pearson test, the Jarque–Bera test and chi-squared test. Power comparisons of these eight tests were obtained via the Monte Carlo simulation of sample data generated from alternative distributions that follow symmetric short-tailed, symmetric long-tailed and asymmetric distributions. Our simulation results show that for symmetric short-tailed distributions, D'Agostino and Shapiro–Wilk tests have better power. For symmetric long-tailed distributions, the power of Jarque–Bera and D'Agostino tests is quite comparable with the Shapiro–Wilk test. As for asymmetric distributions, the Shapiro–Wilk test is the most powerful test followed by the Anderson–Darling test.     

Empirical behavior of some tests for normality

The behavior of a range of tests assessing the normality of a sequence of independent and identically distributed random variables is investigated.An examination of the empirical significance level of the tests is undertaken for different sample sizes. The empirical power associated with these tests is also calculated under some alternative distributions.     

Efficiency of ranked set sampling in tests for normality

In this article, we consider the ranked set sampling (RSS) and investigate seven tests for normality under RSS. Each test is described and then power of each test is obtained by Monte Carlo simulations under various alternatives. Finally, the powers of the tests based on RSS are compared with the powers of the tests based on the simple random sampling and the results are discussed.     

Testing high-dimensional normality based on classical skewness and Kurtosis with a possible small sample size

Abstract

By using the idea of principal component analysis, we propose an approach to applying the classical skewness and kurtosis statistics for detecting univariate normality to testing high-dimensional normality. High-dimensional sample data are projected to the principal component directions on which the classical skewness and kurtosis statistics can be constructed. The theory of spherical distributions is employed to derive the null distributions of the combined statistics constructed from the principal component directions. A Monte Carlo study is carried out to demonstrate the performance of the statistics on controlling type I error rates and a simple power comparison with some existing statistics. The effectiveness of the proposed statistics is illustrated by two real-data examples.     

On tests for normality of experimental error in ridge regression

Considered are tests for normality of the errors in ridge regression. If an intercept is included in the model, it is shown that test statistics based on the empirical distribution function of the ridge residuals have the same limiting distribution as in the one-sample test for normality with estimated mean and variance. The result holds with weak assumptions on the behavior of the independent variables; asymptotic normality of the ridge estimator is not required.     

On tests detecting difference in means and variances simultaneously under normality

In this paper, we propose several tests for detecting difference in means and variances simultaneously between two populations under normality. First of all, we propose a likelihood ratio test. Then we obtain an expression of the likelihood ratio statistic by a product of two functions of random quantities, which can be used to test the two individual partial hypotheses for differences in means and variances. With those individual partial tests, we propose a union-intersection test. Also we consider two optimal tests by combining the p-values of the two individual partial tests. For obtaining null distributions, we apply the permutation principle with the Monte Carlo approach. Then we compare efficiency among the proposed tests with well-known ones through a simulation study. Finally, we discuss some interesting features related to the simultaneous tests and resampling methods as concluding remarks.