Abstract
We discuss an efficient implementation of the iterative proportional scaling procedure in the multivariate Gaussian graphical models. We show that the computational cost can be reduced by localization of the update procedure in each iterative step by using the structure of a decomposable model obtained by triangulation of the graph associated with the model. Some numerical experiments demonstrate the competitive performance of the proposed algorithm.
Acknowledgment
The authors are grateful to two anonymous referees for constructive comments and suggestions, which have led to improvements in the presentation of the article.