# Hyperspherical Manifold for EEG Signals of Epileptic Seizures

Authors: Tahir Ahmad; Vinod Ramachandran
Journal: Journal of Applied Mathematics
 (2012)

---

## 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.


---

## Body

## 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 Section5.Definition 4.1 (σ-algebra).
SuposeX is a set and 𝒜 is a collection of subsets of X. 𝒜 is called a σ-algebra of subsets of X provided it contains the set X and is closed under taking complements (with respect to X), countable unions, and countable intersections [36].Definition 4.2 (set function).
A set function is a functionv 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)
μ(⊘)=0;(4)
ifEn is an increasing sequence in ℛ whose union E is also in ℛ, then μ(E)=lub[μ(En)] [37].Definition 4.4 (measurable space).
A measurable space is a pair(X,𝒫), where 𝒫 is a σ-ring of subsets of X. Elements of 𝒫 are said to be measurable [37].Definition 4.5 (probability space).
A probability space is a measurable space with total measure1, that is if X is a probability space, and μ is the measure defined on X, then μ(X)=1 [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).
LetX be a discrete random variable. A probability mass function for  X  is a function p(x) which assigns a probability to each value of the random variable such that  (1)
p(x)≥0 for all x, and(2)
∑‍p(x)=1 [40].Definition 4.8 (Kolmogorov-Smirnov test).
The Kolmogorov-Smirnov test is designed to testH0:F(x)=G(x),for allx versus H0:F(x)≠G(x)for  at  least  onex, where F(x) and G(x) are distributions associated with two independent groups [41].Definition 4.9 (von Mises-Fisher distribution).
Ad-dimensional random unit vector is said to have a d-variate von Mises-Fisher distribution if its probability density function is given by
(4.1)f(x∣μ,κ)=cd(κ)eκμTx,
where ∥μ∥=1, κ≥0, and d≥2. The normalizing constant cd(κ) is given by
(4.2)cd(κ)=κ(d/2)-1(2π)d/2I(d/2)-1(κ),
where Ir represents the modified Bessel function of the first kind with order r. μ 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 (Euclideann-space).
Euclideann-space, En, is the set of all possible ordered n-tuples of real numbers, that is, En={(x1,x2,…,xn)∣xi∈ℝ} [43].Definition 4.11 (hypersphere).
A hypersphere orn-sphere is given as Sn-1={x∈ℝn∣∥x∥=r2}, where r denotes the radius of the hypersphere. A hypersphere with radius 1 is called a unit hypersphere [44].Definition 4.12 (manifold).
LetS be a set. If there exists a set of coordinate systems 𝒜 for S which satisfies the following 2 conditions.  (1)
Each elementϕ of 𝒜 is a one-to-one mapping from S to some open subset of ℝn;(2)
for allϕ∈𝒜, given any one-to-one mapping ϕ from S to ℝn, ϕ∈𝒜⇔ϕ∘ϕ-1 is a C∞ diffeomorphism.Then (S,𝒜) is called a manifold [45].Theorem 4.13 (hyperspheres as mainfolds).
The hypersphereSn-1 is a manifold [46].The proof of Theorem4.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).
LetX be a nonempty set. A class τ of subsets of X is a topology on X if and only if τ satisfies the following axioms:  (1)
X 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 X together with τ, that is the pair (X,τ) is called a topological space [47].Definition 4.15 (topological subspace).
LetX be a topological space, and let S⊆X be any subset. The subset topology τS on S is defined by
(4.3)τs={U⊆S∣U=S∩VforsomeopensubsetV⊆X}.
The pair (S,τS) is called a topological subspace of X [46].


### 4.3.Lp-Normalization

L
p-normalization is the process of scaling a given vector in p-dimensional space so that its Lp 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 Lp norm of a vector is defined as follows.Definition 4.16 (Lp-norm).
For a vectorx, the Lp-norm is given by
(4.4)‖x‖=(∑i=1∞‍|xi|p)1/p.
A special case of the Lp norm is when p=2, 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 resultingLp norm meets its desired value. Common transformations include the use of trigonometric operators, and division by a common factor.


## 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 Section5.Definition 4.1 (σ-algebra).
SuposeX is a set and 𝒜 is a collection of subsets of X. 𝒜 is called a σ-algebra of subsets of X provided it contains the set X and is closed under taking complements (with respect to X), countable unions, and countable intersections [36].Definition 4.2 (set function).
A set function is a functionv 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)
μ(⊘)=0;(4)
ifEn is an increasing sequence in ℛ whose union E is also in ℛ, then μ(E)=lub[μ(En)] [37].Definition 4.4 (measurable space).
A measurable space is a pair(X,𝒫), where 𝒫 is a σ-ring of subsets of X. Elements of 𝒫 are said to be measurable [37].Definition 4.5 (probability space).
A probability space is a measurable space with total measure1, that is if X is a probability space, and μ is the measure defined on X, then μ(X)=1 [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).
LetX be a discrete random variable. A probability mass function for  X  is a function p(x) which assigns a probability to each value of the random variable such that  (1)
p(x)≥0 for all x, and(2)
∑‍p(x)=1 [40].Definition 4.8 (Kolmogorov-Smirnov test).
The Kolmogorov-Smirnov test is designed to testH0:F(x)=G(x),for allx versus H0:F(x)≠G(x)for  at  least  onex, where F(x) and G(x) are distributions associated with two independent groups [41].Definition 4.9 (von Mises-Fisher distribution).
Ad-dimensional random unit vector is said to have a d-variate von Mises-Fisher distribution if its probability density function is given by
(4.1)f(x∣μ,κ)=cd(κ)eκμTx,
where ∥μ∥=1, κ≥0, and d≥2. The normalizing constant cd(κ) is given by
(4.2)cd(κ)=κ(d/2)-1(2π)d/2I(d/2)-1(κ),
where Ir represents the modified Bessel function of the first kind with order r. μ 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 (Euclideann-space).
Euclideann-space, En, is the set of all possible ordered n-tuples of real numbers, that is, En={(x1,x2,…,xn)∣xi∈ℝ} [43].Definition 4.11 (hypersphere).
A hypersphere orn-sphere is given as Sn-1={x∈ℝn∣∥x∥=r2}, where r denotes the radius of the hypersphere. A hypersphere with radius 1 is called a unit hypersphere [44].Definition 4.12 (manifold).
LetS be a set. If there exists a set of coordinate systems 𝒜 for S which satisfies the following 2 conditions.  (1)
Each elementϕ of 𝒜 is a one-to-one mapping from S to some open subset of ℝn;(2)
for allϕ∈𝒜, given any one-to-one mapping ϕ from S to ℝn, ϕ∈𝒜⇔ϕ∘ϕ-1 is a C∞ diffeomorphism.Then (S,𝒜) is called a manifold [45].Theorem 4.13 (hyperspheres as mainfolds).
The hypersphereSn-1 is a manifold [46].The proof of Theorem4.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).
LetX be a nonempty set. A class τ of subsets of X is a topology on X if and only if τ satisfies the following axioms:  (1)
X 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 X together with τ, that is the pair (X,τ) is called a topological space [47].Definition 4.15 (topological subspace).
LetX be a topological space, and let S⊆X be any subset. The subset topology τS on S is defined by
(4.3)τs={U⊆S∣U=S∩VforsomeopensubsetV⊆X}.
The pair (S,τS) is called a topological subspace of X [46].


## 4.3.Lp-Normalization

L
p-normalization is the process of scaling a given vector in p-dimensional space so that its Lp 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 Lp norm of a vector is defined as follows.Definition 4.16 (Lp-norm).
For a vectorx, the Lp-norm is given by
(4.4)‖x‖=(∑i=1∞‍|xi|p)1/p.
A special case of the Lp norm is when p=2, 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 resultingLp 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 Figure1. 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.Standard electrode positioning as per the International 10–20 system [49].
(a)
(b)
(c)Definition 5.1 (electrode set).
The standard set of all monitored EEG electrodes during an epileptic seizure is called the electrode set,E, and is given by
(5.1)E={e1,e2,…,en},
where ei represents an EEG electrode that is being monitored during an epileptic seizure.
In this study, the electrode set utilized is(5.2)E={FP1,FP2,F3,F4,C3,C4,P3,P4,O1,O2,F7,F8,FT3,T4,T5,T6,A1,A2,FZ,CZ,PZ,T7,T8,P7,P8}.
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,R, and is given by:
(5.3)R={e∈E∣eisareference  electrode}.
In this study, the reference subset utilized is:(5.4)R={T7,T8,P7,P8}.Definition 5.3 (electrode potential set).
In any given patient, the set of all potential differences detected during a seizure at timet is given by
(5.5)Vt={ve∣ve∈R,e∈E∖R},
where ve is the potential difference detected at electrode e at time t in relation to its closest reference electrode. We call Vt the electrode potential set at time t.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 usingQ-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 0.5 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 t=23 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 t=15 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 Section2, 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 Section2, 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 onVt, 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 Vt.Lemma 5.4 (σ-algebra of Vt).
The power set ofVt defines a σ-algebra on the electrode potential set.Proof.
The power set of an electrode potential setVt={vFP1,vFP2,…,vFPZ} is given by
(5.6)P(Vt)={⊘,{vFP1,…,vFPZ},Vt}.
Observe that (1)
⊘∈P(Vt),(2)
Vt∈P(Vt) since Vt={vFP1,vFP2,…,vFPZ}∈P(Vt),(3)
⊘c∈P(Vt) since ⊘c=Vt∖⊘=Vt∈P(Vt),(4)
Vtc∈P(Vt) since Vtc=Vt∖Vt=⊘∈P(Vt),(5)
for anyA,B∈P(Vt), A∪B∈P(Vt) since all elements of P(Vt) are subsets (proper or improper) of Vt, and the union of any number of these subsets cannot exceed Vt,(6)
for any non-emptyA,B∈P(Vt), Ac,Bc∈P(Vt) since Ac=Vt∖A and Bc=Vt∖B are proper subsets of Vt, and thus must belong to the power set of Vt,(7)
for anyA,B∈P(Vt), A∩B=(Ac∪Bc)c by De Morgan's law. Note that Ac and Bc are both in P(Vt) by (6). Applying (5), we have that Ac∪Bc∈P(Vt). Applying (6), we see that A∩B=(Ac∪Bc)c∈P(Vt). Thus, P(Vt) defines a σ-algebra on Vt as per Definition 4.1.Lemma5.4 leads to the following theorem.Theorem 5.5 (measurability of elements ofP(Vt)).
(
V
t
,
P
(
V
t
)
) is a measurable space, and the elements of P(Vt) are measurable.Proof.
By Lemma5.4, P(Vt) is a σ-algebra on Vt, that is, it is a σ-ring that contains Vt. Thus (Vt,P(Vt)) is a measurable space by Definition 4.4, and the measurability of elements of P(Vt) 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 theemulated baseline and jitter, and are defined as follows.Definition 5.6 (emulated baseline).
The emulated baseline,βt, is the average magnitude of electrical potentials detected across all monitored electrodes at a given time t in a seizure and is given by
(5.7)βt=∑i=1n|vi|nforvi∈Vt.Definition 5.7 (jitter).
The jitter,δt, is the total difference in magnitude of detected electrical potentials and the emulated baseline at a given time t in a seizure and is given by
(5.8)δt=∑i=1n||vi|-βt|forvi∈Vt.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 Table1. 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.
t
=
1
t
=
2
t
=
3
t
=
4
t
=
5
Mean
97.40
77.92
75.21
50.70
81.70
Std. Dev
112.47
72.14
61.77
50.16
58.18
K-SZ
1.35
1.03
1.03
1.08
.64
Asymp. Sig. (2-tailed)
.05
.24
.24
.19
.80
t
=
6
t
=
7
t
=
8
t
=
9
t
=
10
Mean
116.53
137.22
144.79
129.05
97.92
Std. Dev
104.45
125.66
106.72
93.47
101.03
K-SZ
.68
1.13
.89
.86
1.18
Asymp. Sig. (2-tailed)
.75
.16
.41
.46
.12Definition 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:
(5.9)μ:P(Vt)⟶Resuchthatμ(A)={0ifn(A)=0,∑i=1n||vi|-βδ|otherwise,
where vi∈A, A∈P(Vt), n(A) is the number of elements in the set A, and ℝe is the set of extended reals.The following lemma verifies that the set function described in Definition5.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 ofVt are mutually exclusive, that is, A∩B=⊘ holds for any A,B∈P(Vt) since no two potential differences can be recorded at one electrode at the same time.  (1)
Pick anyA,B∈P(Vt) such that A∩B=⊘  and consider the following cases:
(a)
Case  I. LetA=B=⊘(5.10)μ(A∪B)=μ(⊘∪⊘)sinceA=B=⊘=μ(⊘)=0bydefinitionofμ=0+0=μ(⊘)+μ(⊘)bydefinitionofμ=μ(A)+μ(B)sinceA=B=⊘.(b)
Case  II. LetA≠⊘,B=⊘(5.11)μ(A∪B)=μ(A∪⊘)sinceB=⊘=μ(A)=μ(A)+0=μ(A)+μ(⊘)bydefinitionofμ=μ(A)+μ(B)sinceB=⊘
If A=⊘,B≠⊘, merely swaping the sets A and B.(c)
Case  III. LetA≠⊘,B≠⊘(5.12)μ(A∪B)=∑i=1n||vi|-βδ|wherevi∈A∪B=||v1|-βδ|+||v2|-βδ|+⋯+||vj|-βδ|+||vj+1|-βδ|+⋯+||vn|-βδ|rearrange and reindex the above so that(5.13)vi∈{Afori=1,2,…,jBfori=j+1,j+2,…,nthen(5.14)μ(A∪B)=∑i=1j‍||vi|-βδ|+∑i=j+1n‍||vi|-βδ|=μ(A)+μ(B).Henceμ is additive for mutually exclusive sets.(2)
Picking anyA∈P(Vt),
(5.15)μ(A)={0ifn(A)=0∑i=1n||vi|-βδ|otherwise.Thusμ is positive since |(|vi|-β)/δ|≥0.(3)
μ(⊘)=0 since n(⊘)=0 and μ(A)=0 for any set A such that n(A)=0 (by definition).(4)
Consider a sequence{A}n of increasing sets in P(Vt), that is, A1⊂A2⊂⋯⊂An. Observe that A2 can be rewritten as A2=(A2-A1)∪A1. Thus,
(5.16)μ(A2)=μ[(A2-A1)∪A1]=μ(A2-A1)+μ(A1)≥μ(A1)since(A2-A1)∩A1=⊘.Similarly,
(5.17)μ(A2)≤μ(A3)μ(A3)≤μ(A4)⋮μ(An-1)≤μ(An).
Hence, the lub of μ({A}n) is μ(An). Denote the union of all sets in {A}n as A. Then,
(5.18)A=⋃i=1nAi=A1∪A2∪⋯∪An=AnsinceA1⊂A2⊂⋯⊂An.
Therefore,μ(A)=μ(An), and, consequenitly, lub[μ({A}n)]=μ(A).
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 ofVt).
The total measure ofVt is 1 with respect to μ.Proof.
The total measure ofVt={v1,v2,…,vn} is given by
(5.19)μ(Vt)=μ({v1,v2,…,vn})=∑i=1n||vi|-βδ|=∑i=1n||vi|-β||δ|=∑i=1n||vi|-β|δsinceδisnonnegative=1δ∑i=1n||vi|-β|Sincethevalueofδisconstant∀vi∈Vt(5.20)=1δ⋅δsinceδ=∑i=1n‍||vi|-β|byDefinition5.7=1
Thus, Vt has a total measure of 1 with respect to μ.With Lemma5.10, the following theorem can be proven.Theorem 5.11 (probability space of EEG signals).
(
V
t
,
μ
) is a probability space.Proof.
By Theorem5.5 and Lemma 5.10, Vt is a measurable space with total measure 1 with respect to μ. Thus, by Definition 4.5, Vt is a probability space.The Delia measureμ operates as a probability mass function (PMF), assigning each element vi∈P(Vt) a probability value μt(vi). 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. Table2 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 (t)
β
t
δ
t
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.88A plot ofμ versus t 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, Figure5 shows that two electrodes display significant changes in measure value throughout the seizure, namely, electrodes FP1 and FP2. 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 0.15 units. Both of the electrodes FP1 and FP2 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 (Figure6) exhibits sharp fluctuations in the potential difference detected at electrode  F8, 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 (Figure7), whose seizure lasts 25 seconds, exhibits a very focused fluctuation in potential difference, namely at electrodes FP1 and O1. Both these electrodes are at exact opposite ends of the left hemisphere of the brain, one being in the anterior region (FP1), and the other at the posterior region (O1). The graph also indicates that the readings at FP1 and O1 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[0,1] 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 thescaled 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
(7.1)Vtμ={μ(A)∣A∈P(Vt),n(A)=1},
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 0<μ<1 and ∑μ∈Vtμ‍μ=1, each scaled potential set represents the discrete probability distribution of EEG signals observed at all electrodes e∈E during any given time frame t 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 Vtμ to such a manifold is accomplished using Lemma 7.2.Lemma 7.2 (spherical compatibility of probability values).
A discrete probability distribution withn outcomes can be represented as a point on a unit hypersphere.Proof.
Consider a discrete probability distributionX with n outcomes, and let μ(xi) denote the probability of each event x∈X with respect to the measure μ. Then Xμ={μ(x1),μ(x2),…,μ(xn)}, and it follows that 0<μ(xi)<1 for each μ(xi), with ∑x∈X‍μ(xi)=1. By letting ρi=μ(xi) for each x∈X,
(7.2)12=1=∑x∈Xμ(xi)=∑x∈Xρi2,
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 lengthn seconds, n discrete electrode potential sets are constructed, each one representing the EEG potential differences detected at all connected electrodes at 1s intervals. Using Definitions 5.8 and 7.1, n scaled potential sets are formed, namely, V1μ,V2μ,…,Vnμ. Suppose there are k electrodes in use. By Lemma 7.2, the positive square root of elements of each scaled potential set Vtμ can be viewed as points on the surface of a k-dimensional unit hypersphere. Thus, by forming n scaled potential sets, the EEG signals detected at k electrodes during an epileptic seizure can be ported to points on a unit hypersphere of dimension k.Definition7.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 k-dimensional hypersphere and is given as
(7.3)Θ:⋃t=1n‍Vtμ⟶S1k-1,
where
(7.4)Θ(Vtμ)=(μe1,μe2,…,μek),μei∈Vtμ refers to the Delia measure value corresponding to the potential difference detected at electrode ei at time t.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,O is the finite subset Θ(Vtμ)⊂S1k-1.With Definition7.5, the following theorem takes form.Theorem 7.6 (topology of EEG signals).
The observed signal set,O, is a topological space.Proof.
Let(ℝk,Γ) be the standard k-dimensional Euclidean space with standard Euclidean topology Γ. Each τ∈Γ is given in the form τx={y∈ℝn∣d(x,y)<ϵ}, where d is the standard Euclidean metric, and ϵ>0. Define the topology for O to be ΓO, where each ψ∈ΓO is in the form ψ=τ∩O. Then,  (1)
⊘∈ΓO since ⊘=⊘∩O for ⊘∈Γ,(2)
O∈ΓO since O=ℝk∩O for ℝk∈Γ,(3)
⋃‍ψ∈ΓO since
(7.5)⋃‍ψ=(τ1∩O)∪(τ2∩O)∪⋯=(⋃i‍τi)∩O=τ∩Osinceτ=⋃i‍τi∈Γ,(4)
⋂i=1k‍ψi∈ΓO since
(7.6)⋂i=1k‍ψi=(τ1∩O)∩(τ2∩O)∩⋯=(⋂i=1k‍τi)∩O=τ∩Osinceτ=⋂i=1k‍τi∈Γ.Thus (O,ΓO) 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 Section5. The subset of seizure signals on the manifold is shown to be a topological space, verifying Ahmad's approach to use topological modelling.



---

*Source: 926358-2012-08-01.xml*