Testing factor–covariate interaction in rank repeated-measures analysis of covariance models

This article derives the asymptotic properties of rank-based tests for the covariate effects in rank repeated-measures analysis of covariance (ANCOVA) models (Fan and Zhang 2017 Fan, C., and D. Zhang. 2017. Rank repeated measures analysis of covariance. Communications in Statistics - Theory and Methods 46:1158–83.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) employing generalized estimating equation (GEE) techniques. One interested application of the proposed tests is to check the validity of the assumption of homogeneous covariate effects in different levels of the factors. Performance of the proposed tests has been confirmed by simulation studies and illustrated using the famous seizure count data. While the article mainly focuses on interaction tests, the scope of the proposed tests includes testing any contrast of the covariate effect such as the null of no overall covariate effect.     

Testing homogeneity in semiparametric mixture case–control models

Parametric and semiparametric mixture models have been widely used in applications from many areas, and it is often of interest to test the homogeneity in these models. However, hypothesis testing is non standard due to the fact that several regularity conditions do not hold under the null hypothesis. We consider a semiparametric mixture case–control model, in the sense that the density ratio of two distributions is assumed to be of an exponential form, while the baseline density is unspecified. This model was first considered by Qin and Liang (2011 Qin, J., Liang, K.Y. (2011). Hypothesis testing in a mixture case–control model. Biometrics 67(1):182–198.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar], biometrics), and they proposed a modified score statistic for testing homogeneity. In this article, we consider alternative testing procedures based on supremum statistics, which could improve power against certain types of alternatives. We demonstrate the connection and comparison among the proposed and existing approaches. In addition, we provide a unified theoretical justification of the supremum test and other existing test statistics from an empirical likelihood perspective. The finite-sample performance of the supremum test statistics was evaluated in simulation studies.     

Testing lack of symmetry in spatial–temporal processes

Symmetry and separability of a covariance function are common assumptions to simplify the modeling effort of spatial–temporal processes. However, many studies in environmental sciences show that real data have complex spatial–temporal dependency structures resulting from lack of symmetry or violation of other standard assumptions of the covariance function. In this study, we propose new formal tests for lack of symmetry by using spectral representations of the spatial–temporal covariance functions of regularly spaced spatial–temporal data. The advantage of the proposed tests is that classical analysis of variance (ANOVA) models can be used for detecting lack of symmetry inherent in spatial–temporal processes. We evaluate the performance of the tests with simulation studies and we apply them to air pollution data.     

Testing for random effect in the Fuller–Battese model

We consider the Fuller–Battese model where random effects are allowed to be from non-normal universes. The asymptotic distribution of the F-statistic in this model is derived as the number of groups tends to infinity (is large) and sample size from any group is either fixed or large. The result is used to establish an approximate test for the significance of the random effect variance component. Robustness of the established approximate test is given.     

Optimal trajectory to the general equilibrium state of the unstable stochastic lattice gas of prey–predator model with pair-approximation

This paper focuses on the study of the optimal control problem of the unstable stochastic lattice gas of prey–predator model with pair-approximation. The Pontryagin minimum procedure is used to derive the optimal behaviour of the proposed system that makes it trajectory closer to its general equilibrium states with respect to a selected performance measure during a specified period of time T. Wide numerical simulations and illustrations are presented in graphs and tables.     

Investigation of age–treatment interaction in the SPACE trial using different statistical approaches

Selection of treatments to fit the specific needs for a certain patient is one major challenge in modern medicine. Personalized treatments rely on established patient–treatment interactions. In recent years, various statistical methods for the identification and estimation of interactions between relevant covariates and treatment were proposed. In this article, different available methods for detection and estimation of a covariate–treatment interaction for a time-to-event outcome, namely the standard Cox regression model assuming a linear interaction, the fractional polynomials approach for interaction, the modified outcome approach, the local partial-likelihood approach, and STEPP (Subpopulation Treatment Effect Pattern Plots) were applied to data from the SPACE trial, a randomized clinical trial comparing stent-protected angioplasty (CAS) to carotid endarterectomy (CEA) in patients with symptomatic stenosis, with the aim to analyse the interaction between age and treatment. Time from primary intervention to the first relevant event (any stroke or death) was considered as outcome parameter. The analyses suggest a qualitative interaction between patient age and treatment indicating a lower risk after treatment with CAS compared to CEA for younger patients, while for elderly patients a lower risk after CEA was observed. Differences in the statistical methods regarding the observed results, applicability, and interpretation are discussed.     

Testing Hypotheses in Nonparametric Models of Production

Data envelopment analysis (DEA) and free disposal hull (FDH) estimators are widely used to estimate efficiency of production. Practitioners use DEA estimators far more frequently than FDH estimators, implicitly assuming that production sets are convex. Moreover, use of the constant returns to scale (CRS) version of the DEA estimator requires an assumption of CRS. Although bootstrap methods have been developed for making inference about the efficiencies of individual units, until now no methods exist for making consistent inference about differences in mean efficiency across groups of producers or for testing hypotheses about model structure such as returns to scale or convexity of the production set. We use central limit theorem results from our previous work to develop additional theoretical results permitting consistent tests of model structure and provide Monte Carlo evidence on the performance of the tests in terms of size and power. In addition, the variable returns to scale version of the DEA estimator is proved to attain the faster convergence rate of the CRS-DEA estimator under CRS. Using a sample of U.S. commercial banks, we test and reject convexity of the production set, calling into question results from numerous banking studies that have imposed convexity assumptions. Supplementary materials for this article are available online.     

Testing for co-directional interactions using union–intersection and intersection–union methods

When interaction terms exist in a two-factor, factorial experiment, the consideration and analysis of main effects are often restricted to those situations where the interaction between factors is not significant. Hinkelman and Kempthorne [4 Hinkelmann, K. and Kempthorne, O. 1994. Design and Analysis of Experiments. Volume 1: Introduction to Experimental Design, New York: Wiley.  [Google Scholar]] softened that stance somewhat and advocate testing main effects when the interaction is deemed co-directional but not anti-directional. A test for the main effects in that situation may be pragmatic to the practitioner and appealing to researchers in other disciplines. Intersection–union and union–intersection methods are examined for assessing the directional nature of significant interactions so that the main effects in a two-factor factorial may be evaluated. The tests suggested are conceptually straightforward and practical and maintain the nominal Type-I error rate. Examples are provided to illustrate the methods.     

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.     

Testing for the partial Koziol–Green model with covariates

In this paper we consider some non-parametric goodness-of-fit statistics for testing the partial Koziol–Green regression model. In this model, the response at a given covariate value is subject to random right censoring by two independent censoring times. One of these censoring times is informative in the sense that its survival function is some power of the survival function of the response. The goodness-of-fit statistics are based on an underlying empirical process for which large sample theory is obtained.     

Testing the difference between two Kolmogorov–Smirnov values in the context of receiver operating characteristic curves

The maximum vertical distance between a receiver operating characteristic (ROC) curve and its chance diagonal is a common measure of effectiveness of the classifier that gives rise to this curve. This measure is known to be equivalent to a two-sample Kolmogorov–Smirnov statistic; so the absolute difference D between two such statistics is often used informally as a measure of difference between the corresponding classifiers. A significance test of D is of great practical interest, but the available Kolmogorov–Smirnov distribution theory precludes easy analytical construction of such a significance test. We, therefore, propose a Monte Carlo procedure for conducting the test, using the binormal model for the underlying ROC curves. We provide Splus/R routines for the computation, tabulate the results for a number of illustrative cases, apply the methods to some practical examples and discuss some implications.     

Inference for the common mean of several Birnbaum–Saunders populations

The Birnbaum–Saunders distribution is a widely used distribution in reliability applications to model failure times. For several samples from possible different Birnbaum–Saunders distributions, if their means can be considered as the same, it is of importance to make inference for the common mean. This paper presents procedures for interval estimation and hypothesis testing for the common mean of several Birnbaum–Saunders populations. The proposed approaches are hybrids between the generalized inference method and the large sample theory. Some simulation results are conducted to present the performance of the proposed approaches. The simulation results indicate that our proposed approaches perform well. Finally, the proposed approaches are applied to analyze a real example on the fatigue life of 6061-T6 aluminum coupons for illustration.     

Asymptotic bayesian inference in some nonstandard cases: Bernstein–von Mises Results and regular bayes' estimators

Asymptotically with probability close to one, the convergence in variation (also in distribution) to the multivariate normal, of the aposteriori density function of a parameter agains an apriori density, viz. the BERNSTEIN–VON MISES results are established when observations are not necessarily indenpendent or identically distributed but satisfy weak regularity assumptions on their joint density function. Regular BAYES' estimators are defined with respect to regular loss functions and a positive apriori density and proved consistent, asymptotically efficient and asymptotically normal. Examples and applications to conjugate families of densities, to inference in MARKOV Chains and other nonstandard cases illustrate results     

Characteristics of some classes of space–time covariance functions

Although a wide list of classes of space–time covariance functions is now available, selecting an appropriate class of models for a variable under study is still difficult and it represents a priority problem with respect to the choice of a particular model of a specified class. Then, knowing the characteristics of various classes of covariances, and their auxiliary functions, and matching those with the characteristics of the empirical space–time covariance surface might be helpful in the selection of a suitable class. In this paper some characteristics, such as behavior at the origin, asymptotic behavior, nonseparability and anisotropy aspects, are studied for some well known classes of covariance models of stationary space–time random fields. Moreover, some important issues related to modeling choices are described and a case study is presented.     

Power Maps in Goodness-of-fit Testing Based on Censored Samples

In goodness-of-fit testing based on Type-II right censored samples, mainly either classical tests are applied to transformed samples, which can be considered as full samples under the null hypothesis, or modifications of the tests are applied to the original samples. Results of an extensive simulation study are presented and discussed in order to compare the different approaches and tests by means of powers against beta alternatives. Depending on the beta parameters, illustrative plots lead to recommended tests in given situations, where there is some prior belief about the size of these parameters in the alternative. Moreover, the effects of adding artificial realizations via random dilation are examined.     

S–S method for stochastically ordered multinomial populations with missing data

Stochastic ordering is a useful concept in order restricted inferences. In this paper, we propose a new estimation technique for the parameters in two multinomial populations under stochastic orderings when missing data are present. In comparison with traditional maximum likelihood estimation method, our new method can guarantee the uniqueness of the maximum of the likelihood function. Furthermore, it does not depend on the choice of initial values for the parameters in contrast to the EM algorithm. Finally, we give the asymptotic distributions of the likelihood ratio statistics based on the new estimation method.     

Testing for the Marshall–Olkin extended form of the Weibull distribution

Marshall–Olkin extended distributions offer a wider range of behaviour than the basic distributions from which they are derived and therefore may find applications in modeling lifetime data, especially within proportional odds models, and elsewhere. The present paper carries out a simulation study of likelihood ratio, Wald and score tests for the parameter that distinguishes the extended distribution from the basic one, for the Weibull and exponential cases, allowing for right censored data. The likelihood ratio test is found to perform better than the others. The test is shown to have sufficient power to detect alternatives that correspond to interesting departures from the basic model and can be useful in modeling.