Population correlation matrices for sampling experiments

A procedure is given for generating correlation matrices which can be used as population correlation matrices for sampling experiments. The algorithm specifies the eigenvalues and randomly selects a correlation matrix from the class of all correlation matrices which possess these same eigenvalues. It is possible to obtain a set of correlation matrices which are indexed by the degree of interdependence among the variables by parameterizing the eigenvalues with a single parameter. An example is the case in which the eigenvalues form a geometric progression. Examples are given and an application to the problem of stopping rules in stepwise regression is discussed. Other applications are also briefly discussed.     

Two tests for monotone failure rate with incomplete data

Two statistics are proposed for testing for the exponential distribution against monotone failure rate alternatives when ran-domly right censored data are available. One of them is a general-ization of the Billmann, Antle and Bain test based on the MLE of the shape parameter of the Weibull distribution. The second has the advantage of being given in closed form. For this test the asymptotic null distribution is given. Consistency of the two tests is proved starting from an expected value inequality characterizing monotone failure rate.     

An empirical study of eight tests of partial correlation coefficients

Partial correlations can be used to statistically control for the effects of unwanted variables.Perhaps the most frequently used test of a partial correlation is the parametric F test,which requires normality of the joint distribution of observations.The possibility that this assumption may not be met in practice suggests a need for procedures that do not require normality.Unfortunately,the statistical literature provides little guidance for choosing other tests when the normalityassumption is not satisfied.Several nonparametric tests of partial correlations are investigated using a computer simulation study.Recommendations are made for selecting certain tests under particular conditions     

Clustering of correlation matrices: A step-wise procedure

An alternative to the criteria proposed by King [9] for clustering correlation matrices is explored. Several simple examples examined by total enumeration employing our criterion indicate that the criterion yields clusters with high within group correlation. A step-wise routine is suggested for problems of realistic size and is applied to the examples presented in King's article. Additional evidence in the form of several Monte Carlo experiments indicates that the routine performs satisfactorily in determining the optimal separation of variables into two groups.     

Generation of random matrices with orthonormal columns and multivariate normal wariates with given sample mean and covariance

In this paper, an algorithm for generating random matrices with orthonormal columns is introduced. As pointed out by a referee, the algorithm is almost identical to Wedderburn's (1975) unpublished method. The method can also be considered as an extension of Stewart's (1980) method, which was designed to generate random orthogonal matrices. It is found outperforming a simple extension of the QR factorization method and that of Heiberger's (1978) method. This paper also demonstrates how the algorithm can be used in generating multivariate normal variates with given sample mean and sample covariance matrix.     

Partial and inverse autocorrelations in portmanteau-type tests for time series

We present a decomposition of the correlation coefficient between x_t and x_t?k into three terms that include the partial and inverse autocorrelations. The first term accounts for the portion of the autocorrelation that is explained by the inner variables {x_t?1 , x_t?2 , …, x _{t? k+1}}, the second one measures the portion explained by the outer variables {x _t+1, x _t+2, } ∪ {x _t?k?1, x _t?k?2,…} and the third term measures the correlation between x _t and x_t?k given all other variables. These terms, squared and summed, can form the basis of three portmanteau-type tests that are able to detect both deviation from white noise and lack of fit of an entertained model. Quantiles of their asymptotic sample distributions are complicated to derive at an adequate level of accuracy, so they are approximated using the Monte Carlo method. A simulation experiment is carried out to investigate significance levels and power of each test, and compare them to the portmanteau test.     

Combined estimators for the correlation coefficient

Three combined estimators for the bivariate normal correlation parameter are considered. The data consist of k independent sample correlation coefficients and it is assumed that the underlying correlation parameters are all equal to ρ. Based upon the joint density function of the sample correlations a combined estimator of ρ is obtained as an approximation to the maximum likelihood solution. Two linearly combined estimators are also considered. One of them is based on Fisher's z-transformation of the sample correlations and the other on an unbiased estimator of ρ. The comparison of these three estimators indicates that the combined (approximate) MLE has a slightly smaller estimated mean squared error relative to the other two combined methods of estimation, but it does so at the expense of a relatively larger bias.     

New results on the convergence of random matrices

This paper extends the previous convergence results in Cerqueti and Costantini (2008) to a more general case using larger normed set of functions. In this regard, the weight-based convergence of the random matrices and their generalized eigenvalues is obtained under less restrictive requirements for the weights.     

Entropy-based tests of uniformity: A Monte Carlo power comparison

In this article, we consider different entropy estimators and propose some entropy-based tests of uniformity. Critical values of the proposed test statistics are obtained by Monte Carlo simulation. Then the power values of the tests for various alternatives and sample sizes are compared. Finally, some recommendations for the application of the proposed tests in practice are presented.     

Monte Carlo comparison of five exponentiality tests using different entropy estimates

The paper studies five entropy tests of exponentiality using five statistics based on different entropy estimates. Critical values for various sample sizes determined by means of Monte Carlo simulations are presented for each of the test statistics. By simulation, we compare the power of these five tests for various alternatives and sample sizes.     

Robust parametric tests of constant conditional correlation in a MGARCH model

This article provides a rigorous asymptotic treatment of new and existing asymptotically valid conditional moment (CM) testing procedures of the constant conditional correlation (CCC) assumption in a multivariate GARCH model. Full and partial quasi maximum likelihood estimation (QMLE) frameworks are considered, as is the robustness of these tests to non-normality. In particular, the asymptotic validity of the LM procedure proposed by Tse (2000 Tse, Y. K. (2000). A test for constant correlations in a multivariate GARCH model. Journal of Econometrics 98 (1):107–127.[Crossref], [Web of Science ®] , [Google Scholar]) is analyzed, and new asymptotically robust versions of this test are proposed for both estimation frameworks. A Monte Carlo study suggests that a robust Tse test procedure exhibits good size and power properties, unlike the original variant which exhibits size distortion under non-normality.     

Power studies of tests for uniformity,II

Results from a power study of six statistics for testing that a sample is from a uniform distribution on the unit interval (0,1) are reported. The test statistics are all well-known and each of them was originally proposed because they should have high power against some alternative distributions. The tests considered are the Pearson probability product test, the Neyman smooth test, the Sukhatme test, the Durbin-Kolmogorov test, the Kuiper test, and the Sherman test. Results are given for each of these tests against each of four classes of alternatives. Also, the most powerful test against each member of the first three alternatives is obtained, and the powers of these tests are given for the same sample sizes as for the six general "omnibus" test statistics. These values constitute a "power envelope" against which all tests can be compared. The Neyman smooth tests with 2nd and 4th degree polynomials are found to have good power and are recommended as general tests for uniformity.     

On population-based simulation for static inference

In this paper we present a review of population-based simulation for static inference problems. Such methods can be described as generating a collection of random variables {X _n } _n=1,…,N in parallel in order to simulate from some target density π (or potentially sequence of target densities). Population-based simulation is important as many challenging sampling problems in applied statistics cannot be dealt with successfully by conventional Markov chain Monte Carlo (MCMC) methods. We summarize population-based MCMC (Geyer, Computing Science and Statistics: The 23rd Symposium on the Interface, pp. 156–163, 1991; Liang and Wong, J. Am. Stat. Assoc. 96, 653–666, 2001) and sequential Monte Carlo samplers (SMC) (Del Moral, Doucet and Jasra, J. Roy. Stat. Soc. Ser. B 68, 411–436, 2006a), providing a comparison of the approaches. We give numerical examples from Bayesian mixture modelling (Richardson and Green, J. Roy. Stat. Soc. Ser. B 59, 731–792, 1997).     

Exploring the Impact of Multivariate GARCH Innovations on Hypothesis Testing for Cointegrating Vectors

This article investigates the impact of multivariate generalized autoregressive conditional heteroskedastic (GARCH) errors on hypothesis testing for cointegrating vectors. The study reviews a cointegrated vector autoregressive model incorporating multivariate GARCH innovations and a regularity condition required for valid asymptotic inferences. Monte Carlo experiments are then conducted on a test statistic for a hypothesis on the cointegrating vectors. The experiments demonstrate that the regularity condition plays a critical role in rendering the hypothesis testing operational. It is also shown that Bartlett-type correction and wild bootstrap are useful in improving the small-sample size and power performance of the test statistic of interest.     

Confidence interval estimation and hypothesis testing for a common correlation coefficient

Several procedures have been proposed for testing 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 is then of interest to estimate and perhaps test hypotheses for the common correlation. In this paper, two versions of five different test statistics are compared via simulation in terms of adequacy of the normal approximation, coverage probabilities of confidence intervals, control of Type I error, and power. The results indicate that two test statistics based on the average of the Fisher z-transforms of the sample correlations should be used in most cases. A statistic based on the sample eigenvalues also gives reasonable results for confidence intervals and lower-tailed tests.     

Bayesian tests for unit root and multiple breaks

A Bayesian approach is considered for identifying sources of nonstationarity for models with a unit root and breaks. Different types of multiple breaks are allowed through crash models, changing growth models, and mixed models. All possible nonstationary models are represented by combinations of zero or nonzero parameters associated with time trends, dummy for breaks, or previous levels, for which Bayesian posterior probabilities are computed. Multiple tests based on Markov chain Monte Carlo procedures are implemented. The proposed method is applied to a real data set, the Korean GDP data set, showing a strong evidence for two breaks rather than the usual unit root or one break.     

On the Value of derivative evaluations and random walk suppression in Markov Chain Monte Carlo algorithms

Two strategies that can potentially improve Markov Chain Monte Carlo algorithms are to use derivative evaluations of the target density, and to suppress random walk behaviour in the chain. The use of one or both of these strategies has been investigated in a few specific applications, but neither is used routinely. We undertake a broader evaluation of these techniques, with a view to assessing their utility for routine use. In addition to comparing different algorithms, we also compare two different ways in which the algorithms can be applied to a multivariate target distribution. Specifically, the univariate version of an algorithm can be applied repeatedly to one-dimensional conditional distributions, or the multivariate version can be applied directly to the target distribution.     

Assessing convergence diagnostic tests for Bayesian Cox regression

The Markov chain Monte Carlo (MCMC) method generates samples from the posterior distribution and uses these samples to approximate expectations of quantities of interest. For the process, researchers have to decide whether the Markov chain has reached the desired posterior distribution. Using convergence diagnostic tests are very important to decide whether the Markov chain has reached the target distribution. Our interest in this study was to compare the performances of convergence diagnostic tests for all parameters of Bayesian Cox regression model with different number of iterations by using a simulation and a real lung cancer dataset.     

Bayesian phylogenetic inference from animal mitochondrial genome arrangements

Summary. The determination of evolutionary relationships is a fundamental problem in evolutionary biology. Genome arrangement data are potentially more informative than deoxyribonucleic acid sequence data for inferring evolutionary relationships between distantly related taxa. We describe a Bayesian framework for phylogenetic inference from mitochondrial genome arrangement data using Markov chain Monte Carlo methods. We apply the method to assess evolutionary relationships between eight animal phyla.