A Bayesian Meta-Modeling Approach for Gaussian Stochastic Process Models Using a Non Informative Prior

In this article, an efficient Bayesian meta-modeling approach is proposed for Gaussian stochastic process models in computer experiments. Different prior densities and particularly, a non informative hyper prior have been employed on the parameters involved in the correlation matrix. And the estimation of related parameters is obtained by the expectation-maximization algorithm. Compared with the recent work of Li and Sudjianto (2005 Li , R. , Sudjianto , A. ( 2005 ). Analysis of computer experiments using penalized likelihood in Kriging models . Technometrics 47 : 111 – 120 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), the proposed approach is not only of higher prediction accuracy but also of lower computational cost, due to the utilization of the non informative prior and the absence of tuning parameters. Experimental results demonstrate that our approach yields state-of-the-art performance.     

The Scale Invariant Wigner Spectrum Estimation of Gaussian Locally Self-Similar Processes

We study locally self-similar processes (LSSPs) in Silverman’s sense. By deriving the minimum mean-square optimal kernel within Cohen’s class counterpart of time–frequency representations, we obtain an optimal estimation for the scale invariant Wigner spectrum (SIWS) of Gaussian LSSPs. The class of estimators is completely characterized in terms of kernels, so the optimal kernel minimizes the mean-square error of the estimation. We obtain the SIWS estimation for two cases: global and local, where in the local case, the kernel is allowed to vary with time and frequency. We also introduce two generalizations of LSSPs: the locally self-similar chirp process and the multicomponent LSSP, and obtain their optimal kernels. Finally, the performance and accuracy of the estimation is studied via simulation.     

Obtaining a Trapezoidal Distribution

Given a most believed value for a quantity together with upper and lower possible deviations from that value, a rectangular distribution might be used to represent state-of-knowledge about the quantity. If the deviations are themselves known by probability distributions, and the value conditioned on the deviations is rectangular, then the marginal distribution of the value is determined by the distributions of the deviations. Here we show under quite general conditions that conversely, given the marginal distribution, the distributions of the deviations are uniquely determined. The case in which the marginal distribution is trapezoidal is studied in some detail.     

Bimodal skew-symmetric normal distribution

ABSTRACT

We introduce a new parsimonious bimodal distribution, referred to as the bimodal skew-symmetric Normal (BSSN) distribution, which is potentially effective in capturing bimodality, excess kurtosis, and skewness. Explicit expressions for the moment-generating function, mean, variance, skewness, and excess kurtosis were derived. The shape properties of the proposed distribution were investigated in regard to skewness, kurtosis, and bimodality. Maximum likelihood estimation was considered and an expression for the observed information matrix was provided. Illustrative examples using medical and financial data as well as simulated data from a mixture of normal distributions were worked.     

Measures of uniformity in experimental designs: A selective overview

ABSTRACT

Discrepancies are measures which are defined as the deviation between the empirical and the theoretical uniform distribution. In this way, discrepancy is a measure of uniformity which provides a way of construction a special kind of space filling designs, namely uniform designs. Several discrepancies have been proposed in recent literature. A brief, selective review of these measures including some construction algorithms are given in this paper. Furthermore, a critical discussion along with some comparisons is provided, as well.     

Parameter estimation by minimizing a probability generating function-based power divergence

Abstract

Generating function-based statistical inference is an attractive approach if the probability (density) function is complicated when compared with the generating function. Here, we propose a parameter estimation method that minimizes a probability generating function (pgf)-based power divergence with a tuning parameter to mitigate the impact of data contamination. The proposed estimator is linked to the M-estimators and hence possesses the properties of consistency and asymptotic normality. In terms of parameter biases and mean squared errors from simulations, the proposed estimation method performs better for smaller value of the tuning parameter as data contamination percentage increases.     

A simulation study of the nearest neighbour analysis method of papadakis

The relative performances of randomised block, balanced lattice squares and Papadakis nearest neighbour analyses were compared on two simulated fields whose soil heterogeneity profiles were generated, one with a few evenly spaced contours and the other with many unevenly spaced contours. Four levels of random error were generated to simulate different proportions of random error and soil heterogeneity. Dummy treatments, corresponding to 7 x 7 and 11 x 11 balanced lattice squares were applied to the fields. The results from simulated experiments showed an interaction of error mean square (EMS) between size of experiment (7 x 7, 11 x 11) and levels of soil heterogeneity in the lattice analyses, but no such interaction in the Papadakis analyses. The Papadakis EMS decreased as random error decreased but at a rate depending on the map andthe ratio of soil heterogeneity to random error.     

Inferences on regression coefficients in a regression model under heteroscedasticity and robustness with respect to departure from normality

In this paper we consider a simple linear regression model under heteroscedasticity and nonnormality. A statistical test for testing the regression coefficient is then derived by assuming normality for the random disturbances and by applying Welch's method. Some Monte Carlo studies are generated for assessing robustness of this test. By combining Tiku's robust procedure with the new test, a robust but more powerful test is developed.     

A new test for the extreme value distribution

In this paper a test statistic which is a modification of the W statistic for testing the goodness of fit for the two paremeter extreme value (smallest element) distribution is proposed. The test statistic Is obtained as the ratio of two linear estimates of the scale parameter. It Is shown that the suggested statistic is computationally simple and has good power properties. Percentage points of the statistic are obtained by performing Monte Carlo experiments. An example is given to illustrate the test procedure.     

On sequential closed testing dose groups with a control

Methods for a sequential test of a dose-response effect in pre-clinical studies are investigated. The objective of the test procedure is to compare several dose groups with a zero-dose control. The sequential testing is conducted within a closed family of one-sided tests. The procedures investigated are based on a monotonicity assumption. These closed procedures strongly control the familywise error rate while providing information about the shape of the dose-responce relationship. Performance of sequential testing procedures are compared via a Monte Carlo simulation study. We illustrae the procedures by application to a real data set.