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

Tue Mar 12 15:24:20 CET 2019

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<mailto: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<mailto:vabelyaeva94 at gmail.com>> wrote:

---------- Forwarded message ---------
From: Valeriya Belyaeva <vabelyaeva94 at gmail.com<mailto:vabelyaeva94 at gmail.com>>
Date: ср, 13 февр. 2019 г. в 16:28
Subject: correlation with behavioural data and cluster statistics
To: <fieldtrip at science.ru.nl<mailto: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

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