Questions on sourceinterpolate after the sourcestatistics

[gao zai]

Thank you very much Robert. I tried your suggestion, and set the
cfg.parameter='stat', still it gets stuck. As you mentioned, it maybe out of
memory. I am now trying to change to a powerful one.


Feng

On Tue, Apr 13, 2010 at 10:59 AM, Robert Oostenveld <
r.oostenveld at fcdonders.ru.nl> wrote:

> Dear Feng,
>
> I suspect that your computer is running out of memory while it is trying to
> interpolate all functional volumes onto the anatomical MRI grid. Instead of
> specifying
>
> cfg.parameter = 'all'
>
> I suggest that you only specify those parameters that you want to have
> interpolated. The negclusterslabelmat volume for example is not one that you
> want to have interpolated.
>
> best regards,
> Robert
>
>
> On 12 Apr 2010, at 15:55, gao zai wrote:
>
>  Dear all,
>>
>> I am now working on the sourcestatistics of the LCMV beamfomer. After
>> finished the volumenormalisation, sourcegrandaverage and sourcestatistics,
>> now I want to plot the t-values to the anatomical MRI. However, when I run
>> the ft_sourceinterpolate (codes see below), I wait for hours and response
>> with the matlab informing that "reslicing and interpolating
>> negclusterslabelmat "
>> --------------------------------------
>> %%statistics on the grandaverage
>> cfg=[];
>> cfg.dim         = gs42.dim;
>> cfg.parameter   = 'nai';
>> cfg.method      = 'montecarlo';
>> cfg.statistic   = 'depsamplesT';
>> cfg.correctm    = 'cluster';
>> cfg.numrandomization = 100;
>> cfg.alpha       = 0.05;
>> cfg.tail        = 0;
>>
>> nsubj=length(gs42.trial);
>> cfg.design(1,:) = [ones(1,nsubj) ones(1,nsubj)*2];
>> cfg.design(2,:) = [1:nsubj 1:nsubj];
>>
>> cfg.ivar        = 1; % row of design matrix that contains independent
>> variable (the conditrions)
>> cfg.uvar        = 2; % row of design matrix that contains subjects
>> number-2 groups
>>
>> stat = sourcestatistics(cfg, gs42,gs50);
>>
>> sMRI = read_mri(fullfile(spm('dir'), 'canonical', 'single_subj_T1.nii'));
>> cfg = [];
>> cfg.downsample = 2;
>> cfg.parameter = 'all';
>> statplot = ft_sourceinterpolate(cfg, stat, sMRI);
>>
>> --------------------------------------------------------------------------------------
>>
>> Does anybody how to deal with it? Thanks a lot in advance.
>>
>> Best,
>> FENG
>> ----------------------------------
>>
>> The aim of this list is to facilitate the discussion between users of the
>> FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and
>> EEG analysis.
>>
>> http://listserv.surfnet.nl/archives/fieldtrip.html
>>
>> http://www.ru.nl/fcdonders/fieldtrip/
>>
>>
> ----------------------------------
> The aim of this list is to facilitate the discussion between users of the
> FieldTrip  toolbox, to share experiences and to discuss new ideas for MEG
> and EEG analysis. See also
> http://listserv.surfnet.nl/archives/fieldtrip.html and
> http://www.ru.nl/neuroimaging/fieldtrip.
>

----------------------------------
The aim of this list is to facilitate the discussion between users of the FieldTrip  toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html and http://www.ru.nl/neuroimaging/fieldtrip.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20100413/3579ecf8/attachment-0002.html>