fNIRS data acquisition and analyses
Hemodynamic response signals were acquired using the NIRSport2 System (NIRx Medical technologies LLC, Berlin, Germany) operating at two wavelengths (760 and 850 nm) with a sampling frequency of 6.78 Hz. In line with previous studies (The fNIRS raw data were analyzed using the NIRS-KIT software which can be used for both resting-state and task-based fNIRS data analyses (We also employed the representational similarity analysis (RSA) technique to investigate whether the received treatment of one participant (test subject, leave-one-out method) could be correctly classified based on their similarity of activation patterns to the mean pattern of subjects receiving the same treatment (group models except for the test subject each time; 