To assess the relationship between functional connectivity of NAcc and delayed-reward discounting gradient, we performed a voxel-wise correlation of the Fisher Z transformed correlation coefficients with ln(k) using an in-house software maintained by Neuro-Imaging Laboratory (Washington University). Monte Carlo simulation was implemented to account for multiple comparisons (threshold of 53 contiguous voxels with a Z-value > 2.25).
