Li-Dan Kuang1, Qiu-Hua Lin2, Xiao-Feng Gong1, Fengyu Cong3, Jing Sui4, Vince D Calhoun5. 1. School of Information and Communication Engineering, Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116024, China. 2. School of Information and Communication Engineering, Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116024, China. Electronic address: qhlin@dlut.edu.cn. 3. Department of Biomedical Engineering, Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116024, China; Department of Mathematical Information Technology, University of Jyvaskyla, Finland. 4. Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, China. 5. The Mind Research Network, Albuquerque, NM 87106, USA; Department of Electrical and Computer Engineering, University of New Mexico, Albuquerque, NM 87131, USA.
Abstract
BACKGROUND: Canonical polyadic decomposition (CPD) may face a local optimal problem when analyzing multi-subject fMRI data with inter-subject variability. Beckmann and Smith proposed a tensor PICA approach that incorporated an independence constraint to the spatial modality by combining CPD with ICA, and alleviated the problem of inter-subject spatial map (SM) variability. NEW METHOD: This study extends tensor PICA to incorporate additional inter-subject time course (TC) variability and to connect CPD and ICA in a new way. Assuming multiple subjects share common TCs but with different time delays, we accommodate subject-dependent TC delays into the CP model based on the idea of shift-invariant CP (SCP). We use ICA as an initialization step to provide the aggregating mixing matrix for shift-invariant CPD to estimate shared TCs with subject-dependent delays and intensities. We then estimate shared SMs using a least-squares fit post shift-invariant CPD. RESULTS: Using simulated fMRI data as well as actual fMRI data we demonstrate that the proposed approach improves the estimates of the shared SMs and TCs, and the subject-dependent TC delays and intensities. The default mode component illustrates larger TC delays than the task-related component. COMPARISON WITH EXISTING METHOD(S): The proposed approach shows improvements over tensor PICA in particular when TC delays are large, and also outperforms SCP with SM orthogonality constraint and SCP with ICA-based SM initialization. CONCLUSIONS: TCs with subject-dependent delays conform to the true situation of multi-subject fMRI data. The proposed approach is suitable for decomposing multi-subject fMRI data with large inter-subject temporal and spatial variability.
BACKGROUND: Canonical polyadic decomposition (CPD) may face a local optimal problem when analyzing multi-subject fMRI data with inter-subject variability. Beckmann and Smith proposed a tensor PICA approach that incorporated an independence constraint to the spatial modality by combining CPD with ICA, and alleviated the problem of inter-subject spatial map (SM) variability. NEW METHOD: This study extends tensor PICA to incorporate additional inter-subject time course (TC) variability and to connect CPD and ICA in a new way. Assuming multiple subjects share common TCs but with different time delays, we accommodate subject-dependent TC delays into the CP model based on the idea of shift-invariant CP (SCP). We use ICA as an initialization step to provide the aggregating mixing matrix for shift-invariant CPD to estimate shared TCs with subject-dependent delays and intensities. We then estimate shared SMs using a least-squares fit post shift-invariant CPD. RESULTS: Using simulated fMRI data as well as actual fMRI data we demonstrate that the proposed approach improves the estimates of the shared SMs and TCs, and the subject-dependent TC delays and intensities. The default mode component illustrates larger TC delays than the task-related component. COMPARISON WITH EXISTING METHOD(S): The proposed approach shows improvements over tensor PICA in particular when TC delays are large, and also outperforms SCP with SM orthogonality constraint and SCP with ICA-based SM initialization. CONCLUSIONS:TCs with subject-dependent delays conform to the true situation of multi-subject fMRI data. The proposed approach is suitable for decomposing multi-subject fMRI data with large inter-subject temporal and spatial variability.