Abstract
OBJECTIVE: Identification of networks from resting brain signals is an important step in understanding the dynamics of spontaneous brain activity. We approach this problem using a tensor-based model. METHODS: We develop a rank-recursive scalable and robust sequential canonical polyadic decomposition (SRSCPD) framework to decompose a tensor into several rank-1 components. Robustness and scalability are achieved using a warm start for each rank based on the results from the previous rank. RESULTS: In simulations we show that SRSCPD consistently outperforms the multi-start alternating least square (ALS) algorithm over a range of ranks and signal-to-noise ratios (SNRs), with lower computation cost. When applying SRSCPD to resting in-vivo stereotactic EEG (SEEG) data from two subjects with epilepsy, we found components corresponding to default mode and motor networks in both subjects. These components were also highly consistent within subject between two sessions recorded several hours apart. Similar components were not obtained using the conventional ALS algorithm. CONCLUSION: Consistent brain networks and their dynamic behaviors were identified from resting SEEG data using SRSCPD. SIGNIFICANCE: SRSCPD is scalable to large datasets and therefore a promising tool for identification of brain networks in long recordings from single subjects.