[FieldTrip] single-trial phase at stimulus onset using Hilbert transform
Tzvetan Popov
tzvetan.popov at uni-konstanz.de
Thu Jun 18 15:00:43 CEST 2020
Using filters is always a “problem”. A symmetric filter will at least not misplace your signal of interest. My reply referred to the “weirdness at the right edge”.
Best
tzvetan
> Am 18.06.2020 um 14:44 schrieb george michael <gmichael_ at hotmail.com>:
>
> 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> on behalf of Tzvetan Popov <tzvetan.popov at uni-konstanz.de>
> Sent: Thursday, June 18, 2020 14:27
> To: FieldTrip discussion list <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://mailman.science.ru.nl/mailman/listinfo/fieldtrip>
> https://doi.org/10.1371/journal.pcbi.1002202 <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/aa7afcf3/attachment.htm>
More information about the fieldtrip
mailing list