[Rd] Cholesky update/downdate

From: Yves Deville <deville.yves_at_alpestat.com>
Date: Thu, 29 Dec 2011 15:05:28 +0100

Dear R-devel members,

I am looking for a fast Cholesky update/downdate. The matrix A being symmetric positive definite (n, n) and factorized as A = L %*% t(L), the goal is to factor the new matrix A +- C %*% t(C) where C is (n, r). For instance, C is 1-column when adding/removing an observation in a linear regression. Of special interest is the case where A is sparse.  

Looking at the 'Matrix' package (help and source code), it seems that the CHOLMOD library shipped with 'Matrix' allows this, but is not (yet?) interfaced in 'Matrix', where the 'update' method for Cholesky decomposition objects seems limited to a new matrix A + m*I with a scalar (diagonal) modification.

If this is true: are there plans to implement such up/downdates?  



Yves Deville, statistical consultant, France.

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Thu 29 Dec 2011 - 14:09:15 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 29 Dec 2011 - 23:00:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive