<div dir="ltr">Dear Pollet,<br><br>> I don’t know where I’ve got it wrong. Is there an obvious error I have made?<br><br>I can't answer to this question. I also want someone to take a look at my groupSIFT code to guarantee there is no error.<br><br>> I wonder if the way EEGLAB calculates things when producing scalp maps is different from the script I have, and if you might know of a way that I could check this? Can I access the data that is used to create the scalp maps?<div><br></div><div>Yes, you can use a debug mode in Matlab. Open one of the functions being used to generate the figure, put a red circle on the left column of the edit window, then run the code. The process will stop there, so you can trace the process line by line. It's a basic technique so if you haven't tried it yet, I recommend you challenge it this time!</div><div><br></div><div>Makoto</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 6, 2017 at 3:14 AM, Pollet S.C. <span dir="ltr"><<a href="mailto:sebastien.pollet@soton.ac.uk" target="_blank">sebastien.pollet@soton.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="EN-GB" link="#0563C1" vlink="#954F72">
<div class="m_4686952079170469228WordSection1">
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have 20 epoched datasets from 10 different subjects – 2 datasets per subject as there were 2 different conditions. I created averaged scalp maps (for all 10 subjects together, for all channels) in the alpha range (8-12Hz) using the ‘Study/Precompute
 channel measures’ and ‘Study/Plot channel measures’ functions. I used 10 EEGLAB ‘.set’ files for one condition and another 10 for a different condition.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">When I ran the ‘Precompute channel measures routine’, I used these parameters: 'cycles', [3 0.5], 'baseline', [-3500 -1500], 'freqs', [2 60], 'alpha', 0.05, 'ntimesout', 400, 'padratio',1.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I then used ‘Plot channel measures’, using ‘8 12’ (8-12Hz) as the frequency parameter and did this repeatedly for 100ms time windows to see what is happening from -1000ms before the event to 1500ms after the event.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I then wanted to plot ERSP over time for a select number of channels (as visually identified using the scalp maps – those with strong ERD/ERS ), so I used a function created in MatLab that uses the ‘newtimef’ function, which creates the
 ‘ersp_sig’ variable containing significant ERSP. ‘ersp_sig’ is a 3-element vector – ersp_sig(a,b,c) where ‘a’ is frequency,’b’ is time points (400 per epoch), and ‘c’ is channels (35 in this case).<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Here is the function I used:<u></u><u></u></p>
<p class="MsoNormal">------------------------------<wbr>------------------------------<wbr>-------------------<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue"><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">function</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> [ersp,ersp_sig] = plotERSPTopo(dirname,filename)</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">eegchannels = 1:35;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">cycles = [3 0.5];</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">freqs = [2 60];</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">baseline = [-3500 -1500];</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">tlimits = [-4000 4000];</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">alpha = 0.05;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">mine = -7;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">maxe = 7;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">eeg = pop_loadset(filename,dirname);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">close
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">all</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> i=eegchannels</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    tmpsig = eeg.data(i,:,:);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    tmpsig = tmpsig(:);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    triallength = size(eeg.data,2);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    close
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">all</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    fprintf(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'PROCESSING CHANNEL #%2.0f\n'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">
 ,i);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    [ersp(:,:,i),itc,powbase,<wbr>times,freqs,erspboot,itcboot] =
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">...</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">        newtimef(tmpsig,triallength,<wbr>tlimits,eeg.srate,cycles,
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">...</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'baseline'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,baseline,
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'alpha'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">, alpha,
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'freqs'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">, freqs,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'padratio'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,1,
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">...</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'plotersp'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'on'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'plotitc'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'<wbr>off'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'timesout'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,400,
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">...</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">,</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'title'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,eeg.chanlocs(i).<wbr>labels,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'erspmax'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,7);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    erspimage{i} = getframe(gcf,[116 80 347 308]);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    </span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">    </span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">
 i1=1:size(ersp,1);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> i2=1:size(ersp,2);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">           
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">if</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> (ersp(i1,i2,i)>erspboot(i1,1)) & (ersp(i1,i2,i)<erspboot(i1,2))<wbr>;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">                ersp_sig(i1,i2,i)=0;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">           
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">else</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">                ersp_sig(i1,i2,i)=ersp(i1,i2,<wbr>i);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">           
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">       
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">------------------------------<wbr>------------------------------<wbr>-------------------<u></u><u></u></p>
<p class="MsoNormal">I then used the following function to plot what I was interested in.<u></u><u></u></p>
<p class="MsoNormal">X is the file with all 10 datasets appended together for one condition and Y is another file with 10 appended datasets for the other condition.<u></u><u></u></p>
<p class="MsoNormal">(These are the same datasets I used to create the scalp maps).<u></u><u></u></p>
<p class="MsoNormal">------------------------------<wbr>------------------------------<wbr>-------------------<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">function</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> erd_alpha_AB10 = erd_alpha_AB10(X,Y</span><span lang="FR-CA" style="font-size:10.0pt;font-family:"Courier New";color:black">)</span><span lang="FR-CA" style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="FR-CA" style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">% Plot, for AB10, ERD% for selected channels</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%Load file for left movements</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ersp_alpha_left_sig = load(X,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'ersp_sig'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ersp_alpha_left_sig = ersp_alpha_left_sig.ersp_sig;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%Load file for right movements</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ersp_alpha_right_sig = load(Y,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'ersp_sig'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ersp_alpha_right_sig = ersp_alpha_right_sig.ersp_sig;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%Convert to ERD/ERS%</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">erd_alpha_left_sig=((10.^((<wbr>ersp_alpha_left_sig)/10))-1)*<wbr>100;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">erd_alpha_right_sig=((10.^((<wbr>ersp_alpha_right_sig)/10))-1)*<wbr>100;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%Calculate mean values for the alpha range<u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%(alpha range (8-12Hz) is actually indices 7-11 as frequency range starts at 2Hz)</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue"><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> j=1:400;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> k=1:35;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">            erd_alpha_left_sig_mean(j,k)= (erd_alpha_left_sig(7,j,k)+<wbr>erd_alpha_left_sig(8,j,k)+erd_<wbr>alpha_left_sig(9,j,k)+erd_<wbr>alpha_left_sig(10,j,k)+erd_<wbr>alpha_left_sig(11,j,k))/5;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> j=1:400;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">for</span><span style="font-size:10.0pt;font-family:"Courier New";color:black"> k=1:35;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">            erd_alpha_right_sig_mean(j,k)= (erd_alpha_right_sig(7,j,k)+<wbr>erd_alpha_right_sig(8,j,k)+<wbr>erd_alpha_right_sig(9,j,k)+<wbr>erd_alpha_right_sig(10,j,k)+<wbr>erd_alpha_right_sig(11,j,k))/<wbr>5;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">   
</span><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">% Plots</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">t=[-4000:20:3999];</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">% Left</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch1_left=12;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%FC2</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch2_left=27;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%CP4</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch3_left=23;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%CP3</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">figure;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">title([</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'Alpha ERD% for all AB10 subjects for left sided movements'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">]);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">xlabel(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'Time, in ms (relative to movement onset at 0)'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ylabel(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'ERD/ERS %'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">)</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">axis([-4000 4000 -100 100]);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">hold
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">on</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_left_sig_<wbr>mean(:,ch1_left));</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_left_sig_<wbr>mean(:,ch2_left),</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'r'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_left_sig_<wbr>mean(:,ch3_left),</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'g'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">legend(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'FC2'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'CP4'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'CP3'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'<wbr>Location'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'NorthWest'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">% Right</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch1_right=11;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%FCz</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch2_right=23;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%CP3</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ch3_right=27;
</span><span style="font-size:10.0pt;font-family:"Courier New";color:forestgreen">%CP4</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">figure;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">title([</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'Alpha ERD% for all AB10 subjects for right sided movements'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">]);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">xlabel(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'Time, in ms (relative to movement onset at 0)'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">ylabel(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'ERD/ERS %'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">)</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">axis([-4000 4000 -100 100]);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">hold
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">on</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">;</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_right_sig_<wbr>mean(:,ch1_right));</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_right_sig_<wbr>mean(:,ch2_right),</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'r'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">plot(t,erd_alpha_right_sig_<wbr>mean(:,ch3_right),</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'g'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:black">legend(</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'FCz'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'CP3'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'CP4'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'<wbr>Location'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:10.0pt;font-family:"Courier New";color:#a020f0">'NorthWest'</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:12.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Courier New";color:blue">end</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">------------------------------<wbr>------------------------------<wbr>-------------------<u></u><u></u></p>
<p class="MsoNormal">When I look at the plots, they don’t quite match what I see on the scalp maps. For example when I plot FC2, CP3 and CP4 together, my plots show that CP4 is activated before FC2 and a with a greater value, which is not what I see on the
 scalp maps – FC2 has a dark blue ‘patch’ (strong ERD) whilst CP4 is green – so I was expecting to see this reflected on the plot.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="m_4686952079170469228MsoPlainText">I don’t know where I’ve got it wrong. Is there an obvious error I have made? I wonder if the way EEGLAB calculates things when producing scalp maps is different from the script I have, and if you might know of a way that I could check
 this? Can I access the data that is used to create the scalp maps?<u></u><u></u></p>
<p class="m_4686952079170469228MsoPlainText"><u></u> <u></u></p>
<p class="MsoNormal">Thanks for any advice on this.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Seb<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Sebastien Pollet<u></u><u></u></p>
<p class="MsoNormal">University of Southampton<u></u><u></u></p>
<p class="MsoNormal"><a href="mailto:sebastien.pollet@soton.ac.uk" target="_blank">sebastien.pollet@soton.ac.uk</a>
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>

<br>______________________________<wbr>_________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/<wbr>eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">eeglablist-unsubscribe@sccn.<wbr>ucsd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu">eeglablist-request@sccn.ucsd.<wbr>edu</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego<br></div></div>
</div>