Hello.I have changed the code as you showed ,but there was still a problem(Error in FUN(array(newX[, i], d.call, dn.call), ...) :
object 'Mask' not found).and as you can see the data can be loaded done,
***** Summary information of data structure *****
2 subjects : chen duan
20 response values
5 levels for factor taste : bitter salty sour sweet umami
2 levels for factor session : one two
20 levels for factor Tstat : chen_11_21_final_REML_MNI.nii[1] chen_11_21_final_REML_MNI.nii[10] chen_11_21_final_REML_MNI.nii[13] chen_11_21_final_REML_MNI.nii[4] chen_11_21_final_REML_MNI.nii[7] chen_11_26_final_REML_MNI.nii[1] chen_11_26_final_REML_MNI.nii[10] chen_11_26_final_REML_MNI.nii[13] chen_11_26_final_REML_MNI.nii[4] chen_11_26_final_REML_MNI.nii[7] duan_01_01_final_REML_MNI.nii[1] duan_01_01_final_REML_MNI.nii[10] duan_01_01_final_REML_MNI.nii[13] duan_01_01_final_REML_MNI.nii[4] duan_01_01_final_REML_MNI.nii[7] duan_01_02_final_REML_MNI.nii[1] duan_01_02_final_REML_MNI.nii[10] duan_01_02_final_REML_MNI.nii[13] duan_01_02_final_REML_MNI.nii[4] duan_01_02_final_REML_MNI.nii[7]

Contingency tables of subject distributions among the categorical variables:

***** End of data structure information *****
++++++++++++++++++++++++++++++++++++++++++++++++++++

Reading input files now...

Reading input files for effect estimates: Done!

Reading input files for tStat: Done!

Error in FUN(array(newX[, i], d.call, dn.call), ...) :
object 'Mask' not found
Calls: array -> apply -> FUN
Execution halted