<div dir="ltr">Dear Hendrik,<div><br></div>> Using 3 cycles at lowest frequency to 12.8 at highest.<br>> Okay, so this means that the wavelet has 3 cycles at my lowest estimated frequency and 12.8 cycles at the highest. If my frequency resolution based on the 5th output argument from newtimef() shows me 46 frequencies with the highest one being 50, then the 12.8 cycles correspond to the 50 Hz correct?<br><br>Yes, correct.<br> <br>>Generating 200 time points (-715.0 to 1713.0 ms)<br>Finding closest points for time variable<br>Time values for time/freq decomposition is not perfectly uniformly distributed<br>The window size used is 285 samples (570 ms) wide.<br>> I assume this is what you mean with the sliding window? Does this mean the length of each wavelet in time is 570 ms?<br><br>Yes. And it slides 200 times to generates 200 time points in this case. The sliding width is 1713-(-715) divided by 200. It makes sense, doesn't it?<br><br>> Estimating 46 linear-spaced frequencies from 5.9 Hz to 50.0 Hz.<br>> And I was wondering if this means that my time-freq decomposition is based on 46 wavelets? I am just curious since we use wavelets for electromyography processing in our lab and they are typically a lot wider space so I just want to make sure that I get this right. I have 46 wavelets to cover up to a frequency range of 50 Hz?<div class="gmail_extra"><br></div><div class="gmail_extra">Yes, you generated 46 wavelet daughters covering 5.9 to 50.0 Hz. These wavelets have both temporal and frequency distribution, so they will nicely overlap in both domains.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Makoto</div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 6, 2015 at 12:34 PM, Hendrik Enders <span dir="ltr"><<a href="mailto:endersh@ucalgary.ca" target="_blank">endersh@ucalgary.ca</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 lang="EN-CA" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Hi Makoto,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">As always, thanks for your help and your quick reply. I ran the code on some data of my subjects to obtain the information. I copy the command line
output just to check if I draw the right conclusions.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Computing Event-Related Spectral Perturbation (ERSP) and Inter-Trial Phase Coherence (ITC) images based on 222 trials of 1500 frames sampled at 500
Hz.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Each trial contains samples from -1000 ms before to 1998 ms after the timelocking event.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Image frequency direction: normal<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">
<b><i>Makes all sense.<u></u><u></u></i></b></span></p>
<p class="MsoNormal"><b><i><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></i></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Using 3 cycles at lowest frequency to 12.8 at highest.<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><b><i><span style="font-size:11pt;font-family:Calibri,sans-serif">Okay, so this means that the wavelet has 3 cycles at my lowest estimated frequency and 12.8 cycles at the highest.
If my frequency resolution based on the 5<sup>th</sup> output argument from newtimef() shows me 46 frequencies with the highest one being 50, then the 12.8 cycles correspond to the 50 Hz correct?<u></u><u></u></span></i></b></p>
<p class="MsoNormal" style="margin-left:36pt"><b><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Generating 200 time points (-715.0 to 1713.0 ms)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Finding closest points for time variable<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Time values for time/freq decomposition is not perfectly uniformly distributed<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">The window size used is 285 samples (570 ms) wide.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">
<b><i>I assume this is what you mean with the sliding window? Does this mean the length of each wavelet in time is 570 ms?<u></u><u></u></i></b></span></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Estimating 46 linear-spaced frequencies from 5.9 Hz to 50.0 Hz.<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><b><i><span style="font-size:11pt;font-family:Calibri,sans-serif">And I was wondering if this means that my time-freq decomposition is based on 46 wavelets? I am just curious since
we use wavelets for electromyography processing in our lab and they are typically a lot wider space so I just want to make sure that I get this right. I have 46 wavelets to cover up to a frequency range of 50 Hz?<u></u><u></u></span></i></b></p>
<p class="MsoNormal" style="margin-left:36pt"><b><i><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></i></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Thanks for your help. As always, much appreciated.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif">Hendrik<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(46,116,181)">Hendrik Enders</span></b><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(46,116,181)"><u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">PhD Candidate<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">NSERC Vanier scholar<u></u><u></u></span></b></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">Faculty of Kinesiology<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">University of Calgary<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">2500 University Dr NW<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">Calgary, AB, Canada, T2N 1N4<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">T: <a href="tel:%28403%29%20220-2413" value="+14032202413" target="_blank">(403) 220-2413</a><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(77,77,77)">E:
<a href="mailto:henders@kin.ucalgary.ca" target="_blank"><span style="color:rgb(77,77,77)">endersh@ucalgary.ca</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif"> Makoto Miyakoshi [mailto:<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>]
<br>
<b>Sent:</b> September-02-15 12:40 PM<br>
<b>To:</b> Hendrik Enders <<a href="mailto:endersh@ucalgary.ca" target="_blank">endersh@ucalgary.ca</a>><br>
<b>Cc:</b> <a href="mailto:eeglablist@sccn.ucsd.edu" target="_blank">eeglablist@sccn.ucsd.edu</a><br>
<b>Subject:</b> Re: [Eeglablist] wavelet properties<u></u><u></u></span></p><div><div class="h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Dear Hendrik,<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">When you run newtimef(),<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"> >> [ersp,itc,powbase,times,freqs,erspboot,itcboot] = ...<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> newtimef(data, frames, epochlim, srate, cycles,...<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> 'key1',value1, 'key2',value2, ... );<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">The 5th output 'freqs' is the list of center frequencies. I think it uses 'freqresol' output from dftfilt3.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">> The second point that I need to find out is the duration in the time domain for each wavelet. Can anybody help how I can generate that information to provide a table with the wavelet information?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">When you run newtimef() it shows the length of the sliding window. You can report the length there.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">In EEGLAB's newtimef, by default the number of cycles increases linearly as the central frequency increases, which may be tricky to explain. The purpose of doing this is to adjust the ratio of time/frequency resolution along with the frequency
axis (if you use the fixed cycle numbers the output image looks vertically stretched i.e. too much temporal resolution with too less freq resolution at high freqs)<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Makoto<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Aug 31, 2015 at 5:05 PM, Hendrik Enders <<a href="mailto:endersh@ucalgary.ca" target="_blank">endersh@ucalgary.ca</a>> wrote:<u></u><u></u></p>
<blockquote style="border-style:none none none solid;border-left-color:rgb(204,204,204);border-left-width:1pt;padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">I have finished my time-frequency analysis of my EEG data. Some people have asked me for some very specific information regarding the wavelet properties and I am unsure how to retrieve
this information. I have a sampling frequency of 500 Hz and I used the input [3 0.5] for cycles and a pad ratio of 2.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">How can I calculate the center frequency of each individual wavelet that was generated between 3 – 50 Hz?<u></u><u></u></p>
<p class="MsoNormal">I have tried the following:<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="font-size:10pt;font-family:'Courier New';color:black">freqs = [3 50];</span><u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="font-size:10pt;font-family:'Courier New';color:black">cycles = [3 0.5];</span><u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="font-size:10pt;font-family:'Courier New';color:black">srate = 500;</span><u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="font-size:10pt;font-family:'Courier New';color:black">[wavelet,cycles,freqresol,timeresol] = dftfilt3( freqs, cycles, srate);</span><u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="font-size:10pt;font-family:'Courier New';color:black"> </span><u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
This only returns two cells with wavelet information, time resolution and frequency resolution. How can I use the function to obtain center frequencies?<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
The second point that I need to find out is the duration in the time domain for each wavelet. Can anybody help how I can generate that information to provide a table with the wavelet information?<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
Thank you very much.<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">
<span style="color:rgb(136,136,136)">Hendrik<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(136,136,136)"> <u></u><u></u></span></p>
</div>
</div>
<p class="MsoNormal"><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" target="_blank">
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" target="_blank">
eeglablist-request@sccn.ucsd.edu</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">-- <u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Makoto Miyakoshi<br>
Swartz Center for Computational Neuroscience<br>
Institute for Neural Computation, University of California San Diego<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div></div></div>
</div>
</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>