[FieldTrip] One-sample monte-carlo permutation statistics in Fieldtrip?

Julian Keil julian.keil at gmail.com
Thu May 2 12:31:12 CEST 2019


Hi Lars,

we recently faced a similar problem (see here: https://www.nature.com/articles/s41598-019-42380-x <https://www.nature.com/articles/s41598-019-42380-x>). Our solution was to create „dummy“ data based on the actual data and use these in the comparison.
In short, this relates to Eric Maris’ statement about the null hypothesis: In our case, the H0 was „the regression weights are independent of the stimulus category“, so we test our empirically found regression weights against random regression weights. Please note that we don’t have a formal proof that his is a valid approach.

Unfortunately, I don’t think this is directly applicable to your situation, as it sounds like you only have one category, but maybe there is a way for you to come up with dummy data? 
Again, please note that there might be problems with our approach we didn’t consider.

Best,

Julian


> Am 02.05.2019 um 12:05 schrieb Lars Costers <larscosters at gmail.com>:
> 
> Hi all,
> 
> I was wondering whether there are any one-sample permutation statistics to test for a difference from zero implemented in Fieldtrip, ideally to work with monte-carlo permutations? Didn’t find any standard of 'statfun’ options in the documentation.
> 
> My goal is to do maxstat correction on ERF MEG data over voxel space.
> cfg.latency = [-0.2 0.8];
> cfg.parameter = 'avg';
> cfg.method = 'montecarlo';
> cfg.numrandomization = 2000;
> cfg.correctm = 'max'; % MaxStat correction
> cfg.design    = ones(1,nsub);
> cfg.ivar = 1;  % independent variable
> cfg.statistic = ???;
> [stat] = ft_timelockstatistics (cfg, timelock{:,1});
> 
>  I read the discussions about one-sample cluster-based permutation tests (e.g. https://mailman.science.ru.nl/pipermail/fieldtrip/2018-August/012314.html <https://mailman.science.ru.nl/pipermail/fieldtrip/2018-August/012314.html>). However, for me it's not possible to permute the baseline and post stimulus condition because I have expectation reactions in the baseline. 
> Anybody could suggest me how I could validly test whether my ERF is signficantly different from zero at every time-frame and electrode?  
> 
> Thanks, 
> 
> Lars
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> https://doi.org/10.1371/journal.pcbi.1002202

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20190502/911e259a/attachment-0002.html>
-------------- next part --------------
_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://doi.org/10.1371/journal.pcbi.1002202


More information about the fieldtrip mailing list