Sparse principal component analysis in medical shape modeling

Abstract

Principal component analysis (PCA) is a widely used tool in medical image analysis for data reduction, model building, and data understanding and exploration. While PCA is a holistic approach where each new variable is a linear combination of all original variables, sparse PCA (SPCA) aims at producing easily interpreted models through sparse loadings, i.e. each new variable is a linear combination of a subset of the original variables. One of the aims of using SPCA is the possible separation of the results into isolated and easily identifiable effects.This article introduces SPCA for shape analysis in medicine. Results for three different data sets are given in relation to standard PCA and sparse PCA by simple thresholding of small loadings. Focus is on a recent algorithm for computing sparse principal components, but a review of other approaches is supplied as well. The SPCA algorithm has been implemented using Matlab and is available for download.The general behavior of the algorithm is investigated, and strengths and weaknesses are discussed. The original report on the SPCA algorithm argues that the ordering of modes is not an issue. We disagree on this point and propose several approaches to establish sensible orderings. A method that orders modes by decreasing variance and maximizes the sum of variances for all modes is presented and investigated in detail.

Selecting signature optical emission spectroscopy variables using sparse principal component analysis

“…LARS-EN is a modified implementation of LARS designed to fit the elastic net framework [8]. In particular, it includes a special algorithm, called soft-thresholding, that is scalable to high dimensional data sets [12].…”

Deformable model fitting with a mixture of local experts

“…Objects that are typically modeled in this way include the human face [3,19] and organs in medical image analysis [22,24]. Numerous representations and fitting strategies have been proposed for these objects, most of which can be categorized based on their representations as being either holistic or patch-based.…”

A nonlinear principal component decomposition

“…Factor analysis is used for forecasting (Stock and Watson (1999)) and for Engel curves construction in demand analysis (Lewbel (1991)). More broadly, applications can be found in many fields of statistics (Loève (1978)) and include medical imaging (Sjöstrand, Stegmann, and Larsen (2006)), data compression (Wallace (1991)) and even search engines (Brin and Page (1998)). …”