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.     

A class of rectangle-screened multivariate normal distributions and its applications

A screening problem is tackled by proposing a parametric class of distributions designed to match the behavior of the partially observed screened data. This class is obtained from the nontruncated marginal of the rectangle-truncated multivariate normal distributions. Motivations for the screened distribution as well as some of the basic properties, such as its characteristic function, are presented. These allow us a detailed exploration of other important properties that include closure property in linear transformation, in marginal and conditional operations, and in a mixture operation as well as the first two moments and some sampling distributions. Various applications of these results to the statistical modelling and data analysis are also provided.     

Bayesian multivariate normal analysis under the extended reflected normal loss function

This paper considers the Bayesian analysis of the multivariate normal distribution under a new and bounded loss function, based on a reflection of the multivariate normal density function. The Bayes estimators of the mean vector can be derived for an arbitrary prior distribution of [d]. When the covariance matrix has an inverted Wishart prior density, a Bayes estimator of[d] is obtained under a bounded loss function, based on the entropy loss. Finally the admissibility of all linear estimators c[d]+ d for the mean vector is considered     

Using mathematical programming to compute singlular multivariate normal probablities

This paper describes two new, mathematical programming-based approaches for evaluating general, one- and two-sidedp-variate normal probabilities where the variance-covariance matrix (of arbitrary structure) is singular with rankr(r<pand r and p can be of unlimited dimensions. In both cases, principal components are used to transform the original, ill-definedp-dimensional integral into a well-definedrdimensional integral over a convex polyhedron. The first algorithm that is presented uses linear programming coupled with a Gauss-Legendre quadrature scheme to compute this integral, while the second algorithm uses multi-parametric programming techniques in order to significantly reduce the number of optimization problems that need to be solved. The application of the algorithms is demonstrated and aspects of computational performance are discussed through a number of examples, ranging from a practical problem that arises in chemical engineering to larger, numerical examples.     

A note on the inequalities for tail probability of the multivariate normal distribution

In this paper inequalities given by Harkness Godambe (1976) for the rail probabilities of the multivariate normal distribution in the equicorrelated case are improved by using the properties of the characteristic roots of a matrix and of the convex function.     

The inverse problem of multivariate and matrix-variate skew normal distributions

In this paper, we prove that the joint distribution of random vectors Z ₁ and Z ₂ and the distribution of Z ₂ are skew normal provided that Z ₁ is skew normally distributed and Z ₂ conditioning on Z ₁ is distributed as closed skew normal. Also, we extend the main results to the matrix variate case.     

Wavelet threshold based on Stein's unbiased risk estimators of restricted location parameter in multivariate normal

In this paper, the problem of estimating the mean vector under non-negative constraints on location vector of the multivariate normal distribution is investigated. The value of the wavelet threshold based on Stein''s unbiased risk estimators is calculated for the shrinkage estimator in restricted parameter space. We suppose that covariance matrix is unknown and we find the dominant class of shrinkage estimators under Balance loss function. The performance evaluation of the proposed class of estimators is checked through a simulation study by using risk and average mean square error values.     

A Unified Approach to the Approximation of Multivariate Densities

Approximation of a density by another density is considered in the case of different dimensionalities of the distributions. The results have been derived by inverting expansions of characteristic functions with the help of matrix techniques. The approximations obtained are all functions of cumulant differences and derivatives of the approximating density. The multivariate Edgeworth expansion follows from the results as a special case. Furthermore, the density functions of the trace and eigenvalues of the sample covariance matrix are approximated by the multivariate normal density and a numerical example is given     

Bivariate,multivariate, and matrix variate normal characterizations: A brief survey II

The paper entitled “Bivariate and Multivariate Normal Characterizations: A Brief Survey,” by Hamedani, which was published in 1992, covered the published characterizations of bivariate and multivariate normal (MVN) distributions from 1941 to 1991. The present work is a follow-up to the 1991/1992 survey which includes not only characterizations of the bivariate and MVN distributions, but also characterizations of the matrix variate normal distribution, which have appeared from 1991/1992 to the present.     

Analytical expression for the integrated squared density partial derivative of a multivariate normal mixture distribution

This paper describes the derivation of the analytical expression for the integrated squared density partial derivative (ISDPD) in a multivariate normal mixture model. The analytical expression of the ISDPD is derived for arbitrary dimensions with partial derivative orders up to four. Although the value of the ISDPD can be obtained by using the common numerical integration via mathematical software (such as Maple, Mathematica, Matlab, etc), it suffers from the limitation of computation time when the dimension or the number of mixture components of the considered multivariate normal mixture model are large. Moreover, numerical comparison indicates the benefits of speed offered by our proposed analytical expression are far superior to the numerical integration performed by Maple. With this analytical expression, the ISDPD can apace be calculated with no assistance of numerical integration.     

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.     

Prediction from a normal model using a generalized inverse Gaussian prior

In this paper, we derive prediction distribution of future response(s) from the normal distribution assuming a generalized inverse Gaussian (GIG) prior density for the variance. The GIG includes as special cases the inverse Gaussian, the inverted chi-squared and gamma distributions. The results lead to Bessel-type prediction distributions which is in contrast with the Student-t distributions usually obtained using the inverted chi-squared prior density for the variance. Further, the general structure of GIG provides us with new flexible prediction distributions which include as special cases most of the earlier results obtained under normal-inverted chi-squared or vague priors.     

Improved hypothesis testing in a general multivariate elliptical model

This paper investigates improved testing inferences under a general multivariate elliptical regression model. The model is very flexible in terms of the specification of the mean vector and the dispersion matrix, and of the choice of the error distribution. The error terms are allowed to follow a multivariate distribution in the class of the elliptical distributions, which has the multivariate normal and Student-t distributions as special cases. We obtain Skovgaard's adjusted likelihood ratio (LR) statistics and Barndorff-Nielsen's adjusted signed LR statistics and we compare the methods through simulations. The simulations suggest that the proposed tests display superior finite sample behaviour as compared to the standard tests. Two applications are presented in order to illustrate the methods.     

A matching prior for the product of normal means based on the modified profile likelihood

In this paper, we develop a matching prior for the product of means in several normal distributions with unrestricted means and unknown variances. For this problem, properly assigning priors for the product of normal means has been issued because of the presence of nuisance parameters. Matching priors, which are priors matching the posterior probabilities of certain regions with their frequentist coverage probabilities, are commonly used but difficult to derive in this problem. We developed the first order probability matching priors for this problem; however, the developed matching priors are unproper. Thus, we apply an alternative method and derive a matching prior based on a modification of the profile likelihood. Simulation studies show that the derived matching prior performs better than the uniform prior and Jeffreys’ prior in meeting the target coverage probabilities, and meets well the target coverage probabilities even for the small sample sizes. In addition, to evaluate the validity of the proposed matching prior, Bayesian credible interval for the product of normal means using the matching prior is compared to Bayesian credible intervals using the uniform prior and Jeffrey’s prior, and the confidence interval using the method of Yfantis and Flatman.     

A generalized two-piece skew normal distribution and some of its properties

In this paper, we develop a generalized version of the two-piece skew normal distribution of Kim [On a class of two-piece skew-normal distributions, Statistics 39(6) (2005), pp. 537–553] and derive explicit expressions for its distribution function and characteristic function and discuss some of its important properties. Further estimation of the parameters of the generalized distribution is carried out.     

Computing some bayes estimates for the mean of a normal distribution

For sampling from a normal population with unknown mean, two families of prior densities for the mean are discussed. The corresponding posterior densities are found. A data analyst may choose a prior from these families to represent prior beliefs and then compute the corresponding Bayes estimator, using the techniques discussed.