|
|  |
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: inverting large, sparse matrices
"Pierre Azoulay" <[email protected]> asks:
I am looking into the best strategy (and software package) to invert a
sparse positive definite 25,000 by 25,000 matrix. Is that something
feasible at all? I really have no idea. Does one need a supercomputer
to do this? Or would Stata/MP and 16 gigs of RAM do on a server do?
Would stata and matlab serve me equally well? Are there references to
look at? Are cholinv() and invsym() the right mata functions to look
at? Will these functions make use of the fact that the matrix is
sparse?
If I am not mistaken, neither -cholinv- nor -invsum- take into account the
fact that the matrix is sparse. I think I understood from the literature
that for large sparse system a preconditioned conjugate gradient algorithm
(see for example Abowd, Creecy and Kramarz 2002) is a good solution, but I
know of no implementation in Stata.
As far as I remember, the problem treated in Abowd, Creecy and Kramarz 2002
involves a higher dimension than 25,000 by 25,000.
Abowd J., Creecy R. and Kramarz F. (2002): "Computing Person and Firm
Effects Using Linked Longitudinal Employer-Employee Data', Technical Paper
No. TP-2002-06, U.S. Census Bureau.
Thomas
-------------------------------------------------------------------------
Thomas Cornelissen
Institute of Empirical Economic Research
Leibniz Universit�t Hannover, Germany
[email protected]
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
