<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">hi roni,<div><br></div><div><div><blockquote type="cite"><div dir="ltr">1. How come when I use multitaper with 1 cycle I see a range starting at 5 hz, but when I use morelet wavelet with 1 cycle, the range starts at 10hz?<br></div></blockquote><div><br></div><div>i'm not sure. what do you mean with "range": the lowest frequency where you get a spectral estimation? just looking at your plots the difference doesn't seem dramatic: the multitaper seems to start somewhere at 6-7, the wavelet at 5. you resolution is set to 1 Hz.</div><br><blockquote type="cite"><div dir="ltr">2. Can I conclude that when the trial is rather short multitaper will be more accurate, or more sensitive to changes?<br></div></blockquote><div><br></div><div>without being an expert: with the multitaper you have more possibilities to do "fine-tuning" (i.e. spectral and temporal resolution can be specified separately). for the have a spectral smoothing of +/- foi*.4. for your wavelet it's +/- foi/1 (your width which for wavelets fixes the temporal AND spectral resolution). this means your frequency resolution in case of the wavelet is *terrible*!!</div><div><br></div><div>you should take a look here:</div><div><a href="http://fieldtrip.fcdonders.nl/tutorial/timefrequencyanalysis">http://fieldtrip.fcdonders.nl/tutorial/timefrequencyanalysis</a></div><div><br></div><div>in general it seems that for doing "good" time-frequency analysis, your trial lengths are far too short, particularly if you are interested in lower frequency activity. you should be more generous when chopping your data into trials (i usually do +/- 2s or so, do TF-analysis and then home-in into interpretable periods). you can see that in your case you hardly have values in the baseline (particularly <20 Hz).</div><div><br></div><div>when you do that (i.e. longer epochs) and use sensible m-factors (~5-7; or just use the fieldtrip default), then you will definitely get a completely different picture. </div><div><br></div><div>i'm not religious about wavelets and multitapers. with proper use both give good and usually very similar results. as said multitapers offer more possibilities for fine-tuning, but more freedom in choice is sometimes a heavy burden :-)</div><br><blockquote type="cite"><div dir="ltr">
3. How come the freq range in the multitaper changes depending on the baseline settings? It started at 5hz when cfg.baseline = [-0.2 -0.05], but started at 10hz when cfg.baseline = [-0.2 -0.1];<br></div></blockquote><div><br></div><div>again: take a look at the pictures. you basically have NaN's everywhere between -.2 and -.1. longer epochs are your friend.</div><div><br></div><div>good luck,</div><div>nathan</div><br><blockquote type="cite"><div dir="ltr"><br>Thanks a lot!<br>Roni <br>
<br> <br><br><div class="gmail_quote">On 3 May 2010 18:36, Roemer van der Meij <span dir="ltr"><<a href="mailto:r.vandermeij@donders.ru.nl">r.vandermeij@donders.ru.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Roni,<br>
<br>
I haven't looked in detail at all the other settings you use, but you apparently use a different number of cycles for the wavelets in both methods. For the multitapers, you use 2 cycles per wavelet (which are then convolved with the set of dpss tapers), whereas you use 1 cycle for your Morlet wavelets (which are then convolved with a hanning taper (or something similar)). This will make the time-window for spectral estimation twice as long for the multitapers than for your Morlet wavelets. You could try setting cfg.width = 2 for 'wltconvol'.<br>
<br>
Hope it helps a bit,<br>
<br>
Roemer<div><div></div><div class="h5"><br>
<br>
<br>
<br>
On 5/3/2010 5:08 PM, Roni Tibon wrote:<br>
</div></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="h5">
Dear all,<br>
<br>
I've been trying to analyze the same data using both methods, but I'm<br>
getting very different results, depending on the method I'm using.<br>
<br>
Enclosed is a picture of the results of the analysis in 3 channels, using<br>
both methods. In the left column I performed a multitaper analysis for each<br>
channel, using the following script:<br>
<br>
cfg = [];<br>
cfg.output = 'pow';<br>
cfg.channel = 'A195';<br>
cfg.method = 'mtmconvol';<br>
cfg.foi = 1:1:50;<br>
cfg.t_ftimwin = 2./cfg.foi;<br>
cfg.tapsmofrq = 0.4 *cfg.foi;<br>
cfg.toi = -0.2:0.05:1;<br>
cfg.pad = 'maxperlen';<br>
f12a_a195 = ft_freqanalysis(cfg, data);<br>
<br>
In the right one, I performed wavelet analysis using:<br>
<br>
cfg = [];<br>
cfg.channel = 'A195';<br>
cfg.method = 'wltconvol';<br>
cfg.output = 'pow';<br>
cfg.foi = 1:1:50;<br>
cfg.width = 1;<br>
cfg.toi = -0.2:0.05:1;<br>
cfg.keeptrials = 'yes';<br>
f12b_a195=freqanalysis(cfg, data);<br>
<br>
As you can see, channel A106 looks about the same using both methods.<br>
However, for channels A176 and A195, the multitaper analysis revealed<br>
activation in frequencies ranging from ~15-25, which was not seen in the<br>
wavelet analysis.<br>
<br>
I used a different program (msi) to divide the trial in two (0-500ms and<br>
500-1000ms) and compare activation in this frequency range between the 2<br>
halves, and indeed- there was more activation in the second half, as would<br>
be suggested by the multitaper.<br>
<br>
Can it be that when using wavelet analysis data is lost like this?<br>
Am I doing something wrong with the wavelet?<br>
<br>
Thanks,<br>
Roni<br>
<br>
<br>
<br></div></div>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
<br>
</blockquote>
<br>
-- <br>
Roemer van der Meij MSc<br>
Scientific Programmer& Data-Analyst<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 3612631<br>
E-mail: <a href="mailto:r.vandermeij@donders.ru.nl" target="_blank">r.vandermeij@donders.ru.nl</a><br>
<br>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
</blockquote></div><br></div><p>----------------------------------</p><p>The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis.</p><p> <a href="http://listserv.surfnet.nl/archives/fieldtrip.html">http://listserv.surfnet.nl/archives/fieldtrip.html</a></p><p> <a href="http://www.ru.nl/fcdonders/fieldtrip/">http://www.ru.nl/fcdonders/fieldtrip/</a></p>
</blockquote></div><br></div><br><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><div><div><div><div><div><div><div><div><b>--------------------------------------------</b></div><div>Dr. Nathan Weisz</div><div><br></div><div>OBOB-Lab</div><div>University of Konstanz</div><div>Department of Psychology</div><div>P.O. Box D23</div><div>78457 Konstanz</div><div>Germany</div><div><br></div><div>Tel: ++49 - (0)7531 - 88 45 84</div><div>Email: <a href="mailto:nathan.weisz@uni-konstanz.de">nathan.weisz@uni-konstanz.de</a></div><div>Homepage: <a href="http://www.uni-konstanz.de/obob">http://www.uni-konstanz.de/obob</a></div><div><br></div><div><div>"Nothing shocks me. I'm a scientist." (Indiana Jones)</div></div></div></div></div></div></div></div></div></div></div></div></span></span>
</div>
<br></body></html><p>----------------------------------</p>
<p>The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis.</p>
<p> http://listserv.surfnet.nl/archives/fieldtrip.html</p>
<p> http://www.ru.nl/fcdonders/fieldtrip/</p>