[Eeglablist] How can I get the power of each frequency band?

[Makoto Miyakoshi]

Dear Shan,

> I'm not very proficient in Matlab programming...

This could be a concern but you can try!

%%%%%%%%%%%%%%%%%%%
% for your epoched data, channel 2
[spectra,freqs] = spectopo(EEG.data(2,:,:), 0, EEG.srate);

% delta=1-4, theta=4-8, alpha=8-13, beta=13-30, gamma=30-80
deltaIdx = find(freqs>1 & freqs<4);
thetaIdx = find(freqs>4 & freqs<8);
alphaIdx = find(freqs>8 & freqs<13);
betaIdx  = find(freqs>13 & freqs<30);
gammaIdx = find(freqs>30 & freqs<80);

% compute absolute power
deltaPower = mean(10.^(spectra(deltaIdx)/10));
thetaPower = mean(10.^(spectra(thetaIdx)/10));
alphaPower = mean(10.^(spectra(alphaIdx)/10));
betaPower  = mean(10.^(spectra(betaIdx)/10));
gammaPower = mean(10.^(spectra(gammaIdx)/10));
%%%%%%%%%%%%%%%%%

Makoto


On Tue, Apr 14, 2015 at 6:23 AM, #SHAN XIN# <SHAN0055 at e.ntu.edu.sg> wrote:

>  Hi EEGLab team,
>
> 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...thank you so much!!
>
>  Best wishes,
>  SHAN Xin
>
>
>
> _______________________________________________
> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
> To unsubscribe, send an empty email to
> eeglablist-unsubscribe at sccn.ucsd.edu
> For digest mode, send an email with the subject "set digest mime" to
> eeglablist-request at sccn.ucsd.edu
>



-- 
Makoto Miyakoshi
Swartz Center for Computational Neuroscience
Institute for Neural Computation, University of California San Diego
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20150415/33c2cb22/attachment.html>