  We present a non-parametric joint estimation method for fMRI task activation
values and the hemodynamic response function (HRF). The HRF is modeled as a
Gaussian process, making continuous evaluation possible for jittered paradigms
and providing a variance estimate at each point.
