Modeling Covariate Effects in Group Independent Component Analysis with Applications to Functional Magnetic Resonance Imaging