Dear Donald,
Yes, I have already considered a similar approach = specifiy the models with the volumes included, take the time course from xX.X (which I want need for building some time modulators - I'm aware of the script from Wager lab, but I prefer it this way), delete the first few volumes and then add this as text file for the final specification/estimation together with the other regressors (and I have dummy regressors for fast head motion anyway).
Then again, as there are already these 32 dummy time bins in the code it should be simple to adjust it. It's just that I don't understand why the original file goes with 32 dummy bins (and accordingly +33 to have non-zero value for the first bin after trial onset) for ton but +128 in
sf = sparse((k*T + 128),size(u,2));
I understand what sparse as such does, but I don't understand why there's a need to add 128 and/or why this exact value was chosen.
Best
Helmut
|