[FieldTrip] ft_sourcestatistics

[jan-mathijs schoffelen]

Hi Max,

At the moment we are designing some major code changes to deal with source level data, and we need to ask you for some patience in this regard.
For now, you could try to move the call to ft_checkdata to outside ft_sourcestatistics, and change the powdimord into 'pos_rpt'; 
Does it work if you subsequently call ft_sourcestatistics (with cfg.implementation = 'new')? This new implementation has never fully matured and may be instable. It has been found to work in certain situations, but it has never been fully tested on all possible processing histories that users may throw at it. Hence the urgent need to clean up the code. 

Best,

Jan-Mathijs



On Mar 7, 2012, at 6:46 PM, Stenner, Max-Philipp wrote:

> Dear fieldtrip users,
> 
> 
> 
> I am encountering a problem with ft_sourcestatistics when using output from ft_sourcegrandaverage. My parameter of interest is a logratio: log(mean(moments(time of interest))/mean(moments(baseline))) which I simply named 'pow' to be compatible with ft functions. Apart from this, nothing else deviates from how fieldtrip handles the data/nothing else was changed 'manually'.
> 
> 
> 
> The problem arises from the fact that the existence of 10 different subjects in the grandaverage (in the usual 1x10 struct field 'trial') isn't represented in the 'powdimord' field of the data after ft_checkdata(varargin{i}, 'datatype', {'source', 'volume'}, 'feedback', 'no', 'inside', 'index', 'sourcerepresentation', 'new'); in the new implementation of ft_sourcestatistics:
> 
> 
> 
> varargin{1}
> 
> ans =
> 
> 
> 
> (before ft_checkdata:)
> 
> 
> 
>        pos: [8196x3 double]
>        var: [1x1 struct]
>     dimord: 'voxel'
>      trial: [1x10 struct]
>     inside: [8196x1 double]
>    outside: [0x1 double]
>         df: [8196x1 double]
>        cfg: [1x1 struct]
> 
> 
> 
> (after ft_checkdata:)
> 
> 
> 
>          pos: [8196x3 double]
>          var: [1x1 struct]
>       dimord: 'pos'
>       inside: [8196x1 double]
>      outside: [0x1 double]
>           df: [8196x1 double]
>          cfg: [1x1 struct]
>          pow: [8196x10 double]
>    powdimord: 'pos'
> 
> 
> 
> This leads to a (wrong) transpose of the design (originally row1 = uvar, row2 = ivar [with two conditions for a depsamples T test]) and ultimately to the message
> 
> 'The data must contain at least two units (usually subjects).'
> 
> 
> 
> Any help would be much appreciated, thanks very much!
> 
> Max
> 
> (fieldtrip version: fieldtrip-lite-20111130 on a Windows7 PC)
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip

Jan-Mathijs Schoffelen, MD PhD 

Donders Institute for Brain, Cognition and Behaviour, 
Centre for Cognitive Neuroimaging,
Radboud University Nijmegen, The Netherlands

Max Planck Institute for Psycholinguistics,
Nijmegen, The Netherlands

J.Schoffelen at donders.ru.nl
Telephone: +31-24-3614793

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20120308/0147cacd/attachment-0002.html>