THE OPTIMAL CLASSIFICATION RULE FOR EXPONENTIAL POPULATIONS

Let π₁ and π₂ be two exponential populations with location parameters α and β respectively. It is shown that the probability of misclassification for the optimal rule, R_E, depends on α and β only through their ratio α/β=r, and that the smaller the ratio r (i.e. the smaller α is compared to β), the greater the superiority of the optimal rule R_E over the commonly used rule R_N. When α and β are unknown, the sample-based version R_E(s) of R_E exhibits the same pattern of superiority over the sample-based R_N(s) of R_N.     

SELECTING THE t BEST POPULATIONS — SCALE PARAMETER CASE

Selection from k independent populations of the t (< k) populations with the smallest scale parameters has been considered under the Indifference Zone approach by Bechhofer k Sobel (1954). The same problem has been considered under the Subset Selection approach by Gupta & Sobel (1962a) for the normal variances case and by Carroll, Gupta & Huang (1975) for the more general case of stochastically increasing distributions. This paper uses the Subset Selection approach to place confidence bounds on the probability of selecting all “good” populations, or only “good” populations, for the Case of scale parameters, where a “good” population is defined to have one of the t smallest scale parameters. This is an extension of the location parameter results obtained by Bofinger & Mengersen (1986). Special results are obtained for the case of selecting normal populations based on variances and the necessary tables are presented.     

COMPARISONS AND SELECTION OF TWO-PARAMETER EXPONENTIAL POPULATIONS

For two-parameter exponential populations with the same scale parameter (known or unknown) comparisons are made between the location parameters. This is done by constructing confidence intervals, which can then be used for selection procedures. Comparisons are made with a control, and with the (unknown) “best” or “worst” population. Emphasis is laid on finding approximations to the confidence so that calculations are simple and tables are not necessary. (Since we consider unequal sample sizes, tables for exact values would need to be extensive.)     

The Coverage Probability of Confidence Intervals in One‐Way Analysis of Covariance after Two F Tests

Volume 3 of Analysis of Messy Data by Milliken & Johnson (2002) provides detailed recommendations about sequential model development for the analysis of covariance. In his review of this volume, Koehler (2002) asks whether users should be concerned about the effect of this sequential model development on the coverage probabilities of confidence intervals for comparing treatments. We present a general methodology for the examination of these coverage probabilities in the context of the two‐stage model selection procedure that uses two F tests and is proposed in Chapter 2 of Milliken & Johnson (2002). We apply this methodology to an illustrative example from this volume and show that these coverage probabilities are typically very far below nominal. Our conclusion is that users should be very concerned about the coverage probabilities of confidence intervals for comparing treatments constructed after this two‐stage model selection procedure.     

BARTLETT-TYPE CORRECTIONS FOR TWO-PARAMETER EXPONENTIAL FAMILY MODELS

ABSTRACT

Cordeiro and Ferrari[1] Cordeiro, G.M. and Ferrari, S.L.P. 1991. A Modified Score Test Statistic Having Chi-Squared Distribution to Order n^?1. Biometrika, 78: 573–582. [Web of Science ®] , [Google Scholar] obtained a Bartlett-type correction to the score statistic that is given by a polynomial of second degree in the statistic itself with coefficients that depend on cumulants of log–likelihood derivatives. Although the corrected statistic has good size properties, it is not always a monotone transformation of the original statistic. Recently, some monotone transformations of the score statistic have been proposed as alternatives to the polynomial transformation. In this paper we derive simple formulae for various modified score statistics for testing a scalar parameter in two-parameter exponential models which do not require knowledge of cumulants. The formulae are readily applicable to cover many important and commonly used distributions and involve only trivial operations on certain functions and their derivatives.     

CHANCE OF MISCLASSIFICATION IN THE CASE OF SEVERAL NORMAL POPULATIONS

When a new observation is to be classified into one of several multivariate normal populations with different means and the same covariance matrix, by Rao's method of scoring, the chance of misclassification is expressed as a multiple integral. This paper gives a practical method of obtaining reasonable approximations to this integral by using tables prepared by Gibbons, Olkin & Sobel (1977) for a different task.     

A NOTE ON SAMPLE QUANTILE PROCESSES OF FINITE POPULATIONS

A representation is established for sample quantile processes of simple random samples drawn without replacement from finite populations. The result is then used to obtain weak limits of properly normalized quantile processes and their linear functions.     

BAYES EMPIRICAL BAYES APPROACH TO ESTIMATION OF THE FAILURE RATE IN EXPONENTIAL DISTRIBUTION

ABSTRACT

In the empirical Bayes (EB) decision problem consisting of squared error estimation of the failure rate in exponential distribution, a prior Λ is placed on the gamma family of prior distributions to produce Bayes EB estimators which are admissible. A subclass of such estimators is shown to be asymptotically optimal (a.o.). The results of a Monte Carlo study are presented to demonstrate the a.o. property of the Bayes EB estimators.     

A NEW TEST FOR EQUALITY OF VARIANCES FOR k NORMAL POPULATIONS

ABSTRACT

A simple test based on Gini's mean difference is proposed to test the hypothesis of equality of population variances. Using 2000 replicated samples and empirical distributions, we show that the test compares favourably with Bartlett's and Levene's test for the normal population. Also, it is more powerful than Bartlett's and Levene's tests for some alternative hypotheses for some non-normal distributions and more robust than the other two tests for large sample sizes under some alternative hypotheses. We also give an approximate distribution to the test statistic to enable one to calculate the nominal levels and P-values.     

Estimating Average Worth of the Selected Subset from Two-Parameter Exponential Populations

ABSTRACT

Suppose independent random samples are available from k(k ≥ 2) exponential populations ∏₁,…,∏ _k with a common location θ and scale parameters σ₁,…,σ _k , respectively. Let X _i and Y _i denote the minimum and the mean, respectively, of the ith sample, and further let X = min{X ₁,…, X _k } and T _i  = Y _i  ? X; i = 1,…, k. For selecting a nonempty subset of {∏₁,…,∏ _k } containing the best population (the one associated with max{σ₁,…,σ _k }), we use the decision rule which selects ∏ _i if T _i  ≥ c max{T ₁,…,T _k }, i = 1,…, k. Here 0 < c ≤ 1 is chosen so that the probability of including the best population in the selected subset is at least P* (1/k ≤ P* < 1), a pre-assigned level. The problem is to estimate the average worth W of the selected subset, the arithmetic average of means of selected populations. In this article, we derive the uniformly minimum variance unbiased estimator (UMVUE) of W. The bias and risk function of the UMVUE are compared numerically with those of analogs of the best affine equivariant estimator (BAEE) and the maximum likelihood estimator (MLE).     

PREDICTION FOR EXPONENTIAL LIFETIMES BASED ON STEP-STRESS TESTING

ABSTRACT

This paper presents methods for constructing prediction limits for a step-stress model in accelerated life testing. An exponential life distribution with a mean that is a log-linear function of stress, and a cumulative exposure model are assumed. Two prediction problems are discussed. One concerns the prediction of the life at a design stress, and the other concerns the prediction of a future life during the step-stress testing. Both predictions require the knowledge of some model parameters. When estimates for the model parameters are available, a calibration method based on simulations is proposed for correcting the prediction intervals (regions) obtained by treating the parameter estimates as the true parameter values. Finally, a numerical example is given to illustrate the prediction procedure.     

The two-sample t test versus satterthwaite's approximate f test

A comparison between the two-sample t test and Satterthwaite's approximate F test is made, assuming the choice between these two tests is based on a preliminary test on the variances. Exact formulas for the sizes and powers of the tests are derived. Sizes and powers are then calculated and compared for several situations.     

PMSE performance of two different types of preliminary test estimators under a multivariate t error term

Abstract

In this paper, assuming that the error terms follow a multivariate t distribution, we derive the exact formula for the predictive mean squared error (PMSE) of two different types of pretest estimators. It is shown analytically that one of the pretest estimator dominates the SR estimator if a critical value of the pretest is chosen appropriately. Also, we compare the PMSE of the pretest estimators with the MMSE, AMMSE, SR and PSR estimators by numerical evaluations. Our results show that the pretest estimators dominate the OLS estimator for all combinations when the degrees of freedom is not more than 5.     

THE USE OF AGGREGATE DATA TO ESTIMATE GOMPERTZ‐TYPE OLD‐AGE MORTALITY IN HETEROGENEOUS POPULATIONS

We consider two related aspects of the study of old‐age mortality. One is the estimation of a parameterized hazard function from grouped data, and the other is its possible deceleration at extreme old age owing to heterogeneity described by a mixture of distinct sub‐populations. The first is treated by half of a logistic transform, which is known to be free of discretization bias at older ages, and also preserves the increasing slope of the log hazard in the Gompertz case. It is assumed that data are available in the form published by official statistical agencies, that is, as aggregated frequencies in discrete time. Local polynomial modelling and weighted least squares are applied to cause‐of‐death mortality counts. The second, related, problem is to discover what conditions are necessary for population mortality to exhibit deceleration for a mixture of Gompertz sub‐populations. The general problem remains open but, in the case of three groups, we demonstrate that heterogeneity may be such that it is possible for a population to show decelerating mortality and then return to a Gompertz‐like increase at a later age. This implies that there are situations, depending on the extent of heterogeneity, in which there is at least one age interval in which the hazard function decreases before increasing again.     

A Note on the Arithmetic—Geometric—Harmonic Mean Inequalities

The familiar inequalities relating the arithmetic, geometric, and harmonic means are derived as corollaries to likelihood ratio tests.     

ON SELECTING VARIABLES AND ASSESSING THEIR PERFORMANCE IN LINEAR DISCRIMINANT ANALYSIS

Linear discriminant analysis between two populations is considered in this paper. Error rate is reviewed as a criterion for selection of variables, and a stepwise procedure is outlined that selects variables on the basis of empirical estimates of error. Problems with assessment of the selected variables are highlighted. A leave-one-out method is proposed for estimating the true error rate of the selected variables, or alternatively of the selection procedure itself. Monte Carlo simulations, of multivariate binary as well as multivariate normal data, demonstrate the feasibility of the proposed method and indicate its much greater accuracy relative to that of other available methods.     

Confidence intervals centred on bootstrap smoothed estimators

Bootstrap smoothed (bagged) parameter estimators have been proposed as an improvement on estimators found after preliminary data‐based model selection. A result of Efron in 2014 is a very convenient and widely applicable formula for a delta method approximation to the standard deviation of the bootstrap smoothed estimator. This approximation provides an easily computed guide to the accuracy of this estimator. In addition, Efron considered a confidence interval centred on the bootstrap smoothed estimator, with width proportional to the estimate of this approximation to the standard deviation. We evaluate this confidence interval in the scenario of two nested linear regression models, the full model and a simpler model, and a preliminary test of the null hypothesis that the simpler model is correct. We derive computationally convenient expressions for the ideal bootstrap smoothed estimator and the coverage probability and expected length of this confidence interval. In terms of coverage probability, this confidence interval outperforms the post‐model‐selection confidence interval with the same nominal coverage and based on the same preliminary test. We also compare the performance of the confidence interval centred on the bootstrap smoothed estimator, in terms of expected length, to the usual confidence interval, with the same minimum coverage probability, based on the full model.