<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Francesco,<br>
<br>
I don't entirely understand what you referring to either, but it
sounds at least like a within channel phenomenon? In that case, you
can use cfg.output = 'pow', or 'fourier', which do not compute the
between channel CSD, saving quite a bit memory.<br>
<br>
If you are interested in the total energy over time of your signal,
but only for the 1-150 Hz part, it might be easier to do a bandpass,
and taking the square of your raw signal (energy is usually squared
amplitude I thought?). If that is what your intending though.<br>
<br>
Cheers,<br>
Roemer<br>
<br>
<br>
On 08/11/2011 11:33, Eelke Spaak wrote:
<blockquote
cite="mid:CABPNLUryxd-8Frrb3kb6F+zw9OgOkPXRG0R_1R+Y1jMhknx5UQ@mail.gmail.com"
type="cite">
<pre wrap="">Dear Francesco,
Your question is not entirely clear to me. Do you want to estimate the
timecourse of the total (i.e. non-frequency-resolved) power in your
signal? In that case, frequency analysis is not the way to go;
basically your signal itself contains this information already (you
may want to rectify (abs()) it first, or apply a hilbert transform).
Also, I don't understand the reference to ft_componentanalysis.
Furthermore, the code you paste in your email seems valid to me. (You
could possibly simplify it by using cfg.toi = comp.time{1}; to get an
estimate of power at each time point in your signal. ) What does not
seem to work, in your case?
Best,
Eelke
2011/11/8 <a class="moz-txt-link-rfc2396E" href="mailto:f.dipompeo@unich.it"><f.dipompeo@unich.it></a>:
</pre>
<blockquote type="cite">
<pre wrap="">I would like to know if there is a ft function to calculate the Time Course
of Power of Electric and Magnetic Channels and IC coming from
ft_componentanalysis in the case of continuous resting state data.
I tried to do it using ft_freqanalysis with 'mtmconvol' as method and a
"toi" large as my whole acquisition time window (the script is below).
Time Course of Power should be obtained by summing up the results over all
the frequencies.
Thank you very much,
Francesco
-------------------------------------------------------
cfg=[];
cfg.method = 'mtmconvol';
cfg.output = 'powandcsd';
cfg.taper = 'hanning';
cfg.trials = 'all';
cfg.foi = 1:150;
cfg.toi = [comp.time{1}(1,1):0.02:comp.time{1}(1,end)];
cfg.t_ftimwin = 4./cfg.foi;
[freq_comp_tr] = ft_freqanalysis(cfg,comp);
_______________________________________________
fieldtrip mailing list
<a class="moz-txt-link-abbreviated" href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a>
<a class="moz-txt-link-freetext" href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a>
</pre>
</blockquote>
<pre wrap="">_______________________________________________
fieldtrip mailing list
<a class="moz-txt-link-abbreviated" href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a>
<a class="moz-txt-link-freetext" href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a>
</pre>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<font size="2"><font color="darkblue"><font face="calibri">Roemer
van der Meij M.Sc.<br>
PhD student<br>
Donders Institute for Brain, Cognition and Behaviour<br>
Centre for Cognition<br>
P.O. Box 9104<br>
6500 HE Nijmegen<br>
The Netherlands<br>
Tel: +31(0)24 3655932<br>
E-mail: <a class="moz-txt-link-abbreviated" href="mailto:r.vandermeij@donders.ru.nl">r.vandermeij@donders.ru.nl</a><br>
</font></font></font></div>
</body>
</html>