[FieldTrip] Fwd: correlation with behavioural data and cluster statistics

Arjen Stolk a.stolk8 at gmail.com
Sun Mar 17 07:12:29 CET 2019 

Dear Juan,

Your code looks fine to me. Indeed, also with ft_statfun_indepsamplesregrT,
the behavioral vector goes into the design matrix.

All p-values computed with the montecarlo method correspond to one tail
and, thus, you would need to correct for this one-sided test approach if
you intended to do a two-sided test. Either cfg.correcttail = 'alpha' or
'prob' will achieve this, as further explained here
<http://www.fieldtriptoolbox.org/faq/why_should_i_use_the_cfg.correcttail_option_when_using_statistics_montecarlo/>
.

Have a great weekend,
Arjen

On Tue, Mar 12, 2019 at 7:35 AM Juan Lei <juan.lei at hotmail.com> wrote:

> Dear Arjen,
>
> if I use ft_statfun_indepsamplesregrT, are all the Parameters the same as correlationT
> method ? Below are all my parameters used. Are they correct? Another
> question is whether I need the "cfg.correcttail  =  'alpha'" for the
> indepsamplesregrT methoed?
>
>    cfg=[];
>    cfg.channel          = {'*'};​
>    cfg.frequency        = [15 35];
>    cfg.avgoverfreq      ='no';
>    cfg.latency          = [-0.7 0];​
>    cfg.avgoverchan      = 'no';​
>    cfg.avgovertime      = 'no';​
>    cfg.parameter        = 'powspctrm';​
>    cfg.method           = 'montecarlo';​
>    cfg.alpha            = 0.05;
>    cfg.numrandomization = 1000;​
>    cfg.correctm         = 'cluster';​
>    cfg.clusteralpha     = 0.05;​
>    cfg.clusterstatistic = 'maxsum';​
>    cfg.correcttail      =  'alpha';
>    cfg.tail             = 0;​
>    cfg.clustertail      = 0;​
>    cfg.minnbchan        = 2;​
>    cfg.neighbours       = neighbours;​
>    cfg.statistic        ='ft_statfun_indepsamplesregrT';
>    cfg.ivar             = 1;​
>    cfg.design = [];​
>    cfg.design(1,:) = xa;​% xa is my behavioral vector. Each subject has
> one behavioral value.
>    freq_stat = ft_freqstatistics(cfg,GA_RI);% GA_RI is my time frequency
> data averaged across group but keeping individuals.
>
>
> Best regards,
>
> Juan Lei
>
>
>
>
> ------------------------------
> *From:* fieldtrip <fieldtrip-bounces at science.ru.nl> on behalf of Valeriya
> Belyaeva <vabelyaeva94 at gmail.com>
> *Sent:* Tuesday, February 19, 2019 7:11 PM
> *To:* FieldTrip discussion list
> *Subject:* Re: [FieldTrip] Fwd: correlation with behavioural data and
> cluster statistics
>
> Thank you a lot, Arjen.
>
> Best regards,
> Valeriya
>
> вт, 19 февр. 2019 г. в 18:44, Arjen Stolk <a.stolk8 at gmail.com>:
>
> Hi Valeriya,
>
> That looks to be correct use of the correlationT method. The behavioral
> data goes into the design matrix, as you did.
>
> Best,
> Arjen
>
> On Feb 19, 2019, at 6:35 AM, Valeriya Belyaeva <vabelyaeva94 at gmail.com>
> wrote:
>
>
>
> ---------- Forwarded message ---------
> From: *Valeriya Belyaeva* <vabelyaeva94 at gmail.com>
> Date: ср, 13 февр. 2019 г. в 16:28
> Subject: correlation with behavioural data and cluster statistics
> To: <fieldtrip at science.ru.nl>
>
>
> Dear fieldtrip developers,
>
> I have a question about correlation of neuronal and behavioural data. You
> have discussed this question here:
> http://www.fieldtriptoolbox.org/faq/how_can_i_test_for_correlations_between_neuronal_data_and_quantitative_stimulus_and_behavioural_variables/
> <https://eur02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.fieldtriptoolbox.org%2Ffaq%2Fhow_can_i_test_for_correlations_between_neuronal_data_and_quantitative_stimulus_and_behavioural_variables%2F&data=02%7C01%7C%7Cf5787d0a00b04a9a6a5a08d696962ca2%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636861969065513175&sdata=PlfWV7rf8x9WKS373TuqhMfhfqyE0C4RNTwGnbJP54M%3D&reserved=0>
>
> I have taken the code, which was given as an example and added
> cluster-based correction method. My goal is to perform correlation between
> neuronal and behavioural data and then to account for multiple comparisons
> with cluster statistics.
>
> Could you, please, tell me, whether this script does it?
>
> Here is the script:
>
> % example script
>
> cfg = [];
>
> cfg.statistic        = 'ft_statfun_correlationT';
>
> cfg.method           = 'montecarlo';
>
> cfg.numrandomization = 1000;
>
> cfg.correctm         = 'cluster';
>
>
> cfg.clusteralpha     = 0.025;
>
> cfg_neighb.method = 'triangulation';
>
> cfg_neighb.layout = 'acticap_62';
>
> cfg.neighbours = ft_prepare_neighbours(cfg_neighb, data_brain{:});
>
>
> n1 = 42;    % n1 is the number of subjects
>
> design(1,1:n1)       = behav_data(:,1); % mean reaction time by subject
>
> cfg.design           = design;
>
> cfg.ivar             = 1;
>
> stat = ft_freqstatistics(cfg, data_brain{:});
>
> Thank you a lot in advance!
> --
> Best regards,
> Valeriya Belyaeva
> HSE, Cognitive science and technologies
>
>
>
> --
> Best regards,
> Valeriya Belyaeva
> HSE, Cognitive science and technologies
>
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> <https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7Cf5787d0a00b04a9a6a5a08d696962ca2%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636861969065523187&sdata=r3MoHOSBgBaGepkioATjwMBN5Pw0LKaFlTLMns7O4JU%3D&reserved=0>
> https://doi.org/10.1371/journal.pcbi.1002202
> <https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7Cf5787d0a00b04a9a6a5a08d696962ca2%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636861969065533192&sdata=WgU%2FWOboFOmczUC7XlDwB4xBym2CpNKJTKg4%2FSuEfD8%3D&reserved=0>
>
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> <https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7Cf5787d0a00b04a9a6a5a08d696962ca2%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636861969065543204&sdata=y%2BPI%2FkTrU8qHK724FjFsqmip%2BQX6vPlj5yHEf7YZorc%3D&reserved=0>
> https://doi.org/10.1371/journal.pcbi.1002202
> <https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7Cf5787d0a00b04a9a6a5a08d696962ca2%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636861969065553209&sdata=nlElRbo3ZfgypPKhwVkHm7bfbJVnYlOVZWkvJzo%2BkDg%3D&reserved=0>
>
>
>
> --
> Best regards,
> Valeriya Belyaeva
> HSE, Cognitive science and technologies
>
> _______________________________________________
> 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/20190316/71d5f66c/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