<div dir="ltr">Dear Mubeen,<div><br></div><div>1) My short answer would be:</div><div><br></div><div>a) There are very good "Fieldtip Lectures" on the topic that addresses questions of frequency resolution and time-windows in general and which would be a good place to start. </div><div>There are several people presenting nearly the same slides, so you can take a pick based on flavour, e.g.:</div><div><a href="https://www.youtube.com/watch?v=QLvsa1r1Voc">https://www.youtube.com/watch?v=QLvsa1r1Voc</a><br></div><div><a href="https://www.youtube.com/watch?v=vwPpSglPJTE">https://www.youtube.com/watch?v=vwPpSglPJTE</a><br></div><div><a href="https://www.youtube.com/watch?v=dHTuzMsjVJA">https://www.youtube.com/watch?v=dHTuzMsjVJA</a><br></div><div><br></div><div>b) If you really want to estimate power changes *over time*, using wavelets is one way, and using tapered sliding time-windows using FFT is another ('cfg.method = 'mtmconvol'). The latter allows you even more freedom in how much you want to smooth over frequencies (see lectures for explanation of "multi tapering").</div><div><br></div><div>d) Your question suggests, however, that you would like to capture the 1-4Hz peaks in SWDs *averaged over time*. For this you would not necessarily need a time-resolved frequency estimate, right? E..g. an averaged FFT done on each segment might suffice and be more to the point. Wavelets do not come into play then. Correct me if I'm wrong! Perhaps doing a peak detection on the spikes - and calculating the inter-spike intervals - might be another approach.<br></div><div><br></div><div>c) If you are set on doing a time-resolved
spectral
analysis, then you will probably be able to see
broad-band high-frequency power associated with the SWD peaks. However, to accurately estimate power at lower frequencies, especially as low as 1Hz, you might need time-widows of over a second (e.g. 3 seconds if you want to base it on 3 cycles of 1Hz, which would be a bare minimum). So you might not be able to resolve those low frequencies over time. Again, I think the lecture might clarify this point, and I am not sure you are asking for time-resolved analysis.<br></div><div><div><br></div><div>d) Feel free to experiment with different parameters, this will give you a sense of the data and limits of frequency analyses, and the effect of certain parameters. In the end it also depends on how stationary the spectral content of the signal is, which is sometimes hard to determine a-priori.</div><div><br></div><div>e) Wavelets can be easily used in a way that creates far more data than necessary, e.g. by resolving at every sample, instead of in e.g. steps of tens of milliseconds. So you might want to keep memory usage and CPU time in mind. A bit of experimentation might help you set parameters that are "good enough" but don't overload your computer's resources.</div><div></div></div><div><br></div><div>2a) Baseline normalization of the timecourses is a good idea before frequency analysis:</div><div><a href="https://www.fieldtriptoolbox.org/faq/why_does_my_tfr_look_strange/">https://www.fieldtriptoolbox.org/faq/why_does_my_tfr_look_strange/</a><br></div><div><br></div><div>2b) Baseline normalization of the time-resolved power estimates can only be done <i>after </i>doing the spectral estimates, right? If/when depends on the reason you want to do so though. For exploring and plotting, note that you can do it 'on the fly' in FieldTrip's plotting functions.<br></div><div> </div><div>Anyway, that's the 'short' answer. I might be able to say more if I know more about the goal of the frequency analysis of PWDs (with which I am somewhat familiar). And check the videos!</div><div><br></div><div>Happy fieldtripping,</div><div>Stephen</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Op do 18 feb. 2021 om 07:10 schreef mubeen afzal <<a href="mailto:mubafzal@hotmail.com">mubafzal@hotmail.com</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi all,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I am aiming to calculate average spectral values for segment of EEG containing spike/wave discharges in patients which can last anywhere between 0.5 seconds to 6 seconds. I am planning on morlet wavelet convolution but am unsure what the best parameters would
be in terms of width/cycles/timewin of the morlet configurations. Typically these spike/waves on visual inspection are anywhere between 1hz to more commonly between 3-4 Hz or less often 5 hz frequency discharges. </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<ol>
<li>What would be the best parameters for a wavelet convolution?</li></ol>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
2. is baseline normalization usually done before spectral analysis or can it be done after? </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Any help would be highly appreciated. </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Regards,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Mubeen J </div>
</div>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
<a href="https://doi.org/10.1371/journal.pcbi.1002202" rel="noreferrer" target="_blank">https://doi.org/10.1371/journal.pcbi.1002202</a><br>
</blockquote></div>