Abstract
We propose a method for selecting edges in undirected Gaussian graphical models. Our algorithm takes after our previous work, an extension of Least Angle Regression (LARS), and it is based on the information geometry of dually flat spaces. Non-diagonal elements of the inverse of the covariance matrix, the concentration matrix, play an important role in edge selection. Our iterative method estimates these elements and selects covariance models simultaneously. A sequence of pairs of estimates of the concentration matrix and an independence graph is generated, whose length is the same as the number of non-diagonal elements of the matrix. In our algorithm, the next estimate of the graph is the nearest graph to the latest estimate of the concentration matrix. The next estimate of the concentration matrix is not just the projection of the latest estimate, and it is shrunk to the origin. We describe the algorithm and show results for some datasets. Furthermore, we give some remarks on model identification and prediction.
Similar content being viewed by others
References
Amari, S.: Differential-Geometrical Methods in Statistics. Lecture Notes in Statistics, vol. 28. Springer, Berlin (1985)
Amari, S., Nagaoka, H.: Methods of Information Geometry. Translations of Mathematical Monographs, vol. 191. Oxford University Press, London (2000)
Dempster, A.P.: Covariance selection. Biometrics 28, 157–175 (1972)
Efron, B., Hastie, T., Johnstone, I., Tibshirani, R.: Least angle regression (with discussion). Ann. Stat. 32, 407–499 (2004)
Friedman, J., Hastie, T., Höfling, H., Tibshirani, R.: Pathwise coordinate optimization. Ann. Appl. Stat. 1, 302–332 (2007)
Hastie, T., Tibshirani, R., Friedman, J.: The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd edn. Springer, New York (2009)
Hirose, Y., Komaki, F.: An extension of least angle regression based on the information geometry of dually flat spaces. J. Comput. Graph. Stat. 19, 1007–1023 (2010)
Kass, R., Vos, P.: Geometrical Foundations of Asymptotic Inference. Wiley, New York (1997)
Lauritzen, S.L.: Graphical Models. Clarendon Press, Oxford (1996)
R Development Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna (2012). ISBN 3-900051-07-0. http://www.R-project.org/.
Speed, T.P., Kiiveri, H.T.: Gaussian Markov distributions over finite graphs. Ann. Stat. 14, 138–150 (1986)
Tibshirani, R.: Regression shrinkage and selection via the lasso. J. R. Stat. Soc. B 58, 267–288 (1996)
Wermuth, N., Scheidt, E.: Fitting a covariance selection model to a matrix. Appl. Stat. 26, 88–92 (1977)
Whittaker, J.: Graphical Models in Applied Multivariate Statistics. Wiley, Chichester (1990)
Yuan, M., Lin, Y.: Model selection and estimation in the Gaussian graphical model. Biometrika 94, 19–35 (2007)
Acknowledgements
This research was supported in part by Grant-in-Aid for Scientific Research (23300104, 23650144), by Global COE Program “The research and training center for new development in mathematics,” MEXT, Japan, and by the Aihara Project, the FIRST program from JSPS, initiated by CSTP.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Hirose, Y., Komaki, F. Edge selection based on the geometry of dually flat spaces for Gaussian graphical models. Stat Comput 23, 793–800 (2013). https://doi.org/10.1007/s11222-012-9347-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11222-012-9347-3