<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 "help spectopo" in Matlab, for the variable "spectra":<br>
</p>
<p><br>
</p>
<div>Outputs:</div>
<div>         spectra  = (nchans,nfreqs) power spectra (mean power over epochs), in dB</div>
<div>         freqs    = frequencies of spectra (Hz)<br>
</div>
<div><br>
</div>
<div>If the unit of variable "spectra" is in dB unit already, then when computing the absolute power, do we still need to do 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>)? <br>
</div>
<div><br>
</div>
<div>I know when plotting the figure, the spectrum is plotted in log scale. But if the unit of spectra is in dB, then there's no need to transform, right? <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 <mmiyakoshi@ucsd.edu><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> </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>
> 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] = spectopo(EEG.data(2,:,:), 0, EEG.srate);</div>
<div><br>
</div>
<div>% delta=1-4, theta=4-8, alpha=8-13, beta=13-30, gamma=30-80</div>
<div>deltaIdx = find(freqs>1 & freqs<4);</div>
<div>thetaIdx = find(freqs>4 & freqs<8);</div>
<div>alphaIdx = find(freqs>8 & freqs<13);</div>
<div>betaIdx  = find(freqs>13 & freqs<30);</div>
<div>gammaIdx = find(freqs>30 & freqs<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  = 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">
<<a href="mailto:SHAN0055@e.ntu.edu.sg" target="_blank">SHAN0055@e.ntu.edu.sg</a>></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 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 or script so that I can get the power of individual band (e.g. power of alpha band or beta band) from your computation, instead of just a continuous line in the power spectrum figure?  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 "set digest mime" 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>