Hi,
Did you set this up via the GUI or did you create an fsf (possibly via the GUI) and edit it?
It is easy to make mistakes when editing the fsf files, especially when changing the number of EVs or contrasts. So if you did not set it up via the GUI with this particular number of EVs, contrasts and basis functions, I would do that via the GUI. You could then look at what it creates and see how that is different from the fsf file that you edited.
If you did use the GUI and not edit the fsf file then let us know.
All the best,
Mark
On 25 Jul 2013, at 08:38, Mark Austin <[log in to unmask]> wrote:
> Hi,
>
> We have developed a successful processing pipeline using the command line FEAT and Featquery tools to analysis a simple visual block stimulus. Our endpoint analysis is looking at a trapezoid fit of the activation within a specific ROI and we have realized that the default dual-gamma HRF model is not flexible enough to be sensitive to the variations in activation patterns that we are looking for. We therefore wanted to introduce a more flexible HRF model and it appears the FLOBS is the tool to achieve this.
>
> We have modified our design.fsf file (snippet attached below), selecting the optimal/custom HRF style and using the default basis function file, however, following this change we get the following error in FEAT:
>
> /usr/share/fsl/5.0/bin/feat_model design
> terminate called after throwing an instance of 'NEWMAT::IndexException'
> Aborted
>
> Is this anything that has been seen before? Is there anything we should try? Do you need more information to be able to understand the issue? Any help would be appreciated.
>
> Thanks
>
>
> design.fsf snippet:
>
> # Convolution (EV 1)
> # 0 : None
> # 1 : Gaussian
> # 2 : Gamma
> # 3 : Double-Gamma HRF
> # 4 : Gamma basis functions
> # 5 : Sine basis functions
> # 6 : FIR basis functions
> set fmri(convolve1) 7
>
> # FIR basis functions number (EV 1)
> set fmri(basisfnum1) 3
>
> # Optimal/custom HRF convolution file (EV 1)
> set fmri(bfcustom1) "/usr/share/fsl/5.0/fsl/etc/default_flobs.flobs/hrfbasisfns.txt"
>
> # Orth basis functions wrt each other
> set fmri(basisorth1) 1
>
> # Convolve phase (EV 1)
> set fmri(convolve_phase1) 0
>
> # Apply temporal filtering (EV 1)
> set fmri(tempfilt_yn1) 1
>
> # Add temporal derivative (EV 1)
> set fmri(deriv_yn1) 0
|