arxiv:physics/ v1 [physics.bio-ph] 5 Oct 2004

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/4124v1 [physics.bio-ph] 5 Oct 24 Ryerson University, Toronto, Ontario, Canada 2 Department of Physiology, University of Toronto, Toronto, Ontario, Canada (Dated: March 3, 28) Abstract The α predominant electroencephalographic (EEG) recording of the human brain during eyes open and closed is studied using the zero-crossing time statistics. A model is presented to demonstrate and compare the key characteristics of the brain state. We found the zero-crossing time statistic is more accurate than the power spectral analysis and the detrend fluctuation analysis. Our results indicate different EEG fractal scaling in eyes closed and open for individuals capable of strong α rhythm. PACS numbers: 1

2 Recent studies showed that the electroencephalographic (EEG) measurement of the cortical layer activity of the human brain exhibits power law characteristic 1,2,3. This implies the complex brain dynamics is fractal or scale-free. In biological systems, band-limited feature often coexists in the power law environment. For example, respiration can synchronize with the cardiovascular dynamics and create the narrow-band feature in the 1/f-like power law spectrum of the heart rate 4. The similar situation in brain dynamics is the appearance of α rhythm that creates the narrow-band feature in the 8 12 Hz band of the EEG power law spectrum 1,2,3. Because of its implication in the cognitive process and brain functioning 5, a better characterization of the fractal dynamics in α predominant brain state is important. For low α intensity, the power law background of the EEG is still discernable. However, for strong α intensity, a situation that is capable of by a trained meditator or a Yoga master, the power law background is no longer evident. Fig. 1 shows the EEG from healthy subjects showing moderate and strong α rhythm. It is seen that the EEG power law background becomes ambiguous with increasing α intensity. In addition, the detrend fluctuation analysis 6 generates only spurious scaling (see below). In this work, we propose to use the zero-crossing time of the EEG to study the background fractal scaling in the α predominant brain state. Our idea is based on the observation of intermittent bursts of α period from the EEG record. It implies the possible fractal background in α predominant brain state may be revealed by subtracting the α component. 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. The successive crossing-time-interval (CTI) is obtained by τ i = t i+1 t i. Let C = {τ i }. We define the subset of successive zero-crossing A = {τ k ; τ k I A } C where I A is a real interval. If I A spans a small range, A captures the zero-crossing of the narrow-band dynamics. For example, I A = [1/12, 1/8].5 (SEC) for the α frequency band. Note, the factor.5 is necessary since the wave form crosses x(t) = twice per cycle. The observed intermittent α periods in EEG implies A consists of disjoint subsets A = A 1 A 2 where A i = {τ i1, τ i2, } and i m+1 = i m + 1. The CTI for the fractal process is known to follow a power law distribution 7 : p(τ) τ ν, where p(τ) is the probability density function. For example, ν = h 2 for fractional Brownian motion B h (t) where h is the Hurst exponent. Intuitively, if A can be related to the α component of the EEG, the set difference C\A should contain the fractal crossing that can be characterized by its power law distribution. However, there are difficulties to apply this idea directly in practice. First, the zero-crossing time t i, i = 1,, is only determined 2

3 approximately (by linear interpolation). Hence, A is not exact. Second, there is evidence of power law scaling in the amplitude fluctuation of the α band-passed EEG 8. This wide-band feature suggests the CTI of the α component can span a considerable range outside the 8 12 Hz band. Thus, to define A for the α component, the interval I A should cover a much larger range. The set C\A so defined will capture only those zero-crossings in the neighborhood of very large or small τ. These zero-crossing events are likely located in the fractal period of the signal. To test the above idea, we generated synthetic fractional Brownian motion B h (t) of h =.3 and.8. Based on the reported scaling in α band-passed EEG 8, we focused on the amplitude process A h (t) defined by the absolute value of the Hilbert transformed 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 let I A = [min(τ i ) 5, max(τ i ).5], τ i C, to define the successive zero-crossing in A. Both theoretical ν values are confirmed before and after subtracting the set A (Fig. 2). This should be the case since A h (t) is scale free and no band-limited component exists in A h (t). We next construct the model for α predominant EEG to test our method. It is necessary to clarify that this model was designed to capture only the signal feature of the EEG without making any reference to the physiological process of the α dynamics. Let the synthetic EEG be y(t). We first generate N random numbers m k, k = 1,, N, as the number of samples in N intervals of y(t). So the size of y(t) is m k. The construction of the α component in y(t) is achieved in two steps. First, a pure α tone is built with the central frequency of 1 Hz: x α (t) = sin[2π(1(1+n f ))t](1+n a ) where N f is a Gaussian variable of mean zero and variance σ 2 f and N a is an uniform random variable in [,.5]. The α component of y(t) is obtained after multiplying the amplitude process A h (t): y α (t) = x α (t)(ca h (t) + 1) where C is a constant. Finally, the value of y(t) in the ith interval (t [ i 1 k m k + 1, i k m k ]) is determined from an uniformly distributed random variable V in [,1]: y(t) = y α (t), if V > p, or y(t) = A h (t) if V p, where p is a real number. A h (t) is further re-scaled to the same data range as y α (t) to minimize artificiality when V p. To mimic the α predominant brain state, a small p =.15 is used. Remaining parameters used in the construction are N = 3, C = 1, σ f =.2, h =.8, and m k, an uniform random variable in [, 2] when V > p and in [, 15] when V p, For comparison, we also calculated the spectral density function, S(f) = 3

4 y(t) exp(2iπft)dt 2, and the detrend fluctuation 6 F(l) = 1 N [z(t) z l (t)] N 2 t where z(t) = y(t) ȳ, ȳ being the mean value of y(t), and z l (t) is the linear trend of z(t) in the data window of length l. For a fractal process 6, S(f) f β, F(l) l γ and 2γ = 1+β. Fig. 3a shows a segment of the synthetic y(t). Figs. 3b and 3c show the S(f) and F(l) of y(t), respectively. The S(f) captures the narrow-band feature due to y α (t) and shows an ambiguous power law background. The detrend fluctuation F(l) has a knee located at the α period: l K.1 (SEC). For l < l K, F(l) is an increasing function of l, capturing mainly the detrend fluctuation of y α (t). For l l K, the detrend fluctuation of y α (t) approaches a constant since its local trend z l (t). Hence, for l l K, F(l) is dominated by the random components in y(t): N f, N a and A h (t). A power law property can be ascertained in this range with the estimated scaling exponent γ.5. It implies the uncorrelated random components in y(t) result in the fluctuation of the random-walk. While the fractal in y(t) cannot be clearly revealed from S(f) and F(l), the crossing-time statistic accurately captures the scaling property of the fractal process A.8 (t) in y(t). The CTI of y(t) is shown in Fig. 4a. The densely populated region τ i.5 (SEC) is contributed by y α (t). Fig. 4b shows the p(τ) of y(t), y α (t), A.8 (t) and the p(τ) determined from C\A. It is seen that CTI of y α (t) not only covers its designed frequency range 1(1 + N f ) but also follows a power law distribution over a much wider range due to the fractal component A.8 (t). We let I A = [min(τ i ) 5, max(τ i ).5], τ i C, to define A. The p(τ) s before and after deleting A are qualitatively different (Fig. 4b). After deleting A, the fractal component A.8 (t) is revealed in the power law p(τ) τ ν with the estimated ν very close to the theoretical value ν.8 2 = 1.2. We also found consistent results for other h values, indicating the robustness of the method (Fig. 4c). This example illustrates (a) the CTI of a process showing both fractal and band-limited characteristics can be significantly different from the CTI consisting of only the fractal component, and (b) the fractal component can be captured from C\A with properly defined A. We now apply the crossing time analysis to the EEG from seven healthy subjects in eyes open (EO) and closed (EC) [gender: 4 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 Fz. 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 4

5 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 three subjects). In order to compare 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. The R α of the data lies between low α power range of.5 to the well above-normal α power range of.7; (Fig. 5). Also, R α is always higher in EC than in EO. Three of the seven subjects (subjects 2 4) were able to show well above-normal α rhythm with R α >.45 in EC. The EEG with strongest α intensity has been shown in Fig. 1 and is reproduced in Fig. 3 to compare with the synthetic y(t) generated above. In addition to the similar appearance between the two, the S(f) and F(l) of the real EEG are also close to its artificial counterparts (Figs 3b, 3c): both EEG and y(t) show ambiguous power law background in S(f) and the knee in F(l) at l.1 (SEC) with an estimated scaling exponent γ.5. The CTI distribution of all EEG data exhibits power law p(τ) τ ν, suggesting the fractal dynamics exist in both moderate and α predominant brain states. Qualitative difference is seen in p(τ) before and after deleting the set A (Fig. 6). Of all the subjects, the ν appears to be smaller in magnitude in EC than EO (Fig. 7). According to ν = h 2, this implies EC has a larger Hurst exponent than EO. Moreover, the magnitude of ν appears to be inversely proportional to the α intensity (Fig. 7). This is found in EO versus EC of the entire group and in subjects 1,5,6,7, which show particularly strong α rhythm. In these cases, the large ν is seen to associate with small R α. Finally, differences of ν of individuals in EO and EC also appear to vary proportionally with the α intensity. In summary, the technique to extract scaling property in EEG showing band-limited characteristic is proposed and successfully demonstrated. Using this technique, the scaling in the brain dynamics in α predominant EEG is obtained. Our results suggest that the difference in the fractal dynamics is closely related to individual s capability of strong α rhythm. Further study is necessary to analyze the statistics from larger popuation with strong α capability. Acknowledgment The authors would like to acknowledgment supports from Natural Science and Engineering Research Council of Canada. 5

6 Reference [1] Y. georgelin, L. Poupard, R. Sartene, J.C. Wallet, Eur. Phys. J. B, 12, 33 (1999). [2] P.A. Watters, Int. J. Sys. Sci., 31, 819 (2). [3] R.C. Hwa, T.C. Ferree, Phys. Rev. E., 66, 2191 (22). [4] J.P. Saul, News Physiol. Sci., 5, 32 (199); A. Malliani, M. Pagani, F. Lombardi, S. Cerutti, Circulation, 84, 482 (1991). B.F. Womack, IEEE Trans Biomed. Eng., 18, 399 (1971); D. Cysarz, D. von Bonin, H. Lackner, P. Heusser, M. Moser, H. Bettermann, Am. J. Physiol. Heart Circ. Physiol., 287, H579 (24). [5] C. Neuper, G. Pfurtscheller, Int. J. Psychophy., 43, 41 (21). [6] C.-K. Peng, S. Havlin, H.E. Stanley, A.L. Goldberger, Chaos, 5, 82 (1995). [7] M. Ding, W. Yang, Phys. Rev. E, 52, 27 (1995). [8] K. Linkenkaer-Hansen, V.V. Nikulin, J.M. Palva, K. Kaila, R. J. Ilmoniemi, Eur. J. Neurosci., 19, 23 (24). 6

7 Figure Caption Fig. 1 Typical EEG with moderate (top curve) and strong (bottom curve) α intensity: (a) EEG record, (b) power spectral densities S(f) and (c) detrend fluctuation F(l) of EEG s shown in (a). The dotted line in (b) is located at log(1) 2.3 Hz. In (c), the estimated long-term scaling exponents are 1.1 and.5 for moderate and strong α rhythm, respectively (solid lines are shown with the given slopes). Fig. 2 (a) CTI s of A.3 (t) (top curve) and A.8 (t) (bottom curve). The CTI of A.3 (t) is approximately 1/8 the vertical range of A.8 (t) due to the persistence correlation in the latter. (b) log(p(τ)) versus log(τ) of the CTI for A.3 (t) (top curves) and A.8 (t) (bottom curves) before (in open circle) and after (in cross) deleting the successive zero-crossing set A. Constants are added to log(p(τ)) to show identical scaling. The solid lines are drawn with the theoretical slope 1.2(=.8 2) and 1.7(=.3 2). Fig. 3 Comparison of synthetic EEG y(t) and real EEG x(t): (a) time series, (b) power spectral density S(f) and (c) detrend fluctuation F(l). Results from y(t) is always put at the top and x(t) at the bottom. In (b), dotted line is drawn at log(1) 2.3 Hz and the cross is located at the 6 Hz power line frequency: log(6) 4.1 Hz. In (c), the estimated long-term scaling exponent is.5 for both synthetic y(t) and x(t) (solid lines are given with the slope.5). Fig. 4 (a) CTI of y(t). Densely populated τ.5 sec. is contributed by y α (t) (see text). (b) Power law p(τ) of y(t) (solid line), y α (t) ( ), A.8 (t) ( + ) and power law p(τ) determined from C\A where A = {min(τ i ) 5 τ i max(τ i ).5} (dash line). The straight line has a slope 1.2. The 1-Hz line, τ = log(.5) is shown as the vertical dotted line. (c) Power law p(τ) determined for h =.2,.4,.6,.8 (bottom to top) after deleting the set A = {exp( 5.5) τ k exp( 3.5)}. The straight lines have the slope 1.8, 1.6, 1.4, 1.2, respectively (bottom to top). Fig. 5 (a) R α of the seven subjects: solid (empty) bars correspond to EC (EO). (b) S(f) of subject 4 (sbj4) during EC. (c) S(f) of subject 7 (sbj7) during EC. The dotted lines in (b) and (c) are located at 1 Hz (log(1) 2.3). Fig. 6 (a) CTI of the typical α predominant EEG (same as Fig. 1). Densely populated τ.5 (SEC) is due to the α rhythm. (b) Power law p(τ) before (thin solid line) and after (thick solid line) deleting the set A. The solid line has a slope of

8 Fig. 7 R α versus ν for the seven sujects in EC and EO. The subjects 2, 3, 4 showing strongest α intensity are drawn in solid (EC) and open (EO) squares. Subject 7 showing the weakest α rhythm is drawn in solid (EC) and open (EO) triangles. Subjects 1,5,6 are drawn in solid (EC) and open (EO) circles. 8

9 (a) (b) 6 4 (c) log(s(f)) log(f(l)) 2 1 sec. 2 4 log(f) log(l) Fig. 1

10 (a) 6 (b) log(p(τ)) log(τ) Fig. 2

11 (a) synthetic EEG y(t) 1 sec. real EEG Fig. 3a

12 1 8 (b) 9 8 (c) 6 7 log(s(f)) 4 2 log(f(l)) log(f) log(l) Fig. 3bc

13 .2 (a) τ # (crossing event) Fig 4a

14 7 (b) 5 log(p(τ)) Fig. 4b log(τ)

15 1 (c) 8 log(p(τ)) log(τ) Fig. 4c

16 1.8 (a) EC EO log(s(f)) 5 (b) log(f) R α.4 4 (c).2 log(s(f)) sbj1 sbj2 sbj3 sbj4 sbj5 sbj6 sbj log(f) Fig. 5

17 .2 6 (a) 5 4 (b) τ # (crossing event) Fig. 6

18 ν Fig. 7 R α