The glm log-rank test:general linear modeling of log-rank scores as a method of analysis for survival data

The use of general linear modeling (GLM) procedures based on log-rank scores is proposed for the analysis of survival data and compared to standard survival analysis procedures. For the comparison of two groups, this approach performed similarly to the traditional log-rank test. In the case of more complicated designs - without ties in the survival times - the approach was only marginally less powerful than tests from proportional hazards models, and clearly less powerful than a likelihood ratio test for a fully parametric model; however, with ties in the survival time, the approach proved more powerful than tests from Cox's semi-parametric proportional hazards procedure. The method appears to provide a reasonably powerful alternative for the analysis of survival data, is easily used in complicated study designs, avoids (semi-)parametric assumptions, and is quite computationally easy and inexpensive to employ.     

Testing axial symmetry and separability of lattice processes

Data collected on a rectangular lattice are common in many areas, and models used often make simplifying assumptions. These assumptions include axial symmetry in the spatial process and separability. Some different methods for testing axial symmetry and separability are considered. Using the sample periodogram is shown to provide some simple satisfactory tests of both hypotheses, but tests for separability given axial symmetry have low power for small lattices.     

Testing the monotonicity or convexity of a function using regression splines

Methods for smoothed isotonic or convex regression are useful in many applications. Sometimes the shape assumptions constitute a priori knowledge about the regression function, but often the shape is part of the research question. The authors propose tests for monotonicity and convexity using constrained and unconstrained regression splines. The tests have good large‐sample properties and the small‐sample behaviour is illustrated through simulations. Extensions to the partial linear model and the generalized regression model are presented. The Canadian Journal of Statistics 39: 89–107; 2011 © 2011 Statistical Society of Canada     

A MONTE CARLO COMPARISON OF ELLIPTICAL-THEORY AND OTHER TESTS FOR EQUALITY OF COVARIANCE MATRICES

A Monte Carlo study of the size and power of tests of equality of two covariance matrices is carried out. Tests based upon normality assumptions, elliptical distribution assumptions as well as distribution-free tests are compared. Samples are generated from normal, elliptical and non-elliptical populations. The elliptical-theory tests, in particular, have poor size properties for both elliptical distributions with moderate sample sizes and for non-elliptical distributions.     

BOOTSTRAP TESTS FOR THE ERROR DISTRIBUTION IN LINEAR AND NONPARAMETRIC REGRESSION MODELS

In this paper we investigate several tests for the hypothesis of a parametric form of the error distribution in the common linear and non‐parametric regression model, which are based on empirical processes of residuals. It is well known that tests in this context are not asymptotically distribution‐free and the parametric bootstrap is applied to deal with this problem. The performance of the resulting bootstrap test is investigated from an asymptotic point of view and by means of a simulation study. The results demonstrate that even for moderate sample sizes the parametric bootstrap provides a reliable and easy accessible solution to the problem of goodness‐of‐fit testing of assumptions regarding the error distribution in linear and non‐parametric regression models.     

Estimating and testing generalized linear models under inequality restrictions

We consider maximum likelihood estimation and likelihood ratio tests under inequality restrictions on the parameters. A special case are order restrictions, which may appear for example in connection with effects of an ordinal qualitative covariate. Our estimation approach is based on the principle of sequential quadratic programming, where the restricted estimate is computed iteratively and a quadratic optimization problem under inequality restrictions is solved in each iteration. Testing for inequality restrictions is based on the likelihood ratio principle. Under certain regularity assumptions the likelihood ratio test statistic is asymptotically distributed like a mixture of χ², where the weights are a function of the restrictions and the information matrix. A major problem in theory is that in general there is no unique least favourable point. We present some empirical findings on finite-sample behaviour of tests and apply the methods to examples from credit scoring and dentistry.     

Multivariate Tests of Mean-Variance Efficiency and Spanning With a Large Number of Assets and Time-Varying Covariances

We develop a finite-sample procedure to test the mean-variance efficiency and spanning hypotheses, without imposing any parametric assumptions on the distribution of model disturbances. In so doing, we provide an exact distribution-free method to test uniform linear restrictions in multivariate linear regression models. The framework allows for unknown forms of nonnormalities as well as time-varying conditional variances and covariances among the model disturbances. We derive exact bounds on the null distribution of joint F statistics to deal with the presence of nuisance parameters, and we show how to implement the resulting generalized nonparametric bounds tests with Monte Carlo resampling techniques. In sharp contrast to the usual tests that are not even computable when the number of test assets is too large, the power of the proposed test procedure potentially increases along both the time and cross-sectional dimensions.     

Performance of some multiple comparison tests under heteroscedasticity and dependency

In this article, 18 multiple comparison tests are compared according to powers and type I error measures under some violations of analysis of variance assumptions with a Monte Carlo simulation study. Considered violations of assumptions are heterogeneity in subgroup variances and dependency between subgroups. Various numbers of subgroups and subgroup sizes are considered simultaneously with the violations of assumptions. Simulation results are analysed by using visual inspection, graphical representations, decision-tree and correspondence analyses. Wide inferences are drawn on the behaviour of considered tests with respect to measures used. Some general suggestions are given on which tests should be used or avoided under violations of assumptions.     

Use of a more general model for bioavailability studies

Use of a multivariate model for bioavailability crossover studies is discussed. The model affords a convenient way to obtain from untransformed data exact confidence intervals and general hypothesis tests for ratios that are pertinent for comparing formulations. Previously, confidence intervals for the case of a two period study of two formulations have been presented. This paper treats the case of three or more formulations in a study in which each subject receives each formulation. The model requires weaker assumptions than the assumptions of other models that have been commonly used. In addition to inferences on ratios, it is pointed out how to obtain from the model tests of the hypotheses of equal formulation effects, equal period effects, and equal sequence effects. The use of these concepts for studies in which a formulation is received more than once by subjects isalso discussed.     

Prediction bands for the EDF of a future sample

We develop both nonparametric and parametric methods for obtaining prediction bands for the empirical distribution function (EDF) of a future sample. These methods yield simultaneous prediction intervals for all order statistics of the future sample, and they also correspond to tests for the two-sample problem. The nonparametric prediction bands correspond to the two-sample Kolmogorov-Smirnov test and related nonparametric tests, but the parametric prediction bands correspond to entirely new parametric two-sample tests. The parametric prediction bands tend to outperform the nonparametric bands when the parametric assumptions hold, but they may have true coverage probabilities well below their nominal levels when the parametric assumptions fail. A new computational algorithm is used to obtain critical values in the nonparametric case.     

PERMUTATION TESTS AND OTHER TEST STATISTICS FOR ILL-BEHAVED DATA: EXPERIENCE OF THE NINDS t-PA STROKE TRIAL

Permutation tests are often used to analyze data since they may not require one to make assumptions regarding the form of the distribution to have a random and independent sample selection. We initially considered a permutation test to assess the treatment effect on computed tomography lesion volume in the National Institute of Neurological Disorders and Stroke (NINDS) t-PA Stroke Trial, which has highly skewed data. However, we encountered difficulties in summarizing the permutation test results on the lesion volume. In this paper, we discuss some aspects of permutation tests and illustrate our findings. This experience with the NINDS t-PA Stroke Trial data emphasizes that permutation tests are useful for clinical trials and can be used to validate assumptions of an observed test statistic. The permutation test places fewer restrictions on the underlying distribution but is not always distribution-free or an exact test, especially for ill-behaved data. Quasi-likelihood estimation using the generalized estimating equation (GEE) approach on transformed data seems to be a good choice for analyzing CT lesion data, based on both its corresponding permutation test and its clinical interpretation.     

Bayes estimators of the average failure rates

In the present paper we have proposed a Bayesian approach for making inferences from accelerated life tests which do not require distributional assumptions     

On a new family of shifted logarithmic transformations

Most parametric statistical methods are based on a set of assumptions: normality, linearity and homoscedasticity. Transformation of a metric response is a popular method to meet these assumptions. In particular, transformation of the response of a linear model is a popular method when attempting to satisfy the Gaussian assumptions on the error components in the model. A particular problem with common transformations such as the logarithm or the Box–Cox family is that negative and zero data values cannot be transformed. This paper proposes a new transformation which allows negative and zero data values. The method for estimating the transformation parameter consider an objective criteria based on kurtosis and skewness for achieving normality. Use of the new transformation and the method for estimating the transformation parameter are illustrated with three data sets.     

Robust testing for stationarity of global surface temperature

Surface temperature is a major indicator of climate change. To test for the presence of an upward trend in surface-temperature (global warming), sophisticated statistical methods are typically used which depend on implausible and/or unverifiable assumptions, in particular on the availability of a very large number of measurements. In this paper, the validity of these methods is challenged. It is argued that the available series are simply not long enough to justify the use of methods which are based on asymptotic arguments, because only a small fraction of the information contained in the data is utilizable to distinguish between a trend and natural variability. Thus, a simple frequency-domain test is proposed for the case when all but a very small number of frequencies may be corrupted by transitory fluctuations. Simulations confirm its robustness against short-term autocorrelation. When applied to a global surface-temperature series, significance can be achieved with far fewer frequencies than required by conventional tests.     

Linear bootstrap methods for vector autoregressive moving-average models

We provide the theoretical justification of bootstrapping stationary invertible echelon vector autoregressive moving-average (VARMA) models using linear methods. The asymptotic validity of the bootstrap is established with strong white noise under parametric and nonparametric assumptions. Our methods are practical and useful for building reliable simulation-based inference and forecasting without implementing nonlinear estimation techniques such as ML which is usually burdensome, time demanding or impractical, particularly in big or highly persistent systems. The relevance of our procedures is more pronounced in the context of dynamic simulation-based techniques such as maximized Monte Carlo (MMC) tests [see Dufour J-M. Monte Carlo tests with nuisance parameters: a general approach to finite-sample inference and nonstandard asymptotics in econometrics. J Econom. 2006;133(2):443–477 and Dufour J-M, Jouini T. Finite-sample simulation-based tests in VAR models with applications to Granger causality testing. J Econom. 2006;135(1–2):229–254 for the VAR case]. Simulation evidence shows that, compared with conventional asymptotics, our bootstrap methods have good finite-sample properties in approximating the actual distribution of the studentized echelon VARMA parameter estimates, and in providing echelon parameter confidence sets with satisfactory coverage.     

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.     

Some Deficiencies of Using Moment Generating Functions

Dissatisfaction is expressed with the well-known definitions of partial and multiple correlation. Partial correlation is defined here as the ordinary correlation between two random variables after their linear dependence on the partialed-out variables has been subtracted. This leads to familiar results; it also clarifies the auto-regressivelike assumptions required for the partial autocorrelations to have a cut-off property. Multiple correlation is similarly considered.     

On quadratic estimation of heteroscedastic variances

Beginning with a brief introduction to the general theory the concept of Bayes invariant quadratic unbiased estimators (BAIQUEs) founded by Kleffe and Pingus(1974)is applied to combined samples with a common mean and different variances.Explicite formulas for Baique under these special assumptions are derived.Finally,some numerical comparisons of the variance function of Baiques under different prior distributions are given.     

Nearly exact tests in factorial experiments using the aligned rank transform

A procedure is studied that uses rank-transformed data to perform exact and estimated exact tests, which is an alternative to the commonly used F-ratio test procedure. First, a common parametric test statistic is computed using rank-transformed data, where two methods of ranking-ranks taken for the original observations and ranks taken after aligning the observations-are studied. Significance is then determined using either the exact permutation distribution of the statistic or an estimate of this distribution based on a random sample of all possible permutations. Simulation studies compare the performance of this method with the normal theory parametric F-test and the traditional rank transform procedure. Power and nominal type I error rates are compared under conditions when normal theory assumptions are satisfied, as well as when these assumptions are violated. The method is studied for a two-factor factorial arrangement of treatments in a completely randomized design and for a split-unit experiment. The power of the tests rivals the parametric F-test when normal theory assumptions are satisfied, and is usually superior when normal theory assumptions are not satisfied. Based on the evidence of this study, the exact aligned rank procedure appears to be the overall best choice for performing tests in a general factorial experiment.     

Testing for nonzero impulse responses in vector autoregressive processes

Impulse response functions are often used to investigate the relationships between the components of a VAR (vector autoregressive) process. A hypothesis of particular interest is that a variable does not react to impulses in another variable, i.e., the impulse responses are zero. Two types of tests for such hypotheses are considered. The first type is based on finite-order VAR assumptions and the second allows for possibly infinite-order processes. It is found that both types of tests have to be used cautiously because small sample and asymptotic properties may differ substantially.