A nonparametric test to compare survival distributions with covariate adjustment

Summary.  The analysis of covariance is a technique that is used to improve the power of a k -sample test by adjusting for concomitant variables. If the end point is the time of survival, and some observations are right censored, the score statistic from the Cox proportional hazards model is the method that is most commonly used to test the equality of conditional hazard functions. In many situations, however, the proportional hazards model assumptions are not satisfied. Specifically, the relative risk function is not time invariant or represented as a log-linear function of the covariates. We propose an asymptotically valid k -sample test statistic to compare conditional hazard functions which does not require the assumption of proportional hazards, a parametric specification of the relative risk function or randomization of group assignment. Simulation results indicate that the performance of this statistic is satisfactory. The methodology is demonstrated on a data set in prostate cancer.     

A goodness-of-fit test for generalized error distribution

The Generalized Error Distribution is a widespread flexible family of symmetric probability distribution. Thanks to its properties it is becoming more and more popular in many science fields therefore determining if a sample is drawn from a GED is an important issue that usually is pursued with a graphical approach. In this paper we present a new goodness-of-fit test for GED that shows good performances for detecting non GED distribution when the alternative distribution is either skewed or a mixture. A comparison between well known tests and this new procedure is performed through a simulation study. We have developed a function that performs the analysis described in this paper in the R environment. The computational time required to compute this procedure is negligible.     

Wavelet goodness-of-fit test for dependent data

We introduce a new goodness-of-fit test which can be applied to hypothesis testing about the marginal distribution of dependent data. We derive a new test for the equivalent hypothesis in the space of wavelet coefficients. Such properties of the wavelet transform as orthogonality, localisation and sparsity make the hypothesis testing in wavelet domain easier than in the domain of distribution functions. We propose to test the null hypothesis separately at each wavelet decomposition level to overcome the problem of bi-dimensionality of wavelet indices and to be able to find the frequency where the empirical distribution function differs from the null in case the null hypothesis is rejected. We suggest a test statistic and state its asymptotic distribution under the null and under some of the alternative hypotheses.     

A nonparametric test for current status data with unequal censoring

Nonparametric tests for the comparison of different treatments based on current status data are proposed. For this problem, most methods proposed in the literature require that observation times on all subjects follow the same distribution. In other words, censoring distributions are identical between the treatment groups. In this paper, we focus on the situation where the censoring distributions may be different for subjects in different treatment groups and the test that can take this unequal censoring into account is given. The asymptotic distribution of the test proposed is derived. The method proposed is applied to data arising from a tumorigenicity experiment.     

Error variance function estimation in nonparametric regression models

In this article, a new class of variance function estimators is proposed in the setting of heteroscedastic nonparametric regression models. To obtain a variance function estimator, the main proposal is to smooth the product of the response variable and residuals as opposed to the squared residuals. The asymptotic properties of the proposed methodology are investigated in order to compare its asymptotic behavior with that of the existing methods. The finite sample performance of the proposed estimator is studied through simulation studies. The effect of the curvature of the mean function on its finite sample behavior is also discussed.     

A simple bera-jarque normality test for nonparametric residuals

This paper demonstrates that the Bera-Jarque normality test using residuals from nonparametric series estimates has the usual X²(2) asymptotic distribution. Monte Carlo results shows that the test has good properties for reasonably sized samples.     

A simple bera-jarque normality test for nonparametric residuals

This paper demonstrates that the Bera-Jarque normality test using residuals from nonparametric series estimates has the usual X ²(2) asymptotic distribution. Monte Carlo results shows that the test has good properties for reasonably sized samples.     

A robust and efficient estimation method for nonparametric models with jump points

Nonparametric models with jump points have been considered by many researchers. However, most existing methods based on least squares or likelihood are sensitive when there are outliers or the error distribution is heavy tailed. In this article, a local piecewise-modal method is proposed to estimate the regression function with jump points in nonparametric models, and a piecewise-modal EM algorithm is introduced to estimate the proposed estimator. Under some regular conditions, the large-sample theory is established for the proposed estimators. Several simulations are presented to evaluate the performances of the proposed method, which shows that the proposed estimator is more efficient than the local piecewise-polynomial regression estimator in the presence of outliers or heavy tail error distribution. What is more, the proposed procedure is asymptotically equivalent to the local piecewise-polynomial regression estimator under the assumption that the error distribution is a Gaussian distribution. The proposed method is further illustrated via the sea-level pressures.     

A new goodness-of-fit test for certain bivariate distributions applicable to traffic accidents

T. Cacoullos and H. Papageorgiou [On some bivariate probability models applicable to traffic accidents and fatalities, Int. Stat. Rev. 48 (1980) 345–356] studied a special class of bivariate discrete distributions appropriate for modeling traffic accidents, and fatalities resulting therefrom. The corresponding random variable may be written as $Z={(N,Y)}^{\prime }$, with $Y={\sum }_{j=1}^N{X}_j$, where ${\left\{X_j\right\}}_{j=1}^N$, are independent copies of a (discrete) random variable $X$, and $N$ is independent of ${\left\{X_j\right\}}_{j=1}^N$, and follows a Poisson law. If $X$ follows a Poisson law (resp. Binomial law), the resulting distribution is termed Poisson–Poisson (resp. Poisson–Binomial). L2-type goodness-of-fit statistics are constructed for the ‘general distribution’ of this kind, where $X$ may be an arbitrary discrete nonnegative random variable. The test statistics utilize a simple characterization involving the corresponding probability generating function, and are shown to be consistent. The proposed procedures are shown to perform satisfactorily in simulated data, while their application to accident data leads to positive conclusions regarding the modeling ability of this class of bivariate distributions.     

A cubic smoothing spline based lack of fit test for nonlinear regression models

The use of a statistic based on cubic spline smoothing is considered for testing nonlinear regression models for lack of fit. The statistic is defined to be the Euclidean squared norm of the smoothed residual vector obtained from fitting the nonlinear model, The asymptotic distribution of the statistic is derived under suitable smooth local alternatives and a numerical example is presented.     

An empirical likelihood goodness-of-fit test for time series

Summary. Standard goodness-of-fit tests for a parametric regression model against a series of nonparametric alternatives are based on residuals arising from a fitted model. When a parametric regression model is compared with a nonparametric model, goodness-of-fit testing can be naturally approached by evaluating the likelihood of the parametric model within a nonparametric framework. We employ the empirical likelihood for an α -mixing process to formulate a test statistic that measures the goodness of fit of a parametric regression model. The technique is based on a comparison with kernel smoothing estimators. The empirical likelihood formulation of the test has two attractive features. One is its automatic consideration of the variation that is associated with the nonparametric fit due to empirical likelihood's ability to Studentize internally. The other is that the asymptotic distribution of the test statistic is free of unknown parameters, avoiding plug-in estimation. We apply the test to a discretized diffusion model which has recently been considered in financial market analysis.     

Hypothesis testing in smoothing spline models

Nonparametric regression models are often used to check or suggest a parametric model. Several methods have been proposed to test the hypothesis of a parametric regression function against an alternative smoothing spline model. Some tests such as the locally most powerful (LMP) test by Cox et al. (Cox, D., Koh, E., Wahba, G. and Yandell, B. (1988). Testing the (parametric) null model hypothesis in (semiparametric) partial and generalized spline models. Ann. Stat., 16, 113–119.), the generalized maximum likelihood (GML) ratio test and the generalized cross validation (GCV) test by Wahba (Wahba, G. (1990). Spline models for observational data. CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM.) were developed from the corresponding Bayesian models. Their frequentist properties have not been studied. We conduct simulations to evaluate and compare finite sample performances. Simulation results show that the performances of these tests depend on the shape of the true function. The LMP and GML tests are more powerful for low frequency functions while the GCV test is more powerful for high frequency functions. For all test statistics, distributions under the null hypothesis are complicated. Computationally intensive Monte Carlo methods can be used to calculate null distributions. We also propose approximations to these null distributions and evaluate their performances by simulations.     

An optimization interpretation of integration and back-fitting estimators for separable nonparametric models

We provide an optimization interpretation of both back-fitting and integration estimators for additive nonparametric regression. We find that the integration estimator is a projection with respect to a product measure. We also provide further understanding of the back-fitting method.     

Statistical models for e-learning data

In this paper we analyse a real e-learning dataset derived from the e-learning platform of the University of Pavia. The dataset concerns an online learning environment with in-depth teaching materials. The main focus of this paper is to supply a measure of the relative importance of the exercises (test) at the end of each training unit; to build predictive models of student’s performance and finally to personalize the e-learning platform. The methodology employed is based on nonparametric statistical methods for kernel density estimation and generalized linear models and generalized additive models for predictive purposes.     

A gamma goodness-of-fit test based on characteristic independence of the mean and coefficient of variation

Characterization theorems in probability and statistics are widely appreciated for their role in clarifying the structure of the families of probability distributions. Less well known is the role characterization theorems have as a natural, logical and effective starting point for constructing goodness-of-fit tests. The characteristic independence of the mean and variance and of the mean and the third central moment of a normal sample were used, respectively, by Lin and Mudholkar [1980. A simple test for normality against asymmetric alternatives. Biometrika 67, 455–461] and by Mudholkar et al. [2002a. Independence characterizations and testing normality against skewness-kurtosis alternatives. J. Statist. Plann. Inference 104, 485–501] for developing tests of normality. The characteristic independence of the maximum likelihood estimates of the population parameters was similarly used by Mudholkar et al. [2002b. Independence characterization and inverse Gaussian goodness-of-fit. Sankhya A 63, 362–374] to develop a test of the composite inverse Gaussian hypothesis. The gamma models are extensively used for applied research in the areas of econometrics, engineering and biomedical sciences; but there are few goodness-of-fit tests available to test if the data indeed come from a gamma population. In this paper we employ Hwang and Hu's [1999. On a characterization of the gamma distribution: the independence of the sample mean and the sample coefficient of variation. Ann. Inst. Statist. Math. 51, 749–753] characterization of the gamma population in terms of the independence of sample mean and coefficient of variation for developing such a test. The asymptotic null distribution of the proposed test statistic is obtained and empirically refined for use with samples of moderate size.     

An omnibus test of goodness-of-fit for conditional distributions with applications to regression models

We introduce an omnibus goodness-of-fit test for statistical models for the conditional distribution of a random variable. In particular, this test is useful for assessing whether a regression model fits a data set on all its assumptions. The test is based on a generalization of the Cramér–von Mises statistic and involves a local polynomial estimator of the conditional distribution function. First, the uniform almost sure consistency of this estimator is established. Then, the asymptotic distribution of the test statistic is derived under the null hypothesis and under contiguous alternatives. The extension to the case where unknown parameters appear in the model is developed. A simulation study shows that the test has good power against some common departures encountered in regression models. Moreover, its power is comparable to that of other nonparametric tests designed to examine only specific departures.     

Gini index based goodness-of-fit test for the logistic distribution

To model growth curves in survival analysis and biological studies the logistic distribution has been widely used. In this article, we propose a goodness-of-fit test for the logistic distribution based on an estimate of the Gini index. The exact distribution of the proposed test statistic and also its asymptotic distribution are presented. In order to compute the proposed test statistic, parameters of the logistic distribution are estimated by approximate maximum likelihood estimators (AMLEs), which are simple explicit estimators. Through Monte Carlo simulations, power comparisons of the proposed test with some known competing tests are carried. Finally, an illustrative example is presented and analyzed.     

A simple nonparametric test for bivariate symmetry about a line

Over the years many researchers have dealt with testing the hypotheses of symmetry in univariate and multivariate distributions in the parametric and nonparametric setup. In a multivariate setup, there are several formulations of symmetry, for example, symmetry about an axis, joint symmetry, marginal symmetry, radial symmetry, symmetry about a known point, spherical symmetry, and elliptical symmetry among others. In this paper, for the bivariate case, we formulate a concept of symmetry about a straight line passing through the origin in a plane and accordingly develop a simple nonparametric test for testing the hypothesis of symmetry about a straight line. The proposed test is based on a measure of deviance between observed counts of bivariate samples in suitably defined pairs of sets. The exact null distribution and non-null distribution, for specified classes of alternatives, of the test statistics are obtained. The null distribution is tabulated for sample size from n=5 up to n=30. The null mean, null variance and the asymptotic null distributions of the proposed test statistics are also obtained. The empirical power of the proposed test is evaluated by simulating samples from the suitable class of bivariate distributions. The empirical findings suggest that the test performs reasonably well against various classes of asymmetric bivariate distributions. Further, it is advocated that the basic idea developed in this work can be easily adopted to test the hypotheses of exchangeability of bivariate random variables and also bivariate symmetry about a given axis which have been considered by several authors in the past.     

A method for high-dimensional smoothing

We consider the problem of the computation of smoothed additive functionals, which are some integrals with respect to the joint smoothing distribution. It is a key issue in inference for general state-space models as these quantities appear naturally for maximum likelihood parameter inference. The computation of smoothed additive functionals is very challenging as exact computations are not possible for non-linear non-Gaussian state-space models. It becomes even more difficult when the hidden state lies in a high dimensional space because traditional numerical methods suffer from the curse of dimensionality. We propose a new algorithm to efficiently calculate the smoothed additive functionals in an online manner for a specific family of high-dimensional state-space models in discrete time, which is named the Space–Time Forward Smoothing (STFS) algorithm. The cost of this algorithm is at least $O\left(N^2{d}^{2}T\right)$, which is polynomial in $d$. $T$ and $N$ denote the number of time steps and the number of particles respectively, while $d$ is the dimension of the hidden state space. Its superior performance over other existing methods is illustrated by various simulation studies. Moreover, STFS algorithm is successfully applied to perform Maximum Likelihood estimation for static model parameters both in an online and an offline manner.