After one-way analysis of variance (ANOVA) for response height (p = 0.0469), time-to-peak (p = 0.0026), and FWHM (p = 0.0268) of region-specific HRF separately, we found response height as well as time to peak and FWHM were significantly different across the layers (p < 0.05 FDR corrected) for left orbitofrontal cortex. In addition, the three parameters were distinct across layers for left primary somatosensory cortex as well (p = 0.0035 for response height, p = 0.0014 for time-to-peak, and p = 0.0312 for FWHM). Furthermore, multiple comparison of means in one-way ANOVA was employed for each parameter. Here, for the left orbitofrontal cortex, we found that the response height and FWHM of layer I was significantly (p < 0.05 corrected) larger than for layer VI, and time to peak of layers I, II, and III was significantly (p < 0.05 corrected) larger than layer VI. For the primary somatosensory cortex, we found that the response height and time to peak of layer I were significantly (p < 0.05 corrected) larger than for layers III–VI, and at the same time,
