Thanks ChrisGorgolweski and effigies for the replies. It turns out that we’ll probably just stick with a surface-based analysis after all.
Given that, am I correct in thinking that nipype’s interface to freesurfer’s glm_fit is the only option currently available? about using nipype for a first-level, surface-based glm,? Nistats was listed in Simple first-level GLM in Python, but it seems like it doesn’t yet interface with nipype (and, would it do glms on a gifti image?)