[FieldTrip] Choose minimun time-width of clusters for cluster based permutations on freq_statistics

[Constantino Méndez Bértolo]

Dear Jörn,

Thanks for your response. You are sensibly confused. It is because we
wanted to do the statistics (montecarlo; correctm = cluster) with just
one channel per data/individual. [As was pointed out above:
size(Template.powspctrm) = 9 1 1 320; Template.dimord =
subj_chan_freq_time], and ft_timelockstatistics would ask me for a
neighbourstructure relenteless, while ft_freqstatistics is ok with
cfg.neighbour = [] and so clusters in time. Hence the unorthodoxy.

I did a function that parses the mask of the already computed stat
finding '1s', determining the span and setting zeros on an arbitrary
width-criteria; which serves the purpose but it is not as ellegant as
it would be my own statfun.

Cheers,
Tino


2013/7/29 "Jörn M. Horschig" <jm.horschig at donders.ru.nl>
>
> Dear Constantino,
>
> I'm pretty confused why you want to trick FieldTrip and put timelock data into a frequency structure. You are aware that there is a ft_timelockstatistics function that can cluster in time, aren't you? Anyway, imho cfg.minnbchan defines how many neighbouring samples another samples needs to have, and a timebin can have max. two neighbouring samples, so it does not make much sense here (and I think it really applies to channels only).
>
> The error you describe is because FT wants to compute neighbouring channels, because you did not set cfg.avgoverchan='yes'. If you do that, it will run through smoothly without an error. However, the best way for you would be to write your own statfun by making a copy of the statfun of your choice and modifying stat.stat to be 1 if the clustersize is < 10. That way, you could achieve what you want.
>
> Best,
> Jörn
>
>
> On 7/26/2013 8:55 PM, Constantino Méndez Bértolo wrote:
>>
>> Using cfg.minnbchan = 1, 5 or 40 and cfg.neighbours = [] finds no cluster (not even with minnbchan = 1) were there was a cluster before (without cfg.minnbchan)
>>
>> Using cfg.neighbours = 'triangulation', or 'distance', with or without cfg.minnbchan computes the statistics but fails prior to compute the clusters in the randomizations
>>
>> This is the error (with neighbours) V V V
>>
>> Error in clusterstat>makechanneighbstructmat (line 519)
>> [seld] = match_str(cfg.channel, cfg.neighbours(chan).label);
>>
>> Error in clusterstat (line 59)
>> channeighbstructmat = makechanneighbstructmat(cfg);
>>
>> Error in statistics_montecarlo (line 320)
>> [stat, cfg] = clusterstat(cfg, statrand, statobs,'issource',issource);
>>
>> Error in ft_freqstatistics (line 267)
>> [stat, cfg] = statmethod(cfg, dat, cfg.design);
>>
>>
>> While it seems it is minnbchan what will be helpful, we do not come up with an idea. Maybe, besides tricking fieldtrip to think we have frequency data we should also trick him to believe that points in time are channels? :)
>>
>>
>> El 26 de julio de 2013 19:54, smoratti at psi.ucm.es <smoratti at psi.ucm.es> escribió:
>>>
>>>
>>> I think cfg.minnbchan, applies also to time bins.
>>>
>>> s
>>>
>>> ________________________________________________________
>>> Stephan Moratti, PhD
>>>
>>> see also: http://web.me.com/smoratti/
>>>
>>> Universidad Complutense de Madrid
>>> Facultad de Psicología
>>> Departamento de Psicología Básica I
>>> Campus de Somosaguas
>>> 28223 Pozuelo de Alarcón (Madrid)
>>> Spain
>>>
>>> and
>>>
>>> Center for Biomedical Technology
>>> Laboratory for Cognitive and Computational Neuroscience
>>> Parque Científico y Tecnológico de la Universidad Politecnica de Madrid
>>> Campus Montegancedo
>>> 28223 Pozuelo de Alarcón (Madrid)
>>> Spain
>>>
>>>
>>> email: smoratti at psi.ucm.es
>>> Tel.:    +34 679219982
>>>
>>> El 26/07/2013, a las 18:46, Constantino Méndez Bértolo escribió:
>>>
>>>> cfg.minnbchan,
>>>
>>>
>>>
>>> _______________________________________________
>>> fieldtrip mailing list
>>> fieldtrip at donders.ru.nl
>>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>>
>>
>>
>>
>> --
>> Constantino Méndez-Bértolo
>> Laboratorio de Neurociencia Clínica, Centro de Tecnología Biomédica (CTB)
>>
>> Parque Científico y Tecnológico de la UPM, Campus de Montegancedo
>>
>> 28223 Pozuelo de Alarcón, Madrid, SPAIN
>>
>>
>> [image]
>>
>>
>>
>> _______________________________________________
>> fieldtrip mailing list
>> fieldtrip at donders.ru.nl
>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>
>
>
> --
> Jörn M. Horschig
> PhD Student
> Donders Institute for Brain, Cognition and Behaviour
> Centre for Cognitive Neuroimaging
> Radboud University Nijmegen
> Neuronal Oscillations Group
> FieldTrip Development Team
>
> P.O. Box 9101
> NL-6500 HB Nijmegen
> The Netherlands
>
> Contact:
> E-Mail: jm.horschig at donders.ru.nl
> Tel:    +31-(0)24-36-68493
> Web: http://www.ru.nl/donders
>
> Visiting address:
> Trigon, room 2.30
> Kapittelweg 29
> NL-6525 EN Nijmegen
> The Netherlands
>
>
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip




-- 
Constantino Méndez-Bértolo
Laboratorio de Neurociencia Clínica, Centro de Tecnología Biomédica (CTB)

Parque Científico y Tecnológico de la UPM, Campus de Montegancedo

28223 Pozuelo de Alarcón, Madrid, SPAIN


[image]