Journal
JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS
Volume 26, Issue 3, Pages 709-724Publisher
AMER STATISTICAL ASSOC
DOI: 10.1080/10618600.2017.1279548
Keywords
Generalized linear array models; Multidimensional smoothing; Penalized estimation; Proximal gradient algorithm
Categories
Funding
- Danish Cancer Society
- Danish Strategic Research Council/Innovation Fund Denmark
- Villum Fonden [00013358] Funding Source: researchfish
Ask authors/readers for more resources
Large-scale generalized linear array models (GLAMs) can be challenging to fit. Computation and storage of its tensor product design matrix can be impossible due to time and memory constraints, and previously considered design matrix free algorithms do not scale well with the dimension of the parameter vector. A new design matrix free algorithm is proposed for computing the penalized maximum likelihood estimate for GLAMs, which, in particular, handles nondifferentiable penalty functions. The proposed algorithm is implemented and available via the R package glamlasso. It combines several ideas-previously considered separately-to obtain sparse estimates while at the same time efficiently exploiting the GLAM structure. In this article, the convergence of the algorithm is treated and the performance of its implementation is investigated and compared to that of glmnet on simulated as well as real data. It is shown that the computation time for glamlasso scales favorably with the size of the problem when compared to glmnet. Supplementary materials, in the form of R code, data and visualizations of results, are available online.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available