Hybrid test for the hypothesis of symmetry

In recent years, McWilliams and Tajuddin have proposed new and more powerful non-parametric tests of symmetry for continuous distributions about a known center. In this paper, we propose a simple non-parametric two-stage procedure based on the sign test and a percentile-modified two-sample Wilcoxon test. The small-sample properties of this test, Tajuddin's test, McWilliams' test and a modified runs test of Modarres and Gastwirth are investigated in a Monte Carlo simulation study. The simulations indicate that, for a wide variety of asymmetric alternatives in the lambda family, the hybrid test is more powerful than are existing tests in the literature.     

Testing for Circular Reflective Symmetry about a Known Median Axis

Circular data arise in many contexts, a particularly rich source being animal orientation experiments. Often, in the analysis of such data, a fundamental question of scientific interest is whether the underlying distribution is reflectively symmetric about some specific axis. In this paper, the situation in which the axis of interest is known to be a median axis is considered and a simple, asymptotically distribution- free test for circular reflective symmetry against skew alternatives is developed. The results from a simulation study lead to a testing strategy incorporating the new test and the circular analogue of the modified runs test of Modarres & Gastwirth (1996). The application of the testing strategy is illustrated using circular data arising from two animal orientation experiments.     

A Comparison of the Runs Test for Volatility Forecastability and the LM Test for GARCH Using Aggregated Returns

Christoffersen and Diebold (2000) have introduced a runs test for forecastable volatility in aggregated returns. In this note, we compare the size and power of their runs test and the more conventional LM test for GARCH by Monte Carlo simulation. When the true daily process is GARCH, EGARCH, or stochastic volatility, the LM test has better power than the runs test for the moderate-horizon returns considered by Christoffersen and Diebold. For long-horizon returns, however, the tests have very similar power. We also consider a qualitative threshold GARCH model. For this process, we find that the runs test has greater power than the LM test. Theresults support the use of the runs test with aggregated returns.     

Bayesian analysis of screening data: Application to AIDS in blood donors

The rapid increase in the number of AIDS cases during the 1980s and the spread of the disease from the high-risk groups into the general population has created widespread concern. In particular, assessing the accuracy of the screening tests used to detect antibodies to the HIV (AIDS) virus in donated blood and determining the prevalance of the disease in the population are fundamental statistical problems. Because the prevalence of AIDS varies widely by geographic region and data on the number of infected blood donors are published regularly, Bayesian methods, which utilize prior results and update them as new data become available, are quite useful. In this paper we develop a Bayesian procedure for estimating the prevalence of a rare disease, the sensitivity and specificity of the screening tests, and the predictive value of a positive or negative screening test. We apply the procedure to data on blood donors in the United States and in Canada. Our results augment those described in Gastwirth (1987) using classical methods. Indeed, we show that the inclusion of sound prior knowledge into the statistical analysis does not yield sufficiently precise estimates of the predictive value of a positive test. Hence confirmatory testing is needed to obtain reliable estimates. The emphasis of the Bayesian predictive paradigm on prediction intervals for future data yields a valuable insight. We demonstrate that using them might have detected a decline in the specificity of the most frequently used screening test earlier than it apparently was.     

The difference-sign runs length distribution in testing for serial independence

We investigate the sequence of difference-sign runs length of a time series in the context of non-parametric tests for serial independence. This sequence is, under suitable conditioning, a stationary sequence and we prove that the normalized correlation of two consecutive runs length is small (≈0.0427). We use this result in a test based on the relative entropy of the empirical distribution of the runs length. We investigate the performance of the test in simulated series and test serial independence of cardiac data series in atrial fibrillation.     

An approximate entropy test for randomness

We develop an entropy-based test for randomness of binary time series of finite length. The test uses the frequencies of contiguous blocks of different lengths. A simple condition ib the block lengths and the length of the time series enables one to estimate the entropy rate for the data, and this information is used to develop a statistic to test the hypothesis of randomness. This static measures the deviation of the estimated entropy of the observed data from the theoretical maximum under the randomness hypothesis. This test offers a real alternative to the conventional runs test. Critical percentage points, based on simulations, are provided for testing the hypothesis of randomness. Power calculations using dependent data show that the proposed test has higher power against the runs test for short series, and it is similar to the runs test for long series. The test is applied to two published data sets that wree investigated by others with respect to their randomness.     

Adaptive positive and negative runs test

As a nonparametric randomness test, the positive and negative runs test is widely used in practice due to the simplicity of its procedures. The test can lose efficiency if the alternative distribution is symmetrical at 0.5. In addition, the test can only be applied to test the randomness of a sequence from the uniform distribution. In this paper, we introduce an adaptive positive and negative runs test method to maximize the power function by choosing the optimal cut point. Also, the test is extended to check the randomness of a sequence generated from any other given distributions. Furthermore, we derive the exact distribution and obtain the asymptotical critical values of the proposed test statistics. Compared with the existed test, the efficiency of the proposed adaptive positive and negative runs test is competitive through simulation study.     

Data depth‐based nonparametric scale tests

Liu and Singh (1993, 2006) introduced a depth‐based d‐variate extension of the nonparametric two sample scale test of Siegel and Tukey (1960). Liu and Singh (2006) generalized this depth‐based test for scale homogeneity of k ≥ 2 multivariate populations. Motivated by the work of Gastwirth (1965), we propose k sample percentile modifications of Liu and Singh's proposals. The test statistic is shown to be asymptotically normal when k = 2, and compares favorably with Liu and Singh (2006) if the underlying distributions are either symmetric with light tails or asymmetric. In the case of skewed distributions considered in this paper the power of the proposed tests can attain twice the power of the Liu‐Singh test for d ≥ 1. Finally, in the k‐sample case, it is shown that the asymptotic distribution of the proposed percentile modified Kruskal‐Wallis type test is χ² with k ? 1 degrees of freedom. Power properties of this k‐sample test are similar to those for the proposed two sample one. The Canadian Journal of Statistics 39: 356–369; 2011 © 2011 Statistical Society of Canada     

Modified Wilcoxon–Mann–Whitney Test and Power Against Strong Null

The Wilcoxon–Mann–Whitney (WMW) test is a popular rank-based two-sample testing procedure for the strong null hypothesis that the two samples come from the same distribution. A modified WMW test, the Fligner–Policello (FP) test, has been proposed for comparing the medians of two populations. A fact that may be under-appreciated among some practitioners is that the FP test can also be used to test the strong null like the WMW. In this article, we compare the power of the WMW and FP tests for testing the strong null. Our results show that neither test is uniformly better than the other and that there can be substantial differences in power between the two choices. We propose a new, modified WMW test that combines the WMW and FP tests. Monte Carlo studies show that the combined test has good power compared to either the WMW and FP test. We provide a fast implementation of the proposed test in an open-source software. Supplementary materials for this article are available online.     

Nonparametric Simultaneous Tests for Location and Scale Testing: A Comparison of Several Methods

The two-sample location-scale problem arises in many situations like climate dynamics, bioinformatics, medicine, and finance. To address this problem, the nonparametric approach is considered because in practice, the normal assumption is often not fulfilled or the observations are too few to rely on the central limit theorem, and moreover outliers, heavy tails and skewness may be possible. In these situations, a nonparametric test is generally more robust and powerful than a parametric test. Various nonparametric tests have been proposed for the two-sample location-scale problem. In particular, we consider tests due to Lepage, Cucconi, Podgor-Gastwirth, Neuhäuser, Zhang, and Murakami. So far all these tests have not been compared. Moreover, for the Neuhäuser test and the Murakami test, the power has not been studied in detail. It is the aim of the article to review and compare these tests for the jointly detection of location and scale changes by means of a very detailed simulation study. It is shown that both the Podgor–Gastwirth test and the computationally simpler Cucconi test are preferable. Two actual examples within the medical context are discussed.     

Linear Signed Rank Test for Model Selection

In this article, we consider a linear signed rank test for non-nested distributions in the context of the model selection. Introducing a new test, we show that, it is asymptotically more efficient than the Vuong test and the test statistic based on B statistic introduced by Clarke. However, here, we let the magnitude of the data give a better performance to the test statistic. We have shown that this test is an unbiased one. The results of simulations show that the rank test has the greater statistical power than the Vuong test where the underline distributions is symmetric.     

Unbiased estimation of income inequality

In a previous paper Gastwirth shows that a broad family of measures of inequality can be accurately estimated when the tax data are known in groups (more precisely, when we know the number of returns in each of several class intervals and their corresponding average income). In the present paper we show that some measures of the preceding family can be unbiasedly estimated when the tax data are individually known for a sample from the population. Specifically, we construct unbiased estimators of a particular measure of inequality in the samplings with and without replacement, and in the stratified samplings with and without replacement.     

A parametric bootstrap solution to the MANOVA under heteroscedasticity

In this article, we consider the problem of comparing several multivariate normal mean vectors when the covariance matrices are unknown and arbitrary positive definite matrices. We propose a parametric bootstrap (PB) approach and develop an approximation to the distribution of the PB pivotal quantity for comparing two mean vectors. This approximate test is shown to be the same as the invariant test given in [Krishnamoorthy and Yu, Modified Nel and Van der Merwe test for the multivariate Behrens–Fisher problem, Stat. Probab. Lett. 66 (2004), pp. 161–169] for the multivariate Behrens–Fisher problem. Furthermore, we compare the PB test with two existing invariant tests via Monte Carlo simulation. Our simulation studies show that the PB test controls Type I error rates very satisfactorily, whereas other tests are liberal especially when the number of means to be compared is moderate and/or sample sizes are small. The tests are illustrated using an example.     

Kolmogorov–Smirnov,Fluctuation, and Z g Tests for Convergence of Markov Chain Monte Carlo Draws

We examine the sizes and powers of three tests of convergence of Markov Chain Monte Carlo draws: the Kolmogorov–Smirnov test, fluctuation test, and Geweke's test. We show that the sizes and powers are sensitive to the existence of autocorrelation in the draws. We propose a filtered test that is corrected for autocorrelation. We present a numerical illustration using the Federal funds rate.     

On blocks and runs estimators of the extremal index

Given a sample from a stationary sequence of random variables, we study the blocks and runs estimators of the extremal index. Conditions are given for consistency and asymptotic normality of these estimators. We show that moment restrictions assumed by Hsing (Stochast. Process. Appl. 37(1), 117–139; Ann. Statist. 21(4), 2043-2021) may be relaxed if a stronger mixing condition holds. The CLT for the runs estimator seems to be proven for the first time.     

A powerful and interpretable alternative to the Jarque–Bera test of normality based on 2nd-power skewness and kurtosis,using the Rao's score test on the APD family

We introduce the 2nd-power skewness and kurtosis, which are interesting alternatives to the classical Pearson's skewness and kurtosis, called 3rd-power skewness and 4th-power kurtosis in our terminology. We use the sample 2nd-power skewness and kurtosis to build a powerful test of normality. This test can also be derived as Rao's score test on the asymmetric power distribution, which combines the large range of exponential tail behavior provided by the exponential power distribution family with various levels of asymmetry. We find that our test statistic is asymptotically chi-squared distributed. We also propose a modified test statistic, for which we show numerically that the distribution can be approximated for finite sample sizes with very high precision by a chi-square. Similarly, we propose a directional test based on sample 2nd-power kurtosis only, for the situations where the true distribution is known to be symmetric. Our tests are very similar in spirit to the famous Jarque–Bera test, and as such are also locally optimal. They offer the same nice interpretation, with in addition the gold standard power of the regression and correlation tests. An extensive empirical power analysis is performed, which shows that our tests are among the most powerful normality tests. Our test is implemented in an R package called PoweR.     

Generalized run tests for statistical process control

ABSTRACT

In a sequence of elements, a run is defined as a maximal subsequence of like elements. The number of runs or the length of the longest run has been widely used to test the randomness of an ordered sequence. Based on two different sampling methods and two types of test statistics used, run tests can be classified into one of four cases. Numerous researchers have derived the probability distributions in many different ways, treating each case separately. In the paper, we propose a unified approach which is based on recurrence arguments of two mutually exclusive sub-sequences. We also consider the sequence of nominal data that has more than two classes. Thus, the traditional run tests for a binary sequence are special cases of our generalized run tests. We finally show that the generalized run tests can be applied to many quality management areas, such as testing changes in process variation, developing non-parametric multivariate control charts, and comparing the shapes and locations of more than two process distributions.     

A Comparison of the Runs Test for Volatility Forecastability and the LM Test for GARCH Using Aggregated Returns

Christoffersen and Diebold (2000 Christoffersen , P. F. , Diebold , F. X. ( 2000 ). How relevant is volatility forecasting for financial risk management? Review of Economics and Statistics 82 : 12 – 22 .[Crossref] , [Google Scholar]) have introduced a runs test for forecastable volatility in aggregated returns. In this note, we compare the size and power of their runs test and the more conventional LM test for GARCH by Monte Carlo simulation. When the true daily process is GARCH, EGARCH, or stochastic volatility, the LM test has better power than the runs test for the moderate-horizon returns considered by Christoffersen and Diebold. For long-horizon returns, however, the tests have very similar power. We also consider a qualitative threshold GARCH model. For this process, we find that the runs test has greater power than the LM test. Theresults support the use of the runs test with aggregated returns.     

Robust suptest for the genetic association study under genetic model uncertainty

In case–control studies the Cochran–Armitage trend test is powerful for detection of an association between a risk genetic marker and a disease of interest. To apply this test, a score should be assigned to the genotypes based on the genetic model. When the underlying genetic model is unknown, the trend test statistic is quite sensitive to the choice of the score. In this paper, we study the asymptotic property of the robust suptest statistic defined as a supremum of Cochran–Armitage trend test across all scores between 0 and 1. Through numerical studies we show that small to moderate sample size performances of the suptest appear reasonable in terms of type I error control and we compared empirical powers of the suptest to those of three individual Cochran–Armitage trend tests and the maximum of the three Cochran–Armitage trend tests. The use of the suptest is applied to rheumatoid arthritis data from a genome-wide association study.     

Testing for parameter constancy in the time series direction in panel data models

We propose tests for parameter constancy in the time series direction in panel data models. We construct a locally best invariant test based on Tanaka [Time series analysis: nonstationary and noninvertible distribution theory. New York: Wiley; 1996] and an asymptotically point optimal test based on Elliott and Müller [Efficient tests for general persistent time variation in regression coefficients. Rev Econ Stud. 2006;73:907–940]. We derive the limiting distributions of the test statistics as T→∞ while N is fixed, and calculate the critical values by applying numerical integration and response surface regression. Simulation results show that the proposed tests perform well if we apply them appropriately.