<div dir="ltr">Hi all,<div><br></div><div>I am trying to compare different methods to get ERSP data from EEGLab.</div><div><br></div><div><br></div><div>When pre-computing the channel measures for ERSP, I noticed that default parameters are: 'cycles', [3 0.8], 'nfreqs', 100, 'ntimesout', 200</div><div><br></div><div>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).</div><div><br></div><div><br></div><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="font-size:10pt;font-family:"courier new";color:black">clear </span><span style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">all</span><span style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="font-size:10pt;font-family:"courier new";color:black">clc</span><span style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="font-size:10pt;font-family:"courier new";color:black"> </span><span style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="font-family:"courier new";font-size:10pt">subject_list = { 'CA20SD'};</span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="font-family:"courier new";font-size:12pt"> </span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">nsubj = length(subject_list); </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:forestgreen">%
number of subjects </span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:forestgreen"> </span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">home_path = </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'D:\EEG_FGX_EEGLab_3\10.ChannelsInterpolated\'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">;</span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> </span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:black;font-family:"courier new";font-size:10pt"> allersp = zeros([100 200 1
nsubj]);</span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allitc = zeros([ 100 200 1
nsubj]);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allpowbase = zeros([ 1 100
1 nsubj]);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> alltimes = zeros([ 1 200 1
nsubj]);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allfreqs = zeros([ 1 100 1
nsubj]);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:forestgreen;font-family:"courier new";font-size:10pt"> </span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:blue">for</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> subj=1:nsubj</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:black;font-family:"courier new";font-size:10pt"> </span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:forestgreen">% Load dataset</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> fprintf(</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'\n\n\n**** %subj: Loading dataset ****\n\n\n'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, subject_list{subj});</span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> EEG = pop_loadset(</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'filename'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">,
[subject_list{subj} </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'_epochs.set'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">], </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'filepath'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, home_path);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:forestgreen"> </span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:forestgreen;font-family:"courier new";font-size:10pt">% Compute a time-frequency decomposition for every electrode</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:blue">for</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> elec = 77</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">[ersp,itc,powbase,times,freqs,erspboot,itcboot] = pop_newtimef(EEG, 1,
elec, [EEG.xmin EEG.xmax]*1000, </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'cycles'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, [3 0.8], </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'nfreqs'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, 100, </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'ntimesout'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, 200, </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'plotersp'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">,</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'off'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">, </span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'plotitc'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">,</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:rgb(160,32,240)">'off'</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">);</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allersp (:,:,elec,subj) =
ersp;</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allitc (:,:,elec,subj) = itc;</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> allpowbase (:,:,elec,subj) =
powbase;</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> alltimes (:,:,elec,subj) =
times;</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> </span><span style="font-size:10pt;font-family:"courier new";color:black">allfreqs (:,:,elec,subj) = freqs;</span><span style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:blue">end</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">;</span><br></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:blue">end</span><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black">;</span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-CA" style="font-size:10pt;font-family:"courier new";color:black"> </span><span lang="EN-CA" style="font-size:12pt;font-family:"courier new""><span></span></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt">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. </p><p class="MsoNormal" style="margin-bottom:0.0001pt"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt">I tried to look at the <span style="color:rgb(0,0,0);white-space:pre-wrap">std_ersp() function to see if other inputs could help to make calculations similar but I did not find anything. </span></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:rgb(0,0,0);white-space:pre-wrap">Thanks again,</span></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></p><p class="MsoNormal" style="margin-bottom:0.0001pt"><span style="color:rgb(0,0,0);white-space:pre-wrap">Simon</span></p></div>