Fitting very large sparse Gaussian graphical models

详细信息	查看全文 | 推荐本文 |

• 作者：Harri Kiiveri^a ; ^{harri.kiiveri@csiro.au} ; Frank de Hoog^b ; ^{frank.dehoog@csiro.au}
• 关键词：Covariance selection ; Gene networks ; Graphical models ; High dimensional ; Large scale optimisation ; Limited memory quasi-Newton
• 刊名：Computational Statistics and Data Analysis
• 出版年：2012
• 期刊代码：19_01679473
• 类别：cp
• 出版时间：September, 2012
• 卷：56
• 期：9
• 页码：2626-2636
• 文件大小：453 K

摘要

In this paper we consider some methods for the maximum likelihood estimation of sparse Gaussian graphical (covariance selection) models when the number of variables is very large (tens of thousands or more). We present a procedure for determining the pattern of zeros in the model and we discuss the use of limited memory quasi-Newton algorithms and truncated Newton algorithms to fit the model by maximum likelihood. We present efficient ways of computing the gradients and likelihood function values for such models suitable for a desktop computer. For the truncated Newton method we also present an efficient way of computing the action of the Hessian matrix on an arbitrary vector which does not require the computation and storage of the Hessian matrix. The methods are illustrated and compared on simulated data and applied to a real microarray data set.

The limited memory quasi-Newton method is recommended for practical use.