Volume estimation by monte carlo methods *

In estimating a multiple integral, it is known that Monte Carlo methods are more efficient than analytical techniques when the number of dimensions is beyond seven. In general, the sample-mean method is better than the hit-or-miss Monte Carlo method. However, when the volume of a domain in a high-dimensional space is of interest, the hit-or-miss method is usually preferred. It is because of the difficulty in generalizing the sample-mean method for the computation of the volume of a domain. This paper develops a technique to make such a generalization possible. The technique can be interpreted as a volume-preserving transformation procedure. A volume-preserving transformation is first performed to transform the concerned domain into a hypersphere. The volume of the domain is then evaluated by computing the volume of the hypersphere.     

A geometric approach to transdimensional markov chain monte carlo

The authors present theoretical results that show how one can simulate a mixture distribution whose components live in subspaces of different dimension by reformulating the problem in such a way that observations may be drawn from an auxiliary continuous distribution on the largest subspace and then transformed in an appropriate fashion. Motivated by the importance of enlarging the set of available Markov chain Monte Carlo (MCMC) techniques, the authors show how their results can be fruitfully employed in problems such as model selection (or averaging) of nested models, or regeneration of Markov chains for evaluating standard deviations of estimated expectations derived from MCMC simulations.     

A monte carlo study on the robustness of four MANOVA criterion tests

Four MANOVA tests (Wilk's Lambda, Roy's Largest Root Test, the Hotelling-Lawley Trace and the Pillai-Bartlett Trace) were studied when restricted sample data were drawn from normal populations. Robustness was compared by examining bias at critical points and fluctuations in the standard error of the empirical distributions. The Wilk's Lambda statistic was found to be the least affected by the restricted sampling.     

A monte carlo study of robust estimators of location

Andrews et al (1972) carried out an extensive Monte Carlo study of robust estimators of location. Their conclusions were that the hampel and the skipped estimates, as classes, seemed to be preferable to some of the other currently fashionable estimators. The present study extends this work to include estimators not previously examined. The estimators are compared over short-tailed as well as long-tailed alternatives and also over some dependent data generated by first-order autoregressive schemes. The conclusions of the present study are threefold. First, from our limited study, none of the so-called robust estimators are very efficient over short-tailed situations. More work seems to be necessary in this situation. Second, none of the estimators perform very well in dependent data situations, particularly when the correlation is large and positive. This seems to be a rather pressing problem. Finally, for long-tailed alternatives, the hampel estimators and Hogg-type adaptive versions of the hampels are the strongest classes. The adaptive hampels neither uniformly outperform nor are they outperformed by the hampels. However, the superiority in terms of maximum relative efficiency goes to the adaptive hampels. That is, the adaptive hampels, under their worst performance.     

Exact monte carlo tests using several statistics

Concepts of ranking and boundary of multivariate statistics are discussed and applied to the simultaneous use of several test statistics calculated for data and simulated replicates. An example of residual analysis in regression is given using layer ranks and supplementary simulation with a stopping rule.     

A monte carlo study of collinearity in linear simultaneous equation models ∗

Although the collinearity issue has been studied in previous simulation studies with a simultaneous system of equations, alternative estimators to circumvent this problem have received little attention. Monte Carlo techniques are used to examine the performance of several estimators under a squared error loss criterion. In particular, this study considers the Vinod–Ullah ridge-type estimators at the first and/or second stage of 2SLS. Ridge regression in the second stage only of 2SLS but not the first stage only, seems to be a practical alternative to 2SLS, especially in situations of strong collinearity. The OLS estimator and the ordinary ridge regression estimator also yield favorable results in situations of moderate to strong collinearity.     

A monte carlo analysis of two spectral tests of the martingale hypothesis

Summary The size, power, and robustness properties of the Kolmogorov-Smirnov and Cramér-von Mises spectral tests of the martingale (difference) hypothesis are investigated by Monte Carlo methods. The results highlight a marked superiority of the Cramér-von Mises with respect to the Kolmogorov-Smirnov test. The paper also shows that the Cramér-von Mises test is simple to compute, more general and more powerful than other converntionally used tests.     

A monte carlo evaluation of response surface analysis based on paired-comparison data

The use of the logit transformation on paired-comparison data in the weighted least squares analysis of response surfaces for aesthetic qualities of products is discussed. Monte Carlo simulations are employed to investigate the small sample properties of the estimators and test statistics. A secondary objective of the Monte Carlo simulations is the comparison of two transformation procedures. The simulations are of standard-item paired-compar-ison experiments in which ties are not allowed.     

Monte carlo sampling approach to testing nonnested hypothesis: monte carlo results

Alternative ways of using Monte Carlo methods to implement a Cox-type test for separate families of hypotheses are considered. Monte Carlo experiments are designed to compare the finite sample performances of Pesaran and Pesaran's test, a RESET test, and two Monte Carlo hypothesis test procedures. One of the Monte Carlo tests is based on the distribution of the log-likelihood ratio and the other is based on an asymptotically pivotal statistic. The Monte Carlo results provide strong evidence that the size of the Pesaran and Pesaran test is generally incorrect, except for very large sample sizes. The RESET test has lower power than the other tests. The two Monte Carlo tests perform equally well for all sample sizes and are both clearly preferred to the Pesaran and Pesaran test, even in large samples. Since the Monte Carlo test based on the log-likelihood ratio is the simplest to calculate, we recommend using it.     

Markov-chain monte carlo: Some practical implications of theoretical results

We review and discuss some recent progress in the theory of Markov-chain Monte Carlo applications, particularly oriented to applications in statistics. We attempt to assess the relevance of this theory for practical applications.     

A monte carlo study of robust and least squares response surface methods

Response surface methodology is useful for exploring a response over a region of factor space and in searching for extrema. Its generality, makes it applicable to a variety of areas. Classical response surface methodology for a continuous response variable is generally based on least squares fitting. The sensitivity of least squares to outlying observations carries over to the surface procedures. To overcome this sensitivity, we propose response surface methodology based on robust procedures for continuous response variables. This robust methodology is analogous to the methodology based on least squares, while being much less sensitive to outlying observations. The results of a Monte Carlo study comparing it and classical surface methodologies for normal and contaminated normal errors are presented. The results show that as the proportion of contamination increases, the robust methodology correctly identifies a higher proportion of extrema than the least squares methods and that the robust estimates of extrema tend to be closer to the true extrema than the least squares methods.     

A monte carlo investigation of a statistic for a bivariate missing data problem

Testing the equal means hypothesis of a bivariate normal distribution with homoscedastic varlates when the data are incomplete is considered. If the correlational parameter, ρ, is known, the well-known theory of the general linear model is easily employed to construct the likelihood ratio test for the two sided alternative. A statistic, T, for the case of ρ unknown is proposed by direct analogy to the likelihood ratio statistic when ρ is known. The null and nonnull distribution of T is investigated by Monte Carlo techniques. It is concluded that T may be compared to the conventional t distribution for testing the null hypothesis and that this procedure results in a substantial increase in power-efficiency over the procedure based on the paired t test which ignores the incomplete data. A Monte Carlo comparison to two statistics proposed by Lin and Stivers (1974) suggests that the test based on T is more conservative than either of their statistics.     

A monte carlo comparison of regression trimmed means

Two proposals are made for constructing adaptive estimators of the parameters in a linear regression model. These estimators are based on regression trimmed means and use an idea of Jaeckel [(1971) Ann Math Statist 42, 1540-1552] and the bootstrap respectively. These adaptive trimmed means as well as some nonadaptive trimmed means are studied by Monte Carlo. A one-step biweight is also included for comparison purposes.     

The relevance of large sample properties of estimators for the errors-in-variables model:a monte carlo study

By means of a Monte Carlo study it is investigated whether moments of the asymptotic distributions of two estimators for the errors-in-variables model are appropriate for employment in small-sample applications.     

Fatal shortcomings of stationary AR (1) exogeneous variables in econometric monte carlo tudies of estimators

The paper first shows that the stationary normal AR(1) process (SNAR1), the most frequently used process for generating exogenous variables in econometric Monte Carlo studies, cannot generate realistic exogenous variables, which are generally trended and similar to those generated by ARIMA (p,d,q) process withd≧1 and positive drift (trend). Then, it illustrates that in the context of AR(1) disturbances,trends in exogenous variables can frequently alter the very ranking of two competing estimators, the ordinary least squares estimator (OLS) and the Cochrane-Orcutt estimators (CO). For three common econometric models—a standard regression model, a dynamic model (i.e., a model with a lagged dependent variable), and a seemingly unrelated regression model, OLS becomes superior in many cases. This is so in spite of the fact that the CO estimator in the study utilizes the true value of the first-order autocorrelation coefficient of the disturbances. The message to be derived from these findings should be ccear. If one accepts the fact that most if not all economic time series are trended, and endorses a proposition that the fundamental if not sole purpose of Monte Carlo studies in econometrics should be to provide useful guidelines to practicing econometricians, then, he must not employ SNARl (nor anyother artificially created nontrended series) as a generator of exogenous variables in a Monte Carlo study, at least in the econometrics of autocorrelated disturbances. Alternative methods of generating stochastic exogenous variables that are trended are suggested in the paper. For almost four decades, the principle of the autoregressive transformation of a regression model with first-order autocorrelated disturbances (the Coestimation priciple) has been taken for granted as a method of correcting for the autocorrelation in the disturbances—be it in the two-stage Cochrane—Orcutt estimator, the iterative Cochrane-Orcutt estimator, or an estimator utilizing nonlinear techniques or search procedures. (Comitting the first observation due to transformation is not considered very crucial in general.) The results of the pertinent Monte Carlo studies appear to justify such a procedure only because most studies have employed SNARl exogenous variables, not trended ones. Thus, Monte Carlo experimenters must be blamed, at least partially, for this prevailining malpractice. It is hoped that they will not commit additional sins by not using realistic data in their future experiments.     

Tests based on weighted rankings in complete blocks: exact distribution and monte carlo simulation

Quade (1972, 1979) proposed a family of nonparametric tests based on weighted within-block rankings, for testing the hypothesis of no treatment effects in a complete randomized blocks layout. In this paper we give a table of the exact null distribution of these tests when the number of treatments is 3, the number of blocks is less than or equal to 14 and the block scores are linear. Moreover, a Monte Carlo study was performed to compare the powers of these tests with parametric and nonparametric competitors     

A monte carlo study of the wald,lm, and lr tests in a heteroscedastic linear model

In this paper, we examine by Monte Carlo experiments the small sample properties of the W (Wald), LM (Lagrange Multiplier) and LR (Likelihood Ratio) tests for equality between sets of coefficients in two linear regressions under heteroscedasticity. The small sample properties of the size-corrected W, LM and LR tests proposed by Rothenberg (1984) are also examined and it is shown that the performances of the size-corrected W and LM tests are very good. Further, we examine the two-stage test which consists of a test for homoscedasticity followed by the Chow (1960) test if homoscedasticity is indicated or one of the W, LM or LR tests if heteroscedasticity should be assumed. It is shown that the pretest does not reduce much the bias in the size when the sizecorrected citical values are used in the W, LM and LR tests.     

Small sample size comparisons of tests for homogeneity of variances by Monte-Carlo

A number of tests are available for testing the equality of several population variances. Some are claimed to be robust. We compared six of those claimed robust procedures by Monte Carlo simulated experiments, particularly for cases of small and unequal sample sizes. Our results show that the jack-knife test compares favorably with the other tests.     

The relative performances of improved ridge estimators and an empirical bayes estimator: some monte carlo results

The relative 'performances of improved ridge estimators and an empirical Bayes estimator are studied by means of Monte Carlo simulations. The empirical Bayes method is seen to perform consistently better in terms of smaller MSE and more accurate empirical coverage than any of the estimators considered here. A bootstrap method is proposed to obtain more reliable estimates of the MSE of ridge esimators. Some theorems on the bootstrap for the ridge estimators are also given and they are used to provide an analytical understanding of the proposed bootstrap procedure. Empirical coverages of the ridge estimators based on the proposed procedure are generally closer to the nominal coverage when compared to their earlier counterparts. In general, except for a few cases, these coverages are still less accurate than the empirical coverages of the empirical Bayes estimator.