Hi Alexandra,
Unfortunately, featquery has not really been designed to work with .feat
directories that do not contain statistics. There's a strong likelihood
that, even when you were successfully creating time series plots using FSL
5.0.4, these errors were still occurring, but being ignored - it is
probably now failing in 5.0.8 due to the addition of more rigorous error
checking.
You have two options:
1. As you suggested, manually transform your mask image to subject space,
and run fslmeants with the transformed mask. This really shouldn't be too
hard - assuming that you have the inverse warp files and transformation
matrices, you just need to use applywarp to transform your mask image, and
then call fslmeants. If you don't have the inverse warps/matrices, you can
create them using the invwarp/convert_xfm commands. It should be easy
enough to wrap these commands in a script, which you can call repeatedly
for different subjects.
2. When you use Feat to preprocess your data, enable statistics and add a
simple regressor EV containing all ones. This shouldn't add too much to the
processing time, and will result in all of the statistic files, that
featquery needs to be present, being created.
Thanks,
Paul
On Wed, Mar 4, 2015 at 3:17 PM, SUBSCRIBE FSL Alexandra DiFeliceantonio <
[log in to unmask]> wrote:
> Hello All,
> I used featquery from FSL 5.0.4 to extract a time series from an ROI in a
> first level feat (no stats). I was using it because I could use the same
> standard space mask and featquery would convert it to subject space for me
> and extract the time series.
> Now, when I try and use the newest version of FSL (5.0.8) I get the
> following error:
>
> can't read "fmri(nstats)": no such element in array
> can't read "fmri(nstats)": no such element in array
> while executing
> "for { set n 0 } { $n < $fmri(nstats) } { incr n 1 } {
> if { $fmri(dostats.$n) == 1 } {
> set statslist "$statslist [ lindex $fmri(statslist) $n ]"
> ..."
> (procedure "featquery_proc" line 31)
> invoked from within
> "featquery_proc .r"
> invoked from within
> ".r.btns.apply invoke"
> ("uplevel" body line 1)
> invoked from within
> "uplevel #0 [list $w invoke]"
> (procedure "tk::ButtonUp" line 22)
> invoked from within
> "tk::ButtonUp .r.btns.apply"
> (command bound to event)
>
> It seems like it is looking for a stats image where previously it worked
> with no stats image. Is there a way around this problem? Is there another
> way I should be extracting the time series (I know fslmeants will do it,
> but then I need to create all of my masks in subject space). Any help, or
> pointing me to previous threads I've missed, would be much appreciated.
> Thanks,
> Alexandra
>
