Some modifications of improved estimators of a normal variance

Consider the problem of estimating a normal variance based on a random sample when the mean is unknown. Scale equivariant estimators which improve upon the best scale and translation equivariant one have been proposed by several authors for various loss functions including quadratic loss. However, at least for quadratic loss function, improvement is not much. Herein, some methods are proposed to construct improving estimators which are not scale equivariant and are expected to be considerably better when the true variance value is close to the specified one. The idea behind the methods is to modify improving equivariant shrinkage estimators, so that the resulting ones shrink little when the usual estimate is less than the specified value and shrink much more otherwise. Sufficient conditions are given for the estimators to dominate the best scale and translation equivariant rule under the quadratic loss and the entropy loss. Further, some results of a Monte Carlo experiment are reported which show the significant improvements by the proposed estimators.     

Stein's idea and minimax admissible estimation of a multivariate normal mean

We consider estimation of a multivariate normal mean vector under sum of squared error loss.We propose a new class of minimax admissible estimator which are generalized Bayes with respect to a prior distribution which is a mixture of a point prior at the origin and a continuous hierarchical type prior. We also study conditions under which these generalized Bayes minimax estimators improve on the James–Stein estimator and on the positive-part James–Stein estimator.     

Asymptotic equivalence of the jackknife and infinitesimal jackknife variance estimators for some smooth statistics

The jackknife variance estimator and the infinitesimal jackknife variance estimator are shown to be asymptotically equivalent if the functional of interest is a smooth function of the mean or a trimmed L-statistic with Hölder continuous weight function.     

Sensitivity of minimaxity and admissibility in the estimation of a positive normal mean

On the problem of estimating a positive normal mean with known variance, it is well known that one minimax admissible estimator is the generalized Bayes one with respect to the non-informative prior measure, the Lebesgue measure, restricted on the positive half-line. When the true variance is misspecified, however, it is shown that this estimator does not always retain minimaxity and admisssibility. In particular, it is almost surely inadmissible in the misspecification case.     

Improved estimators for the GMANOVA problem with application to Monte Carlo simulation

The problem of finding classes of estimators which improve upon the usual (e.g., ML, LS) estimator of the parameter matrix in the GMANOVA model under (matrix) quadratic loss is considered. Classes of improved estimators are obtained via combining integration-by-parts methods for normal and Wishart distributions. Also considered is the application of control variates to achieve better efficiency in multipopulation multivariate simulation studies.     

Minimax confidence bound of the normal mean under an asymmetric loss function

This paper considers a minimax confidence bound of the normal mean under an asymmetric loss function. A minimax confidence bound is obtained for the case that the variance is known or unknown. The admissibility of the minimax confidence bound is also considered for the case of known variance.     

Simultaneous estimation of means of classified normal observations

Simultaneous estimation of normal means is considered for observations which are classified into several groups. In a one-way classification case, it is shown that an adaptive shrinkage estimator dominates a Stein-type estimator which shrinks observations towards individual class averages as Stein's (1966,Festschrift for J. Neyman, (ed. F. N. David), 351–366, Wiley, New York) does, and is minimax even if class sizes are small. Simulation results under quadratic loss show that it is slightly better than Stein's (1966) if between variances are larger than within ones. Further this estimator is shown to improve on Stein's (1966) with respect to the Bayes risk. Our estimator is derived by assuming the means to have a one-way classification structure, consisting of three random terms of grand mean, class mean and residual. This technique can be applied to the case where observations are classified into a two-stage hierarchy.     

A sequence of improvements over the James-Stein estimator

In this article, we consider the problem of estimating a p-variate (p ≥ 3) normal mean vector in a decision-theoretic setup. Using a simple property of the noncentral chi-square distribution, we have produced a sequence of smooth estimators dominating the James-Stein estimator and each improved estimator is better than the previous one. It is also shown by using a technique of [5]. J. Multivariate Anal.36 121–126) that our smooth estimators can be dominated by non-smooth estimators.     

Existence of the uniformly minimum risk equivariant estimators of parameters in a class of normal linear models

In this paper, we study the existence of the uniformly minimum risk equivariant (UMRE) estimators of parameters in a class of normal linear models, which include the normal variance components model, the growth curve model, the extended growth curve model, and the seemingly unrelated regression equations model, and so on. The necessary and sufficient conditions are given for the existence of UMRE estimators of the estimable linear functions of regression coefficients, the covariance matrixV and (trV)^α, where α > 0 is known, in the models under an affine group of transformations for quadratic losses and matrix losses, respectively. Under the (extended) growth curve model and the seemingly unrelated regression equations model, the conclusions given in literature for estimating regression coefficients can be derived by applying the general results in this paper, and the sufficient conditions for non-existence of UMRE estimators ofV and tr(V) are expanded to be necessary and sufficient conditions. In addition, the necessary and sufficient conditions that there exist UMRE estimators of parameters in the variance components model are obtained for the first time.     

Minimax Bayes estimators of a multivariate normal mean

In three or more dimensions it is well known that the usual point estimator for the mean of a multivariate normal distribution is minimax but not admissible with respect to squared Euclidean distance loss. This paper gives sufficient conditions on the prior distribution under which the Bayes estimator has strictly lower risk than the usual estimator. Examples are given for which the posterior density is useful in the formation of confidence sets.     

Approximated Bayes and empirical Bayes confidence intervals—The known variance case

In this paper hierarchical Bayes and empirical Bayes results are used to obtain confidence intervals of the population means in the case of real problems. This is achieved by approximating the posterior distribution with a Pearson distribution. In the first example hierarchical Bayes confidence intervals for the Efron and Morris (1975, J. Amer. Statist. Assoc., 70, 311–319) baseball data are obtained. The same methods are used in the second example to obtain confidence intervals of treatment effects as well as the difference between treatment effects in an analysis of variance experiment. In the third example hierarchical Bayes intervals of treatment effects are obtained and compared with normal approximations in the unequal variance case.Financially supported by the CSIR and the University of the Orange Free State, Central Research Fund.     

On confidence sequences for the mean vector of a multivariate normal distribution

Let X₁, X₂,… be idd random vectors with a multivariate normal distribution N(μ, Σ). A sequence of subsets {R_n(a₁, a₂,…, a_n), n ≥ m} of the space of μ is said to be a (1 − α)-level sequence of confidence sets for μ if P(μ R_n(X₁, X₂,…, X_n) for every n ≥ m) ≥ 1 − α. In this note we use the ideas of Robbins Ann. Math. Statist. 41 (1970) to construct confidence sequences for the mean vector μ when Σ is either known or unknown. The constructed sequence R_n(X₁, X₂, …, X_n) depends on Mahalanobis' or Hotelling's according as Σ is known or unknown. Confidence sequences for the vector-valued parameter in the general linear model are also given.     

Admissible minimax estimators of a mean vector of scale mixtures of multivariate normal distributions

The problem of estimating a mean vector of scale mixtures of multivariate normal distributions with the quadratic loss function is considered. For a certain class of these distributions, which includes at least multivariate-t distributions, admissible minimax estimators are given.     

Minimax estimation of a multivariate normal mean under polynomial loss

Let X be an observation from a p-variate (p ≥ 3) normal random vector with unknown mean vector θ and known covariance matrix

. The problem of improving upon the usual estimator of θ, δ⁰(X) = X, is considered. An approach is developed which can lead to improved estimators, δ, for loss functions which are polynomials in the coordinates of (δ ? θ). As an example of this approach, the loss L(δ, θ) = |δ ? θ|⁴ is considered, and estimators are developed which are significantly better than δ⁰. When

is the identity matrix, these estimators are of the form $\text{δ(X) = (1 ? (}\text{b}\text{(d + |X|}{}^2\text{)}\text{))X}$.     

Minimax estimation of a multivariate normal mean under arbitrary quadratic loss

Let X be a p-variate (p ≥ 3) vector normally distributed with mean θ and known covariance matrix

. It is desired to estimate θ under the quadratic loss (δ ? θ)^tQ(δ ? θ), where Q is a known positive definite matrix. A broad class of minimax estimators for θ is developed.     

Entropy loss and risk of improved estimators for the generalized variance and precision

Let the distributions of X(p×r) and S(p×p) be N(, I _r) and W _p(n, ) respectively and let them be independent. The risk of the improved estimator for || or {ei329-1} based on X and S under entropy loss (=d/|| –log(d/||)–1 or d||–log(d||)–1) is evaluated in terms of incomplete beta function of matrix argument and its derivative. Numerical comparison for the reduction of risk over the best affine equivariant estimator is given.Dedicated to Professor Yukihiro Kodama on his 60th birthday.     

Minimax and admissible minimax estimators of the mean of a multivariate normal distribution for unknown covariance matrix

Let X be a p-variate (p ≥ 3) vector normally distributed with mean μ and covariance Σ, and let A be a p × p random matrix distributed independent of X, according to the Wishart distribution W(n, Σ). For estimating μ, we consider estimators of the form δ = δ(X, A). We obtain families of Bayes, minimax and admissible minimax estimators with respect to the quadratic loss function (δ ? μ)′ Σ^?1(δ ? μ) where Σ is unknown. This paper extends previous results of the author [1], given for the case in which the covariance matrix of the distribution is of the form σ²I, where σ is known.     

On the possibilistic mean value and variance of multiplication of fuzzy numbers

In this paper, we introduce the definitions of the possibilistic mean, variance and covariance of multiplication of fuzzy numbers, and show some properties of these definitions. Then, we apply these definitions to build the possibilistic models of portfolio selection under the situations involving uncertainty over the time horizon, by considering the portfolio selection problem from the point of view of possibilistic analysis. Moreover, numerical experiments with real market data indicate that our approach results in better portfolio performance.     

A minimum average risk approach to shrinkage estimators of the normal mean

For the problem of estimating the normal mean based on a random sample X ₁,...,X _n when a prior value ₀ is available, a class of shrinkage estimators % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-qqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xHapdbiqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaubeaeqaje% aWbaGaamOBaaWcbeqdbaGafqiVd0MbaKaaaaGccaqGGaGaaiikaiaa% dUgacaGGPaGaeyypa0Jaam4AaiaacIcadaqfqaqabKqaahaacaqGUb% aaleqaneaacaqGubaaaOGaaiykaiaabccadaqfqaqabKqaahaacaWG% UbaaleqaneaaceqGybGbaebaaaGccaqGGaGaey4kaSIaaeiiaiaacI% cacaaIXaGaaeiiaiabgkHiTiaabccacaWGRbGaaiikamaavababeqc% baCaaiaab6gaaSqab0qaaiaabsfaaaGccaGGPaGaaiykamaavababe% qcbaCaaiaad6gaaSqab0qaaiabeY7aTbaaaaa!5615!\[\mathop {\hat \mu }\nolimits_n {\rm{ }}(k) = k(\mathop {\rm{T}}\nolimits_{\rm{n}} ){\rm{ }}\mathop {{\rm{\bar X}}}\nolimits_n {\rm{ }} + {\rm{ }}(1{\rm{ }} - {\rm{ }}k(\mathop {\rm{T}}\nolimits_{\rm{n}} ))\mathop \mu \nolimits_n \] is considered, where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-qqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xHapdbiqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaubeaeqaje% aWbaGaamOBaaWcbeqcdawaaiaadsfaaaGccaqGGaGaaeypaiaabcca% caWGUbWaaWbaaSqabeaacaaIXaGaai4laiaaikdaaaGccaGGOaWaa0% aaaeaacaWGybaaamaaBaaajeaWbaGaamOBaaWcbeaakiaabccacqGH% sislcaqGGaWaaubeaeqajeaWbaGaaGimaaWcbeqdbaGaaeiVdaaaki% aacMcacaqGGaGaae4laiabeccaGiabeo8aZbaa!4C33!\[\mathop T\nolimits_n {\rm{ = }}n^{1/2} (\overline X _n {\rm{ }} - {\rm{ }}\mathop {\rm{\mu }}\nolimits_0 ){\rm{ /}} \sigma \] and k is a weight function. For certain choices of k, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-qqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xHapdbiqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaubeaeqaje% aWbaGaamOBaaWcbeqdbaGafqiVd0MbaKaaaaGccaqGGaGaaiikaiaa% dUgacaGGPaaaaa!3CEE!\[\mathop {\hat \mu }\nolimits_n {\rm{ }}(k)\] coincides with previously studied preliminary test and shrinkage estimators. We consider choosing k from a natural non-parametric family of weight functions so as to minimize average risk relative to a specified prior p. We study how, by varying p, the MSE efficiency (relative to \-X) properties of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-qqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xHapdbiqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaubeaeqaje% aWbaGaamOBaaWcbeqdbaGafqiVd0MbaKaaaaGccaqGGaGaaiikaiaa% dUgacaGGPaaaaa!3CEE!\[\mathop {\hat \mu }\nolimits_n {\rm{ }}(k)\] can be controlled. In the process, a certain robustness property of the usual family of posterior mean estimators, corresponding to the conjugate normal priors, is observed.