Locally regularized spatiotemporal modeling and model comparison for functional MRI

In this work we treat fMRI data analysis as a spatiotemporal system identification problem and address issues of model formulation, estimation, and model comparison. We present a new model that includes a physiologically based hemodynamic response and an empirically derived low-frequency noise model. We introduce an estimation method employing spatial regularization that improves the precision of spatially varying noise estimates. We call the algorithm locally regularized spatiotemporal (LRST) modeling. We develop a new model selection criterion and compare our model to the SPM-GLM method. Our findings suggest that our method offers a better approach to identifying appropriate statistical models for fMRI studies.