A Matrix Variate Closed Skew-Normal Distribution with Applications to Stochastic Frontier Analysis

In this article, we introduce the matrix extension of the closed skew-normal distribution and give two constructions for it: a marginal one and another based on hidden truncation. Important basic properties of the distribution are presented such as its closure under linear transformation and moment generating function. We also give distributional results for quadratic forms involving random matrices distributed according to two particular cases of it. Using an additive construction, we derive a submodel which can be employed to describe the compound error structure of a very general multivariate stochastic frontier model. Finally, we consider the skew-elliptical extension of the proposed distribution.     

The Use of Noncentral F Approximations for Calculation of Power and Sample Size

Approximations to the noncentral F distribution yield surprisingly accurate results for power and sample size problems arising from linear hypotheses about normal random variables. The approximations are easy to use with a desk (or hand-held) calculator that computes cumulative F probabilities. These approximations are particularly advantageous for testing the hypothesis that differences among the means are small against the alternative that the differences are large.     

On minimaxity of the normal precision matrix estimator of Krishnamoorthy and Gupta

We consider the orthogonally invariant estimation problem of the inverse of the scale matrix of Wishart distribution using Stein's loss (entropy loss). In this problem Krishnamoorthy and Gupta [2] Krishnamoorthy, K. and Gupta, A. K. (1989). Improved minimax estimation of a normal precision matrix. Canad. J. Statist., 17: 91–102. [Crossref], [Web of Science ®] , [Google Scholar] proposed an estimator and showed its good performance in a Monte Carlo simulation. They conjectured their estimator is minimax. Perron [3] Perron, F. (1997). On a conjecture of Krishnamoorthy and Gupta. J. Multivariate Anal., 62: 110–120.  [Google Scholar] proved its minimaxity for p?=?2. In this paper we prove it for p?=?3 by using a new method.     

Testing the Equality of the Noncentrality Parameters of Two Noncentral t-Distributions with Identical Degrees of Freedom

The problem of testing the equality of the noncentrality parameters of two noncentral t-distributions with identical degrees of freedom is considered, which arises from the comparison of two signal-to-noise ratios for simple linear regression models. A test procedure is derived that is guaranteed to maintain Type I error while having only minimal amounts of conservativeness, and comparisons are made with several other approaches to this problem based on variance stabilizing transformations. The new procedure derived in this article is shown to have good properties and will be useful for practitioners.     

基于贝叶斯Wishart波动模型的原油市场与股市动态相依性研究

针对时变相关系数矩阵在多变量随机波动模型的估计问题,构建了贝叶斯动态相关Wishart波动模型。在CC-MSV模型的基础上,设置精度矩阵服从Wishart分布,使得模型的相关系数矩阵具有时变特征。通过模型的统计结构分析,选择参数先验分布,设计相应的Gibbs-MTM-ARMS混合算法,据此估计模型参数;并利用上证综合指数、标普500指数与原油期货价格数据进行实证分析。研究结果表明：模型能够有效地刻画原油市场与股票市场的动态相依性;金融危机期间,股票市场与原油市场的相关性较强,并且难以判断正负方向;金融危机后,中国股票市场与原油市场呈现极微弱的相关性,而美国股票市场与原油市场的正相关性较为明显。     

Computing and approximating multivariate chi-square probabilities

We consider computational methods for evaluating and approximating multivariate chi-square probabilities in cases where the pertaining correlation matrix or blocks thereof have a low factorial representation. To this end, techniques from matrix factorization and probability theory are applied. We outline a variety of statistical applications of multivariate chi-square distributions and provide a system of MATLAB programs implementing the proposed algorithms. Computer simulations demonstrate the accuracy and the computational efficiency of our methods in comparison with Monte Carlo approximations, and a real data example from statistical genetics illustrates their usage in practice.     

A spatial Bayesian semiparametric mixture model for positive definite matrices with applications in diffusion tensor imaging

Studies on diffusion tensor imaging (DTI) quantify the diffusion of water molecules in a brain voxel using an estimated 3 × 3 symmetric positive definite (p.d.) diffusion tensor matrix. Due to the challenges associated with modelling matrix‐variate responses, the voxel‐level DTI data are usually summarized by univariate quantities, such as fractional anisotropy. This approach leads to evident loss of information. Furthermore, DTI analyses often ignore the spatial association among neighbouring voxels, leading to imprecise estimates. Although the spatial modelling literature is rich, modelling spatially dependent p.d. matrices is challenging. To mitigate these issues, we propose a matrix‐variate Bayesian semiparametric mixture model, where the p.d. matrices are distributed as a mixture of inverse Wishart distributions, with the spatial dependence captured by a Markov model for the mixture component labels. Related Bayesian computing is facilitated by conjugacy results and use of the double Metropolis–Hastings algorithm. Our simulation study shows that the proposed method is more powerful than competing non‐spatial methods. We also apply our method to investigate the effect of cocaine use on brain microstructure. By extending spatial statistics to matrix‐variate data, we contribute to providing a novel and computationally tractable inferential tool for DTI analysis.     

Spatial prediction and temporal backcasting for environmental fields having monotone data patterns

The authors develop a methodology for predicting unobserved values in a conditionally lognormal random spatial field like those commonly encountered in environmental risk analysis. These unobserved values are of two types. The first come from spatial locations where the field has never been monitored, the second, from currently monitored sites which have been only recently installed. Thus the monitoring data exhibit a monotone pattern, resembling a staircase whose highest step comes from the oldest monitoring sites. The authors propose a hierarchical Bayesian approach using the lognormal sampling distribution, in conjunction with a conjugate generalized Wishart distribution. This prior distribution allows different degrees of freedom to be fitted for individual steps, taking into account the differential amounts of information available from sites at the different steps in the staircase. The resulting hierarchical model is a predictive distribution for the unobserved values of the field. The method is demonstrated by application to the ambient ozone field for the southwestern region of British Columbia.