Hi Roland,
it seems you end up with some complex hyperparameter estimate... This is
probably caught in the new SPM version right after spm_reml. But why they
get complex is another question. It seems that ReML is converging very
slowly... Try increasing the number of covariance components slowly, e.g.
by assuming independence and allow for unequal variance only for one of
your two factors.
Volkmar
On Mon, 20 Aug 2007, Roland Marcus Rutschmann wrote:
> Hi,
>
> i have a problem while estimating a VBM within subject design. There are 4
> points in time and 16 subjects. I choose 2 Factors (Subject, time) and use
> the flexible factorial method, all scans in one go (with factor matrix) and 2
> Main effects. The design matrix looks as expected.
>
> As long as I choose both factors to be independent and equal I can estimate
> the model fine both with the newest spm5 updates and the version before that.
> But If I (as originally intended) leave any of the factors to allow for
> unequallity and thus estimating hyperparameters, I get errors.
>
> The latest SPM-Version gives me an error after estimating the hyperparameters.
> (see printout below).
>
> The version before that estimates the model fine but I cannot get results
> because SPM.xX.nKX becomes complex (which I think is strange).
>
> If anybody has an idea I'd be thankfull and will of course provide all
> necessary information.
>
> -----------error latest Version
> Running "Model estimation"
>
> SPM5: spm_spm (v$Rev: 783 $) 14:57:45 - 20/08/2007
> ========================================================================
> Initialising parameters : ...done
> Plane 181/181, block 1/1 : ...done
> Temporal non-sphericity (over voxels) : ...REML estimation
> ReML Iteration : 1 ...2.244799e+01
> ReML Iteration : 2 ...1.215567e+04
> ReML Iteration : 3 ...4.844410e+21
> ...
> ReML Iteration : 127 ...1.112480e+00
> ReML Iteration : 128 ...6.819108e-01
> Error running job:
> Output argument "V" (and maybe others) not assigned during call to
> "/misc/spm/spm5_new/spm_reml.m (spm_reml)".
> In file "/misc/spm/spm5_new/spm_spm.m", function "spm_spm" at line 901.
> In file "/misc/spm/spm5_new/spm_config_fmri_est.m", function "run_est" at line
> 394.
> In file "/misc/spm/spm5_new/spm_jobman.m", function "run_struct1" at line
> 1540.
> In file "/misc/spm/spm5_new/spm_jobman.m", function "run_struct1" at line
> 1548.
> In file "/misc/spm/spm5_new/spm_jobman.m", function "run_struct1" at line
> 1548.
> In file "/misc/spm/spm5_new/spm_jobman.m", function "run_struct" at line 1492.
> ??? Output argument "V" (and maybe others) not assigned during call to
> "/misc/spm/spm5_new/spm_reml.m (spm_reml)".
>
> ??? Error while evaluating uicontrol Callback
>
>
> ------------error second to last version
> Running "Model estimation"
>
> SPM5: spm_spm (v$Rev: 714 $) 16:08:31 - 20/08/2007
> ========================================================================
> Done
> --------------------------
> Done.
>>> --------------------------
> Running "Factorial design specification"
> Mapping files : ...done
> Design configuration : ...done
> Saving SPM configuration : ...SPM.mat saved
> Design reporting : ...done
> Done
> --------------------------
> Running "Model estimation"
>
> SPM5: spm_spm (v$Rev: 714 $) 16:09:19 - 20/08/2007
> ========================================================================
> Initialising parameters : ...done
> Plane 181/181, block 1/1 : ...done
> Temporal non-sphericity (over voxels) : ...REML estimation
> ReML Iteration : 1 ...2.245747e+01
> ReML Iteration : 2 ...8.954866e+02
> ReML Iteration : 3 ...2.400890e+01
> ReML Iteration : 4 ...6.112560e+00
> ReML Iteration : 5 ...4.550228e+01
> ReML Iteration : 6 ...2.271020e+01
> ReML Iteration : 7 ...2.198245e+01
> ReML Iteration : 8 ...-4.817806e+01
>
> SPM5: spm_spm (v$Rev: 714 $) 16:11:50 - 20/08/2007
> ========================================================================
> Initialising parameters : ...done
> Output images : ...initialised
> Plane 181/181, block 1/1 : ...done
> Spatial non-sphericity (over scans) : ...writing resels/voxel image
> Saving results : ...done
> Completed : 16:29:03 - 20/08/2007
> ...use the results section for assessment
>
> Done
> --------------------------
> Done.
>>>
> SPM5: spm_results_ui (v$Rev: 656 $) 16:38:32 - 20/08/2007
> ========================================================================
> ??? Error using ==> image
> Error using ==> image
> Image CData can not be complex..
>
> Error in ==> spm_conman at 960
> hDesMtxIm = image((varargin{2}.nKX+1)*32);
>
> Error in ==> spm_conman at 929
> spm_conman('ImDesMtx',xX,H.hDesMtxAx) %-Depict
> DesMtx
>
> Error in ==> spm_conman at 812
> [F,cF] = spm_conman('Initialise','on',SPM,STATmode,n,Prompt,Mcstr,OK2chg);
>
> Error in ==> spm_getSPM at 278
> [Ic,xCon] = spm_conman(SPM,'T&F',Inf,...
>
> Error in ==> spm_results_ui at 264
> [SPM,xSPM] = spm_getSPM;
>
> ??? Error using ==> [hReg,xSPM,SPM] = spm_results_ui;
> Error using ==> image
> Error using ==> image
> Image CData can not be complex..
>
> ??? Error while evaluating uicontrol Callback
>
>
>
>
> Regards and thanks,
>
> Roland
>
>
>
--
Volkmar Glauche
-
Department of Neurology [log in to unmask]
Universitaetsklinikum Freiburg Phone 49(0)761-270-5331
Breisacher Str. 64 Fax 49(0)761-270-5416
79106 Freiburg http://fbi.uniklinik-freiburg.de/
|