Dear Dr. Stephan,
We are wondering about 2 points regarding your nice paper
Comparing hemodynamic models with DCM, NeuroImage 38, 2007.
1) In SPM8, the data fed to spm_gx_hdm (or spm_gx_dcm or spm_gx_fmri)
is the Bold signal S (called Y in the code), but we believe it should
be Delta S/S0.
Is this a mistake in the implementation of the measure model in SPM8?
Namely that the calculation of Delta S/S0 (with an estimation of what
S0 is) is missing
from the code?
2) In Obata et al. 2004, NeuroImage 21, 144, in the Appendix, equation
A12: there is
a factor of v following k2, but this factor is missing from the SPM8
also in various equations in your article (equation 2 and 13). This
makes no difference
to the linear version, but the nonlinear versions appear to be
different. Please could
you let us know if this omission was intentional or not?
Many thanks for your help,
Philippe Pouliot, Ph.D.
Montreal Heart Institute
Ecole Polytechnique Montreal