arxiv:physics/ v3 [physics.bio-ph] 19 Jan 2005

Transcription

1 APS Scaling Property in the α Predominant EEG D.C. Lin 1, A. Sharif 1,2, H.C. Kwan 2 1 Department of Mechanical and Industrial Engineering, arxiv:physics/4124v3 [physics.bio-ph] 19 Jan 25 Ryerson University, Toronto, Ontario, Canada 2 Department of Physiology, University of Toronto, Toronto, Ontario, Canada (Dated: March 29, 28) Abstract The α predominant electroencephalographic (EEG) recording of the human brain during eyesclosed and -open is studied using the zero-crossing time statistics. We found evidence of fractal characteristics which is otherwise ambiguous from the power spectrum or detrend fluctuation analysis. Our results indicate a reverse relationship between the degree of fractal fluctuation and the α rhythm intensity. Implications to the α brain state are discussed. PACS numbers: 1

2 I. INTRODUCTION The fluctuation in biological signals has been known to exhibit fractal characteristics 1. Band-limited components can often coexist in such a fractal environment. They can sometimes contribute to a large portion of the signal power that the potential fractal background can no longer be clearly identified. This is found in the electroencephalographic measurement (EEG) of the cortical layer activity of the human brain 2. In normal condition, EEG is known to consist of spontaneously generated band-limited oscillations in characteristic frequency band and a broad-band background fluctuation. The band-limited oscillations are presumed to result from synchronized synaptic activity of large numbers of neurons. The EEG broad-band fluctuation, whose origin and purpose remain unclear, is characterized by the power law power spectrum, suggesting fractal dynamics in the cortical layer activity 3,4. Normal individuals in wakefulness and eyes-closed can generate a very specific oscillation in the 8 12 Hz frequency band and most pronounced in the occipital area 1. This particular characteristic, known as the α rhythm, has been associated with the resting or idling state of the cortex 2. Its origin has been a subject of intense study, particulary, based on its signal characteristics 5 and actual physiology such as its metabolic and vascular correlates 6. The study of the fractal background in α rhythm did not receive as much attention until recently. Watters 3 and Hwa and Ferree 4 applied window-variance type of approach, known as the detrend fluctuation analysis (DFA) 7, and found fractal property from the EEG record showing some α activity. DFA is an appealing approach since the detrend fluctuation of the band-limited oscillation vanishes in the large time scale. However, a trained meditator or Yoga master is able to achieve concentrated relaxation 8 that the EEG signal power in a similar setting can shift in large proproportion to the 8 12 Hz α frequency band. Most interestingly, the 1/f-like power law spectrum found in normal individuals can no longer be observed. Fig. 1 shows the EEG with moderate and strong α component from a normal subject and a Yoga master, respectively. It is evident that the power law spectrum becomes ambiguous with increasing α intensity. In particular, the power law trend completely diminishes in the low frequency region of the α predominant EEG. This apparently inverse relationship between the fractal and rhythmic signal powers implies either the fractal fluctuation is buried in the α oscillation, and can thus no longer be detected from the amplitude characteristics of the EEG, or there is simply no long-range correlated fluctuation in the α predominant brain 2

3 state. Indeed, DFA on the integrated α predominant EEG reveals a power law exponent of.5, indicating uncorrelated white noise process in the background fluctuation (Fig. 1c). Although there is no theory to assert a fractal EEG in the α predominant environment, the EEG power law spectrum is clearly discernable in the individual showing moderate α activity. It is thus plausible that the α dynamics may somewhat be generated independently in the fractal environment of the brain dynamics. The purpose of this work is to use the zero-crossing property of α predominant EEG to examine its potential fractal background characteristics. The posterior α rhythm is a very interesting phenomenon of fundamental and practical importance; not only because of the inverse relationship between its fractal and rhythmic signal powers (Fig. 1), but also of its implication in cognitive processing and brain functioning 9, and the increasing evidence showing its potential link to the other autonomic function, such as the cardiovascular regulation 8,1. The zero-crossing of EEG has been used primarily for characterizing event-related frequency information 12. It is an appealing approach since amplitude fluctuation due to EEG unrelated factors, such as motion artefects, will have minimum effect on the result. For a fractal process, the zero-crossing actually captures sufficient detail to characterize its scaling property 11. Despite these apparent advantages, the use of EEG zero-crossing to study its fractal background appears scarce. Watters and Martins used the EEG zero-crossing to construct an EEG walk and applied DFA to estimate its scaling property 13. However, they did not consider EEG with predominant rhythmic component and the connection between the EEG walk and the original EEG remains illusive. Our work is different in that we analyze directly the zero-crossing of EEG and we focus on its fractal background with a coexisting band-limited component of predominant signal power. We will show, using artificial data, that a fractal process coexisting with predominant band-limiting oscillation can be effectively characterized by deleting successive zero-crossing intervals. We will then apply the method to six healthy adults with one highly experienced Yoga practitioner whose EEG was shown in Fig. 1. This paper is organized as follows. In the next section, the zero-crossing time method is introduced. In section III, different synthetic data are generated to test the proposed approach. In section IV, EEG from four healthy subjects are analyzed and their power law background fluctuation is characterized. Conclusion is given in the last section. 3

4 II. ZERO-CROSSING TIME ANALYSIS Let EEG be x(t). The zero-crossing time is the level set: {t i, x(t i ) = } where the index i registers the order of the zero crossing event. In practice, {t i } is obtained by linear interpolation and then used to define the set of crossing-time-interval (CTI) C = {τ i = t i+1 t i }. Zero-crossing of a stochastic process is a surprisingly difficult problem; see, e.g., Ref. 13. For fractal processes, thanks to its self-similarity, the CTI is known to follow a power law distribution 1 : p(τ) τ ν, where p(τ) is the probability density function (PDF). For example, ν = h 2 for the fractional Brownian motion B h (t) where h is the Hurst exponent. However, a coexisting band-limited process can result in a concentrated distribution and destruct the pattern of the power law. Closer examination of the α predominant EEG reveals periods of oscillation in the α frequency band interspersed with short patches of random fluctuation (Fig. 1). If fractal exists in α predominant EEG, it can be captured in C\A α where A α denotes the CTI of the α oscillation. However, the fractal CTI can also attain similar values. It is therefore not possible to obtain A α based solely on the value of CTI. In this work, we focus on the subset of large CTI fluctuation since it is a contradicting property to the band-limited (α) oscillation. Specifically, consider A = A 1 A 2, where A 1 = {τ u τ i τ l } {τ u τ i+1 τ l }, A 2 = {τ u τ i τ l } {τ u τ i 1 τ l } (1) and τ l, τ u are the thresholds that define the large CTI fluctuation. The elements in A represent continuous zero-crossing in [τ l, τ u ]. If τ u 1, τ l 1, A α is simply a subset of A. Hence, C\A contains the CTI s from the fractal component. Obviously, A also contains the CTI s from the fractal process. In principle, its CTI s are deleted uniformly over a significant range in A and thus will not introduce bias to the underlying power law of the PDF (see below). III. ARTIFICIAL EXAMPLES To test the above idea, we generated synthetic fractional Brownian motion (fbm) B h (t) of h =.3 and.8. Based on the reported fractal scaling in α band-passed EEG 15, we focused on the amplitude process A h (t) defined by the absolute value of the Hilbert transform of 4

5 B h (t). Note that A h (t) inherits the same scaling characteristics from B h (t). Hence, the ν exponents are 1.7 and 1.2 for h =.3 and.8, respectively. We followed the experimental setting (below) to use a sampling rate of 25 Hz to generate the synthetic data. To define A, τ l = exp( 5) and τ u = exp( 1) are used. They are determined from the range of CTI values of A h (t): τ u.8 max(c) and τ l 1/f NQ where f NQ is the Nyquist frequency. Figure 2 shows the PDF estimate of C\A. It is seen that both theoretical ν values are verified before and after deleting A (Fig. 2). This should be the case since no band-limited component exists in A h (t). It also confirms that deleting the set A of successive CTI does not affect the power law PDF for a pure fractal signal. To examine the influence of a band-limited component on the statistics of CTI, A h (t) is replaced by a narrow-band process, x α (t), in randomly selected time intervals of variable length; i.e., A h (t) x α (t) = M(t)N(t) (2) where M = 1 + A h (t) models the fractal amplitude modulation and N(t) is a narrow-band process with the central frequency at 1 Hz. The modulation function M is added in (2) based on the recent work showing fractal scaling in the amplitude process of the α bandpassed EEG 15. The narrow-band process N(t) is a sine wave of Gaussian amplitude X and frequency F: i.e., X = N(1, σ X ) and F = N(1, σ F ). To simulate the predominance of the band-limited oscillation in a fractal time series, the probability of an intervals being selected for x α (t) is four times of those for A h (t). In addition, the length of the interval assigned for x α (t) is at least three times shorter than those for A h (t). The synthetic data so constructed will hereafter be denoted by y(t) and is shown in Fig. 3a. The sets of CTI before and after deleting A are shown in Figs. 3b, 3c and their PDF s in Fig. 3d based on based on τ l = exp( 5), τ u = exp( 1). It is evident that the fractal characteristics of A h (t) is well captured after deleting the successive zero-crossing A (Fig. 3d). Note that we kept the time unit in all subsequent figures so as to make easy reference to the underlying narrow-band oscillation. One advantage of using CTI to examine the potential fractal property in biological signals is its guard against amplitude-sensitive, but EEG unrelated, artefacts. Clearly, any signal clipping from the instrument will not affect the zero-crossing statistics. Occasional muscle twitch or eyes blinking, the so-called movement artefacts, can result in low frequency drift of the EEG and affect its amplitude reading. However, the level crossing statistics of the fractal process should remain the same. We used the same signal y(t) from above (with the 5

6 1 Hz narrow-band component) and added a sinusoidal drift to demonstrate the robustness of the proposed approach: y 1 (t) = y(t) + µ sin(2πf s t) (3) for real µ, f z. Figure 4a shows the PDF p(τ) of the CTI of y 1 generated by f s = 1/6. Significant variation of p(τ) is seen to result as a function of µ. But the same power law exponent of the p(τ) can be estimated after deleting the set A using τ l = exp( 5.5) and τ u = exp( 2) (Fig. 4b). The variation in f s can create complicated zero-crossing pattern. The proposed method is effective if the artefacts contribute to the harmonics close to the central frequency of the narrow-band oscillation. To demonstrate, f s was set at 5, 7.5, 15 Hz and µ =.6. Figures 4c and 4d shows the p(τ) before and after deleting the set A using τ l exp( 5.2) and τ u = exp( 2). The correct power law exponent was again obtained. IV. CTI ANALYSIS OF α PREDOMINANT EEG We now apply the crossing time analysis to the EEG from six healthy subjects in eyes open (EO) and closed (EC) [gender: 3 males, 3 females; age: 21 3 (mean 24) year-old]. Surface scalp electrodes were attached according to the 1-2 international system at O1, O2 with reference to Cz. For EO, subjects were asked to direct their gaze at certain part of a shielded room to minimize eye movements. For EC, no specific instruction was given to the subjects other than to relax and have their eyes closed. Data recording lasted for five minutes. The EEG signal was first band-passed from.1 to 7 Hz and then digitized at 25 Hz (first four subjects) and 5 Hz (last two subjects). In order to measure the α intensity of the subjects, we used the ratio of EEG signal power in the 8 12 Hz band to the entire frequency range: R α = 12 8 S(f)df/ S(f)df. Figure 5 shows the R α range from the lowest α signal power of R α.5 to the highest α signal power of R α.7 (whose EEG has been shown in Fig. 1). Also, R α is always higher in EC than in EO. Three of the six subjects were able to generate predominant α rhythm with large R α (>.45) measure in EC. The CTI distribution of all EEG data exhibits power law p(τ) τ ν, indicating the fractal dynamics continues to exist in both moderate and α predominant brain states. Numerical results are summarized in Fig. 6. For α predominant EEG, qualitatively different p(τ) s were found before and after deleting A using τ l exp( 5) and τ u = exp( 2) (Fig. 6a). For 6

7 EEG of small R α measure, p(τ) remains almost the same after deleting A (Fig. 6b); see also Fig. 2. Three major observations can be made from our results. (A) An inverse relation between the magnitude ν and the α intensity is observed (Fig. 6c). This means, ν tends to be smaller for subjects showing large α signal power. (B) EC appears to have a smaller ν value in α predominant EEG. If we make a tentative comparison to the result ν = h 2, this implies EC has a larger Hurst exponent or more persistent fluctuation when strong α rhythm is generated. (C) The difference between ν in EC and EO, ν = ν EC ν EO, is proportional to the difference in the corresponding R α measure, R α = R α,ec R α,eo : i.e., ν R α. V. CONCLUSION In this work, EEG zero-crossing is used to characterize the fractal property in α predominant EEG. We demonstrated that deleting successive zero-crossing events, A defined in (1), can effectively eliminate the band-limited oscillation that coexists with the fractal fluctuation. Although some fractal crossing events are deleted, the power law distribution of the underlying fractal process has been consistently verified in synthetic data with added artefacts. While fractal cannot be ascertained in the α predominant EEG using amplitudebased approach, such as the power spectral density or DFA, the zero-crossing statistics confirm a power law distribution after deleting successive zero-crossing events. This implies the mechanism responsible for the EEG fractal continues to be active. Although the objective of this paper is the use of zero-crossing property to extract the fractal property of EEG, we made a number of interesting observations from the results. The fractal property found in the α predominant brain state may not be surprising based on the unsuccessful efforts to locate a single α generator using neuro-imaging techniques 6, which supports the assumption of a distributed α source, as well as a global fractal origin that underlies the background EEG fluctuation. But our data further found that this fractal background is a function of the strength of the α rhythm (Fig. 6). Specifically, the stronger is the α component, the more persistent is the fractal fluctuation. This is consistent to the finding by Stam and de Bruin 16 where a decrease in the DFA scaling exponent was observed in the α band from eyes open (no task) to eyes closed. However, these authors only analyzed the band-passed EEG and the rhythmic components were not predominant. 7

8 If rhythmic oscillation is considered as a metabolically dominant process 6,17, where the α predominant brain state is considered more energetically economical, the change of the EEG scaling property suggests energy efficiency has a signature in its underlying fractal characteristics. Moreover, it is observed that the ability to reach such efficiency is reflected in the ability to reach α predomiance (Fig. 6d). Futher work on larger population size and different physiological settings are necessary to provide more detailed statistics of the findings reported in this work. Acknowledgment The authors would like to acknowledgment supports from Natural Science and Engineering Research Council of Canada. Reference [1] see, for example, Fluctuations and noise in biological, biophysical and biomedical systems II, Ed. D. Abbott et al., SPIE Proceedings, Vol. 5467, 24. [2] H. Berger, Über das elektroenkephalogramm des menschen, Arch. Psychiat. Nervenkr. 87, (1929); ED. Adrain and BH. Matthews, The Berger rhythm: potential changes from the occipital lobes in man, Brain, 57, (1934). [3] PA. Watters, Time-invariant EEG power laws, International Journal of Systems Science, 31, (2). [4] R. Hwa and T. Ferree, Scaling properties of fluctuations in the human electroencephalogram, Phys. Rev. E, 66, 2191 (22); [5] T. Inouye, K. Shinosaki, A. Yagasaki, A. Shimizu, A spatial distribution of generators of alpha activity, Electroencephalogr. Clin. Neurophysiol., 63, (1986); PL. Nunez, BM. Wingerier, RB. Silberstein, Spatiotemporal structure of human alpha rhythms: theory, microcurrent sources, multiscale measurements and global binding of local networks, Hum. Brain Mapp., 13, (21). [6] M. Moosmann, P. Ritter, I. Krastel, A. Brink, S. Thees, F. Blankenburg, B. Taskin, H. Obrig, A. Villringer, Correlates of alpha rhythm in functional magnetic resonance imaging and near infrared spectroscopy, Neuroimage, 2, (23); RI Goldman, JM. Stern, J. Engel Jr., MS. Cohen, Simultaneous EEG and fmri of the alpha rhythm, Neuroreport, 13, (22). [7] C.-K. Peng, S. Havlin, H.E. Stanley, A.L. Goldberger, Chaos, 5, 82 (1995). 8

9 [8] T. Murata, T. Takahashi, T. Hamada, M. Omori, H. Kosaka, H. Yoshida, Y. Wada, Individual trait anxiety levels characterizing the properties of Zen meditation, Neuropsychobiol., 5, (24). [9] C. Neuper, G. Pfurtscheller, Int. J. Psychophy., 43, 41 (21). [1] see; e.g., H. Otzenberger, C. Simon, C. Gronfier, G. Brandenberger, Temporal relationship between dynamic heart rate variability and electroencephalographic activity during sleep in man, Neurosci. Lett., 229, (1997); J. ehrhart, M. Toussaint, C. Simon, C. Gronfier, R. Luthringer, G. Brandenberger, Alpha activity and cardiac correlates: three types of relationships during nocturnal sleep, 111, (2). [11] M. Ding, W. Yang, Phys. Rev. E, 52, 27 (1995). [12] B. Elizabeth et al., EEG Delta activitiy during undisturbed sleep in squirrel monkey Sleep research online, 3, , 2; RM. Rangayyan, Biomedical Signal analysis, Wiley-IEEE Press, 21. [13] PA. Watters and F. Martin, A method for estimating long-range power law correlations from the electroencephalogram, Biol. Psych., 66, (24). [14] SO. Rice, Mathematical analysis of random noise, Bell Syst. Tech. J., 24, (1945); B. Derrida, V. Hakim, R. Zeitak, Persistent spins in the linear diffusion approximation of phase ordering and zeros of stationary gaussian processes, Phys. Rev. Lett. 77, (1996). [15] K. Linkenkaer-Hansen, V.V. Nikulin, J.M. Palva, K. Kaila, R. J. Ilmoniemi, Eur. J. Neurosci., 19, 23 (24). [16] CJ. Stam, EA. de Bruin, Human Brain Mapping, 22, (24). [17] D. Chawla, ED. Lumer, KJ. Friston, The relationship between synchronization among neuronal populations and their mean activity levels, Neural Comput. 11, (1999). 9

10 Figure Caption Fig. 1 Typical EEG with moderate (top) and strong (bottom) α intensity: (a) EEG record, (b) power spectral density functions. The solid line marks the frequency 1 Hz (log(1) 2.3). (c) DFA results on the EEG with moderate and strong α rhythm are characterized by the power law exponent (solid line) of.98 and.5, respectively. Fig. 2 (a) An example of the CTI of A h (t), h =.3. (b) log(p(τ)) versus log(τ) of the CTI for A.3 (t) (top) and A.8 (t) (bottom) before (in open circle) and after (in cross) deleting the set of successive zero-crossing CTI A. Axes are arbitrary. The solid lines are drawn with the theoretical slope 1.2(=.8 2) and 1.7(=.3 2). The filled circles describe the zero-crossing PDF of a gaussian white noise, where no power law can be claimed; see Ref. 13. Fig. 3 (a) A segment of synthetic EEG y(t). (b) A segment of the set C of y(t). (c) The set C\A where A is defined by (1). Note the horizontal lines mark the levels of τ u and τ l. (d) log(p(τ)) versus log(τ) before (connected dots) and after (solid line) subtracting A. The solid line has the theoretic slope, -1.8, of the fractal component in y(t): A.2 (t). The vertical line marks the 1 Hz oscillation (τ = 1/2, log(τ) 2.9 since the wave form crosses the zero axis twice per cycle). Fig. 4 log(p(τ)) versus log(τ) of y 1 (t); see (3) in text. (a) µ =.8,.5,.2, (top to bottom) and f s = 1/6 Hz. (b) log(p(τ)) versus log(τ) after subtracting A. (c) f s = 5, 7.5, 1, 15 Hz (top to bottom) and µ =.6. (d) log(p(τ)) versus log(τ) after subtracting A. Solid lines in (a) (d) have theoretical slope (-1.8). Fig. 5 The R α measure for the six subjects. Fig. 6 log(p(τ)) versus log(τ) for subjects in eyes-closed showing (a) predominant α rhythm (sub. 4 in Fig. 5) and (b) moderate α rhythm (sub. 6 in Fig. 5). The open (close) circle corresponds to p(τ) before (after) subtracting A and the solid lines are the regression lines. (c) The inverse relationship between ν and R α. (d) The relationship ν R α. In (c), (d), the subject index (used in Fig. 5) is given next to the symbol. 1

11 (a) 1 (b) log(s(f)) 2 mv 1 sec log(f) Figure 1

12 8 slope=.98,.5 7 (c) 6 5 log(f(τ)) log(τ)

13 2 15 (a) (b) τ 1 log(p(τ)) # of crossing events 1 4 log(τ) Figure 2

14 4 (a) y(t) t.1 (b).1 (c) τ i.5 τ i i, crossing event i, crossing event

15 5 4 (d) 3 2 log(p(τ)) Figure 3d log(τ)

16 log(p(τ)) log(p(τ)) (a) (b) log(τ) log(τ) Figure 4

17 1 1 8 (c) 8 (d) 6 6 log(p(τ)) 4 2 log(p(τ)) log(τ) log(τ)

18 1 EO.8 EC.6 R α subject # Figure 5

19 log(p(τ)) (a) log(p(τ)) (b) ν o 6c 6o 2o 3o 4o 1c 5o c 2c 3c 4c Figure 6 log(τ) log(τ) R α Figure 6

20 .8.7 (d) ν R α Figure 6d