Bayesian multivariate normal analysis with a wishart prior

This paper considers the Bayesian analysis of the multivariate normal distribution when its covariance matrix has a Wishart prior density under the assumption of a multivariate quadratic loss function. New flexible marginal posterior distributions of the mean μ and of the covariance matrix Σ are developed and univariate cases with graphical representations are given.     

On Spiring's normal loss function

Loss functions express the loss to society, incurred through the use of a product, in monetary units. Underlying this concept is the notion that any deviation from target of any product characteristic implies a degradation in the product performance and hence a loss. Spiring (1993), in response to criticisms of the quadratic loss function, developed the reflected normal loss function, which is based on the normal density function. We give some modifications of these loss functions to simplify their application and provide a framework for the reflected normal loss function that accomodates a broader class of symmetric loss situations. These modifications also facilitate the unification of both of these loss functions and their comparison through expected loss. Finally, we give a simple method for determing the parameters of the modified reflected normal loss function based on loss information for multiple values of the product characteristic, and an example to illustrate the flexibility of the proposed model and the determination of its parameters.     

Estimation of a multivariate normal covariance matrix under a certain structure

We consider the estimation of Σ of the p-dimensional normal distribution N_p (0, Σ) when Σ?=?θ₀ I_p ?+?θ₁ aa′, where a is an unknown p-dimensional normalized vector and θ₀?>?0, θ₁?≥?0 are also unknown. First, we derive the restricted maximum likelihood (REML) estimator. Second, we propose a new estimator, which dominates the REML estimator with respect to Stein's loss function. Finally, we carry out Monte Carlo simulation to investigate the magnitude of the new estimator's superiority.     

The reflected normal loss function

The use of loss functions in quality assurance has grown steadily with the introduction of Taguchi's philosophy. The quadratic loss function has been used by decision-theoretic statisticians and economists for many years. Taguchi uses a modified form of the quadratic loss function to demonstrate the need to consider proximity to the target while assessing quality. Several authors have suggested that the traditional quadratic loss function is inadequate for assessing quality and quality improvement. A new, easily understood loss function, based on a reflection of the normal density function, is presented, and some associated statistical properties discussed.     

On the Admissibility of Linear Estimators in a Multivariate Normal Distribution Under LINEX Loss Function

Consider an estimation problem of a linear combination of population means in a multivariate normal distribution under LINEX loss function. Necessary and sufficient conditions for linear estimators to be admissible are given. Further, it is shown that the result is an extension of the quadratic loss case as well as the univariate normal case.     

Bayesian analysis in multivariate regression models with conjugate priors

In this paper, we consider the full rank multivariate regression model with matrix elliptically contoured distributed errors. We formulate a conjugate prior distribution for matrix elliptical models and derive the posterior distributions of mean and scale matrices. In the sequel, some characteristics of regression matrix parameters are also proposed.     

Simultaneous estimation of the multivariate normal mean under balanced loss function

This paper considers simultaneous estimation of multivariate normal mean vector using Zellner's(1994) balanced loss function which is defined as follows:

where 0 < w < 1 and for i = 1,…,p and j = 1,…,n, X_ij is distributed as normal with mean θi and variance 1. It is shown that the sample mean, X, is admissible when p <3. For p ≥3, we obtain that James-Stein type estimator which has uniformly smaller risk than that of sample mean X.     

Linearly admissible estimators of mean vector with respect to balanced loss function in multivariate statistics

The authors discuss the bias of the estimate of the variance of the overall effect synthesized from individual studies by using the variance weighted method. This bias is proven to be negative. Furthermore, the conditions, the likelihood of underestimation and the bias from this conventional estimate are studied based on the assumption that the estimates of the effect are subject to normal distribution with common mean. The likelihood of underestimation is very high (e.g. it is greater than 85% when the sample sizes in two combined studies are less than 120). The alternative less biased estimates for the cases with and without the homogeneity of the variances are given in order to adjust for the sample size and the variation of the population variance. In addition, the sample size weight method is suggested if the consistence of the sample variances is violated Finally, a real example is presented to show the difference by using the above three estimate methods.     

E-Bayesian and hierarchical Bayesian estimations for the system reliability parameter based on asymmetric loss function

This paper explores properties of the E-Bayesian and hierarchical Bayesian estimations of the system reliability parameter. E-Bayesian estimation and hierarchical Bayesian estimation of Pascal distribution's parameter under two loss function, LINEX loss function and entropy loss function can be found. We obtained limits of that the E-Bayesian estimation and hierarchical Bayesian estimation are equal. A Monte Carlo simulation is used to compare performances of the two methods.     

Statistical applications of the multivariate skew normal distribution

Azzalini and Dalla Valle have recently discussed the multivariate skew normal distribution which extends the class of normal distributions by the addition of a shape parameter. The first part of the present paper examines further probabilistic properties of the distribution, with special emphasis on aspects of statistical relevance. Inferential and other statistical issues are discussed in the following part, with applications to some multivariate statistics problems, illustrated by numerical examples. Finally, a further extension is described which introduces a skewing factor of an elliptical density.     

Linearly admissible estimators of the common mean parameter under balanced loss function

Admissibility of linear estimators of the common mean parameter is investigated in the context of a linear model under balanced loss function. Sufficient and necessary conditions for linear estimators to be admissible in classes of homogeneous and non homogeneous linear estimators are obtained, respectively.     

a bayesian analysis of the multivariate mixed-linear model

The Bayesian analysis of the multivariate mixed linear model is considered. The exact posterior distribution for the fixed effects matrix and the error covariance matrix are obtained. The exact posterior means and variances of the Bayesian estimators for the covariance matrices of random effects are also derived. These posterior moments are computed without constrained optimization and numerical integration. The calculations are feasible for arbitrary models. Reasonable approximations for the posterior distributions for the covariance matrices associated with the random effects are obtained also. Results are illustrated with a numerical example.     

Conjugate analysis of multivariate normal data with incomplete observations

The authors discuss prior distributions that are conjugate to the multivariate normal likelihood when some of the observations are incomplete. They present a general class of priors for incorporating information about unidentified parameters in the covariance matrix. They analyze the special case of monotone patterns of missing data, providing an explicit recursive form for the posterior distribution resulting from a conjugate prior distribution. They develop an importance sampling and a Gibbs sampling approach to sample from a general posterior distribution and compare the two methods.     

The performance of the adaptive optimal estimator under the extended balanced loss function

The adaptive optimal estimator of Farebrother (1975) is discussed by many authors, but the goodness of fitted model criterion that is used to investigate the performance of estimators is quite often ignored. Shalabh, Toutenburg, and Heumann (2009) proposed the extended balanced loss function in which the mean squared error and the Zellner's balanced loss function are just special cases of it. In this paper, we discuss the performance of the adaptive optimal estimator of Farebrother (1975) under the extended balanced loss function. Moreover, a Monte Carlo simulation experiment is conducted to examine the performance of the estimator in finite samples.     

Admissibility for linear estimators of regression coefficient in a general Gauss–Markoff model under balanced loss function

A generalization of Zellner's balanced loss function is proposed according to unified theory of least squares under a general Gauss–Markoff model. Admissibility of linear estimators is investigated under the balanced loss function. And necessary and sufficient conditions that linear estimators are admissible in a class of homogeneous and nonhomogeneous linear estimators are obtained, respectively.     

A procedure to select a ML-II prior in a multivariate normal case

Berger (1985) derived a procedure to select a maximum likelihood II prior distribution. In this paper a method is suggested to construct such a prior distribution from a multivariate ε-contamination class of distributions. The method is illustrated by the conetruction of a ML-II prior in the multivariate normal case.     

New approximation of the normal distribution function

Two new and simple expressions, one for 0 ≤ x ≤ 2 and another for x > 2, for the normal distribution function, are developed which can be easily computed on desk calculators. They are also comparable in accuracy to the one developed by Patry and Keller (1964).     

Bayesian estimation of the dispersion matrix of a multivariate normal distribution

The estimation of the dispersion matrix of a multivariate normal distribution with zero mean on the basis of a random sample is discussed from a Bayesian view. An inverted-Wishart distribu- tion for the dispersion is taken, with its defining matrix of intraclass form. Some consistency properties are described. The posterior distribution is found and its mode investigated as a possible estimate in preference to that of maximum likelihood     

Stein-rule estimation under an extended balanced loss function

This paper extends the balanced loss function to a more general setup. The ordinary least squares estimator (OLSE) and Stein-rule estimator (SRE) are exposed to this general loss function with quadratic loss structure in a linear regression model. Their risks are derived when the disturbances in the linear regression model are not necessarily normally distributed. The dominance of OLSE and SRE over each other and the effect of departure from normality assumption of disturbances on the risk property are studied.