1. Introduction
When MEG researchers measure an MEG time series in a well magnetically shielded environment and view their measurements on a real-time monitor, they see an amplitude- and frequency-modulated rhythmic wave and make no judgement about whether the rhythmic wave is a signal or noise. MEG-detectable magnetic fields by themselves reflect macroscopic electrophysiological collective phenomena in the brain, which are generated by a large number (ten thousand or more) of spatially (including directionally) and temporally localized neurons. That is, even rhythmic MEGs are already results of not a random spatiotemporal synchronization. Phase is more important than amplitude when it comes to the synchronization of rhythmic MEGs.
1.1. Phase-related phenomena
1.1.1. Event-related desynchronization (ERD) and event-related synchronization (ERS)
Rhythmic MEGs are distinguished by frequency bands: e.g., the delta band (1-4 Hz), theta band (4-8 Hz), alpha band (8-13 Hz), beta band (13-26 Hz), and gamma band (26- Hz). There are occasionally slight differences in the above frequency ranges (e.g., separation of beta and gamma bands) and naming (e.g., mu and tau bands). In the arousal state, the alpha-band MEG is largest in amplitude. As is well-known, it increases when subjects are relaxed and has their eyes closed, and hence, it might relate to an inactive state. Alpha-band suppression is observed in event-related experiments; i.e., the alpha-band MEG decreases when subjects receive stimuli and/or executes motor tasks. During this suppression, gamma- and/or beta-band MEGs sometimes increase in amplitude. Decreases and increases in MEG activity relative to that in the resting state might be caused by populations of synchronizing neurons and are called event-related desynchronizations (ERDs) and synchronizations (ERSs), respectively. ERD is a kind of the suppression of MEG in a specific frequency band, and the ERD of the alpha band is often observed in MEG recordings. On the other hand, ERS is not as distinctly observed because of its short-lived and poor time-locked characteristics. The word ERS seems to be used in context with ERD. In fact, appearances of ERD and ERS would be much more complicated (Lopes da Silva, 2006, Pfurtscheller, 2006, Herrmann et al., 2005, Başar, 1998). ERD and ERS are often mentioned in the literature of EEG-based brain-computer interfaces (BCIs) (Dornhege et al., 2007). In particular, it has been found that well-trained BCI users can use imagery to control ERD/ERS.
The event-related field (ERF) is the MEG version of the event-related potential (ERP), the only difference being the sort of measurement. Of course, there is an intrinsic difference between EEG and MEG; i.e., only MEG can detect intracellular current flows. Many stimuli are typically presented to subjects in ERF experiments. MEG data clip out according to the triggers of stimuli.
In what follows, we denote
Equation (1), or conventional averaging, is related to the supposition that the stimulus-locked additive activity is the most important object of analysis, and this supposition has remained implicit in most of the ERP and ERF literature. Additive activities are similar to impulse responses in the system theory, and they are thought to appear only as a result of stimuli and/or motor executions. In fact, conventional averages are experiment-dependent and (mostly) subject-independent; most of the rhythmic activities are averaged out. It is thus natural to think that conventional averages would have important information and rhythmic activities would not. But is this really the case? For instance, do conventional averages converge to additive activities when the number of epochs becomes large? A much more elementary question arises as to whether additive activities actually exist or not. The conventional averages (1) are just mean values in the mathematical sense. Moreover, ‘stimulus-locked’ in the sense of conventional averaging means ‘phase-locked’ to stimuli, and additive activities are often referred as evoked activities.
ERD and ERS researchers are interested in rhythmic activities. Therefore, they often initially calculate the power-spectrograms of MEG epochs,
where WT{ } means a continuous wavelet transform, the mother wavelet of which is the complex Morlet (modulated Gaussian) (Tallon-Baudry et al., 1999). The wavelet transform can be replaced by other time-frequency methods (e.g., the Wigner distribution or the short-time Fourier transform) and the Hilbert transform (frequency-wise). The power calculation of eq. (2) is a simple but nonlinear operation to eliminate stimulus-unlocked phase differences. Therefore, power-spectrogram averaging (2) can reflect ERD and ERS. Note that power-spectrogram averages are usually frequency-wise normalized by the mean power for a pre-trigger duration to enhance the ERD (ERS) changes. Power-spectrogram averaging indeed has played an important role in making researchers notice that a portion of the rhythmic activities are actually event-related. However, power-spectrogram averages may also contain additive activities, if they actually exist. Moreover, if the phases of ERFs are also event-related, another important piece of information regarding the phases in the brain processing is unfortunately lost. ‘Stimulus-locked’ in power-spectrogram averaging means only ‘time-locked to stimuli’ (not care for ‘phase-locked’), and time-locked rhythmic activities are often referred to as induced activities.
Averaging methods for analysing event-related phases have been also proposed. The phase-locking factor (PLF) for a single channel phase synchronization analysis is expressed as
where j is the imaginary unit (Tallon-Baudry et al., 1996). PLF is the absolute value of the average of a specific-frequency
So far, we have explained three averaging methods (1)-(3) for a single channel and their related technical terms. In fact, their applicability to the problem at hand may still be somewhat confusing. In particular, the term ‘phase-locked activities’ includes in its meaning both additive activities and phase-locked rhythmic activities. Accordingly, the conventional averaging (1) picks up both of these activities, and PLF (3) also picks them up, even if the number of epochs is sufficiently large. Moreover, phase-locked activities are just subsets of time-locked activities, and hence, these two activities cannot be separated either. Power-spectrogram averages (2) include not only the above-mentioned activities but also phase-unlocked rhythmic activities.
1.1.2. Long distance synchronization and phase resetting
Why is the phase of ERF modulated if it is actually modulated? It could be that an oscillator is unilaterally drawn to phase synchronization with another oscillator or two are bilaterally drawn to their halfway phase point or a third oscillator’s phase. Long distance synchronization refers to phase synchronization between distant brain areas (Rodriguez et al., 1999). The phase synchronization of two channels would be meaningful when the distance between two phase-synchronized brain areas is large in the sense of the MEG pattern (not only their locations but also their orientations are relevant). The phase-locking value (PLV) was proposed (Lachaux et al., 1999) as a way to analyze cross-channel (between two channels
PLV is sometimes referred as the phase-locking index (PLI). PLV is apparently a cross-channel version of PLF and is also distributed from zero to one. Unlike PLF, the significance of PLV can be statistically tested with a surrogate method as follows. i) Make surrogate sets of phases by random permutation of
The single-trial phase-locking index (SPLI) For a cross-channel single epoch analysis is (Lachaux et al., 2000)
where
Phase resetting is another way of picturing phase modulation. Long-distance synchronization is based on the cross-channel phase relationship. Phase resetting refers to the phenomenon wherein the phase of an oscillator jumps to a specific value in a short period of time. Since the specific value is not determined by another running oscillator, phase resetting is not a cross-channel phenomenon. If stimuli produce phase resetting, nonzero activities will appear in the conventional averages. Probably because of this, phase-resetting has been discussed in the context of whether activities are additive or not. For instance, i) visual ERFs were divided into two bins by post-stimulus alpha amplitude and ERFs were averaged for each bin. ii) the amplitudes of the averages inherited the post-stimulus alpha-band amplitudes (Makeig et al., 2002). This simple fact could provide strong evidence for the existence of phase resetting.
There is a counterexample. The phase preservation index (PPI) is expressed by
for evaluating how much the phase of a specific pre-stimulus time
Finally, we would like to point out that i) indices (3)-(6) discard the amplitude terms, ii) they are not MEG activities, which can be further analyzed with subsequent methods, but rather indices, and iii) the cross-channel indices (4) and (5) are undirectional.
1.2. Simulation: Phase modulation
This section describes a simple simulation performed to investigate the influence of only phase-locked rhythmic activities on the three averages (1)-(3).
A 10-Hz stable sinusoid with an initial phase
One might be interested in the conventional average of
Let us now investigate whether there is any difference between using the cosine and sine in the simulation setup of the initial
2. Epoch filter
ERF is a function of time, channel, and epoch (Fig. 2). There are huge variations for time, or temporal, filters (e.g., bandpass filter, wavelet transform, and Hilbert filter) and for channel, or spatial, filters (e.g., beamformer, Laplacian filter, principal component analysis (PCA), and independent component analysis (ICA)). These filters have their own benefits. Now let us turn to the question of epochs. As shown in the figure, we can consider a third filter for the epoch axis.
2.1. Concept
Let us begin with an interpretation of the three averaging methods (1)-(3) as possible epoch filters. PLF (3) is significantly different from conventional averaging (1) and power-spectrogram averaging (2). The most conspicuous difference from an epoch filter viewpoint is that PLF treats epochs as a series; i.e., instantaneous phases give epochs their own order. The instantaneous phases also serve as an order for a series for indices (4)-(6). PLF is a weighted average of ones (in the absolute value sense), to which the ERF amplitudes are normalized at all time sampling points. The weight is derived from the instantaneous phase of ERF. More precisely, it is the filter kernel of PLF (3), exp(j
The most essential reason for the temporal and spatial filters is that time and channel have own orders based on physical properties. Since time and channel have orders, they are not always assigned the same weights (filter kernels). Epoch filters would be possible if epochs can be given an order, even if this order does not have an unequivocal physical meaning. PLF has built a small bridge to epoch filters.
Filters, in general, have pass- and stop-bands. Although bands might fit in frequency very well (namely, frequency band), bands should be treated here as just subsets. Filters are designed by taking one of the two bands into account. Since the three averaging methods (1)-(3) are not originally designed as epoch filters, they do not handle pass- and stop-bands very well. There is no appropriate band concept corresponding to time-locked and phase-locked activities for epoch filters.
The question thus is can ERFs have bands? Rhythmic activities, which are sometimes referred to as ongoing spontaneous activities, always exist and are amplitude- and/or frequency-modulated. Since phase modulation can be done in place of frequency modulation, it may be fruitful to separate the rhythmic activities into amplitude-modulated (AM) and phase-modulated (PM) components. In addition to these components, additive activities (AD), which are template-like activities corresponding to individual inputs (stimuli) and outputs (motor executions), may exist. Accordingly, let us treat AD, AM, and PM as bands, and let one purpose of epoch filters be to pass one of these bands and to stop the others.
Another purpose of epoch filters is illustrated by the simulation described in section 1.2. In that simulation, the phase of
Let us introduce a quantity that gives epochs an order. For the above-mentioned two purposes, phase is surely important. Therefore, as the first attempt to design an epoch filter, let us use the instantaneous phase to give an order to epochs, and let us discuss independence of the instantaneous phase as an aside. An epoch is defined as the integer number of stimuli/motor executions that make up an ERF recording. Therefore, an epoch is definitely an independent variable. Although instantaneous-phase-sorted epochs are, however, distinguished by the original numbers, they differ by the other independent variables: time and channel.
Finally, we have a perhaps subjective feeling about filters. We wonder if the inputs and outputs of most of the existing filters might belong to the same category. This means the outputs of epoch filters should not be indices and should be spatiotemporal-structure-preserved ERFs. In other words, the outputs of epoch filters can be used as inputs to spatial and temporal filters.
2.2. Preliminary: Temporal filtering
Before explaining the methodologies of epoch filtering, we should mention frequency-selective temporal filtering.
The instantaneous phase will play an important role in ordering epoch filters. To determine the instantaneous phase, we often perform a discrete Hilbert transform (Oppenheim & Schafer, 2007), which retains the real part of the sequence as is and only generates the imaginary part of the sequence under the constraint that the total complex sequence consists of only positive frequency terms; namely, the total complex sequence is analytic. The discrete Hilbert transform consists of i) a discrete Fourier transform, ii) doubling of positive frequency terms (except the DC and Nyquist frequency terms) and zero-padding of negative frequency terms, and iii) an inverse discrete Fourier transform. Therefore, we usually embed temporal filtering in the Hilbert transform by padding zeros to stop-frequency-band terms. Even in the case that only temporal filtering is necessary and an analytic sequence is not, we should perform the pair-wise forward and inverse discrete Fourier transforms and halfway stop-band zero padding for consistency. The need for consistency relates to the fact that the discrete Fourier transform treats a sequence as periodic. In other words, we should not cross-process two sequences, one of which is generated by a convolution (not a periodic procedure) and the other of which is generated by a Hilbert transform (a periodic procedure). Instead of using a Hilbert transform, one can perform a wavelet transform and convolution-based temporal filtering and maintain consistency. However, in so doing, one may encounter difficulty in obtaining a wide-frequency-band-limited analytic sequence. In addition, the wavelet transform with the complex Morlet also generates an analytic sequence unless one takes care about the negative frequency terms for the complex Morlet.
2.3. Phase-interpolated averaging
Phase-interpolated averaging, which is an epoch filter, has been proposed for separately extracting AD, AM, and PM and for reducing noise in the averaging sense (Matani et al., 2011). It is originated from the phase-series analysis for echography (Matani et al., 2006).
2.3.1. Methodology
Time and channel, which are independent variables of ERF, are discrete samples of continuous variables in a computer; needless to say about time, channels are actually input coils and hence are spatial samples. On the other hand, an epoch is ordered according to the instantaneous phase. Phase is continuous, and therefore, an epoch can be virtually treated as a sample. However, phase sampling intrinsically has unavoidable jitter. Another characteristic is that phase has a period of 2π, or more precisely, has an uncertainty of integral multiples of 2π. Phase-interpolated averaging consists of two stages: I) interpolation of ERF by eliminating phase sampling jitter and II) a discrete Fourier transform of the interpolated ERF over the epoch.
Stage I) Interpolation of ERFLet us recall ERF, and suppose it is first frequency-band
where
The pair of +
Since phase is intrinsically periodic,
The superscript ‘G’ means grid-sampling.
The periodic version of the sampling theory is expressed as
where
is the Dirichlet kernel of the 2
Therefore, the jittered and unjittered (grid-on) phase-sample interface can be formulated as
where
The superscript ‘J’ means jittered-sampling. Note that both
The interpolated ERF is then discrete Fourier transformed. However, so far, we need only the 0th and 1st terms:
They are respectively called the 0th and 1st phase-interpolated averages.
As will be shown in the next simulation, the 0th and 1st phase-interpolated averages are to separately extract the AD activity and the AM component, respectively. The condition number of Ξ
Finally, we would like to stress that phase-interpolated averaging is performed for each time sampling point and each channel (not for cross-channel).
2.3.2. Simulation: Separation ability
In the simulation in section 1.2, only PM occurred in the middle duration (0.1<
The condition number monotonically increases with the number of interpolated epochs, 2
The above-mentioned stability depends on not only the instantaneous phase bias but also the amplitude of AD. This simulation was of a difficult case in which the amplitudes of AD and the rhythmic activities (including AM) were comparable. As mentioned in section 3, also as the common experience, AD (if it actually exists) is usually smaller than the rhythmic activities. If this is not the case, no one tries to perform conventional averaging. Thus, we can say that actual phase interpolation is not so difficult.
One may wonder about strange surges at the ends of the 1st phase-interpolated averages and the condition numbers. This is a transient phenomenon caused by the discrete Hilbert transform, which treats the objective series as periodic. To avoid it, we can set the time duration of interest to the center before the analysis and cut off both ends after the analysis. (We would have done so if this article were for journal papers.) Although we think this treatment is sufficient, such transients can in any event be avoided with Hamming-windowing.
Consequently, phase-interpolated averaging works as an epoch filter to pass only one of AD, AM, and PM bands and to stop the others.
2.4. Phase-compensated averaging
Phase-compensated averaging, which is an epoch filter, has been proposed for visualizing directional phase synchronization between channels and for reducing noise in the averaging sense (Matani et al., 2010).
2.4.1. Methodology
Phase-compensated averaging consists of filter kernel design and actual filtering.
Regarding the filter kernel design, let a narrow-frequency-band
The filter kernel of phase-compensated averaging is
Regarding the actual filtering, we can use a wide-frequency-band
The phase-compensated averaging is expressed as
Thus,
Phase-compensated averages are complex valued. They should be represented by their absolute values and phases (if necessary) rather than their real and imaginary parts. This is because they usually appear as modulated complex sinusoids having a single center frequency (they look like a synchronization burst), and therefore, the absolute value representations contain the envelopes of the phase-compensated averages, which are the most useful for characterization.
Regarding the significance test, we shall fundamentally adhere to the statistical test of PLV (Lachaux et al., 1999). A surrogate method is used for evaluating the significance, as follows. i) Surrogates are generated by random permutation of
Finally, we would like to point out a number of methodological differences between phase-interpolated averaging and phase-compensated averaging: i) The filter kernel of the former is calculated for each time sampling point, whereas that of the latter is calculated for only a reference time sampling point
2.4.2. Simulation: Directional synchronization
We will again refer to the ERFs
Phase-compensated averaging of the given ERF is performed for reference times
The significance levels (
The auto-averages of channel 1 are always a constant value because channel 1 has an independent oscillation (dashed lines in right column, Fig. 4). Note that the phases of channel 1 over
These results show that phase-compensated averaging works as an epoch filter to pass a phase-synchronized band and to stop the other bands.
3. Actual ERF analysis with epoch filters
We shall explain the practical issues affecting the use of epoch filters by illustrating an analysis of the ERF of an actual semantic priming experiment (Ihara et al., 2007). In typical semantic priming experiments, a pair of words (prime and target) is sequentially presented to a subject and a task (e.g., button-pressing) is required to be performed only on the target word. Reaction time differs depending on whether or not there is a semantic relationship between the prime and target words. ERF sometimes is recorded during priming experiments as a physiological index. Now, although the objective of this priming experiment originally addressed a somewhat complicated issue, or selection of a meaning of a target word having semantic ambiguities through a context with the paired prime, it included a typical semantic priming experiment as a control experiment. Hence, for the sake of this discussion, we will illustrate only the epoch filtering of the ERF from the typical portion.
3.1. Semantic priming ERF
The main topic of this chapter is to introduce epoch filters. Here, we will briefly explain the procedures of a sematic priming experiment and the reconstruction of the signal source activities. Such a reconstruction is fraught with problems, due to the ill-posedness of the MEG inverse problem. Although the problems affect the overall performance of epoch filtering, we would believe that this matter should be dealt with as a separate problem from the one at hand.
All the prime and target words are Japanese nouns written with two morphograms (kanji) and three to five syllabograms (kana), respectively (left-panel, Fig. 5). The target words are either semantically related (related condition) or unrelated (unrelated condition) to the prime words. One hundred word pairs for each condition are visually presented in random order with a stimulus onset asynchrony of 1000 ms; the duration of each word is 300 ms. We ask subjects to judge whether the prime and target words are semantically related or not and to press a button after the delayed cue.
MEG in a frequency bandwidth of 0.1-200 Hz is recorded with a sampling frequency of 678 Hz using a 148-channel whole-head system (Magnes 2500WH, BTi). ERFs occasioned upon miss-pressed buttons and eye-blink artifacts are excluded from each condition. From our experience, eye-blinks often cause phase-resets, and therefore, these artifacts must be excluded when analyzing phase-related phenomena.
Dipole activities are tentatively estimated by means of the selective minimum-norm (SMN) method, which obtains minimum L1-norm solutions (Matsuura and Okabe, 1996). The SMN is applied to the conventional averages for a duration of 0-700 ms (target onset: 0 s) for all conditions. All dipoles, the amplitude of which exceeded the noise level of five times the root mean squares of the magnetic fields for a pre-trigger duration of −0.2-0 s, are gathered, neglecting the neighbors within 20 mm (Fujimaki et al., 2002). Thus, approximately 60 dipoles, which span most of the conventional average space, are selected.
The current distribution forms a vector field, and thereby, each discretized site has three orientations. The target time series of the epoch filters should have its related orientation fixed, e.g., parallel to the axons of pyramidal cells in the site. A simple way to fix the orientation with only ERF is to perform a PCA of the three dipole activities for each site. The projection vector corresponding to the first principal component is used for the orientation-fixed dipoles; i.e., three orientations of each site are projected onto the eigenvector corresponding to the largest eigenvalue of the covariance matrix of the three dipole activities in the site.
Up to this point in the experiment, only the dipoles to be fitted have been selected. In the methodological sense, a leadfield matrix, which is based on a real head shape with ASA (ANT Software B.V.), is determined. The epoch-wise (orientation-fixed) dipole activities, which are objects of the epoch filters, are subsequently estimated with least squares fitting.
We shall describe epoch filtering for a typical subject. Fifty-one orientation-fixed dipoles are used for building the leadfield matrix and 12 dipoles (Dipoles 1-12) out of the 51 are selected. Note that these 12 dipole activities showed conditional differences (Ihara et al., 2007). Each of these 12 dipoles belongs to one of six regions of interest (ROIs) (right-side, Fig. 5); the left anterior inferior frontal cortex (LaIFC: Dipoles 1 and 2), the left posterior inferior frontal cortex (LpIFC: Dipole 3), the anterior medial temporal lobe (LaMedT: Dipole 4), the left anterior middle and inferior temporal areas (LaT: Dipoles 5-7), and the left posterior superior temporal and the inferior parietal areas (LpST1: Dipoles 8-11, LpST2: Dipole 12).
The dipole activity version of the ERF is epoch filtered from the next subsection.
3.2. Phase-interpolated averages
First, we will show how phase interpolation works. As a typical example, we will choose the source activity of LaT (Dipole 6) for the related condition. The actual dipole activities for the first seven epochs show rhythmic activities, and the conventional average of the all 98 epochs seem to run through the approximate center (top row, Fig. 6). It is of interest to determine whether the wave (spike is not clear for LaT) of the conventional average come from AD or PM. As shown in the simulation in section 1.2, even only PM can generate nonzero activity in the conventional average. While the instantaneous phases at 0 ms are almost uniformly distributed, those at 200 and 400 ms are clearly biased (second row, Fig. 6). Therefore, the dipole activities seem to at least partially come from PM. The phase interpolation is performed for
Next, we perform phase-interpolated averaging (
The 1st phase-interpolated averages of the alpha-band-limited LpIFC and LpST1 dipole activities show alpha-band ERD (top-row, Fig 8). On the other hand, those of the gamma-band-limited LpIFC and LpST1 dipole activities show intermittent amplitude increments (bottom-row, Fig 8). These results might be gamma-band ERS. Since gamma-band ERS is generally not much more time-locked than alpha-band ERD, phase-interpolated averaging only vaguely visualizes gamma-band ERS. Thus, although 1st phase-interpolated averages of the wide-frequency-band limited ERF give almost no information, those of the narrow-frequency-band limited ERF are at least useful for extracting alpha-band ERD.
3.3. Phase-compensated averages
Unlike phase-interpolated averaging, phase-compensated averaging generates many auto- and cross-averages. Therefore, an overview of the all averages is extremely important and should be performed first. Furthermore, non-significant averages, in the surrogate test sense, should be masked in the overview. For this purpose, we would recommend using image-like matrix (
We cannot deny that the above-mentioned procedure is somewhat paradoxical. Our approach ends up finding a small but significant difference in almost the same two averages. Only the significance of the differences may be important; the significances of the individuals are not so important. On the other hand, we would like to narrow down the number of significant averages. It might be difficult for researchers, in particular, those having a hypothesis bias, to fairly compare a huge number of data on the same basis. Individual researchers should thus decide whether the auto- and cross-averages for
Next, let us investigate the phase-compensated averages (
First, regarding the LaIFC and LaT pair, the auto- and cross-averages (top-left four panels, Fig. 10) and their differences (bottom-left four panels, Fig. 10) for both conditions significantly appear around the reference time (vertical gray lines), but the cross-averages for the unrelated condition (thick lines) are very small. Therefore, the phase-synchronization between LaIFC and LaT might occur especially for the related condition (thin lines) although the local phase-synchronizations within LaIFC and LaT might occur for whatever the condition.
Let us discuss the related condition by inspecting the panels in a row-wise fashion. Suppose the current time is the reference time. The auto-average of LaIFC reaches a significant peak at the current time (no temporal lag), and the cross-average of LaIFC|LaT does so before the current time. Thus, the current phase of LaIFC is up-to-date whereas that of LaT is older. On the other hand, the auto-average of LaT reaches a significant peak at the reference time, and the inverted cross-average (LaT|LaIFC) does so after the current time. Thus, the current phase of LaT is up-to-date whereas that of LaIFC is younger. This row-wise comparison indicates that the phase of LaIFC precedes that of LaT. Therefore, LaT might be drawn to phase synchronization with LaIFC around the reference time. Such temporal comparisons of the peaks are relatively easy to perform when the auto- and cross-averages are all single-peaked.
We should point out that only the auto-average of LaIFC for the unrelated condition consists of significant intermittent bursts, which accordingly yields a significant conditional difference (top-left, bottom-left four panels, Fig. 10). For instance, if one is interested in which dipoles are drawn to phase synchronization with Dipole 1 (at LaIFC) for unrelated condition, the cross-averages between Dipole 1 and Dipole 8 should be subsequently evaluated by consulting the corresponding column of Dipole 1 (Unrelated panels, Fig. 9).
Regarding the LpIFC and LpST1 pair, the auto- and cross-averages (top-right four panels, Fig. 10) and their differences (bottom-right four panels, Fig. 10) for both conditions significantly appear around the reference time (vertical gray lines). There is a consistent tendency in their differences that the auto- and cross-averages of LpIFC and LpST1 for the related condition are greater than those for the unrelated condition just before the reference time and the amplitude relation is inverted just after the reference time (bottom-right four panels, Fig. 10).
Let us examine the results for the related condition in more detail. LpIFC might be drawn to phase synchronization with LpST1 for the same reason as the LaIFC and LaT pair. However, for the unrelated condition, the cross-average of LpST1|LpIFC is double-peaked across the reference time whereas the auto-average of LpST1 has a single peak at the reference time. Therefore, phase-drawing between LpIFC and LpST1 might be bidirectional for the unrelated condition. If one is interested in such double-peaked cross averages, one may wish to investigate the cross averages between Dipole 3 (at LaIFC) and Dipole 12 by consulting the corresponding columns of Dipole 3 (Unrelated and
3.4. Discussion
We will briefly discuss the actual ERF analyses using the two epoch filters.
In the typical sematic priming ERF, the conventional average for the unrelated condition is larger (in absolute value) than that of the related condition after a certain latency these two averages branch off. In fact, this original study concluded so (Ihara et al, 2007). However, the same is not necessarily true for the phase-compensated averages (Fig. 10). The phase-interpolated averages indicate that the conventional averages might consist of AD and PM and these two might cooperatively act when they have the same polarity (Fig. 7). If there is a gap between the conditional difference between the conventional averages and the phase-compensated averages, the activity that fills in the gap might be PM. In addition, directional relationships between ROIs can be discussed with these two epoch filters. This is the most conspicuously different point from the original study.
Are these two epoch filters able to analyze phase-related phenomena? Phase-interpolated averaging can be used to visualize ERD and ERS without additive activities (Fig. 8). Long-distance directional synchronization between LaIFC and LaT and between LpIFC and LpST1 can be visualized by phase-compensated averaging. In fact, these long-distance directional synchronizations likely happen since LaIFC and LaT are connected by the uncinate fasciculus and LpIFC and LpST1 by the arcuate fasciculus. In contrast, we can only speculate about phase-resetting. When the phase is reset, by which instantaneous phase jumps to a specific phase regardless of the phase just before phase-resetting, the condition number of phase-interpolated averaging could be very large at that moment. In fact, there are spike-like shapes in the condition numbers (Fig. 8). However, it must be noted that the spikes always appear if
4. Conclusion
In this chapter, we explained the methodological and practical aspects of two epoch filters. The epoch filters are made possible by giving epochs an order. Since the order is of the instantaneous phases, these filters can be used to analyze phase-related phenomena. If another order is introduced, a different epoch filter can be designed and a different concept can be created accordingly. We hope that new epoch filters besides the ones presented here will be designed.
References
- 1.
Başar E. 1998 Springer,978-3-54064-345-6 Berlin, Germany - 2.
Dornhege G. del Millán R. Hinterberger J. Mc Farland T. D. J. Müller K. R. . 2007 The MIT Press,978-0-26204-244-4 Cambridge, US - 3.
Fujimaki N. Hayakawa T. Nielsen M. Knosche T. R. Miyauchi S. 2002 An fMRI-constrained MEG Source Analysis with Procedures for Dividing and Grouping Activation . ,17 1 324 343 ,1053-8119 - 4.
Herrmann C. S. Grigutsch M. Busch N. A. 2005 EEG Oscillations and Wavelet Analysis, In: Handy T.C. (Ed),229 259 , The MIT Press,978-0-26208-333-1 Cambridge, US - 5.
Ihara A. Hayakawa T. Wei Q. Munetsuna S. Fujimaki N. 2007 Lexical Access and Selection of Contextually Appropriate Meaning for Ambiguous Words . ,38 3 576 588 ,1053-8119 - 6.
Lachaux J. P. Rodriguez E. Le Van Quyen M. Lutz A. Martinerie J. Varela R. J. 2000 Studying Single-trials of Phase Synchronous Activity in the Brain ,10 10 2429 2439 ,0218-1274 - 7.
Lachaux J. P. Rodriguez E. Martinerie J. Varela F. J. 1999 Measuring Phase Synchrony in Brain Sginals.8 4 194 208 ,1065-9471 - 8.
Lopes da Silva F.H. 2006 Event-related Neural Activities: What about Phase ?, In: C. Neuper and W. Klimesch, (Eds.),3 27 , Elsevier,978-0-44452-183-5 Amsterdam, The Netherlands - 9.
Makeig S. Westerfield Jung. M. T. P. Enghoff S. Townsend J. Courchesne E. Sejnowski T. J. 2002 Dynamic Brain Sources of Visual Evoked Responses . ,295 690 694 ,0036-8075 - 10.
Matani A. Shigeno T. 2006 Phase Series Echography with Prior Waveform Distortion for Evaluating Posterior Waveform Distortion, 53 11 2019 2025 ,0885-3010 - 11.
Matani A. Naruse Y. Terazono Y. Fujimaki N. Murata T. 2011 Phase-interpolated Averaging for Analyzing Electroencephalography and Magnetoencephalography Epochs .58 1 71 80 ,0018-9294 - 12.
Matani A. Naruse Y. Terazono Y. Iwasaki T. Fujimaki N. Murata T. 2010 Phase-compensated Averaging for Analyzing Electroencephalography and Magnetoencephalography Epochs .57 5 1117 1123 ,0018-9294 - 13.
Matsuura K. Okabe Y. 1995 Selective Minimum-norm Solution of the Biomagnetic Inverse Problem .42 6 608 EOF 615 EOF ,0018-9294 - 14.
Mazaheri A. Jensen O. 2006 Posterior α Activity Is Not Phase-reset by Visual Stimuli.103 8 2948 2952 ,0027-8424 - 15.
Oppenheim A. V. Schafer R. W. 2007 Prentice Hall,978-0-13198-842-2 Upper Saddle River, US - 16.
Pfurtscheller G. 2006 The Cortical Activation Model (CAM) , In: C. Neuper and W. Klimesch, (Eds.),19 27 , Elsevier,978-0-44452-183-5 Amsterdam, The Netherlands - 17.
Rodriguez E. George N. Lachaux J. P. Renault B. Varela F. J. 1999 Perception’s Shadow: Long-distance Synchronization of Human Brain Activity. ,397 430 433 ,0028-0836 - 18.
Tallon-Baudry C. Bertrand O. 1999 Oscillatoiry Gamma Activity in Humans and its Role inObject Representation.3 4 151 162 ,1364-6613 - 19.
Tallon-Baudry C. Bertrand O. Delpuech C. Pernier J. 1996 Stimulus Specificity of Phase-locked and Non-phase-locked 40 Hz Visual Responses in Human. 16 13 4240 4249 ,0306-4522