Predictive modeling involving brain morphological features and other covariates is of paramount interest in such heterogeneous mental disorders as PTSD. We
propose a comprehensive shape analysis framework representing brain substructures, such as the hippocampus,
amygdala, and putamen, as parameterized surfaces and
quantifying their shape differences using an elastic shape
metric. Under this metric, we compute shape summaries
(mean, covariance, PCA) of subcortical data and represent individual shapes by their principal scores under a shape PCA basis. These representations are rich enough to allow visualizations of full 3D structures and help understand localized changes. Subsequently, we use these PCs, the auxiliary exposure variables, and their interactions for regression modeling and prediction. We apply our method
to data from the Grady Trauma Project (GTP), where the
goal is to predict clinical measures of PTSD. The framework seamlessly integrates accurate morphological features and
other clinical covariates to yield superior predictive performance when modeling PTSD outcomes. This approach
reveals considerably greater predictive power under the
elastic shape analysis than the current approaches and
helps identify local deformations in brain shapes associated with PTSD severity.