Abstract:Summary Many studies collect functional data from multiple subjects that have both multilevel and multivariate structures. An example of such data comes from popular neuroscience experiments where participants’ brain activity is recorded using modalities such as electroencephalography and summarized as power within multiple time-varying frequency bands within multiple electrodes, or brain regions. Summarizing the joint variation across multiple frequency bands for both whole-brain variability between subjects, as well as location–variation within subjects, can help to explain neural reactions to stimuli. This article introduces a novel approach to conducting interpretable principal components analysis on multilevel multivariate functional data that decomposes total variation into subject-level and replicate-within-subject-level (i.e., electrode-level) variation and provides interpretable components that can be both sparse among variates (e.g., frequency bands) and have localized support over time within each frequency band. Smoothness is achieved through a roughness penalty, while sparsity and localization of components are achieved by solving an innovative rank-one based convex optimization problem with block Frobenius and matrix $L_1$-norm-based penalties. The method is used to analyze data from a study to better understand reactions to emotional information in individuals with histories of trauma and the symptom of dissociation, revealing new neurophysiological insights into how subject- and electrode-level brain activity are associated with these phenomena. Supplementary materials for this article are available online.

Localized Sparse Principal Component Analysis of Multivariate Time Series in Frequency Domain

Principal Component Analysis using Frequency Components of Multivariate Time Series

Interpretable Principal Components Analysis for Multilevel Multivariate Functional Data, with Application to EEG Experiments

Interpretable principal component analysis for multilevel multivariate functional data

Sparse Functional Principal Component Analysis in High Dimensions

Sparse Principal Component Analysis With Preserved Sparsity Pattern

Sparse and Functional Principal Components Analysis

Sparse principal component analysis for high-dimensional stationary time series

Classification of High-dimensional Time Series in Spectral Domain using Explainable Features

Sparse Principal Component Analysis via Variable Projection

Multilevel Functional Principal Component Analysis for High-Dimensional Data

Structured Principal Component Analysis Model With Variable Correlation Constraint

Adaptive Wavelet Domain Principal Component Analysis for Nonstationary Time Series

Multivariate Functional Principal Component Analysis for Data Observed on Different (Dimensional) Domains

Adaptive Bayesian Spectral Analysis of High-Dimensional Nonstationary Time Series

Modeling Spectral Properties in Stationary Processes of Varying Dimensions with Applications to Brain Local Field Potential Signals

Sparse Principal Component Analysis for High Dimensional Vector Autoregressive Models

Frequency Band Analysis of Nonstationary Multivariate Time Series

Sparse Principal Component Analysis

Multilevel hybrid principal components analysis for region-referenced functional electroencephalography data