A comparison of estimators of a common correlation coefficient

It is often of interest to test the hypothesis that all off-diagonal elements of the correlation matrix of a multivariate normal distribution are equal. If the hypothesis of equal correlation can be accepted, it then may be of interest to estimate the common correlation coefficient. In this paper, four estimators of the common correlation are compared in terms of bias, variance, mean squared error, adequacy of the normal approximation, and ease of calculation. The average sample correlation is seen to be comparable to the other estimators and is recommended here since it is the easiest to calculate. The estimators are compared using simulation.     

Fixed-width sequential confidence interval for the correlation coefficient of a bivariate normal distribution

A sequential confidence interval of fixed width 2d d > 0, is constructed for the correlation coefficient of a bivariate normal distribution. It is shown that the coverage probability is approximately equal to a preassigned number γ, 0 < γ < as d → 0.     

Confidence interval estimation for the population coefficient of variation using ranked set sampling: a simulation study

In this paper, an evaluation of the performance of several confidence interval estimators of the population coefficient of variation (τ) using ranked set sampling compared to simple random sampling is performed. Two performance measures are used to assess the confidence intervals for τ, namely: width and coverage probabilities. Simulated data were generated from normal, log-normal, skew normal, Gamma, and Weibull distributions with specified population parameters so that the same values of τ are obtained for each distribution, with sample sizes n=15, 20, 25, 50, 100. A real data example representing birth weight of 189 newborns is used for illustration and performance comparison.     

A robust test for testing the correlation coefficient

A test based on Tiku's MML (modified maximum likelihood) estimators is developed for testing that the population correlation coefficient is zero. The test is compared with various other tests and shown to have good Type I error robustness and power for numerous symmetric and skew bivariate populations.     

Confidence intervals for the correlation from a bivariate normal

Improved confidence intervals are given for the correlation coefficient of the bivariate normal distribution. These are based on Cornish–Fisher expansions for the distribution, density and quantiles of the sample correlation.     

Confidence interval estimation for a linear contrast in intraclass correlation coefficients under unequal family sizes for several populations

Confidence intervals [based on F-distribution and (Z) standard normal distribution] for a linear contrast in intraclass correlation coefficients under unequal family sizes for several populations based on several independent multinormal samples have been proposed. It has been found that the confidence interval based on F-distribution consistently and reliably produced better results in terms of shorter average length of the interval than the confidence interval based on standard normal distribution for various combinations of intraclass correlation coefficient values. The coverage probability of the interval based on F-distribution is competitive with the coverage probability of the interval based on standard normal distribution. The interval based on F-distribution can be used for both small sample and large sample situations. An example with real life data has been presented.     

On the distribution of the correlation coefficient when sampling from a mixture of two bivariate normal densities: Robustness and the effect of outliers

The distribution of the sample correlation coefficient is derived when the population is a mixture of two bivariate normal distributions with zero mean but different covariances and mixing proportions 1 - λ and λ respectively; λ will be called the proportion of contamination. The test of ρ = 0 based on Student's t, Fisher's z, arcsine, or Ruben's transformation is shown numerically to be nonrobust when λ, the proportion of contamination, lies between 0.05 and 0.50 and the contaminated population has 9 times the variance of the standard (bivariate normal) population. These tests are also sensitive to the presence of outliers.     

A minimum variance adaptive technique for parameter estimation and hypothesis testing

This report presents numerical results of an approach for parameter estimation and hypothesis testing that does not rely on specific assumptions about the underlying distribution of errors in the measured data. This approach combines robust estimation procedures, the bootstrap method for estimation of parameter uncertainties, permutation techniques for hypothesis testing, and adaptive approaches to estimation in order to obtain the minimum variance estimator or test statistic (within a predefined class) for the data under consideration. The technique produces efficient estimators of central tendency and powerful test statistics, even for small sample sizes. (Portions of this work have been presented in preliminary form (Turkheimer et al., 1996)).     

Parameter expansion for sampling a correlation matrix: an efficient GPX-RPMH algorithm

Sampling the correlation matrix (R) plays an important role in statistical inference for correlated models. There are two main constraints on a correlation matrix: positive definiteness and fixed diagonal elements. These constraints make sampling R difficult. In this paper, an efficient generalized parameter expanded re-parametrization and Metropolis-Hastings (GPX-RPMH) algorithm for sampling a correlation matrix is proposed. Drawing all components of R simultaneously from its full conditional distribution is realized by first drawing a covariance matrix from the derived parameter expanded candidate density (PXCD), and then translating it back to a correlation matrix and accepting it according to a Metropolis-Hastings (M-H) acceptance rate. The mixing rate in the M-H step can be adjusted through a class of tuning parameters embedded in the generalized candidate prior (GCP), which is chosen for R to derive the PXCD. This algorithm is illustrated using multivariate regression (MVR) models and a simulation study shows that the performance of the GPX-RPMH algorithm is more efficient than that of other methods.     

A monte carlo study of the LN statistic for the multivariate nonparametric median and rank sum tests for two populations

Small sample tables are not available for the multisample multivariate rank sum test (MMRST) or the multisample multivariate median test (MMMT) L_N statistic. Consequently, the statistic usually is compared to its asymptotic Chi-square value. To investigate the appropriateness of this procedure a Monte Carlo study is used to measure both significance level and relative power for a variety of multivariate dispersion structures.     

Computational approach test for inference about several correlation coefficients: Equality and common

In this article, we consider inference about the correlation coefficients of several bivariate normal distributions. We first propose computational approach tests for testing the equality of the correlation coefficients. In fact, these approaches are parametric bootstrap tests, and simulation studies show that they perform very satisfactory, and the actual sizes of these tests are better than other existing approaches. We also present a computational approach test and a parametric bootstrap confidence interval for inference about the parameter of common correlation coefficient. At the end, all the approaches are illustrated using two real examples.     

Confidence intervals for a two-parameter exponential distribution: One- and two-sample problems

The problems of interval estimating the mean, quantiles, and survival probability in a two-parameter exponential distribution are addressed. Distribution function of a pivotal quantity whose percentiles can be used to construct confidence limits for the mean and quantiles is derived. A simple approximate method of finding confidence intervals for the difference between two means and for the difference between two location parameters is also proposed. Monte Carlo evaluation studies indicate that the approximate confidence intervals are accurate even for small samples. The methods are illustrated using two examples.     

Effects of nonnormality on the spart for the correlation coefficient: bivariate edgeworth series distribution

The effects of sampling from the bivariate Edgeworth series distribution (BVESD) on the sequential probability ratio test (SPRT) for the correlation coefficient are assessed. The values of the average sample number (ASN) and the operating characteristic (OC) are determined by simulation from such populations. The robustness of the SPRT to this type of nonnormality is demonstrated.     

Robust estimation and testing in one-way ANOVA for Type II censored samples: skew normal error terms

Censoring can be occurred in many statistical analyses in the framework of experimental design. In this study, we estimate the model parameters in one-way ANOVA under Type II censoring. We assume that the distribution of the error terms is Azzalini's skew normal. We use Tiku's modified maximum likelihood (MML) methodology which is a modified version of the well-known maximum likelihood (ML) in the estimation procedure. Unlike ML methodology, MML methodology is non-iterative and gives explicit estimators of the model parameters. We also propose new test statistics based on the proposed estimators. The performances of the proposed estimators and the test statistics based on them are compared with the corresponding normal theory results via Monte Carlo simulation study. A real life data is analysed to show the implementation of the methodology presented in this paper at the end of the study.     

Bayesian rank-based hypothesis testing for the rank sum test,the signed rank test,and Spearman's ρ

Bayesian inference for rank-order problems is frustrated by the absence of an explicit likelihood function. This hurdle can be overcome by assuming a latent normal representation that is consistent with the ordinal information in the data: the observed ranks are conceptualized as an impoverished reflection of an underlying continuous scale, and inference concerns the parameters that govern the latent representation. We apply this generic data-augmentation method to obtain Bayes factors for three popular rank-based tests: the rank sum test, the signed rank test, and Spearman''s ρs.     

Locally optimal test for simultaneously testing the mean and variance of a normal distribution

A generalization of the locally most powerful unbiased (LMPU) test for the single parameter case to the k-parameter case was proposed by SenGupta and Vermeire (1986). In particular we defined a locally most mean power unbiased (LMMPU) test based on the mean curvature of the power hypersurface. Compared to the type C tests of Neyman and Pearson and the type D tests (Isaacson, 1951), LMMPU tests possess better theoretical properties and enjoy ease of construction of critical regions. In this paper we present an interesting example of a two-parameter univariate normal population for which Isaacson (1951, p. 233) was unsuccessful in finding a type D test. For the case of one observation, we prove that no Type D region exists but the LMMPU test is obtained - it is an example of a test with singular Hessian matrix for its power but is nevertheless a strictly locally unbiased (LU) test.     

The relationship between the improvement on the point estimation and the improvement on the interval estimation for the disturbance variance in a linear regression model

In this paper the relationship between the improvement on the point estimation and the improvement on the interval estimation for the disturbance variance in a linear regression model is discussed It is shown that substituting the Stein-type estimatoi for the usual estimatoi in the confidence interval leads to the improvement on the interval estimation The equal-tailed and the shoitest unbiased intervals are dealt with Some appealing relationship is also found in the unbiased case.