Bayes estimators of the reliability of the logistic distribution

Bayes estimators of the reliability function of the logistic distribution are obtained using the methods of Lindley (1980) and Tierney & Kadane (1986). Squared-error and log-odds squared-error loss functions are used. A numerical example is presented. Comparisons are made between these two procedures, based on a Monte Carlo simulation study.     

Improved Empirical Bayes Estimation in Group Testing Procedure for Small Proportions

Group testing has been long recognized as an efficient method to classify all the experimental units into two mutually exclusive categories: defective or not defective. In recent years, more attention has been brought to the estimation of the population prevalence rate p of a disease, or of some property, using group testing. In this article, we propose two scaled squared-error loss functions, which improve the Bayesian approach to estimating p in terms of minimizing the mean squared error (MSE) of the Bayes estimators of p for small p. We show that the new estimators are preferred over the estimator from the usual squared-error loss function and the maximum likelihood estimator (MLE) for small p.     

Estimation of scale parameters in mixture distributions

Simultaneous estimation of scale parameters is considered in mixture distributions under squared-error loss. A general class of estimators is obtained which dominates the componentwise best multiple estimators and the moment estimators. As special cases, improved estimators are obtained for the multivariate t-distribution and the p-variate Lomax distribution.     

Equivariant estimation of a normal mean using a normal concomitant variable for covariance adjustment

Equivariant point estimators of one component of a bivariate normal mean vector are considered when the second component is known. Equivariant point estimators are characterized and compared in terms of their risk functions with respect to a normalized squared-error loss function. Specific point estimators that dominate the usual estimator when the squared correlation coefficient is sufficiently large are provided.     

Minimax estimation of a lower-bounded scale parameter of a gamma distribution for scale-invariant squared-error loss

Let X have a gamma distribution with known shape parameter θr;aL and unknown scale parameter θ. Suppose it is known that θ ≥ a for some known a > 0. An admissible minimax estimator for scale-invariant squared-error loss is presented. This estimator is the pointwise limit of a sequence of Bayes estimators. Further, the class of truncated linear estimators C = {θ_ρ|θ_ρ(x) = max(a, ρ), ρ > 0} is studied. It is shown that each θ_ρ is inadmissible and that exactly one of them is minimax. Finally, it is shown that Katz's [Ann. Math. Statist., 32, 136–142 (1961)] estimator of θ is not minimax for our loss function. Some further properties of and comparisons among these estimators are also presented.     

On the admissibility of the maximum-likelihood estimator of the binomial variance

This paper addresses the admissibility of the maximum-likelihood estimator (MLE) of the variance of a binomial distribution with parameters n and p under squared-error loss. We show that the MLE is admissible for n ≤ 5 and inadmissible for n≥ 6.     

Bayesian inference based on a jointly type-II censored sample from two exponential populations

In this paper, based on a jointly type-II censored sample from two exponential populations, the Bayesian inference for the two unknown parameters are developed with the use of squared-error, linear-exponential and general entropy loss functions. The problem of predicting the future failure times, both point and interval prediction, based on the observed joint type-II censored data, is also addressed from a Bayesian viewpoint. A Monte Carlo simulation study is conducted to compare the Bayesian estimators with the maximum likelihood estimator developed by Balakrishnan and Rasouli [Exact likelihood inference for two exponential populations under joint type-II censoring. Comput Stat Data Anal. 2008;52:2725–2738]. Finally, a numerical example is utilized for the purpose of illustration.     

Non-parametric Estimation of the Death Rate in Branching Diffusions

We consider finite systems of diffusing particles in with branching and immigration. Branching of particles occurs at position dependent rate. Under ergodicity assumptions, we estimate the position-dependent branching rate based on the observation of the particle process over a time interval [0, t ]. Asymptotics are taken as t  → ∞. We introduce a kernel-type procedure and discuss its asymptotic properties with the help of the local time for the particle configuration. We compute the minimax rate of convergence in squared-error loss over a range of Hölder classes and show that our estimator is asymptotically optimal.     

Comments on “a generalization of fisher's exact test in pxq contingency tables using more concordant relations”

The purpose of this note is to criticize Nguyen (1985) for his account of the literature on the generalization of Fisher's exact test and to point out parallels with existing algorithms of the algorithm proposed by Nguyen. Subsequently we will briefly raise some questions on the methodology proposed by Nguyen.

Nguyen (1985) suggests that all literature on exact testing prior to Nguyen & Sampson (1985) is based on the “more probable” relation or Exact Probability Test (EPT) as a test statistic. This is not correct. Yates (1934 - Pearson's X²), Lewontin & Felsenstein (1965 - X²), Agresti & Wackerly (1977 - X², Kendall's tau, Kruskal & Goodman's gamma), Klotz (1966 - Wilcoxon), Klotz & Teng (1977 - Kruskall & Wallis' H), Larntz (1978 - X², loglike-lihood-ratio statistic G², Freeman & Tukey statistic), and several others have investigated exact tests with other statistics than the EPT. In fact, Bennett & Nakamura (1963) are incorrectly cited as they investigated both X² and G², rather than EPT. Also, Freeman & Halton (1951) are incorrectly cited for they generalized Fisher's exact test to pxq tables and not 2xq tables as stated. And they are even predated by Yates (1934) who extended the test to 2×3 tables.     

Bayes estimation of Gompertz distribution parameters and acceleration factor under partially accelerated life tests with type-I censoring

In this paper, the Bayesian approach is applied to the estimation problem in the case of step stress partially accelerated life tests with two stress levels and type-I censoring. Gompertz distribution is considered as a lifetime model. The posterior means and posterior variances are derived using the squared-error loss function. The Bayes estimates cannot be obtained in explicit forms. Approximate Bayes estimates are computed using the method of Lindley [D.V. Lindley, Approximate Bayesian methods, Trabajos Estadistica 31 (1980), pp. 223–237]. The advantage of this proposed method is shown. The approximate Bayes estimates obtained under the assumption of non-informative priors are compared with their maximum likelihood counterparts using Monte Carlo simulation.     

Selection of suitable prior for the Bayesian mixture of a class of lifetime distributions under type-I censored datasets

This paper explores the study on mixture of a class of probability density functions under type-I censoring scheme. In this paper, we mold a heterogeneous population by means of a two-component mixture of the class of probability density functions. The parameters of the class of mixture density functions are estimated and compared using the Bayes estimates under the squared-error and precautionary loss functions. A censored mixture dataset is simulated by probabilistic mixing for the computational purpose considering particular case of the Maxwell distribution. Closed-form expressions for the Bayes estimators along with their posterior risks are derived for censored as well as complete samples. Some stimulating comparisons and properties of the estimates are presented here. A factual dataset has also been for illustration.     

Robust Bayesian analysis of generalized inverted family of distributions

ABSTRACT

In the current study we develop the robust Bayesian inference for the generalized inverted family of distributions (GIFD) under an ε-contamination class of prior distributions for the shape parameter α, with different possibilities of known and unknown scale parameter. We used Type II censoring and Bartholomew sampling scheme (1963) for the following derivations under the squared-error loss function (SELF) and linear exponential (LINEX) loss function : ML-II Bayes estimators of the i) parameters; ii) Reliability function and; iii) Hazard function. We also present simulation study and analysis of a real data set.     

Assessing the contribution of R&;D to total factor productivity—a Bayesian approach to account for heterogeneity and heteroskedasticity

This paper proposes a hierarchical Bayes estimator for a panel data random coefficient model with heteroskedasticity to assess the contribution of R&D capital to total factor productivity. Based on Hall (1993) data for 323 US firms over 1976–1990, we find that there appear to have substantial unobserved heterogeneity and heteroskedasticity across firms and industries that support the use of our Bayes inference procedure. We find much higher returns to R&D capital and a more pronounced downswing for the 1981–1985 period, followed by a more pronounced upswing than those yielded by the conventional feasible generalized least squares estimators or other estimates. The estimated elasticities of R&D capital are 0.062 for 1976–1980, 0.036 for 1981–1985 and 0.081 for 1986–1990, while the estimated elasticities of ordinary capital are much more stable over these periods.     

An extreme null distribution approach to the problem of paired comparisons

David (1963) and Davidson & Farquhar(1976) contain extensive bibliographies of proposed approaches to problems involving paired comparisons. However, each of the discussed methods that is based on a hypothesis test, relies heavily on the assumption that all paired comparisons are made independently. In this paper we eliminate this assumption and develop a new procedure based on an adaptation of a statistic considered by Kendall & Babington Smith (1940). We show that their original test procedure substantially underestimates the true significance level if the comparisons are not made independently. Our modification utilizes the approach developed in Costello & Wolfe (1985) for the problem of agreement between two groups of judges and relies heavily on computer-generated tables.     

Estimation of multiple gamma scale-parameters: bayes estimation subject to uniform domination

Simultaneous estimation of p gamma scale-parameters is considered under squared-error loss. The problem of minimizing, subject to uniform risk domination, the Bayes risk (or more generally the posterior expected loss) against certain conjugate or mixtures of conjugate priors is considered. Rather surprisingly, it is shown that the minimization can be done conditionally, thus avoiding variational arguments. Relative savings loss (and a posterior version thereof) are found, and it is found that in the most favorable situations, Bayesian robustness can be achieved without sacrificing substantial subjective Bayesian gains.     

A remark on estimating the mean of a normal distribution with known coefficient of variation

Let X₁, X₂, …, X_n be iid N(μ, aμ²) (a>0) random variables with an unknown mean μ>0 and known coefficient of variation (CV) √a. The estimation of μ is revisited and it is shown that a modified version of an unbiased estimator of μ [cf. Khan RA. A note on estimating the mean of a normal distribution with known CV. J Am Stat Assoc. 1968;63:1039–1041] is more efficient. A certain linear minimum mean square estimator of Gleser and Healy [Estimating the mean of a normal distribution with known CV. J Am Stat Assoc. 1976;71:977–981] is also modified and improved. These improved estimators are being compared with the maximum likelihood estimator under squared-error loss function. Based on asymptotic consideration, a large sample confidence interval is also mentioned.     

Size and shape analysis of palmar interdigital areas

As pointed out by Hauser & Abraham, dermatoglyphic methodology investigates either shape or size rather than both variables simultaneously. Hauser & Abraham have given a method to analyse simultaneously the shape and size of the third palmar interdigital area. We describe analytically their selection of the shape variables and indicate how their method can fit into a general framework of size and shape analysis. Their data are reanalysed and we describe some advantages of their choice of variables. However, their size and shape variables are highly correlated.     

From Ordinary to Shrinkage Square-Root Estimators

Consider a skewed population. Suppose an intelligent guess could be made about an interval that contains the population mean. There may exist biased estimators with smaller mean squared error than the arithmetic mean within such an interval. This article indicates when it is advisable to shrink the arithmetic mean towards a guessed interval using root estimators. The goal is to obtain an estimator that is better near the average of natural origins. An estimator proposed. This estimator contains the Thompson (1968 Thompson , J. R. ( 1968 ). Accuracy borrowing in the estimation of the mean by shrinkage towards an interval . J. Amer. Statist. Assoc. 63 : 953 – 963 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) ordinary shrinkage estimator, the Jenkins et al. (1973 Jenkins , O. C. , Ringer , L. J. , Hartley , H. O. ( 1973 ). Root estimators . J Amer. Statist. Assoc. 68 : 414 – 419 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) square-root estimator, and the arithmetic sample mean as special cases. The bias and the mean squared error of the proposed more general estimator is compared with the three special cases. Shrinkage coefficients that yield minimum mean squared error estimators are obtained. The proposed estimator is considerably more efficient than the three special cases. This remains true for highly skewed populations. The merits of the proposed shrinkage square-root estimator are supported by the results of numerical and simulation studies.     

Control charts based on linear combinations of order statistics

The last 20 years have seen an increasing emphasis on statistical process control as a practical approach to reducing variability in industrial applications. Control charts are used to detect problems such as outliers or excess variability in subgroup means that may have a special cause. We describe an approach to the computation of control limits for exponentially weighted moving average control charts where the usual statistics in classical charts are replaced by linear combinations of order statistics; in particular, the trimmed mean and Gini's mean difference instead of the mean and range, respectively. Control limits are derived, and simulated average run length experiments show the trimmed control charts to be less influenced by extreme observations than their classical counterparts, and lead to tighter control limits. An example is given that illustrates the benefits of the proposed charts. parameters; see, for example, Hunter (1986) and Montgomery (1996). On the other hand, EWMA charts have been shown to be more efficient than Shewharttype charts in detecting small shifts in the process mean; see, for example, Ng & Case (1989), Crowder (1989), Lucas & Saccucci (1990), Amin & Searcy (1991) and Wetherill & Brown (1991). In fact, the EWMA control chart has become popular for monitoring a process mean; see Hunter (1986) for a good discussion. More recently, EWMA charts have been developed for monitoring process variability;     

Some Modifications for Choosing Ridge Parameters

Standard least square regression can produce estimates having a large mean squares error (MSE) when predictor variables are highly correlated or multicollinear. In this article, we propose four modifications to choose the ridge parameter (K) when multicollinearity exists among the columns of the design matrix. The proposed new estimators are extended versions of that suggested by Khalaf and Shukur (2005 Khalaf , G. , Shukur , G. ( 2005 ). Choosing ridge parameter for regression problems . Commun. Statist. A 34 : 1177 – 1182 . [CSA] [Taylor & Francis Online] , [Google Scholar]). The properties of these estimators are compared with those of Hoerl and Kennard (1970a Hoerl , A. E. , Kennard , R. W. ( 1970a ). Ridge regression: biased estimation for non-orthogonal problems . Tech. . 12 : 55 – 67 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and the OLS using the MSE criterion. All estimators under consideration are evaluated using simulation techniques under certain conditions where a number of factors that may affect their properties have been varied. In addition, it is shown that at least one of the proposed estimators either has a smaller MSE than the others or is the next best otherwise.