On asymptotic non-normality of null distributions of mrpp statistics

Severe departures from normality occur frequently for null distributions of statistics associated with applications of mulLi-response permutation procedures (MRPP) for either small or large finite populations. This paper describes the commonly encountered situation associated with asymptotic non-normality for null distributions of MRPP statistics which does not depend on the underlying multivariate distribution. In addition, this paper establishes the existence of a non-degenerate underlying distribution for which the null distributions of MRPP statistics are asymptotically non-normal for essentially all size structure configurations. It is known that MRPP statistics are symmetric versions of a broader class of statistics, most of which are asymmetric. Because of the non-normality associated with null distributions of MRPP statistics, this paper includes necessary results for inferences based on the exact first three moments of anv statistic in this broader class (analogous to existing results for MRPP statistics).     

A generalized scheffé method of multiple comparisons

Standard resulrs on the extrema of quotients of quadratic forms are extended to the non-negative definite case. The maximum and the set over which it is achieved are characterized explicitly both in terms of generalized inverse matrices and generalized eigenvalues. These results become the basis of Scheffe type multiple comparisons in the usual way. To demonstrate their application to statistics with singular covariance matrices, the method is detailed for Mantel-Haenszel, Breslow, and Cox statistics. An example is presented illustrating a situation where the proposed Scheffe type comparisons may be better than the pairwise method.     

A time-varying measurement error model for age of onset of a psychiatric diagnosis: applied to first depressive episode diagnosed by the Composite International Diagnostic Interview (CIDI)

A substantial degree of uncertainty exists surrounding the reconstruction of events based on memory recall. This form of measurement error affects the performance of structured interviews such as the Composite International Diagnostic Interview (CIDI), an important tool to assess mental health in the community. Measurement error probably explains the discrepancy in estimates between longitudinal studies with repeated assessments (the gold-standard), yielding approximately constant rates of depression, versus cross-sectional studies which often find increasing rates closer in time to the interview. Repeated assessments of current status (or recent history) are more reliable than reconstruction of a person's psychiatric history based on a single interview. In this paper, we demonstrate a method of estimating a time-varying measurement error distribution in the age of onset of an initial depressive episode, as diagnosed by the CIDI, based on an assumption regarding age-specific incidence rates. High-dimensional non-parametric estimation is achieved by the EM-algorithm with smoothing. The method is applied to data from a Norwegian mental health survey in 2000. The measurement error distribution changes dramatically from 1980 to 2000, with increasing variance and greater bias further away in time from the interview. Some influence of the measurement error on already published results is found.     

Reliability of global sensitivity indices

Uncertainty and sensitivity analysis is an essential ingredient of model development and applications. For many uncertainty and sensitivity analysis techniques, sensitivity indices are calculated based on a relatively large sample to measure the importance of parameters in their contributions to uncertainties in model outputs. To statistically compare their importance, it is necessary that uncertainty and sensitivity analysis techniques provide standard errors of estimated sensitivity indices. In this paper, a delta method is used to analytically approximate standard errors of estimated sensitivity indices for a popular sensitivity analysis method, the Fourier amplitude sensitivity test (FAST). Standard errors estimated based on the delta method were compared with those estimated based on 20 sample replicates. We found that the delta method can provide a good approximation for the standard errors of both first-order and higher-order sensitivity indices. Finally, based on the standard error approximation, we also proposed a method to determine a minimum sample size to achieve the desired estimation precision for a specified sensitivity index. The standard error estimation method presented in this paper can make the FAST analysis computationally much more efficient for complex models.     

Two Pitfalls of Using Standard Regression Diagnostics When Both X and Y Have Measurement Error

Modern exploratory data analysis produces models that are not based on physical theory but that are consistent with pictures of the data. When both X and Y have error this can be risky, because important features are hidden. Two examples are given that show that systematic model departures and heteroscedasticity may not be detectable with standard regression diagnostics.     

Optimal non-diagonal-type estimators in linear regression under the prediction error sum of squares criterion

This article considers the notion of the non-diagonal-type estimator (NDTE) under the prediction error sum of squares (PRESS) criterion. First, the optimal NDTE in the PRESS sense is derived theoretically and applied to the cosmetics sales data. Second, we make a further study to extend the NDTE to the general case of the covariance matrix of the model and then give a Bayesian explanation for this extension. Third, two remarks concerned with some potential shortcomings of the NDTE are presented and an alternative solution is provided and illustrated by means of simulations.     

The Method to Identify a Biomarker and to Evaluate Its Efficiency for Survival by Using the Joint Model of the Accelerate Failure Time and Longitudinal Data

In this article, we discuss how to identify longitudinal biomarkers in survival analysis under the accelerated failure time model and also discuss the effectiveness of biomarkers under the accelerated failure time model. Two methods proposed by Shcemper et al. are deployed to measure the efficacy of biomarkers. We use simulations to explore how the factors can influence the power of a score test to detect the association of a longitudinal biomarker and the survival time. These factors include the functional form of the random effects from the longitudinal biomarkers, in the different number of individuals, and time points per individual. The simulations are used to explore how the number of individuals, the number of time points per individual influence the effectiveness of the biomarker to predict survival at the given endpoint under the accelerated failure time model. We illustrate our methods using a prothrombin index as a predictor of survival in liver cirrhosis patients.     

HAZARD RATE ESTIMATION FOR CENSORED DATA BY WAVELET METHODS

ABSTRACT

We study the estimation of a hazard rate function based on censored data by non-linear wavelet method. We provide an asymptotic formula for the mean integrated squared error (MISE) of nonlinear wavelet-based hazard rate estimators under randomly censored data. We show this MISE formula, when the underlying hazard rate function and censoring distribution function are only piecewise smooth, has the same expansion as analogous kernel estimators, a feature not available for the kernel estimators. In addition, we establish an asymptotic normality of the nonlinear wavelet estimator.     

Estimation of the mean of the exponential distribution with an outlying observation

Given a life testing experiment consisting of n items, n-1 of which have the expected life λ while one could have an expected life λ/α with 0 < α < 1 the problem is. to find a mean square error (MSE) minimizing estimation function. The standard estimators for the homogeneous case (α = 1) overestimate the expected life and their MSE tend to infinity when a tends to 0.

Looking at the estimation problem as an insurance (see Anscombe (1960)) two different “testimators” are compared with respect to their MSE, Numerical results show that an estimation function based on the “Epstein-statistic” x_(n)/[xbar] is the best one.