Bootstrap method in ranking and slippage problems

In this paper, the bootstrap method of Efron (1979) is given for a ranking and a slippage problem, where the ranking (or slippage) is with respect to the mean of the distributions. The method is also applied to obtain a confidence interval for the largest mean.     

Inferences on the lognormal mean for complete samples

In many engineering problems it is necessary to draw statistical inferences on the mean of a lognormal distribution based on a complete sample of observations. Statistical demonstration of mean time to repair (MTTR) is one example. Although optimum confidence intervals and hypothesis tests for the lognormal mean have been developed, they are difficult to use, requiring extensive tables and/or a computer. In this paper, simplified conservative methods for calculating confidence intervals or hypothesis tests for the lognormal mean are presented. In this paper, “conservative” refers to confidence intervals (hypothesis tests) whose infimum coverage probability (supremum probability of rejecting the null hypothesis taken over parameter values under the null hypothesis) equals the nominal level. The term “conservative” has obvious implications to confidence intervals (they are “wider” in some sense than their optimum or exact counterparts). Applying the term “conservative” to hypothesis tests should not be confusing if it is remembered that this implies that their equivalent confidence intervals are conservative. No implication of optimality is intended for these conservative procedures. It is emphasized that these are direct statistical inference methods for the lognormal mean, as opposed to the already well-known methods for the parameters of the underlying normal distribution. The method currently employed in MIL-STD-471A for statistical demonstration of MTTR is analyzed and compared to the new method in terms of asymptotic relative efficiency. The new methods are also compared to the optimum methods derived by Land (1971, 1973).     

An Easily Implemented Generalised Inverse Gaussian Generator

Using the ‘ratio’ method an easily implemented algorithm is derived for the generalised inverse Gaussian distribution. Computer timings and efficiency calculations show that the procedure is fast over a wide range of distribution parameter values.     

Testing Treatment-Covariate Interaction in a Nonparametric Analysis of Covariance Model

We consider a recently introduced nonparametric model for Analysis of Covariance and derive an asymptotic test for interaction between covariate and treatment. Furthermore, we suggest data depth techniques to obtain joint confidence regions for the covariate effects in this model.

The finite sample behavior of the asymptotic method is evaluated in simulations. Application of the procedures is illustrated using an epileptic seizures and chemotherapy data set.     

Empirical likelihood confidence bands for distribution functions with missing responses

Distribution function estimation plays a significant role of foundation in statistics since the population distribution is always involved in statistical inference and is usually unknown. In this paper, we consider the estimation of the distribution function of a response variable Y with missing responses in the regression problems. It is proved that the augmented inverse probability weighted estimator converges weakly to a zero mean Gaussian process. A augmented inverse probability weighted empirical log-likelihood function is also defined. It is shown that the empirical log-likelihood converges weakly to the square of a Gaussian process with mean zero and variance one. We apply these results to the construction of Gaussian process approximation based confidence bands and empirical likelihood based confidence bands of the distribution function of Y. A simulation is conducted to evaluate the confidence bands.     

On the simple step-stress model for two-parameter exponential distribution

In this paper, we consider the simple step-stress model for a two-parameter exponential distribution, when both the parameters are unknown and the data are Type-II censored. It is assumed that under two different stress levels, the scale parameter only changes but the location parameter remains unchanged. It is observed that the maximum likelihood estimators do not always exist. We obtain the maximum likelihood estimates of the unknown parameters whenever they exist. We provide the exact conditional distributions of the maximum likelihood estimators of the scale parameters. Since the construction of the exact confidence intervals is very difficult from the conditional distributions, we propose to use the observed Fisher Information matrix for this purpose. We have suggested to use the bootstrap method for constructing confidence intervals. Bayes estimates and associated credible intervals are obtained using the importance sampling technique. Extensive simulations are performed to compare the performances of the different confidence and credible intervals in terms of their coverage percentages and average lengths. The performances of the bootstrap confidence intervals are quite satisfactory even for small sample sizes.     

New High-Speed Computing Devices

“One method of error analysis (not the one we will use) is based upon the principles of mathematical statistics. Unfortunately, statistical methods can only be meaningfully applied when one has large amounts of data for a given system. In many cases … these large quantities of data are not available … then statistical methods are not applicable, and some other methods must be devised.”     

Non-Normal Bivariate Distributions with Normal Marginals

The purpose of this note is to indicate that Fieller's Theorem can be expressed in the matrix formulation of the general linear model. The practical consequence is that one general computer program which can estimate the parameters and test the validity of a pertinent model, can also compute confidence limits for the ratios of any linear combinations of the parameters.     

Confidence Intervals for Gini's Diversity Measure and Shannon's Entropy Using Adjusted Proportions

Abstract

Asymptotic confidence intervals are given for two functions of multinomial outcome probabilities: Gini's diversity measure and Shannon's entropy. “Adjusted” proportions are used in all asymptotic mean and variance formulas, along with a possible logarithmic transformation. Exact confidence coefficients are computed in some cases. Monte Carlo simulation is used in other cases to compare actual coverages to nominal ones. Some recommendations are made.     

Criticism of the predictive distribution

The predictive distribution is a mixture of the original distribution model and is used for predicting a future observation. Therein, the mixing distribution is the posterior distribution of the distribution parameters in the Bayesian inference. The mixture can also be computed for the frequentist inference because the Bayesian posterior distribution has the same meaning as a frequentist confidence interval. I present arguments against the concept of predictive distribution. Examples illustrate these. The most important argument is that the predictive distribution can depend on the parameterization. An improvement of the theory of the predictive distribution is recommended.