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

Arjen Stolk a.stolk8 at gmail.com
Sat Jun 29 00:35:36 CEST 2019


Hi Juan,

Yes. Indeed, you will have to use ft_timelockstatistics
specifying cfg.statistic = 'ft_statfun_indepsamplesregrT'. For more
information on testing for correlations between neural and say behavioral
data, please see here
<http://www.fieldtriptoolbox.org/faq/how_can_i_test_for_correlations_between_neuronal_data_and_quantitative_stimulus_and_behavioural_variables/>
.

Best,
Arjen

On Thu, Jun 27, 2019 at 8:33 AM Juan Lei <juan.lei at hotmail.com> wrote:

> Dear Arjen,
>
> can this ft_statfun_indepsamplesregrT also work for ERP data, i. e. the
> data only have three dimensions (time, channel and subject)?
>
> In this case, of course I have to use
> ft_timelockstatistics(cfg,GAalldevRi) in the end.
>
> Best,
>
> Juan
> ------------------------------
> *From:* fieldtrip <fieldtrip-bounces at science.ru.nl> on behalf of Arjen
> Stolk <a.stolk8 at gmail.com>
> *Sent:* Sunday, March 17, 2019 7:12 AM
> *To:* FieldTrip discussion list
> *Subject:* Re: [FieldTrip] Fwd: correlation with behavioural data and
> cluster statistics
>
> 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
> <https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.fieldtriptoolbox.org%2Ffaq%2Fwhy_should_i_use_the_cfg.correcttail_option_when_using_statistics_montecarlo%2F&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747439566&sdata=cCUJE92o9BrAd35xx36cP0fkXukAUhsNW%2BEt4xMnwCM%3D&reserved=0>
> .
>
> 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://nam04.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%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747449571&sdata=qic0njE42iEFnLzLpjI5pOH3Qixs6BfGYLlsS2yuVPk%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://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747459581&sdata=N8cEFn91nSkR8l0OA%2B6%2BqQYs%2Bcnx9ZJytZMu8uQFDYo%3D&reserved=0>
> https://doi.org/10.1371/journal.pcbi.1002202
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747469586&sdata=dr%2FO6kt9rbRIPu3GpQk1D9%2BVuxMfg5OxPJzNQgoZheQ%3D&reserved=0>
>
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747489602&sdata=q91RtnrucbD%2Fd%2Fn9v1f9EyJihjXCUSbstmU9gYK66Wo%3D&reserved=0>
> https://doi.org/10.1371/journal.pcbi.1002202
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747499607&sdata=q3op8RD6C6E3WFxcQ2I8%2Fg7bkfzX6Vyc5jsGHfWOfPk%3D&reserved=0>
>
>
>
> --
> Best regards,
> Valeriya Belyaeva
> HSE, Cognitive science and technologies
>
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747510970&sdata=wHjEE%2B8%2FqnJYLVLjx7OmC9bBLhhV84h4ZttZtC5pUdo%3D&reserved=0>
> https://doi.org/10.1371/journal.pcbi.1002202
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7C6482c0ed5b854521095608d6aaa1afee%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636884008747519617&sdata=gVJd%2BEAz4hsZKd6v2lDgXG98gsnhELRu%2B9%2BHU80jAFc%3D&reserved=0>
>
> _______________________________________________
> 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/20190628/1d4fd883/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