Journal of Applied Mathematics

Volume 2012 (2012), Article ID 926358, 22 pages

doi:10.1155/2012/926358

Research Article

Hyperspherical Manifold for EEG Signals of Epileptic Seizures

Tahir Ahmad and Vinod Ramachandran

Ibnu Sina Institute for Fundamental Science Studies, Universiti Teknologi Malaysia, 81310 Skudai, Johor, Malaysia

Received 3 February 2012; Revised 11 April 2012; Accepted 21 April 2012

Academic Editor: E. S. Van Vleck

Copyright © 2012 Tahir Ahmad and Vinod Ramachandran. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The mathematical modelling of EEG signals of epileptic seizures presents a challenge as seizure data is erratic, often with no visible trend. Limitations in existing models indicate a need for a generalized model that can be used to analyze seizures without the need for apriori information, whilst minimizing the loss of signal data due to smoothing. This paper utilizes measure theory to design a discrete probability measure that reformats EEG data without altering its geometric structure. An analysis of EEG data from three patients experiencing epileptic seizures is made using the developed measure, resulting in successful identification of increased potential difference in portions of the brain that correspond to physical symptoms demonstrated by the patients. A mapping then is devised to transport the measure data onto the surface of a high-dimensional manifold, enabling the analysis of seizures using directional statistics and manifold theory. The subset of seizure signals on the manifold is shown to be a topological space, verifying Ahmad's approach to use topological modelling.

1. Introduction

Epilepsy is a chronic disorder of the brain that affects people all over the world. It is characterized by recurrent seizures—physical reactions to sudden, usually brief, excessive electrical discharges in a group of brain cells [1].

Seizures are classified into two major categories, partial seizures and generalized seizures. Partial seizures are those in which the clinical or electroencephalographic evidence suggests that the attacks have a localized onset in the brain [2]. This type involves only a part of the cerebral hemisphere at seizure onset and produces symptoms in corresponding parts of the body or disturbances in some related mental functions [3]. On the contrary, generalized seizures are said to occur if the evidence suggests that the attacks were well spread [4]. The diagnosis and treatment of epilepsy is greatly aided by the use of electroencephalography as a monitoring tool.

Electroencephalography (commonly referred to by its abbreviation EEG) is the measurement of electrical activity produced by the firing of neurons in the brain. It operates by recording the fluctuations in the potential difference of electrodes attached to the scalp of a patient, with these fluctuations indicating the presence of neural activity [5]. In practice, computerized systems such as NicoletOne are used to digitize EEG signals before they are subjected to statistical analysis. It is not uncommon for the interpretation of clinical EEG to involve speculation as to the possible locations of the source(s) inside the brain which are responsible for the observed activity on the scalp [4]. The mathematical analysis of EEG signals aids medical professionals by providing a description of the brain activity being observed, thus increasing the understanding of human brain function.

2. Related Work

To date, numerous methods have been employed in the mathematical modelling of epileptic seizures, each with varying objectives and results.

In 1986, Babloyantz and Destexhe used the time series of EEG data during a seizure to show the presence of deterministic dynamics of a complex nature for seizures [6]. By identifying and comparing the chaotic attractors during (ictal) and between (interictal) seizures, they argued that an epileptic seizure is a low-dimensional state, which was further clarified by Stam as a “loss of complexity” [7]. This concept was further strengthened by Iasemidis et al. by computing the decrease of the largest Lyapunov exponent during an epileptic seizure [8], a result which was further corraborated in 1990 by Frank's analysis of the EEG of absence seizures, in which he indicated the existence of an underlying chaotic attractor [9]. However, Theiler's analysis of the same dataset showed that the supposed chaotic properties can be attributed to a noisy limit cycle [10]. Support for this opposing view can be seen in studies by Schiff et al. [11], Friedrich and Uhl [12], Hernndez et al. [13], Le van Quyen et al. [14], and Feucht et al. [15], which leads to the conclusion that the use of time series modelling is highly dependent on the quality of the EEG signals produced. In practice, a noise filter is applied prior to signal amplification to help remove white noise. However, one inherent problem with EEG recordings of seizures is the crossing of its input and output sets—all recorded signals (gross output) are a combination of seizure activity (net output), exogenous noise, and the current source used to generate a potential difference between two electrodes (net input). Hence, to accurately model a seizure, dynamical noise reduction needs to be applied to isolate the signals generated from seizure activity alone, which in turn results in degradation of the source signal.

Nevertheless, modelling of the EEG signals of epileptic seizures endured, with focus shifting to forecasting brainstorms, reasoning for which is the ability to warn and prepare treatment for an impending seizure. The ultimate goal is to produce a closed system where a patient can be connected to an automated device that predicts the onset of seizures and administers the appropriate medication as necessary, a concept that was first outlined by Peters et al. in 2001 [16]. When it came to forecasting brainstorms, nonlinear models took preference over linear models. Studies conducted by Elger and Lehnertz [17], and Martinerie et al. used intercranial EEG to show that seizures can be anticipated, sometimes up to 5 minutes in advance [18]. In 2001, Le van Quyen demonstrated that seizure prediction can be carried out with surface recordings alone [7], jumpstarting research into noninvasive models for seizure prediction. Consequently, various approaches were employed to model seizures, including Lyapunov exponents [19] (shown to be inapplicable for low-dimensional deterministic chaotic systems such as epileptic seizures by Lai et al. in 2003 [20]), correlation integrals and dimensions [21] (also shown to be inapplicable for epileptic seizures [22]), phase clustering [8, 23], and entropy measures [24]. Although some of these models were quite successful, they are very specific in nature, only being able to model a very small subset of seizures, a common flaw in nonlinear models. Additionally, Mormann et al. have highlighted that although most of the studies published in the 1990s and around the turn of the millennium yielded rather promising results, more recent evaluations could not reproduce these optimistic findings [25].

It was only after analysis of Martinerie's earlier work that McSharry et al. discovered that linear measures could also be used to model seizures [26]. The acceptability of linear models was solidified with Kugiumtzis and Larsson concluding that nonlinear methods offered no significant advantage over linear models [27]. This result was followed up in 2001 by Jerger, who compared 7 different linear and nonlinear measures and found that both classes of measures produced similar results [28]. This provided strength to existing linear models, such as Baillet and Garnero's Bayesian-based model [29]. Although Bayesian modelling proved successful, its reliance on apriori information yields a low level of accuracy for patients who have a low occurrence of seizures. Furthermore, additional data cannot be generated by inducing seizures, as the Pavlovian effect on seizure provocation is still undetermined [30]; that is, researchers are unable to identify if repeated seizure provocation results in specific brain cells being conditioned to react to stimuli, thus jeopardizing the integrity of the data obtained.

With the introduction of new mathematical tools over the last decade, more complex models were developed. In 2008, Faust et al. successfully applied Burg autoregressive coefficients in the modelling of epileptic seizures, achieving an accuracy of over 90% in detecting seizures [31]. This result was bettered by Sivasankari and Thanushkodi in 2009 using fast independent component analysis coupled with neural networks [32]. These methods, although highly successful, are of high complexity and are notoriously resource-intensive.

In 2000, Ahmad et al. formulated a fuzzy-based topological model to identify the foci of an epileptic seizure [33]. The model, called Fuzzy Topographic Topological Mapping (FTTM), involves inducing a topology on the magnetic field of magnetoencephalographic (MEG) recordings and utilizes fuzzy techniques to estimate the location of the epileptic foci [34]. The epileptic foci here refer to points in the brain from which the seizures are assumed to originate. This method was further refined in 2008 to include EEG input [3]. Apart from estimating the epileptic foci, the FTTM method is able to index key events during the progression of seizures, as demonstrated in Idris's 2010 paper [35]. Despite its successes, the model has drawn some criticism due to the loss of information incurred during the fuzzification/defuzzification process. By using the same set of patient data as Ahmad, this research aims to provide a foundation for the geometric modelling of statistical data of seizures without the use of fuzzy theory.

3. Contributions

The contributions of this paper are as follows.

(1)A discrete probability measure called the Delia measure is developed. Usage of this measure overcomes the need for apriori information without compromising the structural pattern of the inputted data. (2)Digitized raw EEG data is processed using the developed measure. Subsequent analysis shows success in highlighting specific areas of the brain which show a large change in potential difference. In the case of three patients who are known to demonstrate physical symptoms, the model successfully identifies increased potential differences in areas of the brain that correspond to the observed symptoms. (3)A mapping from the set of seizure measure values to a unit hypersphere is defined. This results in a manifold for EEG signals of seizures, thus enabling further analysis of seizures using directional statistics and manifold theory. (4)The set of EEG signals on the manifold is shown to constitute a topological space, thus verifying Ahmad's use of topological modelling.

4. Preliminaries

In this section, some existing definitions and/or results that will be used in this paper are introduced.

4.1. Probability and Measure

The first part of this paper will be the development of a probability measure, purpose-suited to measuring EEG signals detected during epileptic seizures. The given definitions below are primarily sourced from measure theory, motivation of which is given in Section 5.

Definition 4.1 (-algebra). Supose is a set and is a collection of subsets of . is called a -algebra of subsets of provided it contains the set and is closed under taking complements (with respect to ), countable unions, and countable intersections [36].

Definition 4.2 (set function). A set function is a function whose domain of definition is a nonempty class of sets [37].

Definition 4.3 (measure). A measure is defined as a set function whose domain of definition is a ring , such that

(1) is additive; (2) is positive; (3); (4)if is an increasing sequence in whose union is also in , then [37].

Definition 4.4 (measurable space). A measurable space is a pair , where is a -ring of subsets of . Elements of are said to be measurable [37].

Definition 4.5 (probability space). A probability space is a measurable space with total measure , that is if is a probability space, and is the measure defined on , then [38].

Definition 4.6 (discrete random variable). A discrete random variable is a numerical quantity that, in some experiment that involves some degree of randomness, takes one value from some discrete set of possible values [39].

Definition 4.7 (probability mass function). Let be a discrete random variable. A probability mass function for is a function which assigns a probability to each value of the random variable such that

(1) for all , and (2) [40].

Definition 4.8 (Kolmogorov-Smirnov test). The Kolmogorov-Smirnov test is designed to test versus , where and are distributions associated with two independent groups [41].

Definition 4.9 (von Mises-Fisher distribution). A -dimensional random unit vector is said to have a -variate von Mises-Fisher distribution if its probability density function is given by

where , , and . The normalizing constant is given by

where represents the modified Bessel function of the first kind with order . is known as the mean vector, and is called the concentration parameter [42].

4.2. Topology and Manifolds

This section covers some basics on topology and manifold theory that will be used to map the EEG signals to a manifold.

Definition 4.10 (Euclidean -space). Euclidean -space, , is the set of all possible ordered -tuples of real numbers, that is, [43].

Definition 4.11 (hypersphere). A hypersphere or -sphere is given as , where denotes the radius of the hypersphere. A hypersphere with radius is called a unit hypersphere [44].

Definition 4.12 (manifold). Let be a set. If there exists a set of coordinate systems for which satisfies the following 2 conditions.

(1)Each element of is a one-to-one mapping from to some open subset of ; (2)for all , given any one-to-one mapping from to , is a diffeomorphism.

Then is called a manifold [45].

Theorem 4.13 (hyperspheres as mainfolds). The hypersphere is a manifold [46].

The proof of Theorem 4.13 is trivial, and is hereby skipped. For the complete proof, refer to Lee's treatment of topological manifolds [46].

Definition 4.14 (topological space). Let be a nonempty set. A class of subsets of is a topology on if and only if satisfies the following axioms:

(1) and belong to ; (2)The union of any number of sets in belongs to ; (3)The intersection of any two sets in belongs to .

The members of are called -open sets, or simply open sets, and together with , that is the pair is called a topological space [47].

Definition 4.15 (topological subspace). Let be a topological space, and let be any subset. The subset topology on is defined by

The pair is called a topological subspace of [46].

4.3. -Normalization

-normalization is the process of scaling a given vector in -dimensional space so that its norm meets a specific value. It is commonly used in data mining applications such as text categorization and gene expression analysis, where the data possesses a high dimension and is inherently directional in nature [42]. The norm of a vector is defined as follows.

Definition 4.16 (-norm). For a vector , the -norm is given by

A special case of the norm is when , where the norm is called the Hilbert norm [48].

There is no unique method of carrying out the normalization, as any suitable transformation can be used, provided the resulting norm meets its desired value. Common transformations include the use of trigonometric operators, and division by a common factor.

5. Measure Construction

During an epileptic seizure, EEG electrodes are attached to the patient's scalp and the resulting potential differences are recorded. To standardize the locations from which the potential differences were recorded, a system of electrode placement was introduced, called the International 10–20 System [49].

The seizure data used in this research is that utilized by Ahmad in his work, where 25 (21 active and 4 reference) electrodes are attached to a patient's head in accordance with the International 10–20 System. The standard electode positioning for this system is illustrated in Figure 1. When a seizure occurs, the data is recorded and trimmed to correspond to the period of ictal. To mathematically identify the EEG signals observed by their respective electrodes, the set of electrodes used in this research is introduced as the electrode set.

(a)

(b)

(c)

Figure 1: Standard electrode positioning as per the International 10–20 system [49].

Definition 5.1 (electrode set). The standard set of all monitored EEG electrodes during an epileptic seizure is called the electrode set, , and is given by

where represents an EEG electrode that is being monitored during an epileptic seizure.

In this study, the electrode set utilized is

Elements of the electrode set serve as indices to indicate the electrode at which a particular set of electrical potentials is detected.

Definition 5.2 (reference subset). The subset of the electrode set that contains only reference electrodes is called the reference subset, , and is given by:

In this study, the reference subset utilized is:

Definition 5.3 (electrode potential set). In any given patient, the set of all potential differences detected during a seizure at time is given by

where is the potential difference detected at electrode at time in relation to its closest reference electrode. We call the electrode potential set at time .

Prior to developing a custom measure, the data was tested for conformity to existing distributions. The data was analyzed twice, first using the electrodes as the independent variable, and again using one-second time frames as the independent variable. Taking electrodes as the independent variable allow a seizure to be viewed as a breakdown of potential differences detected across monitored electrodes at any fixed one-second time frame. Alternatively, taking time frames as the independent variable results in a view of the progression of a seizure in terms of the potential differences detected at each electrode over a given period of time. Visual inspection using Q-Q plots show that although some electrodes display characteristics that are close to that of normal and log-normal distributions, the same does not hold for all electrodes. As for the one-second time frame case, similar behavior is observable, as can be seen from Figures 2 and 3. Furthermore, Kolmogorov-Smirnov (K-S) testing at a tolerance level of reveals that the data does not entirely conform to that of normal, uniform, nor exponential distributions.

Figure 2: Q-Q of potential differences detected at time  s during a seizure. The data collected at this time frame conforms to that of a power law distribution.

Figure 3: Q-Q of potential differences detected at time  s during the same seizure. Note that the data no longer displays normal tendencies.

To compare the progression of one seizure to another, a general probability distribution that applies to all recorded events during the seizure is required. A generalized nonlinear model is most suited to the task, but limited understanding on the factors that influence a seizure results in a poorly designed model. Given that some electrodes display normal tendencies (including the power law distribution), the authors conjecture that the von Mises-Fisher (vMF) distribution is a natural candidate to represent an epileptic seizure. As mentioned in Section 2, Babloyantz and Destexhe describe a seizure as a low-dimensional state, which Stam calls a “loss of complexity.” This loss of complexity can be explained as a loss of structural information in the seizure data, resulting in a vMF distribution to appear somewhat normal in lower dimensions. Furthermore, the use of a hypersphere to describe a seizure gives the data a fixed structure to lie on, thus confining seizure data to a specific surface. This upscale is possible as a point in a lower dimension in topologically equivalent to its counterpart in a higher dimension.

As mentioned in Section 2, Bayesian modelling is highly dependent on apriori information. As a workaround to assigning probability values based on data derived from prior seizures, a simple measure is designed based on the deviation from the average signal strength being read by the EEG apparatus during a seizure. This allows seizures to be modelled using instantaneous data, thus eliminating the need for multiple measurements. The measure serves as a foothold for transferring the seizure data to a hypersphere so the vMF distribution can be applied. A simplified process flow diagram of the route taken by the authors is given in Figure 4.

Figure 4: Pathway to mapping EEG signals to a manifold.

To define a probability measure on , a -algebra is constructed on the set. Lemma 5.4 illustrates that the choice of -algebra to be defined is simple, namely, the power set of .

Lemma 5.4 (-algebra of ). The power set of defines a -algebra on the electrode potential set.

Proof. The power set of an electrode potential set is given by

Observe that

(1), (2) since ,(3) since ,(4) since ,(5)for any , since all elements of are subsets (proper or improper) of , and the union of any number of these subsets cannot exceed , (6)for any non-empty , since and are proper subsets of , and thus must belong to the power set of , (7)for any , by De Morgan's law. Note that and are both in by . Applying , we have that . Applying , we see that .

Thus, defines a -algebra on as per Definition 4.1.

Lemma 5.4 leads to the following theorem.

Theorem 5.5 (measurability of elements of ). is a measurable space, and the elements of are measurable.

Proof. By Lemma 5.4, is a -algebra on , that is, it is a -ring that contains . Thus is a measurable space by Definition 4.4, and the measurability of elements of is implied.

For the purpose of modelling EEG signals of epileptic seizures, inspiration for the measure is drawn from the standard score, which is commonly used in normalizing data sets. As with most measures, the intended measure will have parameters, that is certain variables that characterize the behavior of the measure. These parameters are the emulated baseline and jitter, and are defined as follows.

Definition 5.6 (emulated baseline). The emulated baseline, , is the average magnitude of electrical potentials detected across all monitored electrodes at a given time in a seizure and is given by

Definition 5.7 (jitter). The jitter, , is the total difference in magnitude of detected electrical potentials and the emulated baseline at a given time in a seizure and is given by

The emulated baseline and jitter represent the mean and standard deviation of the data set, respectively, with the change in terminology due to the nature of the data being recorded, that is, electrical signals. A sample of the computed emulated baseline and jitter values for a 10-second seizure is given in Table 1. The resulting measure, the Delia EEG measure, is formally introduced below.

Table 1: Result of K-S testing against a normal distribution for a 10-second seizure.

Mean
97.40
77.92
75.21
50.70
81.70

Std. Dev
112.47
72.14
61.77
50.16
58.18

K-S
1.35
1.03
1.03
1.08
.64

Asymp. Sig. (2-tailed)
.05
.24
.24
.19
.80

Mean
116.53
137.22
144.79
129.05
97.92

Std. Dev
104.45
125.66
106.72
93.47
101.03

K-S
.68
1.13
.89
.86
1.18

Asymp. Sig. (2-tailed)
.75
.16
.41
.46
.12

Definition 5.8 (Delia EEG measure). The Delia EEG measure, , is a measure defined on the power set of an electrode potential set, and is given by the set function:

where , , is the number of elements in the set , and is the set of extended reals.

The following lemma verifies that the set function described in Definition 5.8 satisfies the requirements to be a measure.

Lemma 5.9 (measure property of the Delia set function). When defined on the power set of an electrode potential set, the set function is a measure.

Proof. Note that elements of are mutually exclusive, that is, holds for any since no two potential differences can be recorded at one electrode at the same time.

(1) Pick any such that and consider the following cases:

(a)Case  I. Let

(b)Case  II. Let

 If , merely swaping the sets and . (c)Case  III. Let

 rearrange and reindex the above so that

 then

 Hence is additive for mutually exclusive sets.

(2)Picking any ,

 Thus is positive since . (3) since and for any set such that (by definition). (4)Consider a sequence of increasing sets in , that is, . Observe that can be rewritten as . Thus,

Similarly,

Hence, the lub of is . Denote the union of all sets in as . Then,

Therefore, , and, consequenitly, .

Thus is a measure by Definition 4.3.

The following lemma shows that the Delia EEG measure can be used as a probability measure on the electrode potential set.

Lemma 5.10 (measure of ). The total measure of is 1 with respect to .

Proof. The total measure of is given by

Thus, has a total measure of 1 with respect to .

With Lemma 5.10, the following theorem can be proven.

Theorem 5.11 (probability space of EEG signals). is a probability space.

Proof. By Theorem 5.5 and Lemma 5.10, is a measurable space with total measure 1 with respect to . Thus, by Definition 4.5, is a probability space.

The Delia measure operates as a probability mass function (PMF), assigning each element a probability value . The following section details the application of the measure to data sets of three patients experiencing epileptic seizures.

6. Measure Implementation

The EEG signal data of three patients A, B, and C, suffering from epileptic seizures were modelled using the Delia measure. Emulated baseline and jitter values are first computed from the raw EEG data. Table 2 gives an example of the emulated baseline and jitter values for a 10-second seizure experienced by Patient A.

Table 2: Emulated baseline and jitter values for a 10-second seizure (Patient A).

Time ()

1
98.00806
1586.979

2
78.68421
1036.199

3
76.41594
1040.109

4
51.93148
892.6314

5
78.81386
939.0676

6
120.1871
1672.974

7
143.4501
2120.952

8
150.1537
1754.462

9
132.2427
1533.646

10
99.98561
1576.88

A plot of versus shows the progression of neural activity at any specific electrode in comparison to the rest. Due to the nature of the measure, ratio of signals in one electrode to another is identical to that of the unprocessed signal, that is, the geometric structure of the signal remains unchanged. can be physically interpreted as the rate of seizure activity at any specific electrode relative to all other monitored electrodes. This gives a simplified view of the seizure, with emphasis on electrodes that show the most changes in activity.

As an example, Figure 5 shows that two electrodes display significant changes in measure value throughout the seizure, namely, electrodes and . Core seizure activity takes place in the range of 0.00–0.15, indicating that the majority of the signals detected deviate from the emulated baseline by at most units. Both of the electrodes and are situated near the prefrontal cortex, which accounts for the behavioral changes associated with some epileptic patients [50, 51].

Figure 5: Plot of measure values over time for a 10-second seizure observed in Patient A.

Comparatively, the measure data of Patient B (Figure 6) exhibits sharp fluctuations in the potential difference detected at electrode , which is situated at the premotor cortex of the right hemisphere of the brain, which is the section of the brain that is responsible for sensory guidance of movement and control of proximal and trunk muscles of the body [52]. Video footage taken during the seizure confirms that the patient experienced body spasms during ictal.

Figure 6: Plot of measure values over time for a 15-second seizure observed in Patient B.

On the other hand, patient C (Figure 7), whose seizure lasts 25 seconds, exhibits a very focused fluctuation in potential difference, namely at electrodes and . Both these electrodes are at exact opposite ends of the left hemisphere of the brain, one being in the anterior region (), and the other at the posterior region (). The graph also indicates that the readings at and are inversely proportional to each other, suggesting that the bulk of the electrical activity during Patient C's seizure was moving back and forth between these two electrodes.

Figure 7: Plot of measure values over time for a 25-second seizure observed in Patient C.

However, application of the Delia measure alone is insufficient to describe the progression of seizures, as it only reformats the data to fit in the interval. Based on the earlier assumption made in Section 5, more data regarding the pattern of EEG signals during a seizure may lie in higher dimensions. For this purpose, a mapping is developed to transfer the output of the Delia measure to a manifold.

7. Mapping Measure Data to a Manifold

This section covers the construction of a mapping from the set of Delia measure values to a high-dimensional manifold. To begin with, the set of selected output generated by applying the Delia measure is formalized as the scaled potential set.

Definition 7.1 (scaled potential set). The scaled potential set is the set of values generated by applying the Delia measure on singleton subsets of an electrode potential set, and is given by

Only the singleton subsets of the electrode potential sets are considered to form the scaled potential set. This is to ensure that each measure value in the set corresponds to only one potential difference, detected at a single electrode during the seizure. Since with and , each scaled potential set represents the discrete probability distribution of EEG signals observed at all electrodes during any given time frame during the seizure. Due to the nature of the data, and the intention of utilizing the vMF distribution, the manifold required is a unit hypersphere. Transferring the contents of the to such a manifold is accomplished using Lemma 7.2.

Lemma 7.2 (spherical compatibility of probability values). A discrete probability distribution with outcomes can be represented as a point on a unit hypersphere.

Proof. Consider a discrete probability distribution with outcomes, and let denote the probability of each event with respect to the measure . Then , and it follows that for each , with . By letting for each ,

which is the general form for a unit hypersphere.

Bahlmann noted in 2006 that any probability distribution can be “wrapped” around the circumference of a circle with unit radius [53]. However, this approach was not utilized as it sacrifices information regarding the electrode of origin for each measure value. The following theorem is a direct result of Lemma 7.2.

Theorem 7.3 (EEG Signals as points on a hypersphere). The EEG signals detected during a seizure can be viewed as points on the surface of a unit hypersphere.

Proof. Given an epileptic seizure of length seconds, discrete electrode potential sets are constructed, each one representing the EEG potential differences detected at all connected electrodes at intervals. Using Definitions 5.8 and 7.1, scaled potential sets are formed, namely, . Suppose there are electrodes in use. By Lemma 7.2, the positive square root of elements of each scaled potential set can be viewed as points on the surface of a -dimensional unit hypersphere. Thus, by forming scaled potential sets, the EEG signals detected at electrodes during an epileptic seizure can be ported to points on a unit hypersphere of dimension .

Definition 7.4 below formally defines the mapping from the scaled potential set to a unit hypersphere, whilst Definition 7.5 specifies the subset of the hypersphere that contains EEG recordings.

Definition 7.4 (EEG hypersphere map). The EEG hypersphere map, , maps elements in a scaled potential set to the surface of a -dimensional hypersphere and is given as

where

refers to the Delia measure value corresponding to the potential difference detected at electrode at time .

The EEG hypersphere map, used in conjunction with the Delia measure, allows for the EEG signals of epileptic seizures to be viewed as points on the surface of a unit hypersphere. Using a spherical distribution such as the von Mises-Fisher distribution, further modelling of seizures can be carried out. However, the vMF modelling of seizures is not included in this paper, as the process involves estimating the concentration parameter , which is notoriously lengthy.

Definition 7.5 (the observed signal set). The observed signal set, is the finite subset .

With Definition 7.5, the following theorem takes form.

Theorem 7.6 (topology of EEG signals). The observed signal set, , is a topological space.

Proof. Let be the standard -dimensional Euclidean space with standard Euclidean topology . Each is given in the form , where is the standard Euclidean metric, and . Define the topology for to be , where each is in the form . Then,

(1) since for ,(2) since for ,(3) since

(4) since

Thus is a topological subspace, and hence a topological space.

The above proof that the seizure signals on the manifold constitutes a topological space verifies Ahmad's claim that a topological space can be constructed to house EEG signals of epileptic seizures [33].

8. Conclusion and Discussion

In this paper, measure theory was utilized to design the Delia measure, a discrete probability measure that reformats EEG data without altering its geometric structure. An analysis of EEG data from three patients experiencing epileptic seizures was made using the measure, with the resulting graph confirming the presence of symptoms such as behavioral changes and body spasms during recording. A mapping called the EEG hypersphere map then is devised to transport the measure data onto the surface of a unit hypersphere, thus allowing seizures to be viewed as a geometric structure. The authors suggest that further statistical analysis of seizures be carried out using the von Mises-Fisher distribution, justification for which is given in Section 5. The subset of seizure signals on the manifold is shown to be a topological space, verifying Ahmad's approach to use topological modelling.

Acknowledgments

The authors would like to thank the Ministry of Science, Technology and Innovation for funding this research under the Fundamental Research Grant Scheme (Vot. no. 78315) and the National Science Fellowship (NSF) Grant.

References

WHO, Epilepsy factsheet, 2009

H. Gastaut, “Clinical and electroencephalographical classification of epileptic seizures,” Epilepsia, vol. 11, no. 1, pp. 102–112, 1970.

T. Ahmad, R. A. Fairuz, F. Zakaria, and H. Isa, “Selection of a subset of EEG channels of epileptic patient during seizure using PCA,” in Proceedings of the 7th WSEAS International Conference on Signal Processing, Robotics and Automation, pp. 270–273, World Scientific and Engineering Academy and Society (WSEAS), Stevens Point, Wis, USA, 2008.

F. Zakaria, Dynamic profiling Of EEG data during seizure using fuzzy information space [Ph.D. thesis], Universiti Teknologi Malaysia, 2008.

E. Niedermeyer and F. H. L. Silva, Electroencephalography: Basic Principles, Clinical Applications, and Related Fields, M—Medicine Series, Lippincott Williams & Wilkins, 2005.

A. Babloyantz and A. Destexhe, “Low-dimensional chaos in an instance of epilepsy,” Proceedings of the National Academy of Sciences, vol. 83, no. 10, pp. 3513–3517, 1986.

C. J. Stam, “Nonlinear dynamical analysis of EEG and MEG: review of an emerging field,” Clinical Neurophysiology, vol. 116, no. 10, pp. 2266–2301, 2005.

L. D. Iasemidis, S. J. Chris, H. P. Zaveri, and W. J. Williams, “Phase space topography and the lyapunov exponent of electrocorticograms in partial seizures,” Brain Topography, vol. 2, no. 3, pp. 187–201, 1990.

G. W. Frank, T. Lookman, M. A. H. Nerenberg, C. Essex, J. Lemieux, and W. Blume, “Chaotic time series analyses of epileptic seizures,” Physica D, vol. 46, no. 3, pp. 427–438, 1990.

J. Theiler, “On the evidence for low-dimensional chaos in an epileptic electroencephalogram,” Physics Letters A, vol. 196, no. 5-6, pp. 335–341, 1995.

N. Schiff, J. Victor, A. Canel, and D. Labar, “Characteristic nonlinearities of the 3/s ictal electroencephalogram identified by nonlinear autoregressive analysis,” Biological Cybernetics, vol. 72, no. 6, pp. 519–526, 1995.

R. Friedrich and C. Uhl, “Spatio-temporal analysis of human electroen-cephalograms: petit-mal epilepsy,” Physica D, vol. 98, no. 1, pp. 171–182, 1996.

J. L. Hernndez, P. A. Valds, and P. Vila, “EEG spike and wave modelled by a stochastic limit cycle,” NeuroReport, vol. 7, no. 13, pp. 2246–2250, 1996.

M. Le van Quyen, J. Martinerie, C. Adam, and F. J. Varela, “Unstable periodic orbits in human epileptic activity,” Physical Review E, vol. 56, no. 3, pp. 3401–3411, 1997.

M. Feucht, U. Mller, H. Witte, et al., “Nonlinear dynamics of 3 hz spike- and-wave discharges recorded during typical absence seizures in children,” Cerebral Cortex, vol. 8, no. 6, pp. 524–533, 1998.

T. E. Peters, N. C. Bhavaraju, M. G. Frei, and I. Osorio, “Network system for automated seizure detection and contingent delivery of therapy,” Journal of Clinical Neurophysiology, vol. 18, no. 6, pp. 545–549, 2001.

C. E. Elger and K. Lehnertz, “Seizure prediction by non-linear time series analysis of brain electrical activity,” European Journal of Neuroscience, vol. 10, no. 2, pp. 786–789, 1998.

J. Martinerie, C. Adam, M. Le van Quyen, et al., “Epileptic seizures can be anticipated by non-linear analysis,” Nature Medicine, vol. 4, pp. 1173–1176, 1998.

H. R. Moser, B. Weber, H. G. Wieser, and P. F. Meier, “Electroencephalograms in epilepsy: analysis and seizure prediction within the framework of lyapunov theory,” Physica D, vol. 130, no. 3-4, pp. 291–305, 1999.

Y. C. Lai, M. A. F. Harrison, M. G. Frei, and I. Osorio, “Inability of lyapunov exponents to predict epileptic seizures,” Physical Review Letters, vol. 91, no. 6, Article ID 068102, 4 pages, 2003.

I. Osorio, M. A. Harrison, Y. C. Lai, and M. G. Frei, “Observations on the application of the correlation dimension and correlation integral to the prediction of seizures,” Journal of Clinical Neurophysiology, vol. 18, no. 3, pp. 269–274, 2001.

M. A. F. Harrison, I. Osorio, M. G. Frei, S. Asuri, and Y. C. Lai, “Correlation dimension and integral do not predict epileptic seizures,” Chaos, vol. 15, no. 3, Article ID 33106, 2005.

S. Kalitzin, J. Parra, D. N. Velis, and F. H. Lopes da Silva, “Enhancement of phase clustering in the EEG/MEG gamma frequency band anticipates transitions to paroxysmal epileptiform activity in epileptic patients with known visual sensitivity,” IEEE Transactions on Biomedical Engineering, vol. 49, no. 11, pp. 1279–1286, 2002.

W. van Drongelen, S. Nayak, D. M. Frim, et al., “Seizure anticipation in pediatric epilepsy: use of kolmogorov entropy,” Pediatric Neurology, vol. 29, no. 3, pp. 207–213, 2003.

F. Mormann, R. G. Andrzejak, C. E. Elger, and K. Lehnertz, “Seizure prediction: the long and winding road,” Brain, vol. 130, no. 2, pp. 314–333, 2007.

P. McSharry, T. He, L. Smith, and L. Tarassenko, “Linear and non-linear methods for automatic seizure detection in scalp electro-encephalogram recordings,” Medical and Biological Engineering and Computing, vol. 40, no. 4, pp. 447–461, 2002.

D. Kugiumtzis and P. G. Larsson, “Prediction of epileptic seizures with linear and nonlinear analysis of EEG,” in Chaos in Brain?: Proceedings of the Workshop, K. Lehnertz and C. E. Elger, Eds., pp. 329–332, World Scientific, 1999.

K. K. Jerger, T. I. Netoff, J. T. Francis, et al., “Early seizure detection,” Journal of Clinical Neurophysiology, vol. 18, pp. 259–268, 2001.

S. Baillet and L. Garnero, “A bayesian approach to introducing anatomo-functional priors in the EEG/MEG inverse problem,” IEEE Transactions on Biomedical Engineering, vol. 44, no. 5, pp. 374–385, 1997.

J. R. Stevens, “Electroencephalographic studies of conditional cerebral response in epileptic subjects,” Electroencephalography and Clinical Neurophysiology, vol. 12, no. 2, pp. 431–444, 1960.

O. Faust, R. U. Acharya, A. R. Allen, and C. M. Lin, “Analysis of eeg signals during epileptic and alcoholic states using ar modeling techniques,” Ingénierie et Recherche Biomédicale, vol. 29, no. 1, pp. 44–52, 2008.

N. Sivasankari and K. Thanushkodi, “Automated epileptic seizure detection in EEG signals using fast ICA and neural networks,” International Journal of Soft Computing Applications, vol. 1, no. 2, 2009.

T. Ahmad, R. S. Ahmad, F. Zakaria, and L. L. Yun, “Development of detection model for neuromagnetic fields,” in Proceedings of the BIOMED 2000, pp. 119–121, University of Malaya, September 2000.

T. Ahmad, R. S. Ahmad, W. E. Z. W. Abdul Rahman, L. L. Yun, and F. Zakaria, “Fuzzy topographic topological mapping for localisation simulated multiple current sources of MEG,” Journal of Interdisciplinary Mathematics, vol. 11, pp. 381–393, 2008.

A. Idris, T. Ahmad, and N. Maan, “A novel technique for visualization electrical activities in the brain during epileptic seizure,” in Proceedings of the International Conference on Applied Mathematics and Informatics, pp. 94–99, 2010.

J. M. Franks, Terse Introduction to Lebesgue Integration, Student mathematical library, American Mathematical Society, Providence, RI, USA, 2009.

S. K. Berberian, Measure and Integration, AMS Chelsea Publishing, Providence, RI, USA, 2011.

G. G. Roussas, An Introduction to Measure-Theoretic Probability, Elsevier/Academic Press, Amsterdam, The Netherlands, 2005.

W. J. Ewens and G. R. Grant, Statistical Methods in Bioinformatics: An Introduction, vol. 10 of Statistics for Biology and Health, Springer, New York, NY, USA, 2005.

M. J. Hassett and D. Stewart, Probability for Risk Management, ACTEX academic series, ACTEX Publications, 2006.

R. R. Wilcox, Introduction to Robust Estimation and Hypothesis Testing, Statistical Modeling and Decision Science, Elsevier Science & Technology, 1997.

A. Banerjee, I. S. Dhillon, J. Ghosh, and S. Sra, “Clustering on the unit hypersphere using von Mises-Fisher distributions,” Journal of Machine Learning Research, vol. 6, pp. 1345–1382, 2005.

E. Zakon, B. Lucier, and T. Zakon, Mathematical Analysis, The Zakon Series on Mathematical Analysis, The Trillia Group, 2004.

A. Dold, Lectures on Algebraic Topology, Classics in Mathematics, Springer, Berlin, Germany, 1995.

S. I. Amari, H. Nagaoka, and D. Harada, Methods of Information Geometry, Translations of Mathematical Monographs, American Mathematical Society, Providence, RI, USA, 2007.

J. M. Lee, Introduction to Topological Manifolds, vol. 202 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 2011.

S. Lipschutz, Schaum's Outline of General Topology, McGraw-Hill, 1988, Séerie Schaum.

Y. Eidelman, V. D. Milman, and A. Tsolomitis, Functional Analysis: An Introduction, Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, USA, 2004.

J. Malmivuo and R. Plonsey, Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields, Oxford University Press, 1995.

J. W. B. Douglas, “A neuropsychiatric study in childhood by Michael Rutter, Philip Graham and William Yule, clinics in developmental medicine, volume 35/36. (pp. 272; illustrated; 375.) heinemann medical books: London. 1970,” Psychological Medicine, vol. 1, no. 5, pp. 437–439, 1971.

S. McDermott, S. Mani, and S. Krishnawami, “A population-based analysis of specific behavior problems associated with childhood seizures,” Journal of Epilepsy, vol. 8, no. 2, pp. 110–118, 1995.

L. M. Surhone, M. T. Tennoe, and S. F. Henssonow, Premotor Cortex, VDM Verlag Dr. Mueller AG & Co. Kg, 2010

C. Bahlmann, “Directional features in online handwriting recognition,” Pattern Recognition, vol. 39, pp. 115–125, 2006.