Nonparametric testing for simple regression under progressive censoring with stagering entry and random withdrawal

In the context of (multi-center) clinical trials and life testins problems, a general model incorporating both the staggering entry and random withdrawal and pertaining to a simple regression problem (including the two-sample location problem as a special case) is conceived, and, within this framework, a scheme allowing progressive censoring (continuous monitoring of experimentation from the beginning) is developed along with the proposal for and study of some nonparametric testing procedures, The proposed tests rest on the construction of certain two-dimensional time-parameter stochastic processes from a triangular array of progressively censored linear rank statistics and their weak convergence to appropriate Gaussian functions. Asymptotic properties of these procedures are studied. A computer program pertaining to the numerical computations and practical administrations of these testing procedures is also provided at the end.     

Model Order Estimation of a Multivariable Stochastic Process

ABSTRACT

This article presents a procedure allowing us to estimate the minimal order of a state-space representation, for a multivariable stochastic process, from a sequence of observations. The method proposes a statistical rule for testing the rank of a block Hankel matrix of data, since this rank is related to the order of the process. A new information criterion is then developed and used to decide upon the order of the model. In this article we generalize the Aoki C-test. Using two representative data sets as the basis for a Monte Carlo experiment and real data based on Danish economy, we estimate the order of multivariable stochastic processes.     

Generalised Rank Regression Estimator with Standard Error Adjusted Lasso

One of the standard variable selection procedures in multiple linear regression is to use a penalisation technique in least‐squares (LS) analysis. In this setting, many different types of penalties have been introduced to achieve variable selection. It is well known that LS analysis is sensitive to outliers, and consequently outliers can present serious problems for the classical variable selection procedures. Since rank‐based procedures have desirable robustness properties compared to LS procedures, we propose a rank‐based adaptive lasso‐type penalised regression estimator and a corresponding variable selection procedure for linear regression models. The proposed estimator and variable selection procedure are robust against outliers in both response and predictor space. Furthermore, since rank regression can yield unstable estimators in the presence of multicollinearity, in order to provide inference that is robust against multicollinearity, we adjust the penalty term in the adaptive lasso function by incorporating the standard errors of the rank estimator. The theoretical properties of the proposed procedures are established and their performances are investigated by means of simulations. Finally, the estimator and variable selection procedure are applied to the Plasma Beta‐Carotene Level data set.     

Expectation of rank statistics under setup of stochastic inequalities

The paper presents a general result on the order of means of rank scores under setup of stochastic inequalities between the original observations. The result also describes the order of the ranks themselves. Some illustrative applications in the change-point problem are introduced.     

Testing order restricted hypotheses with proportional hazards

Inferences for survival curves based on right censored continuous or grouped data are studied. Testing homogeneity with an ordered restricted alternative and testing the order restriction as the null hypothesis are considered. Under a proportional hazards model, the ordering on the survival curves corresponds to an ordering on the regression coefficients. Approximate likelihood methods are obtained by applying order restricted procedures to the estimates of the regression coefficients. Ordered analogues to the log rank test which are based on the score statistics are considered also. Chi-bar-squared distributions, which have been studied extensively, are shown to provide reasonable approximations to the null distributions of these tests statistics. Using Monte Carlo techniques, the powers of these two types of tests are compared with those that are available in the literature.     

Weighted Mean Survival Test Statistics: a Class of Distance Tests for Censored Survival Data

A class of test statistics is introduced which is sensitive against the alternative of stochastic ordering in the two-sample censored data problem. The test statistics for evaluating a cumulative weighted difference in survival distributions are developed while taking into account the imbalances in base-line covariates between two groups. This procedure can be used to test the null hypothesis of no treatment effect, especially when base-line hazards cross and prognostic covariates need to be adjusted. The statistics are semiparametric, not rank based, and can be written as integrated weighted differences in estimated survival functions, where these survival estimates are adjusted for covariate imbalances. The asymptotic distribution theory of the tests is developed, yielding test procedures that are shown to be consistent under a fixed alternative. The choice of weight function is discussed and relies on stability and interpretability considerations. An example taken from a clinical trial for acquired immune deficiency syndrome is presented.     

A Monte Carlo comparison of alternative methods of maximum likelihood ranking in racing sports

Applications of maximum likelihood techniques to rank competitors in sports are commonly based on the assumption that each competitor's performance is a function of a deterministic component that represents inherent ability and a stochastic component that the competitor has limited control over. Perhaps based on an appeal to the central limit theorem, the stochastic component of performance has often been assumed to be a normal random variable. However, in the context of a racing sport, this assumption is problematic because the resulting model is the computationally difficult rank-ordered probit. Although a rank-ordered logit is a viable alternative, a Thurstonian paired-comparison model could also be applied. The purpose of this analysis was to compare the performance of the rank-ordered logit and Thurstonian paired-comparison models given the objective of ranking competitors based on ability. Monte Carlo simulations were used to generate race results based on a known ranking of competitors, assign rankings from the results of the two models, and judge performance based on Spearman's rank correlation coefficient. Results suggest that in many applications, a Thurstonian model can outperform a rank-ordered logit if each competitor's performance is normally distributed.     

Nonparametric multiple comparison procedures for unbalanced one-way factorial designs

In this paper, we present several nonparametric multiple comparison (MC) procedures for unbalanced one-way factorial designs. The nonparametric hypotheses are formulated by using normalized distribution functions and the comparisons are carried out on the basis of the relative treatment effects. The proposed test statistics take the form of linear pseudo rank statistics and the asymptotic joint distribution of the pseudo rank statistics for testing treatments versus control satisfies the multivariate totally positive of order two condition irrespective of the correlations among the rank statistics. Therefore, in the context of MCs of treatments versus control, the nonparametric Simes test is validated for the global testing of the intersection hypothesis. For simultaneous testing of individual hypotheses, the nonparametric Hochberg stepup procedure strongly controls the familywise type I error rate asymptotically. With regard to all pairwise comparisons, we generalize various single-step and stagewise procedures to perform comparisons on the relative treatment effects. To further compare with normal theory counterparts, the asymptotic relative efficiencies of the nonparametric MC procedures with respect to the parametric MC procedures are derived under a sequence of Pitman alternatives in a nonparametric location shift model for unbalanced one-way layouts. Monte Carlo simulations are conducted to demonstrate the validity and power of the proposed nonparametric MC procedures.     

Rank tests for independence based on partially right-censored pairs

Linear rank procedures are developed for testing independence with right-censored matched pairs. It is assumed that censoring Is Independent of the random variables under study. The test statistics are derived as score statistics (Hajek and Sidak, 1967) based on the probability of the generalised rank vectors (Prentice, 1978). Applications to survival data analysis are also discussed.     

Goodness-of-fit tests for centralized Wishart processes

Abstract

In this paper we present several goodness-of-fit tests for the centralized Wishart process, a popular matrix-variate time series model used to capture the stochastic properties of realized covariance matrices. The new test procedures are based on the extended Bartlett decomposition derived from the properties of the Wishart distribution and allows to obtain sets of independently and standard normally distributed random variables under the null hypothesis. Several tests for normality and independence are then applied to these variables in order to support or to reject the underlying assumption of a centralized Wishart process. In order to investigate the influence of estimated parameters on the suggested testing procedures in the finite-sample case, a simulation study is conducted. Finally, the new test methods are applied to real data consisting of realized covariance matrices computed for the returns on six assets traded on the New York Stock Exchange.     

Critical values of some nonparametric statistics used for tests of total stochastic independence

The problem of testing for total independence of the variates of a stochastic p(≧3) component vector using rank correlation statistics is considered. Two distribution free statistics are considered, one based on the determinant of the matrix of rank correlation statistics, the second on their sum of squares. Tables of critical values are given for p=3,4 for the cases when (a) ranks, and (b) exponential scores are used to replace the ordered observations within each variate. Some approximations to the critical values are proposed and evaluated.     

Comparing Survival Times for Treatments with Those for a Control Under Proportional Hazards

Inferences for survival curves based on right censored data are studied for situations in which it is believed that the treatments have survival times at least as large as the control or at least as small as the control. Testing homogeneity with the appropriate order restricted alternative and testing the order restriction as the null hypothesis are considered. Under a proportional hazards model, the ordering on the survival curves corresponds to an ordering on the regression coefficients. Approximate likelihood methods, which are obtained by applying order restricted procedures to the estimates of the regression coefficients, and ordered analogues to the log rank test, which are based on the score statistics, are considered. Mau's (1988) test, which does not require proportional hazards, is extended to this ordering on the survival curves. Using Monte Carlo techniques, the type I error rates are found to be close to the nominal level and the powers of these tests are compared. Other order restrictions on the survival curves are discussed briefly.     

SELECTION PROCEDURES FOR SCALE PARAMETERS USING TWO-SAMPLE U-STATISTICS

Let be k independent populations having the same known quantile of order p (0 p 1) and let F(x)=F(x/_i) be the absolutely continuous cumulative distribution function of the ith population indexed by the scale parameter ₁, i = 1,…, k. We propose subset selection procedures based on two-sample U-statistics for selecting a subset of k populations containing the one associated with the smallest scale parameter. These procedures are compared with the subset selection procedures based on two-sample linear rank statistics given by Gill & Mehta (1989) in the sense of Pitman asymptotic relative efficiency, with interesting results.     

Some invariance principles for mized rank statistics and induced order statistics and some application

The structural affinity of mixed rank statistics and linear combinations of functions of concomitants of order statistics (or induced order statistics) is examined here. Some weal as well as strong invariance principles for these statistics are studied. A variety of models (depend on the nature of stochastic dependence of the two variates) is considered and the regularity conditions are tailored for these diverse situations. Some possible applications of these results in some problems of sequential (statistical) inference are also considered.     

P-splines with derivative based penalties and tensor product smoothing of unevenly distributed data

The P-splines of Eilers and Marx (Stat Sci 11:89–121, 1996) combine a B-spline basis with a discrete quadratic penalty on the basis coefficients, to produce a reduced rank spline like smoother. P-splines have three properties that make them very popular as reduced rank smoothers: (i) the basis and the penalty are sparse, enabling efficient computation, especially for Bayesian stochastic simulation; (ii) it is possible to flexibly ‘mix-and-match’ the order of B-spline basis and penalty, rather than the order of penalty controlling the order of the basis as in spline smoothing; (iii) it is very easy to set up the B-spline basis functions and penalties. The discrete penalties are somewhat less interpretable in terms of function shape than the traditional derivative based spline penalties, but tend towards penalties proportional to traditional spline penalties in the limit of large basis size. However part of the point of P-splines is not to use a large basis size. In addition the spline basis functions arise from solving functional optimization problems involving derivative based penalties, so moving to discrete penalties for smoothing may not always be desirable. The purpose of this note is to point out that the three properties of basis-penalty sparsity, mix-and-match penalization and ease of setup are readily obtainable with B-splines subject to derivative based penalization. The penalty setup typically requires a few lines of code, rather than the two lines typically required for P-splines, but this one off disadvantage seems to be the only one associated with using derivative based penalties. As an example application, it is shown how basis-penalty sparsity enables efficient computation with tensor product smoothers of scattered data.     

Adaptive choice of scale tests in flexible two-stage designs with applications in experimental ecology and clinical trials

In this paper, the two-sample scale problem is addressed within the rank framework which does not require to specify the underlying continuous distribution. However, since the power of a rank test depends on the underlying distribution, it would be very useful for the researcher to have some information on it in order to use the possibly most suitable test. A two-stage adaptive design is used with adaptive tests where the data from the first stage are used to compute a selector statistic to select the test statistic for stage 2. More precisely, an adaptive scale test due to Hall and Padmanabhan and its components are considered in one-stage and several adaptive and non-adaptive two-stage procedures. A simulation study shows that the two-stage test with the adaptive choice in the second stage and with Liptak combination, when it is not more powerful than the corresponding one-stage test, shows, however, a quite similar power behavior. The test procedures are illustrated using two ecological applications and a clinical trial.     

Some useful nonparametric tests for ordered alternatives in randomized block experiments

Testing for ordered alternatives in randomized block designs has been a problem of interest for almost three decades (Jonckheere (1954)). Three classes of rank tests have evolved—tests based on “within-blocks” rankings (W-tests), tests based on “ranking after alignment” within blocks (RAA-tests), and tests based on “among-blocks” rankings (A-Tests). This paper focuses on the latter. A simplified version of the Skillings-Wolfe generalized Purl test (1977) is suggested and two very useful A-tests—a generalized Johnson-Mehrotra “Optimal contrast” procedure and a generalized Tryon-Hettmansperger rank test—are developed. These procedures are compared and contrasted with other recent competitors presented by Skllllngs and Wolfe (1978) and by Salama and Quade (1981).     

Tail-Adaptive Location Rank Test for the Generalized Secant Hyperbolic Distribution

The generalized secant hyperbolic distribution (GSHD) was recently introduced as a modeling tool in data analysis. The GSHD is a unimodal distribution that is completely specified by location, scale, and shape parameters. It has also been shown elsewhere that the rank procedures of location are regular, robust, and asymptotically fully efficient. In this article, we study certain tail weight measures for the GSHD and introduce a tail-adaptive rank procedure of location based on those tail weight measures. We investigate the properties of the new adaptive rank procedure and compare it to some conventional estimators.     

Comparisons of tests for multivariate cointegration

This paper compares the small sample properties of different tests for multivariate cointegration like Johansen's trace test, stock &; Watson's common trend test, Phillips &; Ouliaris' principal component test, as well as cointegration rank decisions based on order selection criteria. Under the null hypothesis of non-cointegration we find a slow convergence rate of the test statistics. In bivariate models the Phillips &; Ouliaris test is extremely dependent on the specification and is outperformed by the other procedures. For trivariate processes we find dependence of the power results on the dynamic specification. The lag order is successfully estimated by order selection criteria.     

Vector autoregression and causality: a theoretical overview and simulation study

This paper provides a theoretical overview of Wald tests for Granger causality in levels vector autoregressions (VAR's) and Johansen-type error correction models (ECM's). The theory is based on results in Toda and Phillips (1991a) and allows for stochastic and deterministic trends as well as arbitrary degrees of cointegration. We recommend some operational procedures for conducting Granger causality tests that are based on the Gaussian maximum likelihood estimation of ECM's. These procedures are applicable in the important practical case of testing the causal effects of one variable on another group of variables and vice versa. This paper also investigates the sampling properties of these testing procedures through simulation exercises. Three sequential causality tests in ECM's are compared with conventional causality tests in levels and differences VAR's.