Fitting large-scale structured additive regression models using Krylov subspace methods

Paul Schmidt, Mark Mühlau and Volker Schmid

Computational Statistics & Data Analysis, 2017, vol. 105, issue C, 59-75

Abstract: Fitting regression models can be challenging when regression coefficients are high-dimensional. Especially when large spatial or temporal effects need to be taken into account the limits of computational capacities of normal working stations are reached quickly. The analysis of images with several million pixels, where each pixel value can be seen as an observation on a new spatial location, represent such a situation. A Markov chain Monte Carlo (MCMC) framework for the applied statistician is presented that allows to fit models with millions of parameters with only low to moderate computational requirements. The method combines a modified sampling scheme with novel accomplishments in iterative methods for sparse linear systems. This way a solution is given that eliminates potential computational burdens such as calculating the log-determinant of massive precision matrices and sampling from high-dimensional Gaussian distributions. In an extensive simulation study with models of moderate size it is shown that this approach gives results that are in perfect agreement with state-of-the-art methods for fitting structured additive regression models. Furthermore, the method is applied to two real world examples from the field of medical imaging.

Keywords: Markov chain Monte Carlo; Krylov subspace methods; Lanczos algorithm; Structured additive regression; Gaussian Markov random field; Image analysis (search for similar items in EconPapers)
Date: 2017
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0167947316301700
Full text for ScienceDirect subscribers only.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:csdana:v:105:y:2017:i:c:p:59-75

DOI: 10.1016/j.csda.2016.07.006

Access Statistics for this article

Computational Statistics & Data Analysis is currently edited by S.P. Azen

More articles in Computational Statistics & Data Analysis from Elsevier
Bibliographic data for series maintained by Catherine Liu ().