Hi

sorry for the late reply. You can extract some shape parameters (width, height, time to peak) and use them as (individual or multiple) voxel-wise dependent variables, or use the full shape of the HRF, for example using AFNI's 3dMVM or 3dMSS.