TY - JOUR
T1 - Is first-order vector autoregressive model optimal for fMRI data?
AU - Ting, Chee Ming
AU - Seghouane, Abd Krim
AU - Khalid, Muhammad Usman
AU - Salleh, Sh Hussain
N1 - Copyright:
Copyright 2015 Elsevier B.V., All rights reserved.
PY - 2015/9/21
Y1 - 2015/9/21
N2 - We consider the problem of selecting the optimal orders of vector autoregressive (VAR) models for fMRI data.Many previous studies used model order of one and ignored that it may vary considerably across data sets depending on different data dimensions, subjects, tasks, and experimental designs. In addition, the classical information criteria (IC) used (e.g., the Akaike IC (AIC)) are biased and inappropriate for the high-dimensional fMRI data typically with a small sample size. We examine the mixed results on the optimal VAR orders for fMRI, especially the validity of the order-one hypothesis, by a comprehensive evaluation using different model selection criteria over three typical data types-a resting state, an event-related design, and a block design data set-with varying time series dimensions obtained from distinct functional brain networks. We use a more balanced criterion, Kullback's IC (KIC) based on Kullback's symmetric divergence combining two directed divergences.We also consider the bias-corrected versions (AICc and KICc) to improve VAR model selection in small samples. Simulation results show better small-sample selection performance of the proposed criteria over the classical ones. Both bias-corrected ICs provide more accurate and consistent model order choices than their biased counterparts, which suffer from overfitting, with KICc performing the best. Results on real data show that orders greater than one were selected by all criteria across all data sets for the small to moderate dimensions, particularly from small, specific networks such as the resting-state default mode network and the task-related motor networks, whereas low orders close to one but not necessarily one were chosen for the large dimensions of full-brain networks.
AB - We consider the problem of selecting the optimal orders of vector autoregressive (VAR) models for fMRI data.Many previous studies used model order of one and ignored that it may vary considerably across data sets depending on different data dimensions, subjects, tasks, and experimental designs. In addition, the classical information criteria (IC) used (e.g., the Akaike IC (AIC)) are biased and inappropriate for the high-dimensional fMRI data typically with a small sample size. We examine the mixed results on the optimal VAR orders for fMRI, especially the validity of the order-one hypothesis, by a comprehensive evaluation using different model selection criteria over three typical data types-a resting state, an event-related design, and a block design data set-with varying time series dimensions obtained from distinct functional brain networks. We use a more balanced criterion, Kullback's IC (KIC) based on Kullback's symmetric divergence combining two directed divergences.We also consider the bias-corrected versions (AICc and KICc) to improve VAR model selection in small samples. Simulation results show better small-sample selection performance of the proposed criteria over the classical ones. Both bias-corrected ICs provide more accurate and consistent model order choices than their biased counterparts, which suffer from overfitting, with KICc performing the best. Results on real data show that orders greater than one were selected by all criteria across all data sets for the small to moderate dimensions, particularly from small, specific networks such as the resting-state default mode network and the task-related motor networks, whereas low orders close to one but not necessarily one were chosen for the large dimensions of full-brain networks.
UR - http://www.scopus.com/inward/record.url?scp=84939640163&partnerID=8YFLogxK
U2 - 10.1162/NECO_a_00765
DO - 10.1162/NECO_a_00765
M3 - Article
C2 - 26161816
AN - SCOPUS:84939640163
SN - 0899-7667
VL - 27
SP - 1857
EP - 1871
JO - Neural Computation
JF - Neural Computation
IS - 9
ER -