Thank you for your response!
I’ve also referenced your solution on this topic.
I have concatenated fMRI using Nilearn’s concat images:
concated_clean_fmri = nimg.concat_imgs([clean_rest_fmri, clean_session_fmri], auto_resample=True)
And concatenated design matrices as below:
And I computed contrast with the following contrast vector: (I wanted to do ‘session - rest’ contrast)
Is my approach correct? I am a bit worried about the fMRI concatenation step