Hi Yael,
> In 1st level-
> 1. Is each image contrast received after computing the 1st level, an
> average of the interpolated images in the time interval I have previously
> specified?
>
Right now, you can only look at averages. However, in principle, you can
enter any contrast (over peri-stimulus time). I'll add this to the code
soon.
> 2. Is there a way to do more sophisticated temporal matrixes, for example
> specify wavelets in the 1st level as described in Kiebel&Friston, 2004b ?
>
> 3. Is the 1st level design not to be estimated? I'm asking this since there
> is no SPM file created in 1st level , however according the SPM5 manual,
> Kiebel&Friston paper and the version before update 456, it was possible to
> make single subject inference from the 1st level.
>
For SPM5, these two things didn't make it into the software. The main
reason was that with 1st level design matrices it can easily happen that
you specify design matrices that don't capture all your interesting
effects. This is due to the huge variability of evoked responses over
different types of experiments. You would lose a lot of sensitivity, at
the 2nd level, if you used the 'wrong' kind of first level design matrix.
Note that you still can do a 1st level analysis by averaging over single
trials (i.e., 1st lvl: single trials, 2nd lvl: effects averaged within
trial type).
> 2cnd level:
> 1. Full factor analysis: 4 factors- got a message that it can't run
> this design.
>
This is a slight limitation of the software. However, I'd recommend
breaking your design apart and not attempt to model everything in one
huge model. If possible, you should keep the tests at the 2nd level as
simple as possible, e.g. a 1-sample t-test. The disadvantage of this
approach is that you have, inconveniently, to specify one model for each
hypothesis. The advantage of that approach is that the modelling is
simpler, with respect to variance parameters.
> 2. Full factor with 2-factors-the design was run, but when estimated,
> doesn't produce any parameters (betas, Res etc.).
> 3. When I go to results I get a message: model not estimated would
> you like to estimate it now? I press yes. I define a contrast and get an
> error:
>
> ??? Reference to non-existent field 'VResMS'.
>
> Error in ==> spm_contrasts at 46
> VHp = SPM.VResMS;
>
> Error in ==> spm_getSPM at 430
> SPM = spm_contrasts(SPM, unique([Ic, Im]));
>
> Error in ==> spm_results_ui at 264
> [SPM,xSPM] = spm_getSPM;
>
> ??? Error while evaluating uicontrol Callback.
>
This is strange. Could you send me, by email, the saved job-file of this
analysis?
Thanks for your email, Stefan
|