Checking a Semiparametric Additive Risk Model

McKeague and Sasieni [A partly parametric additive risk model. Biometrika 81 (1994) 501] propose a restriction of Aalen’s additive risk model by the additional hypothesis that some of the covariates have time-independent influence on the intensity of the observed counting process. We introduce goodness-of-fit tests for this semiparametric Aalen model. The asymptotic distribution properties of the test statistics are derived by means of martingale techniques. The tests can be adjusted to detect particular alternatives. As one of the most important alternatives we consider Cox’s proportional hazards model. We present simulation studies and an application to a real data set.     

On the Power of Bootstrapped Specification Tests

Abstract

Decisions based on econometric model estimates may not have the expected effect if the model is misspecified. Thus, specification tests should precede any analysis. Bierens' specification test is consistent and has optimality properties against some local alternatives. A shortcoming is that the test statistic is not distribution free, even asymptotically. This makes the test unfeasible. There have been many suggestions to circumvent this problem, including the use of upper bounds for the critical values. However, these suggestions lead to tests that lose power and optimality against local alternatives. In this paper we show that bootstrap methods allow us to recover power and optimality of Bierens' original test. Bootstrap also provides reliable p-values, which have a central role in Fisher's theory of hypothesis testing. The paper also includes a discussion of the properties of the bootstrap Nonlinear Least Squares Estimator under local alternatives.     

Data Driven Rank Test for Two-Sample Problem

Traditional linear rank tests are known to possess low power for large spectrum of alternatives. In this paper we introduce a new rank test possessing a considerably larger range of sensitivity than linear rank tests. The new test statistic is a sum of squares of some linear rank statistics while the number of summands is chosen via a data-based selection rule. Simulations show that the new test possesses high and stable power in situations when linear rank tests completely break down, while simultaneously it has almost the same power under alternatives which can be detected by standard linear rank tests. Our approach is illustrated by some practical examples. Theoretical support is given by deriving asymptotic null distribution of the test statistic and proving consistency of the new test under essentially any alternative.     

Regions of alternatives with high and low power for goodness-of-fit tests

In the present paper we find finite dimensional spaces W of alternatives with high power for a given class of tests and non-parametric alternatives. On the orthogonal complement of W the power function is flat. These methods can be used to reduce the dimension of interesting alternatives. We sketch a device how to calculate (approximately) an alternative with maximum power of a fixed test on a given ball of certain non-parametric alternatives.

The calculations are done within different asymptotic models specified by signal detection tests. Specific tests are Kolmogorov–Smirnov type tests, integral tests (like the Anderson and Darling test) and Rényi tests for hazard based models. The statistical meaning and interpretation of the spaces of alternatives with high power is discussed. These alternatives belong to least favorable directions of a class of statistical functionals which are linear combinations of quantile functions. For various cases their meaning is explained for parametric submodels, in particular for location alternatives.     

Distribution-free test in Tobit mean regression model

This paper discusses the problem of fitting a parametric model in Tobit mean regression models. The proposed test is based on the supremum of the Khamaladze-type transformation of a partial sum process of calibrated residuals. The asymptotic null distribution of this transformed process is shown to be the same as that of a time-transformed standard Brownian motion. Consistency of this sequence of tests against some fixed alternatives and asymptotic power under some local nonparametric alternatives are also discussed. Simulation studies are conducted to assess the finite sample performance of the proposed test. The power comparison with some existing tests shows some superiority of the proposed test at the chosen alternatives.     

Comparison of goodness of fit tests for the cox proportional hazards model

The proportional hazards regression model of Cox(1972) is widely used in analyzing survival data. We examine several goodness of fit tests for checking the proportionality of hazards in the Cox model with two-sample censored data, and compare the performance of these tests by a simulation study. The strengths and weaknesses of the tests are pointed out. The effects of the extent of random censoring on the size and power are also examined. Results of a simulation study demonstrate that Gill and Schumacher's test is most powerful against a broad range of monotone departures from the proportional hazards assumption, but it may not perform as well fail for alternatives of nonmonotone hazard ratio. For the latter kind of alternatives, Andersen's test may detect patterns of irregular changes in hazards.     

Optimal tests for carcinogenicity in a model with fatal and incidental tumours

In a recent paper Heimann and Neuhaus (Biometrika 88 (2001) 435) studied the combined test of Peto et al. (Long Term and Short Term Screening Assays for Carcinogens: A Critical Appraisal, IARC Monograph on the evaluation of the carcinogenic risk of chemicals to humans, Annex to supplement 2, WHO, Geneva, pp. 311–426) for the two sample testing problem. This test is an ad hoc test adding the log rank test statistic for the random censorship model and the Mantel Haenszel test statistic for the interval censored data model. Heimann and Neuhaus (Biometrika 88 (2001) 435) also introduced an appropriate statistical model combining randomly censored and interval censored data, to handle lethal and incidental tumours jointly.In the present paper we will derive asymptotically optimal tests for this model in a systematic way. The optimal test statistic is decomposed into two orthogonal parts, representing the fatal and incidental components. This will allow to use well known results from the literature for the fatal part, to estimate unknown quantities for the incidental part, and last but not least to compare with Peto's combined test.Moreover, we introduce conditional permutation versions of our tests which are finite sample distribution free under the null hypothesis with equal censoring and are asymptotically equivalent to their unconditional counterparts even under locally unequal censoring. Crucial for these results is a conditional limit theorem for the test statistics under local alternatives which follows from results of Strasser and Weber (Math. Methods Statist. 2 (1999) 220).     

Conditional Score Tests for Heteroscedasticity in the Two-Way Error Components Model

This paper extends the one-way heteroskedasticity score test of Holly and Gardiol (2000, In: Krishnakumar, J, Ronchetti, E (Eds.), Panel Data Econometrics: Future Directions, North-Holland, Amsterdam, pp. 199–211) to two conditional Lagrange Multiplier (LM) tests of heteroskedasticity under contiguous alternatives within the two-way error components model framework. In each case, the derivation of Rao's efficient score statistics for testing heteroskedasticity is first obtained. Then, based on a specific set of assumptions, the asymptotic distribution of the score under contiguous alternatives is established. Finally, the expression for the score test statistic in the presence of heteroskedasticity and related asymptotic local powers of these score test statistics are derived and discussed.     

Moments of nonparametric probability density functions of entropy estimators applied to testing the inverse Gaussian distribution

The inverse Gaussian (IG) distribution is widely used to model data and then it is important to develop efficient goodness of fit tests for this distribution. In this article, we introduce some new test statistics for examining the IG goodness of fit based on correcting moments of nonparametric probability density functions of entropy estimators. These tests are consistent against all alternatives. Critical points and power of the tests are explored by simulation. We show that the proposed tests are more powerful than competitor tests. Finally, the proposed tests are illustrated by real data examples.     

Testing for lack of fit in linear multiresponse models based on exact or near replicates

Three procedures for testing the adequacy of a proposed linear multiresponse regression model against unspecified general alternatives are considered. The model has an error structure with a matrix normal distribution which allows the vector of responses for a particular run to have an unknown covariance matrix while the responses for different runs are uncorrelated. Furthermore, each response variable may be modeled by a separate design matrix. Multivariate statistics corresponding to the classical univariate lack of fit and pure error sums of squares are defined and used to determine the multivariate lack of fit tests. A simulation study was performed to compare the power functions of the test procedures in the case of replication. Generalizations of the tests for the case in which there are no independent replicates on all responses are also presented.     

Testing generalized linear and semiparametric models against smooth alternatives

We propose goodness-of-fit tests for testing generalized linear models and semiparametric regression models against smooth alternatives. The focus is on models having both continous and factorial covariates. As a smooth extension of a parametric or semiparametric model we use generalized varying-coefficient models as proposed by Hastie and Tibshirani. A likelihood ratio statistic is used for testing. Asymptotic expansions allow us to write the estimates as linear smoothers which in turn guarantees simple and fast bootstrapping of the test statistic. The test is shown to have √ n -power, but in contrast with parametric tests it is powerful against smooth alternatives in general.     

Wavelet semi-parametric inference for long memory in volatility in the presence of a trend

Risk of investing in a financial asset is quantified by functionals of squared returns. Discrete time stochastic volatility (SV) models impose a convenient and practically relevant time series dependence structure on the log-squared returns. Different long-term risk characteristics are postulated by short-memory SV and long-memory SV models. It is therefore important to test which of these two alternatives is suitable for a specific asset. Most standard tests are confounded by deterministic trends. This paper introduces a new, wavelet-based, test of the null hypothesis of short versus long memory in volatility which is robust to deterministic trends. In finite samples, the test performs better than currently available tests which are based on the Fourier transform.     

Score test for testing zero-inflated Poisson regression against zero-inflated generalized Poisson alternatives

In several cases, count data often have excessive number of zero outcomes. This zero-inflated phenomenon is a specific cause of overdispersion, and zero-inflated Poisson regression model (ZIP) has been proposed for accommodating zero-inflated data. However, if the data continue to suggest additional overdispersion, zero-inflated negative binomial (ZINB) and zero-inflated generalized Poisson (ZIGP) regression models have been considered as alternatives. This study proposes the score test for testing ZIP regression model against ZIGP alternatives and proves that it is equal to the score test for testing ZIP regression model against ZINB alternatives. The advantage of using the score test over other alternative tests such as likelihood ratio and Wald is that the score test can be used to determine whether a more complex model is appropriate without fitting the more complex model. Applications of the proposed score test on several datasets are also illustrated.     

Comparison of powers and exact bahadur slopes of some tests for dimensionality in manova models

In this paper we consider test of dimensionality in MANOVA model. For this testing problem, the likelihood ratio (=LR) test, Lawley-Hotelling (=LH) type test and Bartlett-Nanda-Pillai (=BNP) type test are often used. We obtain the asymptotic expansions of powers of these tests under the local alternatives. Also Bahadur exact slopes of these tests are obtained. Based on these results, we obtain a unified opinion concerning comparison of LR test, LH type test and BNP type test.     

A Test for Constant Hazard Against a Change-Point Alternative

In this article, we consider the change-point hazard rate model which arises quite commonly in mechanical or biological systems, which experience a high hazard rate early in their lifetime due to infant mortality and then a constant or steady hazard rate after the threshold time. We first derive the corresponding mean residual life function (MRLF) and observe that the MRLF is initially increasing and then constant. Here, we derive a test statistic for exponentiality against Increasing Initially then Constant Mean Residual Life (ICMRL). We also derive the asymptotic distribution of the test statistic and compare the power of the test with other existing tests such as likelihood ratio, Weibull, and Log gamma tests considered in the literature. The test performs quite well as compared to other alternatives studied.     

Testing for Serial Correlation: Generalized Andrews–Ploberger Tests

This paper considers testing the null hypothesis that a times series is uncorrelated when the time series is uncorrelated but statistically dependent. This case is of interest in economic and finance applications. The GARCH(1, 1) model is a leading example of a model that generates serially uncorrelated but statistically dependent data. The tests of serial correlation introduced by Andrews and Ploberger (1996, hereafter AP) are generalized for the purpose of testing the null. The rationale for generalizing the AP tests is that they have attractive properties for cases for which they were originally designed: they are consistent against all nonwhite-noise alternatives and have good all-round power against nonseasonal alternatives compared to several widely used tests in the literature. These properties are inherited by the generalized AP tests.     

Relative Risk Regression for Binary Outcomes: Methods and Recommendations

Relative risks are often considered preferable to odds ratios for quantifying the association between a predictor and a binary outcome. Relative risk regression is an alternative to logistic regression where the parameters are relative risks rather than odds ratios. It uses a log link binomial generalised linear model, or log‐binomial model, which requires parameter constraints to prevent probabilities from exceeding 1. This leads to numerical problems with standard approaches for finding the maximum likelihood estimate (MLE), such as Fisher scoring, and has motivated various non‐MLE approaches. In this paper we discuss the roles of the MLE and its main competitors for relative risk regression. It is argued that reliable alternatives to Fisher scoring mean that numerical issues are no longer a motivation for non‐MLE methods. Nonetheless, non‐MLE methods may be worthwhile for other reasons and we evaluate this possibility for alternatives within a class of quasi‐likelihood methods. The MLE obtained using a reliable computational method is recommended, but this approach requires bootstrapping when estimates are on the parameter space boundary. If convenience is paramount, then quasi‐likelihood estimation can be a good alternative, although parameter constraints may be violated. Sensitivity to model misspecification and outliers is also discussed along with recommendations and priorities for future research.     

Most powerful tests for use in matched pair experiments when data may be censored

Four distribution-free tests are developed for use in matched pair experiments when data may be censored: a bootstrap based on estimates of the median difference, and three rerandomization tests. The latter include a globally almost most powerful (GAMP) test which uses the original data and two modified Gilbert-Gehan tests which use the ranks. Computation time is reduced by using a binary count to generate subsamples and by restricting subsampling to the uncensored pairs. In Monte Carlo simulations against normal alternatives, mixed normal alternatives, and exponential alternatives, the GAMP test is most powerful with light censoring, the rank test is most powerful with heavy censoring. The bootstrap degenerates to the sign test and is least powerful.     

Robust Testing Procedures in Heteroscedastic Linear Models

There exist many studies which treat the robust tests in homoscedastic linear models. However, the robust testing procedure in heteroscedastic linear models has not been examined. In this article, three classes of testing procedures for testing subhypothesis in heteroscedastic linear models are developed. These are Wald-type, score-type, and drop-in dispersion tests. The asymptotic distributions of these tests are obtained under the null hypothesis and contiguous alternatives. For a robustness criterion, the maximum asymptotic bias of the level of the test for distributions in a shrinking contamination neighborhood is used and the most-efficient robust test is derived. Finally, the performance of these tests in small sample is studied by Monte Carlo simulation.     

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.