Skip Navigation


Biometrika Advance Access originally published online on October 9, 2009
Biometrika 2009 96(4):821-834; doi:10.1093/biomet/asp049
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Wang, H.
Right arrow Articles by West, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2009 Biometrika Trust

Article

Bayesian analysis of matrix normal graphical models

Hao Wang and Mike West

Department of Statistical Science, Duke University, Durham, North Carolina 27708, U.S.A. hao{at}stat.duke.edu mike{at}stat.duke.edu

Received for publication 1 February 2008. Revision received 1 February 2009.

We present Bayesian analyses of matrix-variate normal data with conditional independencies induced by graphical model structuring of the characterizing covariance matrix parameters. This framework of matrix normal graphical models includes prior specifications, posterior computation using Markov chain Monte Carlo methods, evaluation of graphical model uncertainty and model structure search. Extensions to matrix-variate time series embed matrix normal graphs in dynamic models. Examples highlight questions of graphical model uncertainty, search and comparison in matrix data contexts. These models may be applied in a number of areas of multivariate analysis, time series and also spatial modelling.

Key Words: Gaussian graphical model • Graphical model search • Hyper-inverse Wishart distribution • Marginal likelihood • Matrix normal model • Matrix-variate dynamic graphical model • Parameter expansion



References

    Besag J. A candidate’s formula: a curious result in Bayesian prediction. Biometrika (1989) 76:183.[Abstract/Free Full Text]

    Carvalho C. M., Massam H., West M. Simulation of hyper-inverse Wishart distributions in graphical models. Biometrika (2007) 94:647–59.[Abstract/Free Full Text]

    Carvalho C. M., West M. Dynamic matrix-variate graphical models. Bayesian Anal. (2007a) 2:69–98.[CrossRef]

    Carvalho C. M., West M. Dynamic matrix-variate graphical models—a synopsis. Bayesian Statistics, VIII—Bernardo J. M., Bayarri M. J., Berger J. O., Dawid A. P., Heckerman D., Smith A. F. M., West M., eds. (2007b) Oxford: Oxford University Press. 585–90.

    Chib S. Marginal likelihood from the Gibbs output. J. Am. Statist. Assoc. (1995) 90:1313–21.[CrossRef][Web of Science]

    Dawid A. P. Some matrix-variate distribution theory: notational considerations and a Bayesian application. Biometrika (1981) 68:265–74.[Abstract/Free Full Text]

    Dawid A. P., Lauritzen S. L. Hyper-Markov laws in the statistical analysis of decomposable graphical models. Ann. Statist. (1993) 21:1272–317.[CrossRef]

    Dobra A., Jones B., Hans C., Nevins J., West M. Sparse graphical models for exploring gene expression data. J. Mult. Anal. (2004) 90:196–212.[CrossRef]

    Dutilleul P. The MLE algorithm for the matrix normal distribution. J. Statist. Comp. Simul. (1999) 64:105–23.[CrossRef]

    Finn J. D. A General Model for Multivariate Analysis (1974) New York: Holt, Rinehart and Winston.

    Galecki A. General class of covariance structures for two or more repeated factors in longitudinal data analysis. Commun. Statist. (1994) A 23:3105–19.[CrossRef]

    Gelman A. Parameterization and Bayesian modeling. J. Am. Statist. Assoc. (2004) 99:537–45.[CrossRef][Web of Science]

    Gelman A. Prior distributions for variance parameters in hierarchical models. Bayesian Anal. (2006) 3:515–34.

    Giudici P. Learning in graphical Gaussian models. Bayesian Statistics, 5—Bernado J. M., Berger J. O., Dawid A. P., Smith A. F. M., eds. (1996) Oxford: Oxford University Press. 621–28.

    Giudici P., Green P. J. Decomposable graphical Gaussian model determination. Biometrika (1999) 86:785–801.[Abstract/Free Full Text]

    Gupta A. K., Nagar D. K. Matrix Variate Distributions (2000) Monographs and Surveys in Pure & Applied Mathematics 104. London: Chapman & Hall.

    Hans C., Dobra A., West M. Shotgun stochastic search in regression with many predictors. J. Am. Statist. Assoc. (2007) 102:507–16.[CrossRef][Web of Science]

    Hobert J. P., Marchev D. A theoretical comparison of the data augmentation, marginal augmentation and PXDA algorithms. Ann. Statist. (2008) 2:532–54.

    Huizenga H. M., de Munck J. C., Waldorp L. J., Grasman R. Spatiotemporal EEG/MEG source analysis based on a parametric noisecovariance model. IEEE Trans. Biomed. Eng. (2002) 49:533–39.[CrossRef][Web of Science][Medline]

    Jones B., Carvalho C. M., Dobra A., Hans C., Carter C., West M. Experiments in stochastic computation for high-dimensional graphical models. Statist. Sci. (2005) 20:388–400.[CrossRef]

    Jordan M., Ghahramani Z., Jaakkola T., Saul L. An introduction to variational methods for graphical models. Mach. Learn. (1999) 37:183–233.[CrossRef]

    Lauritzen S. L. Graphical Models (1996) Oxford: Clarendon Press.

    Liu C., Rubin D. B., Wu Y. N. Parameter expansion to accelerate EM: the PX-EM algorithm. Biometrika (1998) 85:755–70.[Abstract/Free Full Text]

    Liu J. S., Wu Y. N. Parameter expansion for data augmentation. J. Am. Statist. Assoc. (1999) 94:1264–74.[CrossRef][Web of Science]

    Mardia K. V., Goodall C. R. Spatial-temporal analysis of multivariate environmental monitoring data. In: Multivariate Environmental Statistics—Patil G. P., Rao C. R., eds. (1993) Amsterdam: Elsevier. 347–85.

    McCulloch R. E., Polson N. G., Rossi P. E. Bayesian analysis of the multinomial probit model with fully identified parameters. J. Economet. (2000) 99:173–93.[CrossRef]

    Mitchell M. W., Genton M. G., Gumpertz M. L. A likelihood ratio test for separability of covariances. J. Mult. Anal. (2006) 97:1025–43.[CrossRef]

    Naik D. N., Rao S. S. Analysis of multivariate repeated measures data with a Kronecker product structured covariance matrix. J. Appl. Statist. (2001) 29:91–105.

    Pole A., West M., Harrison P. J. Applied Bayesian Forecasting and Time Series Analysis (1994) New York: Chapman-Hall.

    Quintana J. M., West M. Multivariate time series analysis: new techniques applied to international exchange rate data. Statistician (1987) 36:275–81.[CrossRef][Web of Science]

    Roy V., Hobert J. P. Convergence rates and asymptotic standard errors for MCMC algorithms for Bayesian probit regression. J. R. Statist. Soc. (2007) B 69:607–23.[CrossRef]

    Theobald D. L., Wuttke D. S. Empirical Bayes hierarchical models for regularizing maximum likelihood estimation in the matrix Gaussian procrustes problem. Proc. Nat. Acad. Sci. (2006) 103:18521–7.[Abstract/Free Full Text]

    West M., Harrison P. J. Bayesian Forecasting and Dynamic Models (1997) 2nd ed. New York: Springer.

    Whittaker J. Graphical Models in Applied Multivariate Statistics (1990) Chichester, UK: John Wiley and Sons.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Wang, H.
Right arrow Articles by West, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?