On selecting the best treatment in a generalized linear model

The problem of selecting the best treatment is studied under generalized linear models. For certain balanced designs, it is shown that simple rules are Bayes with respect to any non-informative prior on the treatment effects under any monotone invariant loss. When the nuisance parameters such as block effects are assumed to follow a uniform (improper) prior or a normal prior, Bayes rules are obtained for the normal linear model under more suitable balanced designs, keeping the generality of the loss and the generality of the non-informativeness on the prior of the treatment effects. These results are extended to certain types of informative priors on the treatment effects. When the designs are unbalanced, algorithms based on the Gibbs sampler and the Laplace method are provided to compute the Bayes rules.     

Objective priors for hypothesis testing in one-way random effects models

The Bayes factor is a key tool in hypothesis testing. Nevertheless, the important issue of which priors should be used to develop objective Bayes factors remains open. The authors consider this problem in the context of the one-way random effects model. They use concepts such as orthogonality, predictive matching and invariance to justify a specific form of the priors for common parameters and derive the intrinsic and divergence based prior for the new parameter. The authors show that both intrinsic priors or divergence-based priors produce consistent Bayes factors. They illustrate the methods and compare them with other proposals.     

Bayes estimation of a langevin mean direction

A Langevin distribution with two parameters (mean direction and concentration parameter) has been extensively used for modeling and analyzing problems related to directional data. In this article, we examine the estimation problem for the mean direction. Bayes estimators are derived with respect to a conjugate as well as the Jeffreys’ priors. Further in case of unknown concentration parameter, other priors are also chosen. An extensive analysis of risk behavior of Bayes estimators is carried out with the help of simulations.     

Expansion estimation by Bayes rules

In the problem of estimating a location parameter in any symmetric unimodal location parameter model, we demonstrate that Bayes rules with respect to squared error loss can be expanders for some priors that belong to the family of all symmetric priors. That generalizes the results obtained by DasGupta and Rubin for the one dimensional case. We also consider symmetric priors which either have an appropriate point mass at 0 or are unimodal, and prove that under the latter condition all Bayes rules are shrinkers. Results of such nature are important, for example, in wavelet based function estimation and data denoising, where shrinkage of wavelet coefficients is associated with smoothing the data. We illustrate the results using FIAT stock market data.     

On asymptotic optimality of bayes empirical bayes estimators

In an empirical Bayes decision problem, a prior distribution ? is placed on a one-dimensfonal family G of priors G_w, wεΩ, to produce a Bayes empirical Bayes estimator, The asymptotic optimaiity of the Bayes estimator is established when the support of ? is Ω and the marginal distributions H_w have monotone likelihood ratio and continuous Kullback-Leibler information number.     

A note on james-stein and bayes empiricl bayes estimators

In an empirical Bayes decision problem, a simple class of estimators is constructed that dominate the James-Stein

estimator, A prior distribution A is placed on a restricted (normal) class G of priors to produce a Bayes empirical Bayes estimator, The Bayes empirical Bayes estimator is smooth, admissible, and asymptotically optimal. For certain A rate of convergence to minimum Bayes risk is 0(n^-1)uniformly on G. The results of a Monte Carlo study are presented to demonstrate the favorable risk bebhavior of the Bayes estimator In comparison with other competitors including the James-Stein estimator.     

A balanced approach to region estimation with tables for the normal model

Practitioners of statistics are too often guilty of routinely selecting a 95% confidence level in interval estimation and ignoring the sample size and the expected size of the interval. One way to balance coverage and size is to use a loss function in a decision problem. Then either the Bayes risk or usual risk (if a pivotal quantity exists) may be minimized. It is found that some non-Bayes solutions are equivalent to Bayes results based on non-informative priors. The decision theory approach is applied to the mean and standard deviation of the univariate normal model and the mean of the multivariate normal. Tables are presented for critical values, expected size, confidence and sample size.     

Alternative Bayes factors for model selection

Several alternative Bayes factors have been recently proposed in order to solve the problem of the extreme sensitivity of the Bayes factor to the priors of models under comparison. Specifically, the impossibility of using the Bayes factor with standard noninformative priors for model comparison has led to the introduction of new automatic criteria, such as the posterior Bayes factor (Aitkin 1991), the intrinsic Bayes factors (Berger and Pericchi 1996b) and the fractional Bayes factor (O'Hagan 1995). We derive some interesting properties of the fractional Bayes factor that provide justifications for its use additional to the ones given by O'Hagan. We further argue that the use of the fractional Bayes factor, originally introduced to cope with improper priors, is also useful in a robust analysis. Finally, using usual classes of priors, we compare several alternative Bayes factors for the problem of testing the point null hypothesis in the univariate normal model.     

Empirical bayes estimation of binomial parameter with symmetric priors

This paper deals with the problem of estimating the binomial parameter via the nonparametric empirical Bayes approach. This estimation problem has the feature that estimators which are asymptotically optimal in the usual empirical Bayes sense do not exist (Robbins (1958, 1964)), However, as pointed out by Liang (1934) and Gupta and Liang (1988), it is possible to construct asymptotically optimal empirical Bayes estimators if the unknown prior is symmetric about the point 1/2, In this paper, assuming symmetric priors a monotone empirical Bayes estimator is constructed by using the isotonic regression method. This estimator is asymptotically optimal in the usual empirical Bayes sense. The corresponding rate of convergence is investigated and shown to be of order n^-1, where n is the number of past observations at hand.     

Intrinsic Priors for Testing Two Normal Means with Intrinsic Bayes Factors

ABSTRACT

In this article we consider the problem of comparing two normal means with unknown common variance using a Bayesian approach. Conventional Bayes factors with improper non informative priors are not well defined. The intrinsic Bayes factors are used to overcome such a difficulty. We derive intrinsic priors whose Bayes factors are asymptotically equivalent to the corresponding intrinsic Bayes factors. We illustrate our results with numerical examples.     

Bayes estimation and prediction of the two-parameter gamma distribution

In this article, the Bayes estimates of two-parameter gamma distribution are considered. It is well known that the Bayes estimators of the two-parameter gamma distribution do not have compact form. In this paper, it is assumed that the scale parameter has a gamma prior and the shape parameter has any log-concave prior, and they are independently distributed. Under the above priors, we use Gibbs sampling technique to generate samples from the posterior density function. Based on the generated samples, we can compute the Bayes estimates of the unknown parameters and can also construct HPD credible intervals. We also compute the approximate Bayes estimates using Lindley's approximation under the assumption of gamma priors of the shape parameter. Monte Carlo simulations are performed to compare the performances of the Bayes estimators with the classical estimators. One data analysis is performed for illustrative purposes. We further discuss the Bayesian prediction of future observation based on the observed sample and it is seen that the Gibbs sampling technique can be used quite effectively for estimating the posterior predictive density and also for constructing predictive intervals of the order statistics from the future sample.     

On the problem of selecting good populations

The problem of selecting good populations out of k normal populations is considered in a Bayesian framework under exchangeable normal priors and additive loss functions. Some basic approximations to the Bayes rules are discussed. These approximations suggest that some well-known classical rules are "approximate" Bayes rules. Especially, it is shown that Gupta-type rules are extended Bayes with respect to a family of the exchangeable normal priors for any bounded and additive loss function. Furthermore, for a simple loss function, the results of a Monte Carlo comparison of Gupta-type rules and Seal-type rules are presented. They indicate that, in general, Gupta-type rules perform better than Seal-type rules     

Gamma-minimax results for the class of unimodal priors

Olman and Shmundak proved 1985 that in estimating a bounded normal mean under squared error loss the Bayes estimator with respect to the uniform distribution on the parameter interval is gamma-minimax when the parameter interval is sufficiently small and the class of priors consists of all symmetric and unimodal distributions. Recently, one of the authors showed that this result remains valid for quite general families of distributions which satisfy some regularity conditions. In the present paper a generalization to the class of unimodal priors with fixed mode is derived. It is proved that the Bayes estimator with respect to a suitable mixture of two uniform distributions is gamma-minimax for sufficiently small parameter intervals. To that end appropriate characterizations of a saddle point in the corresponding statistical games are established. Some results of a numerical study are presented.     

Admissible estimation in an one parameter nonregular family of absolutely continuous distributions

Consider the problem of estimating under entropy loss an arbitrarily positive, strictly increasing or decreasing parametric function based on a sample of size n in an one parameter noregular family of absolutly continuous distributions with both endpoints of the support depending on a single parameter. We first provide sufficient conditions for the admissibility of generalized Bayes estimator with respect to some specific priors and then treat several examples which illustrate the admissibility of best invariant estimators is some location or scale parameter problems.     

A family of admissible minimax estimators of the mean of a multivariate,normal distribution

Let X has a p-dimensional normal distribution with mean vector θ and identity covariance matrix I. In a compound decision problem consisting of squared-error estimation of θ, Strawderman (1971) placed a Beta (α, 1) prior distribution on a normal class of priors to produce a family of Bayes minimax estimators. We propose an incomplete Gamma(α, β) prior distribution on the same normal class of priors to produce a larger family of Bayes minimax estimators. We present the results of a Monte Carlo study to demonstrate the reduced risk of our estimators in comparison with the Strawderman estimators when θ is away from the zero vector.     

The equivalence of Bayes and robust Bayes estimators for various loss functions

The problem of Bayes and robust Bayes estimation for various bounded and/or symmetric loss functions in a normal model with conjugate and non-informative prior distributions is considered. The prior distribution is not fully specified and covers the conjugate family of priors. It is of interest to know that the Bayes and robust Bayes estimators for symmetric losses are the same as those for the standard square-error loss function.     

Admissible estimation in an one parameter nonregular family of absolutely continuous distributions

Consider the problem of estimating under squared error loss an arbitrarily positive, strictly increasing or decreasing parametric function based on a sample of size n in an one parameter nonregular family of absolutly continuous distributions with both endpoints of the support depending on a single parameter. We first provide sufficient conditions for the admissibility of generalized Bayes estimators with respect to some specific priors and then treat several examples which illustrate the admissibility of best invariant estimators in some location or scale parameter problems.     

Bayes minimax estimation of a bernoulli p in a restricted parameter space

In many estimation problems the parameter of interest is known,a priori, to belong to a proper subspace of the natural parameter space. Although useful in practice this type of additional information can lead to surprising theoretical difficulties. In this paper the problem of minimax estimation of a Bernoulli pwhen pis restricted to a symmetric subinterval of the natural parameter space is considered. For the sample sizes n = 1,2,3, and 4 least favorable priors with finite support are provided and the corresponding Bayes estimators are shown to be minimax. For n = 5 and 6 the usual constant risk minimax estimator is shown to be the Bayes minimax estimator corresponding to a least favorable prior with finite support, provided the restriction on the parameter space is not too tight.     

On the Choice of a Prior for Bayesian D-Optimal Designs for the Logistic Regression Model with a Single Predictor

The Bayesian design approach accounts for uncertainty of the parameter values on which optimal design depends, but Bayesian designs themselves depend on the choice of a prior distribution for the parameter values. This article investigates Bayesian D-optimal designs for two-parameter logistic models, using numerical search. We show three things: (1) a prior with large variance leads to a design that remains highly efficient under other priors, (2) uniform and normal priors lead to equally efficient designs, and (3) designs with four or five equidistant equally weighted design points are highly efficient relative to the Bayesian D-optimal designs.     

On selecting the treatment with the largest probability of survival

Suppose there are k(>= 2) treatments and each treatment is a Bernoulli process with binomial sampling. The problem of selecting a random-sized subset which contains the treatment with the largest survival probability (reliability or probability of success) is considered. Based on the ideas from both classical approaches and general Bayesian statistical decision approach, a new subset selection procedure is proposed to solve this kind of problem in both balanced and unbalanced designs. Comparing with the classical procedures, the proposed procedure has a significantly smaller selected subset. The optimal properties and performance of it were examined. The methods of selecting and fitting the priors and the results of Monte Carlo simulations on selected important cases are also studied.