<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dear Makoto,<br>
</p>
<p><br>
</p>
<p>I have a doubt for the power unit. I read the &quot;help spectopo&quot;&nbsp;in Matlab, for the variable &quot;spectra&quot;:<br>
</p>
<p><br>
</p>
<div>Outputs:</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;spectra &nbsp;= (nchans,nfreqs) power spectra (mean power over epochs), in dB</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;freqs &nbsp; &nbsp;= frequencies of spectra (Hz)<br>
</div>
<div><br>
</div>
<div>If the unit of&nbsp;variable &quot;spectra&quot; is in dB&nbsp;unit already, then when computing the absolute power, do we still need to do&nbsp;the transformation as you suggested (i.e.
<span style="color: rgb(255, 0, 0); font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12.6666669845581px; background-color: rgb(255, 255, 255);">
10.^(spectra(deltaIdx)/10)</span>)?&nbsp;<br>
</div>
<div><br>
</div>
<div>I know when plotting the figure, the spectrum&nbsp;is&nbsp;plotted in log scale. But if the unit of spectra is&nbsp;in dB, then there's no need to transform, right?&nbsp;<br>
</div>
<div><br>
</div>
<div>I don't know if this is correct... Could you look more into it? Thank you so much!!<br>
</div>
<div><br>
</div>
<div>Shan Xin<br>
</div>
<p><br>
</p>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Makoto Miyakoshi &lt;mmiyakoshi@ucsd.edu&gt;<br>
<b>Sent:</b> Thursday, 16 April 2015 11:11 AM<br>
<b>To:</b> #SHAN XIN#<br>
<b>Cc:</b> eeglablist@sccn.ucsd.edu<br>
<b>Subject:</b> Re: [Eeglablist] How can I get the power of each frequency band?</font>
<div>&nbsp;</div>
</div>
<div>
<div dir="ltr">
<div><span style="color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px">Dear Shan,</span></div>
<span style="color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px">
<div><span style="color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px"><br>
</span></div>
&gt; I'm not very proficient in Matlab programming...</span><br>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">This could be a concern but you can try!</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">
<div style="font-size:12.6666669845581px">%%%%%%%%%%%%%%%%%%%</div>
<div style="font-size:12.6666669845581px">
<div>
<div>% for your epoched data, channel 2</div>
<div>[spectra,freqs] =&nbsp;spectopo(EEG.data(2,:,:), 0,&nbsp;EEG.srate);</div>
<div><br>
</div>
<div>%&nbsp;delta=1-4,&nbsp;theta=4-8,&nbsp;alpha=8-13,&nbsp;beta=13-30, gamma=30-80</div>
<div>deltaIdx = find(freqs&gt;1 &amp; freqs&lt;4);</div>
<div>thetaIdx = find(freqs&gt;4 &amp; freqs&lt;8);</div>
<div>alphaIdx&nbsp;= find(freqs&gt;8 &amp; freqs&lt;13);</div>
<div>betaIdx &nbsp;= find(freqs&gt;13 &amp; freqs&lt;30);</div>
<div>gammaIdx = find(freqs&gt;30 &amp; freqs&lt;80);</div>
<div><br>
</div>
<div><font color="#ff0000">% compute absolute power</font></div>
<div><font color="#ff0000">deltaPower = mean(10.^(spectra(deltaIdx)/10));</font></div>
<div><font color="#ff0000">thetaPower = mean(10.^(spectra(thetaIdx)/10));</font></div>
<div><font color="#ff0000">alphaPower = mean(10.^(spectra(alphaIdx)/10));</font></div>
<div><font color="#ff0000">betaPower &nbsp;= mean(10.^(spectra(betaIdx)/10));</font></div>
<div><font color="#ff0000">gammaPower = mean(10.^(spectra(gammaIdx)/10));</font></div>
</div>
<div>%%%%%%%%%%%%%%%%%</div>
<div><br>
</div>
<div>Makoto</div>
</div>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Tue, Apr 14, 2015 at 6:23 AM, #SHAN XIN# <span dir="ltr">
&lt;<a href="mailto:SHAN0055@e.ntu.edu.sg" target="_blank">SHAN0055@e.ntu.edu.sg</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-color:rgb(204,204,204); border-left-style:solid; padding-left:1ex">
<div dir="ltr">
<div style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; background-color:rgb(255,255,255)">
<p><span style="font-size:12pt">Hi EEGLab&nbsp;team,</span></p>
<div dir="ltr" style="color:rgb(33,33,33)">
<div>
<div style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; background-color:rgb(255,255,255)">
<div><br>
Does EEGlab has plugin&nbsp;or&nbsp;script&nbsp;so that I can get&nbsp;the&nbsp;power of individual band (e.g. power of alpha band or beta band) from your computation, instead of just&nbsp;a continuous line in the power spectrum&nbsp;figure?&nbsp; I'm not very proficient in Matlab programming...t<span style="font-size:12pt">hank
 you so much!!</span></div>
<div><br>
</div>
<div>Best wishes,<br>
</div>
<div>SHAN Xin<br>
</div>
<div><br>
</div>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" target="_blank">
http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">
eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
For digest mode, send an email with the subject &quot;set digest mime&quot; to <a href="mailto:eeglablist-request@sccn.ucsd.edu">
eeglablist-request@sccn.ucsd.edu</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="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>
</div>
</div>
</div>
</div>
</body>
</html>