A comparison of bootstrap method and Edgeworth expansion in approximating the distribution of sample variance — one sample and two sample cases

The performance of the bootstrap method and the Edgeworth expansion in approximating the distribution of sample variance are compared when the data are from a non-normal population. Both approximations are very good. so long as the parent population is close to normal.     

ASYMPTOTIC THEORY FOR PRINCIPAL COMPONENT ANALYSIS: NON-NORMAL CASE1

This paper considers the asymptotic distributions of latent roots and vectors in principal components analysis when the parent population is non-normal. It is shown that sufficient of T. W. Anderson's asymptotic theory in the multivariate normal case carries over for some results to be obtained.     

Power comparisons of single-sample winsorized T,trimmed T and student's T

The objective of this study is to evaluate the performance of Winsorized t which was proposed to be used as one alternative to the student's t statistic when the parent population is non-normal. The criterion we use is the power of the test. Comparative performance of the Winsorized t, trimmed t, and Studen$tCs t for normal populations and for Studen$tCs t distribution with 8, 4, and 2 degrees of freedom was investigated.     

Estimation of parameters for rounded data from non-normal distributions

This paper is concerned with how standard estimation procedures perform in terms of eficiency for non-normal rounded data. Previous research has shown that the loss in eficiency due to rounding normal data is small. However, evidence from the non-normal distribution considered in this paper suggests, if rounding is coarse or the distribution is very skewed the loss in efficiency due to rounding can be considerable.     

Contributions to two-sample statistics

When testing the equality of the means from two independent normally distributed populations given that the variances of the two populations are unknown but assumed equal, the classical Student's two-sample t-test is recommended. If the underlying population distributions are normal with unequal and unknown variances, either Welch's t-statistic or Satterthwaite's approximate F test is suggested. However, Welch's procedure is non-robust under most non-normal distributions. There is a variable tolerance level around the strict assumptions of data independence, homogeneity of variances, and identical and normal distributions. Few textbooks offer alternatives when one or more of the underlying assumptions are not defensible. While there are more than a few non-parametric (rank) procedures that provide alternatives to Student's t-test, we restrict this review to the promising alternatives to Student's two-sample t-test in non-normal models.     

The robustness of the systemwise breusch-godfrey autocorrelation test for non-normal distributed error terms

Using Monte Carlo methods, the properties of systemwise generalisations of the Breusch-Godfrey test for autocorrelated errors are studied in situations when the error terms follow either normal or non-normal distributions, and when these errors follow either AR(1) or MA(1) processes. Edgerton and Shukur (1999) studied the properties of the test using normally distributed error terms and when these errors follow an AR(1) process. When the errors follow a non-normal distribution, the performances of the tests deteriorate especially when the tails are very heavy. The performances of the tests become better (as in the case when the errors are generated by the normal distribution) when the errors are less heavy tailed.     

How non-normality affects the quadratic discriminant function

The quadratic discriminant function is commonly used for the two group classification problem when the covariance matrices in the two populations are substantially unequal. This procedure is optimal when both populations are multivariate normal with known means and covariance matrices. This study examined the robustness of the QDF to non-normality. Sampling experiments were conducted to estimate expected actual error rates for the QDF when sampling from a variety of non-normal distributions. Results indicated that the QDF was robust to non-normality except when the distributions were highly skewed, in which case relatively large deviations from optimal were observed. In all cases studied the average probabilities of misclassification were relatively stable while the individual population error rates exhibited considerable variability.     

Multiple imputation using multivariate gh transformations

Multiple imputation has emerged as a popular approach to handling data sets with missing values. For incomplete continuous variables, imputations are usually produced using multivariate normal models. However, this approach might be problematic for variables with a strong non-normal shape, as it would generate imputations incoherent with actual distributions and thus lead to incorrect inferences. For non-normal data, we consider a multivariate extension of Tukey's gh distribution/transformation [38] to accommodate skewness and/or kurtosis and capture the correlation among the variables. We propose an algorithm to fit the incomplete data with the model and generate imputations. We apply the method to a national data set for hospital performance on several standard quality measures, which are highly skewed to the left and substantially correlated with each other. We use Monte Carlo studies to assess the performance of the proposed approach. We discuss possible generalizations and give some advices to practitioners on how to handle non-normal incomplete data.     

Robust Two-Sample Statistics for Testing Equality of Means: a Simulation Study

When testing the equality of the means from two independent normally distributed populations given that the variances of the two populations are unknown but assumed equal, the classical two-sample t-test is recommended. If the underlying population distributions are normal with unequal and unknown variances, either Welch's t-statistic or Satterthwaite's Approximate F-test is suggested. However, Welch's procedure is non-robust under most non-normal distributions. There is a variable tolerance level around the strict assumptions of data independence, homogeneity of variances and normality of the distributions. Few textbooks offer alternatives when one or more of the underlying assumptions are not defensible.     

Improved autoregressive forecasts in the presence of non-normal errors

This paper is concerned with obtaining more accurate point forecasts in the presence of non-normal errors. Specifically, we apply the residual augmented least-squares (RALS) estimator to autoregressive models to utilize the additional moment restrictions embodied in non-normal errors. Monte Carlo experiments are performed to compare our RALS forecasts to forecasts based on the ordinary least-squares estimator and the least absolute deviations (LAD) estimator. We find that the RALS approach provides superior forecasts when the data are skewed. Compared to the LAD forecast, the RALS forecast has smaller mean squared prediction errors in the baseline case with normal errors.     

Robust 2k factorial design with Weibull error distributions

It is well known that the least squares method is optimal only if the error distributions are normally distributed. However, in practice, non-normal distributions are more prevalent. If the error terms have a non-normal distribution, then the efficiency of least squares estimates and tests is very low. In this paper, we consider the 2^k factorial design when the distribution of error terms are Weibull W(p,σ). From the methodology of modified likelihood, we develop robust and efficient estimators for the parameters in 2^k factorial design. F statistics based on modified maximum likelihood estimators (MMLE) for testing the main effects and interaction are defined. They are shown to have high powers and better robustness properties as compared to the normal theory solutions. A real data set is analysed.     

Sensitivity of Bayesian sampling inspection schemes to a non-normal prior distribution

summary In this paper we derive the predictive density function of a future observation under the assumption of Edgeworth-type non-normal prior distribution for the unknown mean of a normal population. Fixed size single sample and sequential sampling inspection plans, in a decisive prediction framework, are examined for their sensitivity to departures from normality of the prior distribution. Numerical illustrations indicate that the decision to market the remaining items of a given lot for a fixed size plan may be sensitive to the presence of skewness or kurtosis in the prior distribution. However, Bayes'decision based on the sequential plan may not change though expected gains may change with variation in the non-normality of the prior distribution.     

Acceptance sampling plans by variables for a class of symmetric distributions

The estimation of percentage defectives using a normal sampling plan will not be appropriate when the assumption of normality is violated. In this paper, we propose a sampling plan based on a more general symmetric family of distributions with the parameters estimated using the modified maximum likelihood (MML) procedures introduced by Tiku and Suresh . This sampling plan works well for most of the symmetric non-normal distributions. Some numerical study has also been carried out to show the superiority of the proposed plan.     

A NEW TEST FOR EQUALITY OF VARIANCES FOR k NORMAL POPULATIONS

ABSTRACT

A simple test based on Gini's mean difference is proposed to test the hypothesis of equality of population variances. Using 2000 replicated samples and empirical distributions, we show that the test compares favourably with Bartlett's and Levene's test for the normal population. Also, it is more powerful than Bartlett's and Levene's tests for some alternative hypotheses for some non-normal distributions and more robust than the other two tests for large sample sizes under some alternative hypotheses. We also give an approximate distribution to the test statistic to enable one to calculate the nominal levels and P-values.     

A Comparison of Approaches to Interval Estimation of Mediated Effects in Linear Regression Models

A novel approach based on the concepts of a generalized pivotal quantity (GPQ) is developed to construct confidence intervals for the mediated effect. Thereafter, its performance is compared with six interval estimation approaches in terms of empirical coverage probability and expected length via simulation and two real examples. The results show that the GPQ-based and bootstrap percentile methods outperform other methods when mediated effects exist in small and medium samples. Moreover, the GPQ-based method exhibits a more stable performance in small and non-normal samples. A discussion on how to choose the best interval estimation method for mediated effects is presented.     

REDUCED DISTANCE MEASURES AND TRANSFORMATIONS IN PROCESSING MULTIVARIATE OUTLIERS

The study of multivariate outliers raises many problems of definition, principle and manipulation. Well-authenticated tests of discordancy exist only for the multivariate normal distribution. Detection of outliers in non-normal distributions involves the adoption of appropriate criteria to represent 'extremeness' of observations in a sample; corresponding tests of discordancy usually require tedious, or even intractable, distributional and computational manipulations. A class of transformations of the data is considered with a view of transferring some of the familiar and desirable features of discordancy tests for normal samples to non-normal situations.     

Median rankit control chart by the quantile approach

It is desirable that the data for a statistical control chart be normally distributed. However, if the data are not normal, then a transformation can be used, e.g. Box-Cox transformations, to produce a suitable control chart. In this paper we will discuss a quantile approach to produce a control chart and to estimate median rankit for various non-normal distributions. We will also provide examples of logistic data to indicate how a quantile approach could be used to construct a control chart for a non-normal distribution using a median rankit.     

JACKKNIFE TESTS FOR HETEROSCEDASTICITY IN THE GENERAL LINEAR MODEL

In this paper, tests based on the Jackknife technique are proposed to test for heteroscedasticity in the linear regression model when the errors are non-normal. These are the Jackknifed Goldfeld-Quandt (GQ), and jack-knife related variations of White (H), Lagrange multiplier (LM), Glejser (GL) and Bickel (B) tests. The power of the proposed tests is compared with that of GQ, H, LM, GL and B tests; and the robustness to the error distribution is analyzed under several heteroscedastic assumptions. The GQ test is by far the best test if the error distribution is close to normal, however, GQ test is not robust against non-normal errors. By applying the jackknife technique to the regression a more robust statistic (GQJRG) is produced but the cost is a loss in power. The GQJRG statistic generally is not M powerful as the Bickel (BlOLS) and Glejser (GLlOLS) statistics.     

Bayesian Analysis of Two-Level Fractional Factorial Experiments with Non-Normal Responses

An intractable issue on screening experiments is to identify significant effects and to select the best model when the number of factors is large, especially for fractional factorial experiments with non-normal responses. In such cases, a three-stage Bayesian approach based on generalized linear models (GLMs) is proposed to identify which effects should be included in the target model where the principles of effect sparsity, hierarchy, and heredity are successfully considered. Three simulation experiments with non-normal responses which follow Poisson, binomial, and gamma distributions are presented to illustrate the performance of the proposed approach.     

Distribution theory and inference for polynomial-normal densities

This paper considers a class of densities formed by taking the product of nonnegative polynomials and normal densities. These densities provide a rich class of distributions that can be used in modelling when faced with non-normal characteristics such as skewness and multimodality. In this paper we address inferential and computational issues arising in the practical implementation of this parametric family in the context of the linear model. Exact results are recorded for the conditional analysis of location-scale models and an importance sampling algorithm is developed for the implementation of a conditional analysis for the general linear model when using polynomial-normal distributions for the error.