An extension of Monte Carlo hypothesis tests

There are many hypothesis testing settings in which one can calculate a “reasonable” test statistic, but in which the null distribution of the statistic is unknown or completely intractable. Fortunately, in many such situations, it is possible to simulate values of the test statistic under the null hypothesis, in which case one can conduct a Monte Carlo test. A difficulty however arises in that Monte Carlo tests, as they are currently structured, are applicable only if ties cannot occur among the values of the test statistics. There is a frequently occurring scenario in which there are lots of ties, namely that in which the null distribution of the test statistic has a (single) point mass. It turns out that one can modify the current form of Monte Carlo tests so as to accommodate such settings. Developing this modification leads to an intriguing identity involving the binomial probability function and its derivatives. In this article, we will briefly explain the modified procedure, discuss simulation studies which demonstrate its efficacy, and provide a proof of the identity referred to above.     

Goodness-of-fit test for exponentiality based on spacings for general progressive Type-II censored data

There have been numerous tests proposed to determine whether or not the exponential model is suitable for a given data set. In this article, we propose a new test statistic based on spacings to test whether the general progressive Type-II censored samples are from exponential distribution. The null distribution of the test statistic is discussed and it could be approximated by the standard normal distribution. Meanwhile, we propose an approximate method for calculating the expectation and variance of samples under null hypothesis and corresponding power function is also given. Then, a simulation study is conducted. We calculate the approximation of the power based on normality and compare the results with those obtained by Monte Carlo simulation under different alternatives with distinct types of hazard function. Results of simulation study disclose that the power properties of this statistic by using Monte Carlo simulation are better for the alternatives with monotone increasing hazard function, and otherwise, normal approximation simulation results are relatively better. Finally, two illustrative examples are presented.     

Is There a Jump in the Transition?

This article develops a statistical test for the presence of a jump in an otherwise smooth transition process. In this testing, the null model is a threshold regression and the alternative model is a smooth transition model. We propose a quasi-Gaussian likelihood ratio statistic and provide its asymptotic distribution, which is defined as the maximum of a two parameter Gaussian process with a nonzero bias term. Asymptotic critical values can be tabulated and depend on the transition function employed. A simulation method to compute empirical critical values is also developed. Finite-sample performance of the test is assessed via Monte Carlo simulations. The test is applied to investigate the dynamics of racial segregation within cities across the United States.     

Book Review: Introducing Monte Carlo Methods with R

In this article we propose a nonparametric test for poolability in large dimensional semiparametric panel data models with cross-section dependence based on the sieve estimation technique. To construct the test statistic, we only need to estimate the model under the alternative. We establish the asymptotic normal distributions of our test statistic under the null hypothesis of poolability and a sequence of local alternatives, and prove the consistency of our test. We also suggest a bootstrap method as an alternative way to obtain the critical values. A small set of Monte Carlo simulations indicate the test performs reasonably well in finite samples.     

On a data-dependent choice of the tuning parameter appearing in certain goodness-of-fit tests

We propose a data-dependent method for choosing the tuning parameter appearing in many recently developed goodness-of-fit test statistics. The new method, based on the bootstrap, is applicable to a class of distributions for which the null distribution of the test statistic is independent of unknown parameters. No data-dependent choice for this parameter exists in the literature; typically, a fixed value for the parameter is chosen which can perform well for some alternatives, but poorly for others. The performance of the new method is investigated by means of a Monte Carlo study, employing three tests for exponentiality. It is found that the Monte Carlo power of these tests, using the data-dependent choice, compares favourably to the maximum achievable power for the tests calculated over a grid of values of the tuning parameter.     

No-cointegration test based on fractional differencing: Some Monte Carlo results

This paper examines the use of the t-statistic in the Geweke–Porter-Hudak regression for the estimation of the fractional differencing parameter as a test for cointegration. The critical values of the test statistic are estimated using Monte Carlo methods. The results confirm that the test will over-reject the null hypothesis of no cointegration if the standard-normal critical values are used. The estimated critical values are generally robust to the nuisance parameters in the autoregressive or moving average specification of the error process of the component time series. Exceptions occur when the dependent variable in the cointegration regression follows an autoregressive process with a large positive parameter or a moving average process with a large negative parameter.     

Robust estimators and robust tests for the slightly contaminated stochastic logistic population models

This article introduces the robust indirect technique for the slightly contaminated stochastic logistic population models. Based on discrete sampled data with a fixed unit of time between two consecutive observations, we not only construct the robust indirect inference generalized method of moments (GMM) estimator for the model parameters, but also propose a likelihood-ratio-type indirect statistic and a robust indirect GMM saddle-point statistic for testing the parameters of interest. In addition, we develop the robust exponential tilting estimator and the robust exponential tilting test to improve their small sample performances. Finally, their finite-sample properties are studied through Monte Carlo experiments.     

GMM tests for the Katz family of distributions

Generalized method of moments (GMM) is used to develop tests for discriminating discrete distributions among the two-parameter family of Katz distributions. Relationships involving moments are exploited to obtain identifying and over-identifying restrictions. The asymptotic relative efficiencies of tests based on GMM are analyzed using the local power approach and the approximate Bahadur efficiency. The paper also gives results of Monte Carlo experiments designed to check the validity of the theoretical findings and to shed light on the small sample properties of the proposed tests. Extensions of the results to compound Poisson alternative hypotheses are discussed.     

A Goodness-of-fit Test for Exponentiality Based on the Memoryless Property

In this investigation a test of goodness of fit for exponentiality is proposed. This procedure applies equally whether the scale and/or the location parameters of the distribution are known or not. The limiting null and non-null distributions of the test statistic are normal under minimal conditions. Monte Carlo critical values for small sample sizes are given and the power of the test is calculated for various alternatives showing that it compares favourably relatively to other more complicated published procedures.     

A simple test of symmetry about an unknown median

A simple test statistic for testing symmetry of a distribution function about an unknown value is presented. The asymptotic distributions under symmetry and asymmetry are derived. Using the normal as a “calibration” distribution, the critical values of the test are calculated by Monte Carlo methods. Comparisons with other tests indicate that this procedure performs well.     

Repeated chi-square testing

A sequentialized version of the x²; goodness of fit test, called repeated x,²; test, is introduced. The form of the asymptotic distribution of the repeated x² test statistic is given under the null hypothesis as well as under local alternatives. For various numbers of cells Monte Carlo results are given for critical values, power and distribution of stopping time. Finally, the perfor-mance of the repeated and the fixed sample x² test are compared.     

A Generalized Distance Multi-Sample Test of Normality With Applications to Process Control

In statistical process control one typically takes periodic small samples. Statistical inferences made from these samples often assume that the samples come from normal distributions with the means and variances possibly changing over time. A multisample test of normality is proposed to test this assumption. The test statistic is the generalized distance between the standardized order statistic vector averaged across the samples and its expected value under normality. The null distribution of the statistic approaches a chi-squared distribution as the number of samples increases. A Monte Carlo study suggests that the test has desirable power properties relative to competing tests.     

Using the Bootstrap to Test for Symmetry Under Unknown Dependence

This article considers tests for symmetry of the one-dimensional marginal distribution of fractionally integrated processes. The tests are implemented by using an autoregressive sieve bootstrap approximation to the null sampling distribution of the relevant test statistics. The sieve bootstrap allows inference on symmetry to be carried out without knowledge of either the memory parameter of the data or of the appropriate norming factor for the test statistic and its asymptotic distribution. The small-sample properties of the proposed method are examined by means of Monte Carlo experiments, and applications to real-world data are also presented.     

Modification of the kolmogorov-smirnov test for the erlang-2 distribution

An adjusted Kolmogorov-Smirnov statistic and critical values are developed for the Erlang-2 probability distribution using data from Monte Carlo simulations. The process used is similar to that of Stephens in the 1970s. The test statistic produced features of compactness and ease of implementation. It is quite accurate for sample sizes as low as ten.     

Asymptotically refined score and GOF tests for inverse Gaussian models

ABSTRACT

The score test and the GOF test for the inverse Gaussian distribution, in particular the latter, are known to have large size distortion and hence unreliable power when referring to the asymptotic critical values. We show in this paper that with the appropriately bootstrapped critical values, these tests become second-order accurate, with size distortion being essentially eliminated and power more reliable. Two major generalizations of the score test are made: one is to allow the data to be right-censored, and the other is to allow the existence of covariate effects. A data mapping method is introduced for the bootstrap to be able to produce censored data that are conformable with the null model. Monte Carlo results clearly favour the proposed bootstrap tests. Real data illustrations are given.     

Asymptotocally nonparametric tests with censored paired data

A class of asymptotically nonparametric test with contains a test proposed by Wei(1980), is considered for testing the equality of two continuous distribution funcitons when paired observations are subject to arbitrary right censorship. It is shown that under the null hypothesis each test statistic converges in distribution to the standard normal random variable. Furthermore. the Monte Carlo simulation results indicate that some tests in this class are more powerful than Wei's test. A generalization to incomplete censored paired data is also included.     

Unit root test for short panels with serially correlated errors

This paper proposes a unit root test for short panels with serially correlated errors. The proposed test is based on the instrumental variables (IVs) and the generalized method of moments (GMM) estimators. An advantage of the new test over other tests is that it allows for an ARMA-type serial correlation. A Monte Carlo simulation shows that the new test has good finite sample properties. Several methods to estimate the lag orders of the ARMA structure are briefly discussed.     

非线性协整的秩检验方法及其响应面函数研究

 本文对非线性协整关系的秩检验方法进行了系统的梳理,运用Monte Carlo模拟给出了不同样本容量的各个秩检验统计量的临界值,并进一步探讨了其响应面函数,给出了各个秩检验统计量临界值的近似计算公式。对中国上证综指与主要发达国家股指关系的秩协整检验表明,与传统线性协整Johansen检验相比,秩协整检验能够检测到更多的线性和非线性协整关系。     

On Variance Components in Semiparametric Mixed Models for Longitudinal Data

Abstract. First, to test the existence of random effects in semiparametric mixed models (SMMs) under only moment conditions on random effects and errors, we propose a very simple and easily implemented non‐parametric test based on a difference between two estimators of the error variance. One test is consistent only under the null and the other can be so under both the null and alternatives. Instead of erroneously solving the non‐standard two‐sided testing problem, as in most papers in the literature, we solve it correctly and prove that the asymptotic distribution of our test statistic is standard normal. This avoids Monte Carlo approximations to obtain p ‐values, as is needed for many existing methods, and the test can detect local alternatives approaching the null at rates up to root n. Second, as the higher moments of the error are necessarily estimated because the standardizing constant involves these quantities, we propose a general method to conveniently estimate any moments of the error. Finally, a simulation study and a real data analysis are conducted to investigate the properties of our procedures.     

An entropy-based test for goodness of fit of the von mises distribution

The maximum entropy characterization of the von Mises distribution on the circle is exploited to derive a consistent goodness of fit test for the von Mises distribution. Monte Carlo simulation results are tabulated giving critical values of the test statistic for various sample sizes and values of the concentration parameter. A power analysis is presented for various alternative hypotheses, comparing this entropy statistic to two other competing goodness of fit statistics. The entropy statistic is shown to compare favorably and may be more attractive, especially considering its ease of computation.