[Eeglablist] Comparison STUDY ERSPDATA and ERSPDATA from script

Simon Recherche srigoulotresearch at gmail.com
Thu Sep 1 11:07:31 PDT 2016


Hi all,

I am trying to compare different methods to get ERSP data from EEGLab.


When pre-computing the channel measures for ERSP, I noticed that default
parameters are: 'cycles', [3 0.8], 'nfreqs', 100, 'ntimesout', 200

I tried to retrieve the same results by adapting Arnaud's script on the
tutorial , for one subject and one electrode so it goes faster (see below).


clear all

clc



subject_list = { 'CA20SD'};



nsubj = length(subject_list); % number of subjects



home_path  = 'D:\EEG_FGX_EEGLab_3\10.ChannelsInterpolated\';



        allersp = zeros([100 200 1 nsubj]);

        allitc = zeros([ 100 200 1 nsubj]);

        allpowbase = zeros([ 1 100 1 nsubj]);

        alltimes = zeros([ 1 200 1 nsubj]);

        allfreqs = zeros([ 1 100 1 nsubj]);



for subj=1:nsubj



  % Load dataset

        fprintf('\n\n\n**** %subj: Loading dataset ****\n\n\n',
subject_list{subj});

        EEG = pop_loadset('filename', [subject_list{subj} '_epochs.set'],
'filepath', home_path);



% Compute a time-frequency decomposition for every electrode

for elec = 77

[ersp,itc,powbase,times,freqs,erspboot,itcboot] = pop_newtimef(EEG, 1,
elec, [EEG.xmin EEG.xmax]*1000, 'cycles', [3 0.8], 'nfreqs', 100,
'ntimesout', 200, 'plotersp','off', 'plotitc','off');


    allersp (:,:,elec,subj) = ersp;

    allitc (:,:,elec,subj) = itc;

    allpowbase (:,:,elec,subj) = powbase;

    alltimes (:,:,elec,subj) = times;

    allfreqs (:,:,elec,subj) = freqs;

end;

end;





When I compare the matrix A = ersp(:,:,1,1) with the figure I get with the
STUDY GUI (by selecting the proper channel and subjects), I dont have the
same results.


I tried to look at the std_ersp() function to see if other inputs could
help to make calculations similar but I did not find anything.



Thanks again,


Simon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20160901/3dafe5dd/attachment.html>


More information about the eeglablist mailing list