Hi Mike,
After reading Chapter 35 on group-level analyses, I’ve been thinking on how to apply strategies 2a and 2b to a within-subject analysis. Since I work with iEEG, the number and coverage of the electrodes varies considerably across subjects. I imagine that I could define time-frequency-region windows based on my hypotheses, but averaging power within a region might not be the best strategy given the higher spatial resolution of iEEG compared to M/EEG. Furthermore, my sample size is considerably small, and therefore I think it makes more sense to focus on within-subject statistics.
For defining TF-window boundaries based on pixels that are statistically significant, while avoiding circular inference, I was wondering if comparing time-frequency samples in all conditions vs. baseline (which you mention is a common practice in fMRI research) is equivalent to your suggestion in the YouTube videos of choosing the window based on the average of all data (across all conditions)? Moreover, if my hypothesis includes more than one channel (within or across brain regions), does it make sense to choose independent TF windows for each channel? Couldn’t the time-frequency peaks vary across channels, similar to how they may vary across subjects? Treating channels as independent would also imply forming clusters across time-frequency only, and I was wondering if one would then have to correct for multiple comparisons across channels?
Finally, for testing power relative to baseline, is there a difference in shuffling the trial labels of the baseline periods and of the activation periods (I think this is the way it’s performed in Fieldtrip) to form the null-hypothesis distribution, compared to the method of temporally shifting the time series (depicted in Fig. 34.2A)? Specifically, are there differences in the null hypotheses of these two methods? Also, when would it be better to use the method (shifting the baseline period) shown in Fig. 34.2B? Is it better suited for experiments in which your ITI is very short (or non-existent)? Do you gain anything by preserving the temporal structure of the time series?
Thanks!
--Gabriel
Hi Mike,
After reading Chapter 35 on group-level analyses, I’ve been thinking on how to apply strategies 2a and 2b to a within-subject analysis. Since I work with iEEG, the number and coverage of the electrodes varies considerably across subjects. I imagine that I could define time-frequency-region windows based on my hypotheses, but averaging power within a region might not be the best strategy given the higher spatial resolution of iEEG compared to M/EEG. Furthermore, my sample size is considerably small, and therefore I think it makes more sense to focus on within-subject statistics.
For defining TF-window boundaries based on pixels that are statistically significant, while avoiding circular inference, I was wondering if comparing time-frequency samples in all conditions vs. baseline (which you mention is a common practice in fMRI research) is equivalent to your suggestion in the YouTube videos of choosing the window based on the average of all data (across all conditions)?
Moreover, if my hypothesis includes more than one channel (within or across brain regions), does it make sense to choose independent TF windows for each channel? Couldn’t the time-frequency peaks vary across channels, similar to how they may vary across subjects?
Treating channels as independent would also imply forming clusters across time-frequency only, and I was wondering if one would then have to correct for multiple comparisons across channels?
Finally, for testing power relative to baseline, is there a difference in shuffling the trial labels of the baseline periods and of the activation periods (I think this is the way it’s performed in Fieldtrip) to form the null-hypothesis distribution, compared to the method of temporally shifting the time series (depicted in Fig. 34.2A)? Specifically, are there differences in the null hypotheses of these two methods? Also, when would it be better to use the method (shifting the baseline period) shown in Fig. 34.2B? Is it better suited for experiments in which your ITI is very short (or non-existent)? Do you gain anything by preserving the temporal structure of the time series?
--
Thanks!
--Gabriel
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/f4c2a741-d2fa-4dee-9008-e21eee80f071%40googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimeseriesdata+unsub...@googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/f4c2a741-d2fa-4dee-9008-e21eee80f071%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/ad47ed4b-8563-47f6-8741-278b21aeea29%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/5444f7cb-6204-48f3-83c4-5bcb2c2898bd%40googlegroups.com.
Hi Gabriel. I don't think that's necessary. If the frequency window is relatively narrow, then the 1/f won't really bias the spectrum. And if the frequency window is wide enough that the 1/f will introduce a bias, then the window is probably too wide ;)That said, baseline normalization has several benefits aside from getting rid of the 1/f issue. For example, baseline normalization also helps to separate background/ongoing activity from the specific task-related modulation. In general, I recommend baseline normalization unless there is a specific reason not to use it.Mike
On Thu, Sep 12, 2019 at 11:38 AM Gabriel Obregon-Henao <gabrielobr...@gmail.com> wrote:
Hey Mike,
Does one need to account for 1/f when averaging power across a TF-window? In the Fieldtrip example of analyzing high-gamma in human ECoG, for example, they multiply the power values at each time-frequency sample by the square of the corresponding frequency prior to averaging across frequencies within the high-gamma range. I’m mainly asking because I’m planning on using Strategy 2a from the book for a between-trials analysis.
Thanks!
-Gabriel
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimeseriesdata+unsub...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/5444f7cb-6204-48f3-83c4-5bcb2c2898bd%40googlegroups.com.
Understood. That makes me wonder why in your toy example of comparing between the first and second half of the experimental trials you omitted the normalization step. Also, if one were to extrapolate from that example to a comparison between two experimental conditions, does it really make sense to run stats on the whole trial? Wouldn't you just want to compare the post-stimulus periods between conditions regardless of whether you perform baseline normalization?
I was also wondering why do you use n-1 degrees of freedom in Matlab's tinv function if you constructed the null distribution using Welch's t-test. Shouldn't we be using the Welch-Satterthwaite equation for estimating the degrees of freedom (or n-2 if we assume equal variances)?
Finally, I've seen that in Fieldtrip when you use a Montercarlo sample, and not the full permutation, they adjust their p-values because the minimum p-value shouldn't be zero but 1/num_permutations. How do we account for this in your code?
Thanks!--Gabriel
On Saturday, September 14, 2019 at 3:51:51 AM UTC-7, Mike X Cohen wrote:
Hi Gabriel. I don't think that's necessary. If the frequency window is relatively narrow, then the 1/f won't really bias the spectrum. And if the frequency window is wide enough that the 1/f will introduce a bias, then the window is probably too wide ;)That said, baseline normalization has several benefits aside from getting rid of the 1/f issue. For example, baseline normalization also helps to separate background/ongoing activity from the specific task-related modulation. In general, I recommend baseline normalization unless there is a specific reason not to use it.Mike
On Thu, Sep 12, 2019 at 11:38 AM Gabriel Obregon-Henao <gabrielobr...@gmail.com> wrote:
Hey Mike,
Does one need to account for 1/f when averaging power across a TF-window? In the Fieldtrip example of analyzing high-gamma in human ECoG, for example, they multiply the power values at each time-frequency sample by the square of the corresponding frequency prior to averaging across frequencies within the high-gamma range. I’m mainly asking because I’m planning on using Strategy 2a from the book for a between-trials analysis.
Thanks!
-Gabriel
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/5444f7cb-6204-48f3-83c4-5bcb2c2898bd%40googlegroups.com.
--
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/7d2f287c-773d-4ecb-9933-a53a88f3ad38%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/934c0eba-965b-495f-b777-ef5220fdb89b%40googlegroups.com.
You received this message because you are subscribed to a topic in the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/analyzingneuraltimeseriesdata/MQXHm06upo4/unsubscribe.
To unsubscribe from this group and all its topics, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/CAABvs%2BvNjwC3_ndXCnpy9tzu1NfCV2yWfwoA20ngThyV8gLW-A%40mail.gmail.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/889DC832-EF19-4BFF-A6EF-A4684A5830B7%40gmail.com.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimeseriesdata+unsub...@googlegroups.com.
Hi Mike,I saw that you migrated the group to a new platform, for which I'll definitely sign up, but I wanted to go back to this thread. Specifically, you recommended that I plot effect sizes across the brain, and I was wondering if you have a good guide reference for doing so (I've only seen people report p-values)?Thanks,--Gabriel
On Monday, August 26, 2019 at 2:42:04 AM UTC-7, Mike X Cohen wrote:
Hi Gabriel. Apologies for the delayed reply. See below.On Thu, Aug 22, 2019, 23:50 Gabriel Obregon-Henao <gabrielobr...@gmail.com> wrote:Hi Mike,
After reading Chapter 35 on group-level analyses, I’ve been thinking on how to apply strategies 2a and 2b to a within-subject analysis. Since I work with iEEG, the number and coverage of the electrodes varies considerably across subjects. I imagine that I could define time-frequency-region windows based on my hypotheses, but averaging power within a region might not be the best strategy given the higher spatial resolution of iEEG compared to M/EEG. Furthermore, my sample size is considerably small, and therefore I think it makes more sense to focus on within-subject statistics.
For defining TF-window boundaries based on pixels that are statistically significant, while avoiding circular inference, I was wondering if comparing time-frequency samples in all conditions vs. baseline (which you mention is a common practice in fMRI research) is equivalent to your suggestion in the YouTube videos of choosing the window based on the average of all data (across all conditions)?
Yes, averaging over conditions first will avoid biased data selection for condition comparisons.Moreover, if my hypothesis includes more than one channel (within or across brain regions), does it make sense to choose independent TF windows for each channel? Couldn’t the time-frequency peaks vary across channels, similar to how they may vary across subjects?
Interesting thought. It depends on the quality of the data. I'd be concerned about the peak-finding algorithm getting caught up by noise.Treating channels as independent would also imply forming clusters across time-frequency only, and I was wondering if one would then have to correct for multiple comparisons across channels?
I guess it depends on how you set it upand what level you will use to make inferences. "electrode" is the sample of the population of all possible electrodes in this group of patients. So you wouldn't need to correct for the number of electrodes. It might be useful to do a lot of qualitative visualizations, for example, showing effects sizes across the brain.
Finally, for testing power relative to baseline, is there a difference in shuffling the trial labels of the baseline periods and of the activation periods (I think this is the way it’s performed in Fieldtrip) to form the null-hypothesis distribution, compared to the method of temporally shifting the time series (depicted in Fig. 34.2A)? Specifically, are there differences in the null hypotheses of these two methods? Also, when would it be better to use the method (shifting the baseline period) shown in Fig. 34.2B? Is it better suited for experiments in which your ITI is very short (or non-existent)? Do you gain anything by preserving the temporal structure of the time series?
The baseline-shifting method is useful when you want to test every TF pixel in the map.
--
Thanks!
--Gabriel
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/f4c2a741-d2fa-4dee-9008-e21eee80f071%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/5b34a8ed-efca-4da1-976d-9aa58b9ccb1e%40googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimeseriesdata+unsub...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/f4c2a741-d2fa-4dee-9008-e21eee80f071%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimeseriesdata+unsub...@googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/f4c2a741-d2fa-4dee-9008-e21eee80f071%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/ad47ed4b-8563-47f6-8741-278b21aeea29%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "AnalyzingNeuralTimeSeriesData" group.
To unsubscribe from this group and stop receiving emails from it, send an email to analyzingneuraltimes...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/analyzingneuraltimeseriesdata/faa6cffa-f1b8-4197-8e37-f795bc7c4a9c%40googlegroups.com.