spam: A Sparse Matrix R Package with Emphasis on MCMC Methods for Gaussian Markov Random Fields
Abstractspam is an R package for sparse matrix algebra with emphasis on a Cholesky factorization of sparse positive definite matrices. The implemantation of spam is based on the competing philosophical maxims to be competitively fast compared to existing tools and to be easy to use, modify and extend. The first is addressed by using fast Fortran routines and the second by assuring S3 and S4 compatibility. One of the features of spam is to exploit the algorithmic steps of the Cholesky factorization and hence to perform only a fraction of the workload when factorizing matrices with the same sparsity structure. Simulations show that exploiting this break-down of the factorization results in a speed-up of about a factor 5 and memory savings of about a factor 10 for large matrices and slightly smaller factors for huge matrices. The article is motivated with Markov chain Monte Carlo methods for Gaussian Markov random fields, but many other statistical applications are mentioned that profit from an efficient Cholesky factorization as well.
Download InfoIf you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.
Bibliographic InfoArticle provided by American Statistical Association in its journal Journal of Statistical Software.
Volume (Year): 36 ()
Issue (Month): i10 ()
Contact details of provider:
Web page: http://www.jstatsoft.org/
You can help add them by filling out this form.
CitEc Project, subscribe to its RSS feed for this item.
- Toshihiro Hirano & Yoshihiro Yajima, 2011. "Covariance Tapering for Prediction of Large Spatial Data Sets in Transformed Random Fields," CIRJE F-Series CIRJE-F-823, CIRJE, Faculty of Economics, University of Tokyo.
- Millo, Giovanni, 2014. "Maximum likelihood estimation of spatially and serially correlated panels with random effects," Computational Statistics & Data Analysis, Elsevier, vol. 71(C), pages 914-933.
- Toshihiro Hirano & Yoshihiro Yajima, 2013. "Covariance tapering for prediction of large spatial data sets in transformed random fields," Annals of the Institute of Statistical Mathematics, Springer, vol. 65(5), pages 913-939, October.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Christopher F. Baum).
If references are entirely missing, you can add them using this form.