Testing a markov chain for independence

A bootstrap procedure is proposed for testing whether an observed Markov chain is actually an independent process, based on the observed transition probability matrix. The results of simulations showing the power and size of the bootstrap test are presented. The asymptotic distribution of the non-unit eigenvalues is given under the null hypothesis.     

Recent and classical goodness-of-fit tests for the Poisson distribution

We give a critical synopsis of classical and recent tests for Poissonity, our emphasis being on procedures which are consistent against general alternatives. Two classes of weighted Cramér–von Mises type test statistics, based on the empirical probability generating function process, are studied in more detail. Both of them generalize already known test statistics by introducing a weighting parameter, thus providing more flexibility with regard to power against specific alternatives. In both cases, we prove convergence in distribution of the statistics under the null hypothesis in the setting of a triangular array of rowwise independent and identically distributed random variables as well as consistency of the corresponding test against general alternatives. Therefore, a sound theoretical basis is provided for the parametric bootstrap procedure, which is applied to obtain critical values in a large-scale simulation study. Each of the tests considered in this study, when implemented via the parametric bootstrap method, maintains a nominal level of significance very closely, even for small sample sizes. The procedures are applied to four well-known data sets.     

Bootstrap methods for the nonparametric assessment of population bioequivalence and similarity of distributions

A completely nonparametric approach to population bioequivalence in crossover trials has been suggested by Munk and Czado (1999). It is based on the Mallows (1972) metric as a nonparametric distance measure which allows the comparison between the entire distribution functions of test and reference formulations. It was shown that a separation between carry-over and period effects is not possible in the nonparametric setting. However when carry-over effects can be excluded, treatment effects can be assessed when period effects are or not. Munk and Czado (1999) proved bootstrap limit laws of the corresponding test statistics because estimation of the limiting variance of the test statistic is very cumbersome. The purpose of this paper is to investigate the small sample behavior of various bootstrap methods and to compare it with the asymptotic test obtained by estimation of the limiting variance. The percentile (PC) and bias correct- ed and accelerated (BCA) bootstrap were compared for multivariate normal and nonnormal populations. From the simulation results presented, the BCA bootstrap is found to be less conservative and provides higher power compared to the PC bootstrap, especially when skewed multivariate populations are present.     

协整参数的自举推断

 针对完全修正最小二乘(full-modified ordinary least square,简称FMOLS）估计方法,给出一种协整参数的自举推断程序,证明零假设下自举统计量与检验统计量具有相同的渐近分布。关于检验功效的研究表明,虽然有约束自举的实际检验水平表现良好,但如果零假设不成立,自举统计量的分布是不确定的,因而其经验分布不能作为检验统计量精确分布的有效估计。实际应用中建议使用无约束自举,因为无论观测数据是否满足零假设,其自举统计量与零假设下检验统计量都具有相同的渐近分布。最后,利用蒙特卡洛模拟对自举推断和渐近推断的有限样本表现进行比较研究。     

A simple test procedure in standardizing the power of Hosmer–Lemeshow test in large data sets

The Hosmer–Lemeshow (H–L) test is a widely used method when assessing the goodness-of-fit of a logistic regression model. However, the H–L test is sensitive to the sample sizes and the number of groups in H–L test. Cautions need to be taken for interpreting an H–L test with a large sample size. In this paper, we propose a simple test procedure to evaluate the model fit of logistic regression model with a large sample size, in which a bootstrap method is used and the test result is determined by the power of H–L test at the target sample size. Simulation studies show that the proposed method can effectively standardize the power of the H–L test under the pre-specified level of type I error. Application to the two datasets illustrates the usefulness of the proposed model.     

On the choice of precedence tests

The best precedence test (BPT) is derived for testing the hypothesis that the lifetimes of two types of items on test have the same distribution. The test has maximum power in the class of the Lehmann type of alternatives F - 1 - (1-G) , A > 1, where F and G are probability distributions of the lifetimes of two types of items on test. This class includes exponential distributions, the Weibull distribution differing only in scale and distributions with proportional hazard rates. Exact power of the BPT is compared with other nonparametrie and parametric tests. The test may terminate before all the lifetimes of the items on test are recorded. In comparing with competing tests of equal size, the power functions are similar but a considerable number of items can be saved and the time on test can be reduced by using the BPT     

A simple algorithm for the adaptation of scores and power behavior of the corresponding rank test

The purpose of this paper is twofold:On one hand we want to give a very simple algorithm for evaluating a special rank estimator of the type given in Behnen, Neuhaus, and Ruymgaart (1983) for the approximate optimal choice of the scores-generating function of a two-sample linear rank test for the general testing problem H_o:F=G versus H₁:F ≤ G, F ≠ G, in order to demonstrate that the corresponding adaptive rank statistic is simple enough for practical applications. On the other hand we prove the asymptotic normality of the adaptive rank statistic under H (leading to approximate critical values) and demonstrate the adaptive behavior of the corresponding rank test by a Monte Carlo power simulation for sample sizes as low as m=10, n=10.     

Bootstrap LR tests of stationarity,common trends and cointegration

This paper considers the likelihood ratio (LR) tests of stationarity, common trends and cointegration for multivariate time series. As the distribution of these tests is not known, a bootstrap version is proposed via a state- space representation. The bootstrap samples are obtained from the Kalman filter innovations under the null hypothesis. Monte Carlo simulations for the Gaussian univariate random walk plus noise model show that the bootstrap LR test achieves higher power for medium-sized deviations from the null hypothesis than a locally optimal and one-sided Lagrange Multiplier (LM) test that has a known asymptotic distribution. The power gains of the bootstrap LR test are significantly larger for testing the hypothesis of common trends and cointegration in multivariate time series, as the alternative asymptotic procedure – obtained as an extension of the LM test of stationarity – does not possess properties of optimality. Finally, it is shown that the (pseudo-)LR tests maintain good size and power properties also for the non-Gaussian series. An empirical illustration is provided.     

Smooth nonparametric quantile estimation under censoring: simulations and bootstrap methods

Based on right-censored data from a lifetime distribution F , a smooth nonparametric estimator of the quantile function Q (p) is given by Qn(p)=h 1jjQn(t)K((t-p)/h)dt, where QR(p) denotes the product-limit quantile function. Extensive Monte Carlo simulations indicate that at fixed p this kernel-type quantile estimator has smaller mean squared error than (L(p) for a range of values of the bandwidth h. A method of selecting an "optimal" bandwidth (in the sense of small estimated mean squared error) based on the bootstrap is investigated yielding results consistent with the simulation study. The bootstrap is also used to obtain interval estimates for Q (p) after determining the optimal value of h.     

Asymptotic properties of the bootstrap unit root test statistic under possibly infinite variance

ABSTRACT

In this article, the unit root test for the AR(1) model is discussed, under the condition that the innovations of the model are in the domain of attraction of the normal law with possibly infinite variances. By using residual bootstrap with sample size m < n (n being the size of the original sample), we bootstrap the least-squares estimator of the autoregressive parameter. Under some mild assumptions, we prove that the null distribution of the unit root test statistic based on the least-square estimator of the autoregressive parameter can be approximated by using residual bootstrap.     

A new generalized two-sided class of distributions with an emphasis on two-sided generalized normal distribution

The ordinary-G class of distributions is defined to have the cumulative distribution function (cdf) as the value of the cdf of the ordinary distribution F whose range is the unit interval at G, that is, F(G), and it generalizes the ordinary distribution. In this work, we consider the standard two-sided power distribution to define other classes like the beta-G and the Kumaraswamy-G classes. We extend the idea of two-sidedness to other ordinary distributions like normal. After studying the basic properties of the new class in general setting, we consider the two-sided generalized normal distribution with maximum likelihood estimation procedure.     

The Size and Power of Bootstrap and Bartlett-Corrected Tests of Hypotheses on the Cointegrating Vectors

In this paper we compare Bartlett-corrected, bootstrap, and fast double bootstrap tests on maximum likelihood estimates of cointegration parameters. The key result is that both the bootstrap and the Bartlett-corrected tests must be based on the unrestricted estimates of the cointegrating vectors: procedures based on the restricted estimates have almost no power. The small sample size bias of the asymptotic test appears so severe as to advise strongly against its use with the sample sizes commonly available; the fast double bootstrap test minimizes size bias, while the Bartlett-corrected test is somehow more powerful.     

Percentage points of a weighted Kolmogorov-Smirnov statistic

Let F and G be the cumulative distribution functions corresponding to two independent random variables. Define the shift function, Δ(x), by F(x)=G(x+Δ(x)). Doksum and Sievers (1976) compared two confidence bands for Δ(x). The confidence band they found to be best requires the percentage points of a weighted form of the Kolmogorov-Smirnov statistic. The goal in this paper is to supply a table of some of the exact percentage points.     

Testing the rate ratio under inverse sampling based on gradient statistic

Inverse sampling is widely applied in studies with dichotomous outcomes, especially when the subjects arrive sequentially or the response of interest is difficult to obtain. In this paper, we investigate the rate ratio test problem under inverse sampling based on gradient statistic with the asymptotic method and parametric bootstrap technique. The gradient statistic has many advantages, for example, it is simple to calculate and competitive with Wald-type, score and likelihood ratio tests in terms of local power. Numerical studies are carried out to evaluate the performance of our gradient test and the existing tests, namely Wald-type, score and likelihood ratio tests. The simulation results suggest that the gradient test based on the parametric bootstrap method has excellent type I error control and large powers even in small sample design. Two real examples, from a heart disease study and a drug comparison study, are applied to illustrate our methods.     

The resampling of entropies with the application of biodiversity

This paper discusses the bootstrap test of entropies. Since the comparison of entropies is of prime interest in applied fields, finding an appropriate way to carry out such a comparison is of utmost importance. This paper presents how resampling should be performed to obtain an accurate p-value. Although the test using a pair-wise bootstrap confidence interval (CI) has already been dealt with in few works, here the bootstrap tests are studied because it may demand quite a different resampling algorithm compared with the CI. Moreover, the multiple test is studied. The proposed tests appear to yield several appreciable advantages. The easy implementation and the power of the proposed test can be considered as advantages. Here the entropy of the discrete variable is studied. The proposed tests are examined using Monte Carlo investigations and also evaluated using various distributions.     

An asymptotic sequential test based on confidence sequences

In the context of a translation parameter family of distributions F₀(x) = F(x-θ) an asymptotic sequential test of H₀: θ ≤ -△ versus H₁: θ ≥ △ developed. The test is based on confidence sequences. In the special case where F is a specified normal distribution the proposed test is uniformly at least as efficient (in the sense of Rechanter (1960)) as the Wald sequention probibilty ratio test.     

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.     

An approximation to the exact distribution of the wilcoxon-mann-whitney rank sum test statistic

An approximation to the exact distribution of the Wilcoxon rank sum test (Mann-Whitney U-test) and the Siegel-Tukey test based on a linear combination of the two-sample t-test applied to ranks and the normal approximation is compared with the usual normal approximation. The normal approximation results in a conservative test in the tails while the linear combination of the test statistics provides a test that has a very high percentage of agreement with tables of the exact distribution. Sample sizes 3≤m, n≤50 were considered.     

Direct estimation of the percentile 'p-value' for the one-sample median test

In this paper we outline and illustrate an easy-to-use inference procedure for directly calculating the approximate bootstrap percentile-type p-value for the one-sample median test, i.e. we calculate the bootstrap p -value without resampling, by using a fractional order statistics based approach. The method parallels earlier work on fractionalorder-statistics-based non-parametric bootstrap percentile-type confidence intervals for quantiles. Monte Carlo simulation studies are performed, which illustrate that the fractional-order-statistics-based approach to the one-sample median test has accurate type I error control for small samples over a wide range of distributions; is easy to calculate; and is preferable to the sign test in terms of type I error control and power. Furthermore, the fractional-order-statistics-based median test is easily generalized to testing that any quantile has some hypothesized value; for example, tests for the upper or lower quartile may be performed using the same framework.     

Testing the Equality of Multivariate Distributions Using the Bootstrap and Integrated Empirical Processes

We introduce some projected integrated empirical processes for testing the equality of two multivariate distributions. The bootstrap is used for determining the approximate critical values. We show that the bootstrap test is consistent. A number-theoretic method is used for efficient computation of the bootstrap critical values. Some simulation results are also given.