[FieldTrip] single-trial phase at stimulus onset using Hilbert transform

[Schoffelen, J.M. (Jan Mathijs)]

Dear Georgios,

Note that the default filter direction implemented for the ‘firws’ filter is onepass-zerophase. I recommend to use this (thus not to specify ‘onepass’), but in addition to widen up (as per Tzvetan’s suggestion) the epoch to avoid filter edge effects in your time window of interest.
Note that you ask for a very narrow passband on very short data segments. This is in principle not really meaningful, since narrow filters require long filterorders (i.e. they need to integrate the original signal in a long window). As long as they are integrating from the left only, the phase estimates won’t be affected by the stimulus induced transient.


Best wishes,

Jan-Mathijs


On 18 Jun 2020, at 14:44, george michael <gmichael_ at hotmail.com<mailto:gmichael_ at hotmail.com>> wrote:

Dear Tzvetan,

thanks you for the quick reply.
Wouldn't be a problem to use "onepass-zerophase" or "twopass" given that I am interested in phase at stimulus onset?
Don't both filter-types smear poststimulus activity into the prestimulus interval?

Best,
Georgios
________________________________
From: fieldtrip <fieldtrip-bounces at science.ru.nl<mailto:fieldtrip-bounces at science.ru.nl>> on behalf of Tzvetan Popov <tzvetan.popov at uni-konstanz.de<mailto:tzvetan.popov at uni-konstanz.de>>
Sent: Thursday, June 18, 2020 14:27
To: FieldTrip discussion list <fieldtrip at science.ru.nl<mailto:fieldtrip at science.ru.nl>>
Subject: Re: [FieldTrip] single-trial phase at stimulus onset using Hilbert transform

Dear Georgios,

the reason is your cfg.bpfiltdir = ‘onepass’. This should be zero phase or twopass filter. Yet, you will still get filter artefacts at the edges. I suggest you should re-segment your data, e.g. -1000 to 500ms, filter and then re-segment to -500 to 0 again.
Good luck,
tzvetan


Am 18.06.2020 um 14:14 schrieb george michael <gmichael_ at hotmail.com<mailto:gmichael_ at hotmail.com>>:

Dear Fieldtrip Community,

I am trying to extract the alpha phase at stimulus onset from single trials.
My approach is to :
1) Band-pass filter the trial segments from -500ms to 0 (±1Hz around individual alpha frequency).
I use the filter 'firws' and filter direction 'onepass' to avoid smearing poststimulus activity into the prestimulus interval.
2) Compute Hilbert transform to extract the instantaneous phase of the bandpass signal.
3) Select the phase at the last sample (t=0) as the one representing phase at stimulus onset.

 cfg=[];
cfg.channel = 'O1';
cfg.bpfilter='yes';
cfg.bpfreq=[foi-1 foi+1];
cfg.bpfilttype= 'firws';
cfg.bpfiltdir= 'onepass';
cfg.hilbert='angle';
phase = ft_preprocessing(cfg , trials_prestim);

 for tr=1:length(phase .trial)
       phase_onset(tr,:)=phase .trial{tr}(:,end); % select phase at last sample
 end

The problem is that the band-bass signal looks quite weird at the right edge, towards t=0s (see attached figure).

Do you have any idea what would be the reason for this problem?
And which would be the optimal way to extract phase at stimulus onset?

Many thanks,
Georgios
-----------------------
Georgios Michail, MSc.
Department of Psychiatry and Psychotherapy
Charite-Universitätsmedizin Berlin



<bandpass_alpha_singlesubj.png>_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=02%7C01%7C%7C75cdf37ebb8c4629c90d08d813830d14%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637280801062055379&sdata=wX2%2FX5Q7YJ9T3jgil%2B0nWa9rdZ535Axf%2FZEqeMPHhq4%3D&reserved=0>
https://doi.org/10.1371/journal.pcbi.1002202<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=02%7C01%7C%7C75cdf37ebb8c4629c90d08d813830d14%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637280801062055379&sdata=YN7KvEKe8UwFVkTeeQV%2F0B2axKX8LtPgf50ehSehzSA%3D&reserved=0>

_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://doi.org/10.1371/journal.pcbi.1002202

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200618/7cee4b42/attachment.htm>