Non-conjugate prior distribution assessment for multivariate normal sampling

Elicitation methods are proposed for quantifying expert opinion about a multivariate normal sampling model. The natural conjugate prior family imposes a relationship between the mean vector and the covariance matrix that can portray an expert's opinion poorly. Instead we assume that opinions about the mean and the covariance are independent and suggest innovative forms of question which enable the expert to quantify separately his or her opinion about each of these parameters. Prior opinion about the mean vector is modelled by a multivariate normal distribution and about the covariance matrix by both an inverse Wishart distribution and a generalized inverse-Wishart (GIW) distribution. To construct the latter, results are developed that give insight into the GIW parameters and their interrelationships. Certain of the elicitation methods exploit unconditional assessments as fully as possible, since these can reflect an expert's beliefs more accurately than conditional assessments. Methods are illustrated through an example.     

Second-order properties of improved two-stage procedure for a multivariate normal distribution

We consider the classical two problems: (I) fixed-size confidence region estimation problem, and (II) bounded risk point estimation problem, for the mean vector μ in the N_p(μ, ∑) population when ∑ is unknown and positive definite. Healy's (1956) solution given by a two-stage procedure for the prob¬lem (I) is improved in terms of the sample size. The second-order properties of the solution to each problem are verified under the assumptions that the maximum latent root which is simple of ∑ is bounded below for (I) and tr(∑) is bounded below for (II), by respective known and positive number.     

Inequalities for tail probabilities for the multivariate normal distribution

Inequalities for tail probabilities of the multivariate normal distribution are obtained, as a generalization of those given by Feller (1966). Upper and lower bounds are given in the equi-correlated case. For an arbitrary correlation matrix R, an upper bound is obtained, using a result of Slepian (1962) which asserts that certain multivariate normal probabilities are a non-decreasing function of correlations.     

On estimating a function of the mean of a multivariate normal distribution

The unique minimum variance of unbiased estimator is obtained for analysis functions of the mean of a multivariate normal distribution with either unknown covariance matrix or with covariance matrix of the form σ²v where σ² is unknown.     

Empirical bayes estimation of the mean in a multivariate normal distribution

We consider the problem of estimating the mean vector of a multivariate normal distribution under a variety of assumed structures among the parameters of the sampling and prior distributions. We adopt a pragmatic approach. We adopt distributional familites, assess hyperparmeters, and adopt patterned mean and coveariance structures when it is relatively simple to do so; alternatively, we use the sample data to estimate hyperparameters of prior distributions when assessment is a formidable task; such as the task of assessing parameters of multidimensional problems. James-Stein-like estimators are found to result. In some cases, we've been abl to show that the estimators proposed uniformly dominate the MLE's when measured with respect to quadratic loss functions.     

Unbiased estimation of the common mean of a multivariate normal distribution

The problem of unbiased estimation of the common mean of a multivariate normal population is considered. An unbiased estimator is proposed which has a smaller variance than the usual estimator over a large part of the parameter space.     

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     

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.     

On the mvue for the distribution in the multivariate normal regression

In the multivariate normal regression setting, the estimability of a distribution is studied generalizing earlier results for the univariate case. The MVUE of an estimable distribution is obtained.     

Asymptotic distributions of the sphericity test in a complex multivariate normal distribution

In this paper, asymptotic expansions of the null and non-null distributions of the sphericity test criterion in the case of a complex multivariate normal distribution are obtained for the first time in terms of beta distributions. In the null case, it is found that the accuracy of the approximation by taking the first term alone in the asymptotic series is sufficient for practical purposes. In fact for p - 2. the asymptotic expansion reduces to the first term which is also the exact distribution in this case. Applications of the results to the area of inferences on multivariate time series are also given.     

On the exact joint distribution of a linear combination of order statistics and their concomitants in an exchangeable multivariate normal distribution

In this work we prove that for an exchangeable multivariate normal distribution the joint distribution of a linear combination of order statistics and a linear combination of their concomitants together with an auxiliary variable is skew normal. We also investigate some special cases, thus extending the results of Olkin and Viana (J Am Stat Assoc 90:1373–1379, 1995), Loperfido (Test 17:370–380, 2008a) and Sheikhi and Jamalizadeh (Paper 52:885–892, 2011).     

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.     

On a multivariate Lindley distribution

This article proposes a multivariate extension of the generalized Lindley distribution introduced by Abouammoh et al. (2015 Abouammoh, A.M., Alshangiti Arwa, M., Ragab, I.E. (2015). Reliability estimation of the generalized Lindley distributions. J. Stat. Comput. Simul. 85(17):3662–3678.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). The proposed model is based on a probabilistic construction in which several Lindley random variables are connected by a common shock. Many statistical properties of this new distribution are explored. In particular, explicit forms of product moments, moment-generating function, and conditional moments are derived. Explicit expressions of moment-based estimators of the underlying parameters of the proposed model are established. Estimation using the maximum likelihood procedure is also investigated. Simulations attesting to the quality of the proposed estimators are presented. Application of the proposed model in reliability is also discussed.     

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.     

Prior viability assessment for Bayesian analysis

We address the problem of determining whether the cost of a proposed Bayesian analysis is likely to be justified by the potential benefit. A method is described for identifying the likely order of magnitude benefits from the analysis, and this approach is applied to an example concerning trading on the sugar market.