Comparisons of Software Packages for Generalized Linear Multilevel Models

We review five software packages that can fit a generalized linear mixed model for data with more than a two-level structure and a multiple number of independent variables. These five packages are MLn, MLwiN, SAS Proc Mixed (Glimmix Macro), HLM, and VARCL. We first discuss the features of each of the five packages. These features include data input and data management, statistical model capabilities, output, and user friendliness. We then compare their performance on several simulated datasets.     

A Statistical Test of Change‐Point in Mean that Almost Surely Has Zero Error Probabilities

In this paper we develop a non‐conventional statistical test for the change‐point in a mean model by making use of an almost‐sure (a.s.) convergence (or strong convergence) result that we obtain, in respect of the difference between the sums of squared residuals under the null and alternative hypotheses. We prove that both types of error probabilities of the new test converge to zero almost surely when the sample size goes to infinity. This result does not hold for any conventional statistical test where the type I error probability, i.e. the significance level or the size, is prescribed at a low but non‐zero level (e.g. 0.05). The test developed is easy to use in practice, and is ready to be generalised to other change‐point models provided that the relevant almost‐sure convergence results are available. We also provide a simulation study in the paper to compare the new and conventional tests under different data scenarios. The results obtained are consistent with our asymptotic study. In addition we provide least squares estimators of those parameters used in the change‐point test together with their almost‐sure convergence properties.     

Nonparametric predictive comparison of lifetime data under progressive censoring

In reliability and lifetime testing, comparison of two groups of data is a common problem. It is often attractive, or even necessary, to make a quick and efficient decision in order to save time and costs. This paper presents a nonparametric predictive inference (NPI) approach to compare two groups, say X and Y, when one (or both) is (are) progressively censored. NPI can easily be applied to different types of progressive censoring schemes. NPI is a statistical approach based on few assumptions, with inferences strongly based on data and with uncertainty quantified via lower and upper probabilities. These inferences consider the event that the lifetime of a future unit from Y is greater than the lifetime of a future unit from X.     

Spatio-temporal analysis of a plant disease in a non-uniform crop: a Monte Carlo approach

Identification of the type of disease pattern and spread in a field is critical in epidemiological investigations of plant diseases. For example, an aggregation pattern of infected plants suggests that, at the time of observation, the pathogen is spreading from a proximal source. Conversely, a random pattern suggests a lack of spread from a proximal source. Most of the existing methods of spatial pattern analysis work with only one variety of plant at each location and with uniform genetic disease susceptibility across the field. Pecan orchards, used in this study, and other orchard crops are usually composed of different varieties with different levels of susceptibility to disease. A new measure is suggested to characterize the spatio-temporal transmission patterns of disease; a Monte Carlo test procedure is proposed to test whether the transmission of disease is random or aggregated. In addition, we propose a mixed-transmission model, which allows us to quantify the degree of aggregation effect.     

Nonstationarity and Level Shifts With an Application to Purchasing Power Parity

This study considers testing for a unit root in a time series characterized by a structural change in its mean. The analysis is in the spirit of Perron (1990a), who showed that the existence of such a shift in a stationary time series biases the usual tests for a unit root toward nonrejection. The approach is, however, different given that we suppose the date of the change to be unknown. The statistic of interest is then the minimal t statistic over all possible breakpoints in regressions similar to those proposed by Perron (1990a). Other related statistics are also discussed. We derive and tabulate the asymptotic distributions of interest. Most of the emphasis, however, is given to the tabulation of finite-sample critical values using simulation experiments. Particular attention is given to the effect, on the finite-sample critical values, of various procedures to select the appropriate order of the estimated autoregressions. We apply the tests to analyze the issue of purchasing power parity between the United States and the United Kingdom and also between the United States and Finland, whose real exchange rates are characterized by apparent shifts in level when using particular price indexes.     

SEMI-PARAMETRIC ANALYSIS OF COVARIANCE UNDER DEPENDENCE CONDITIONS WITHIN EACH GROUP

Consider the problem of covariance analysis based on regression models whose regression function is the sum of a linear and a non-parametric component. We propose a parametric and a non-parametric statistical test to compare the effects of the linear and non-parametric components, respectively, on the response variable in   L ≥ 2  groups. Serially correlated errors within each group are allowed. The first (second) test compares the differences between the estimates of the parametric (non-parametric) components of each group by means of a Mahalanobis  ( L ²)  distance. The asymptotic distribution of each statistic under the null hypothesis is obtained. A modest simulation study and an application to a real data set illustrate our methodology.     

Bias correction for the regression-based LM fractional integration test

This paper examines the finite-sample behavior of the Lagrange Multiplier (LM) test for fractional integration proposed by Breitung and Hassler (J. Econom. 110:167–185, 2002). We find by extensive Monte Carlo simulations that size distortions can be quite large in small samples. These are caused by a finite-sample bias towards the alternative. Analytic expressions for this bias are derived, based on which the test can easily be corrected.     

Exact Goodness‐of‐Fit Testing for the Ising Model

The Ising model is one of the simplest and most famous models of interacting systems. It was originally proposed to model ferromagnetic interactions in statistical physics and is now widely used to model spatial processes in many areas such as ecology, sociology, and genetics, usually without testing its goodness of fit. Here, we propose various test statistics and an exact goodness‐of‐fit test for the finite‐lattice Ising model. The theory of Markov bases has been developed in algebraic statistics for exact goodness‐of‐fit testing using a Monte Carlo approach. However, finding a Markov basis is often computationally intractable. Thus, we develop a Monte Carlo method for exact goodness‐of‐fit testing for the Ising model that avoids computing a Markov basis and also leads to a better connectivity of the Markov chain and hence to a faster convergence. We show how this method can be applied to analyze the spatial organization of receptors on the cell membrane.     

On hypothesis testing for Poisson processes: Regular case

ABSTRACT

We consider the problem of hypothesis testing in the situation when the firsthypothesis is simple and the second one is local one-sided composite. We describe the choice of the thresholds and the power functions of the Score Function test, of the General Likelihood Ratio test, of the Wald test, and of two Bayes tests in the situation when the intensity function of the observed inhomogeneous Poisson process is smooth with respect to the parameter. It is shown that almost all these tests are asymptotically uniformly most powerful. The results of numerical simulations are presented.     

The ideas of conditional confidence in the simplest setting

A recently developed framework for comparing the properties of various conditional procedures is studied in detail in the setting of testing between two simple hypotheses, where the ideas are most transparent. In that setting, possible goodness criteria are considered, and illustrations are given. The conditional confidence methodology, unlike the Bayes, fiducial, and likelihood techniques, presents a measure of conclusiveness which has frequentist interpretability; and, unlike traditional Neyman-Pearson procedures, the measure is highly data-dependent.