Open Access

EEGgui: a program used to detect electroencephalogram anomalies after traumatic brain injury

  • Justin Sick1,
  • Eric Bray1,
  • Amade Bregy1, 2,
  • W Dalton Dietrich1, 2,
  • Helen M Bramlett1, 2, 4 and
  • Thomas Sick3Email author
Source Code for Biology and Medicine20138:12

https://doi.org/10.1186/1751-0473-8-12

Received: 29 October 2012

Accepted: 23 April 2013

Published: 21 May 2013

Abstract

Background

Identifying and quantifying pathological changes in brain electrical activity is important for investigations of brain injury and neurological disease. An example is the development of epilepsy, a secondary consequence of traumatic brain injury. While certain epileptiform events can be identified visually from electroencephalographic (EEG) or electrocorticographic (ECoG) records, quantification of these pathological events has proved to be more difficult. In this study we developed MATLAB-based software that would assist detection of pathological brain electrical activity following traumatic brain injury (TBI) and present our MATLAB code used for the analysis of the ECoG.

Methods

Software was developed using MATLAB(™) and features of the open access EEGLAB. EEGgui is a graphical user interface in the MATLAB programming platform that allows scientists who are not proficient in computer programming to perform a number of elaborate analyses on ECoG signals. The different analyses include Power Spectral Density (PSD), Short Time Fourier analysis and Spectral Entropy (SE). ECoG records used for demonstration of this software were derived from rats that had undergone traumatic brain injury one year earlier.

Results

The software provided in this report provides a graphical user interface for displaying ECoG activity and calculating normalized power density using fast fourier transform of the major brain wave frequencies (Delta, Theta, Alpha, Beta1, Beta2 and Gamma). The software further detects events in which power density for these frequency bands exceeds normal ECoG by more than 4 standard deviations. We found that epileptic events could be identified and distinguished from a variety of ECoG phenomena associated with normal changes in behavior. We further found that analysis of spectral entropy was less effective in distinguishing epileptic from normal changes in ECoG activity.

Conclusion

The software presented here was a successful modification of EEGLAB in the Matlab environment that allows detection of epileptiform ECoG signals in animals after TBI. The code allows import of large EEG or ECoG data records as standard text files and uses fast fourier transform as a basis for detection of abnormal events. The software can also be used to monitor injury-induced changes in spectral entropy if required. We hope that the software will be useful for other investigators in the field of traumatic brain injury and will stimulate future advances of quantitative analysis of brain electrical activity after neurological injury or disease.

Background

The development of epilepsy is an important secondary consequence of traumatic brain injury [13]. Standard methods for detecting and quantifying epilepsy in an animal model include electrographic detection of epileptiform events in the brain and behavioral assessments using observer scoring systems [4]. Currently, there are a number of techniques used to quantify different phenomena associated with changes in ECoG signals. Based on a review of the current literature, the most popular methods include Power Spectral Density (PSD), Short Time Fourier analysis and Spectral Entropy (SE). For our study, we built a graphical user interface (GUI) in the MATLAB programming platform that allows the user to analyze ECoG signals using the three methods mentioned above. In this paper we present our methods and MATLAB code used for the analysis of ECoG signals.

The basic PSD of a signal x(ti) describes how the power is distributed with respect to frequency. It is calculated by multiplying the Fast Fourier Transform (FFT) by the complex conjugate of the Fourier component.
P f i = X f i × X * f i
Where
X f i = t i x t i e 2 πi f i t i

and X * (f i ) is the complex conjugate of the Fourier component X (f i ) [5].

The PSD is very effective for describing which frequency components of a signal are most prominent. Unfortunately, there is no time information when signals are processed this way. For a signal that is changing over time, this method lacks the ability to identify when specific frequency components are present in the signal.

The Short Time Fourier Transform (STFT) advances through the signal in epochs. Each epoch is transformed individually which allows for frequency analysis in time domain. This is essential for linking specific events seen in the ECoG signal to changes in specific frequency components.

Different frequency components are associated with unique brain functions or pathologies. Delta rhythms are normally defined to be between frequencies 0.5-3.5Hz and are characteristic of deep sleep stages. Frequencies between 3.5-7.5Hz are called theta rhythms and are enhanced by the limbic system. Alpha rhythms (7.5-12.5Hz) are most prominent during drowsiness or coma, beta rhythms (12.5-30Hz) are enhanced with alertness and gamma rhythms (30-60Hz) are linked to common perceptual information [6]. We refer to neurophysiology literature for more information on the different frequency bands since this is not the focus of this manuscript.

Spectral Entropy describes the order of a signal in terms of power distributed across a range of frequencies. In this context, entropy describes the irregularity, complexity, or unpredictability characteristics of a signal [5]. The process for calculating SE from a signal x(ti) is given by [5] and is as follows.

First, the power spectrum is normalized so that the sum of the normalized power spectrum over the selected frequency range is equal to one.
f i = f 2 f 2 Pn f i = C n f i f 1 f 2 P f i = 1
SE is then calculated as the sum,
S f 1 f 2 = f i = f 1 f 2 P n f i log 1 P n f i
The SE is then normalized to a range between 1 (maximum irregularity) and 0 (complete regularity). This is done by dividing by the factor log of the total number of frequency components in the selected frequency range, N[f1,f2].
S N = S f 1 , f 2 log N f 1 , f 2

Methods

Animals and ECoG recordings

All procedures involving the use of live animals for this study were reviewed and approved by the Institutional Animals Care and Use Committee at the University of Miami School of Medicine. In this study we monitored ECoG activity and behavior (by video recording) simultaneously in awake freely moving rats 1 year after moderate (2.0 ATM) parasagittal fluid percussion brain injury (mTBI) or sham injury [2]. Twenty four hours prior to recording, animals were anesthetized with isofluorane (1.5%) and silver-silver choride electrodes were implanted over the dura mater through 1 mm burrholes in the skull. The electrodes were fixed in place with dental acrylic. Twenty four hours after recovery from surgery the electrodes were attached to a custom designed first stage amplifier (gain X 100, 0.1 Hz- 1 KHz bandpass filtered) and the headstage amplifier was tethered to a second stage amplifier by means of a flexible cable. Signals were amplified again and further filtered (final gain X 500, bandpass 0.1Hz – 500 Hz). Analog signals were digitized and stored using AD Instruments Powerlab™ and LabChart software.

EEG seizure analysis

The overall goal of this analysis is to detect abnormal events in the ECoG recordings that are periodic and not seen or seen less frequently in recordings from naïve animals. This was accomplished by: 1) normalizing power density to reduce variability that might occur between recordings due to differences in signal quality (ECoG amplitude) and 2) by detecting events over time defined as power density values that exceeded those found in a predefined segment of “normal” ECoG activity. Our EEG analysis software (EEGgui) was developed by modifying EEGLAB open access software running on the MATLAB platform. EEGgui offers researchers who are not proficient in computer programming to run analysis on ECoG signals using the MATLAB programming platform. The flow of data import, processing and figure production is shown in Figure 1.
Figure 1

Visual outline of the steps taken by EEGgui to analyze an ECoG record.

Import data

All files are imported as text files which allows for signals to be analyzed independent of the recording software. Records should have two dimensions, the voltage data and a time component. Our import function advances through the record looking for discontinuities in the time array based on the sampling frequency (Fs), then fills in missing data with NaN (not a number) values to keep the time component synched with the voltage data in signals that have segments removed due to artifacts. This is especially important for studies tracking changes in ECoG over time after a stimulus. The function (timelock.m) is provided in Additional file 1 along with detailed comments.

One constraint is that the record must be continuous. If the time component of the imported file does not advance in one direction throughout the record, an error will occur. The most common cause of this error when recording with LabChart is if the record is stopped and restarted, resulting in the time returning to 0 at some point during the record. Exampales of imported ECoG data are shown in Figure 2.
Figure 2

Screen shot of the raw data figure returned from EEgui. A. Initial Figure launched with raw data. The user can use the built in zoom function to find a region of awake ECoG. B. Zoomed in to region of interest (red area from A). The x-axis is in total data points, not time, for this figure. We suggest selecting 1 minute of ECoG. From our notes used during observation, we know the animal was awake and freely moving until the 5 minute mark which is equivalent to the 600,000th point, recording at 2000 samples per second. Therefore, for 1 minute of behavior, we select a start point of 300,000 and endpoint of 420,000 (red area in B).

Power normalization

One issue to overcome when analyzing ECoG in frequency domain is the dependence of power on signal strength. It is known that signal strength can vary depending on the distance of the electrodes from the source potential (brain tissue) and the conductivity of the tissue. While the amplitude of the recorded signal might not appear to vary significantly among animals, small changes in amplitude result in larger changes in power as power is a function of the square of the signal amplitude. Therefore, to minimize this potential problem we first select a region of the signal that is associated with awake, exploratory behavior of rats. We then average the power for a frequency range outside the range of interest (e.g. 100Hz-200Hz) and divide all powers by this value. For users with recording systems having more limited bandpass, the software must be modified to normalize power using a frequency band within recording limits (a description of this modification is provided with the source code.

The normalization procedure operates under two debatable assumptions. The first is that all differences in signal amplitude at the normalization frequency will be reflected by proportional changes in amplitude at the frequencies of interest. For example, in the case where signal amplitude (in microvolts) depends on the distance from the recording electrodes to the signal source (brain) we assume that the difference in amplitude of the signal at the normalization frequency will be proportional to the difference in amplitude at the frequencies of interest. In other words, assume a signal contains two frequencies f 1 and f 2 . The measured amplitudes will be, f 1 (x) and f 2 (x) for a distance x from the signal source. If the signal source is then moved to a distance x + Δx, the frequencies will have amplitudes of f 1 (x + Δx) and f 2 (x + Δx). Under our assumption,
f 1 x f 2 x = f 1 x + Δx f 2 x + Δx

Our second assumption is that ECoG at 100-200Hz measured during periods of awake, exploratory behavior is minimally affected by TBI. This will be discussed further below.

Power spectrum density

Our PSDs are calculated using Welch’s method. A detailed explanation of Welch’s method can be found in the MATLAB help. The pwelch function is available in the MATLAB signal processing toolbox, if the signal processing toolbox is not available, Lanspeary © 2006 developed a function that is available via the Free Software Foundation and is provided in our code.

We allow the user to select a region of interest on a plot of the signal to evaluate. This is a simple and efficient way to deduce the frequency components of particular segments of a signal corresponding to different behavior. An example of power spectral analysis from a data sample is shown in Figure 3.
Figure 3

Selected region for PSD analysis from Figure 2a. User can use the built in zoom function in MATLAB to narrow in on a region of interest in the raw signal record. Click “PSD,” then select region and this figure will appear which includes the raw data record (top trace) and the power density versus frequency graph (bottom trace).

Event detection

The STFT is used in part as our “event detection.” As mentioned earlier, the STFT is a Fourier transform performed on short time epochs of a signal. Editable variables in this part of the analysis are the length of each epoch and standard deviations. We define the frequency bands delta, theta, alpha, beta I, beta II and gamma and identify events where the experimental power deviates from control power within a frequency band.

Before events can be detected, the average power and standard deviations for each frequency band must be calculated from the control group to be compared against. Similar to our normalization, we use a 1 minute segment of awake, exploratory behavior. In the code we provide our calculated values. Included in the Additional file 1 is a detailed explanation on how to store new values if the user wishes to include their own recordings.

The user then inputs the number of standard deviations the experimental group power must deviate from the control power necessary to trigger an event. The program uses this information to proceed through the signal in steps, calculating the mean power for each frequency band and identifying events with dots. A figure will appear with the raw signal and power plots with events marked by dots. An example is shown in Figure 4.
Figure 4

First 10 minutes of recording after injection of pentylene tetrazol (PTZ). The two high amplitude bursts near 4 and 7 min are seizures. These are detected as events with power increases in all frequency band. In contrast, following the second seizure ECoG events associated with sleep are detected in lower frequency bands.

After each signal is analyzed, a variable is saved to the MATLAB workspace that holds the event data for each frequency band. The number of events for each recording can then be compared between all groups. For statistical analysis outside of MATLAB, event data can be copied to Microsoft EXCEL or a familiar software package.

Spectral entropy

During periods of awake, exploratory behavior, the entropy will be greatest because of the activity which is present across all frequency bands. During synchronous epileptic discharges, the majority of power is expressed at fewer frequencies, therefore the entropy will decrease. The spectral entropy is normalized to be independent of power, therefore we do not have to normalize each individual animal, instead we ratio it to our segment of normal behavior from naïve animals. No values much greater than one should be detected because of the assumption that normal behavior has the greatest entropy. Values of less than one represent segments of the signal where the SE has decreased. An example is shown in Figure 5 using a longer segment of the same data from which Figure 4 was derived.
Figure 5

Spectral Entropy (lower panel) plotted along with the filtered ECoG signal (top trace) and event detection as described above. The trace shows two seizure events followed by sleeping. The animal’s behavior was confirmed by visual observation and video analysis.

Before SE is calculated, the signal is low pass filtered at 60Hz. The filter can be found in the Additional file 1. This allows removal of electromyographic (EMG) contamination of the SE calculation during seizures. Powers in EMG frequencies diminish the changes in SE by spreading power across a wide range of frequencies that we are not studying. The NaN values used in the loading process are replaced with zeros for the filtering calculation then converted back to NaN before calculating SE.

Results

We were able to distinguish between a variety of behaviors represented in the ECoG signal such as awake and exploratory behavior, muscle artifact, sleep and different seizure categories. Examples of ECoG traces and detected events in the different power bands are shown in Figure 6.
Figure 6

Examples of ECoG recordings from an animal during awake exploratory behavior, sleeping, seizure activity, and chewing (electromyographic (EMG) artifact). Significant changes in spectral power for 2sec intervals are indicated by dots for the respective frequencies. Note that no significant events are detected for ECoG collected during awake behavior (as expected since this is the comparator signal) or during chewing artifact since EMG frequencies are outside of our selected power bands. Significant events are detected at slow frequencies during sleep and across all frequencies during seizures.

The ECoG during awake periods is of small amplitude. No peaks in the PSD are seen, no events in any power band are recorded and the SE is at its greatest during awake, exploratory behavior.

Muscle artifact had no affect on our analysis. Large bursts in the signal could be seen when the animal was chewing and grinding teeth. The power in EMG is expressed at frequencies higher than the power in ECoG and therefore not seen on PSD comparisons. This was confirmed using our event detection where no events were recorded using any of the frequency bands analyzed. Changes in SE are also not seen because the signal is filtered above 60Hz.

Sleeping was identified as a single peak in the PSD below 10Hz. The delta, theta, and alpha bands most often triggered events with the beta I and beta II bands occasionally triggering events during periods of sleep. The gamma band power does not trigger an event during sleep. SE is decreased during periods of sleep.

Seizures showed a massive increase in power across the frequency spectrum. During periods of epileptic activity, events were recorded in each frequency band. Spectral Entropy decreased dramatically.

Discussion

In this paper, we have described a program called EEGgui that was developed in the MATLAB programming platform. The purpose of the program is to allow users to quantify changes in brain activity compared to a pre-defined behavioral and electrographic state (in our case the awake state). ECoG records were used to investigate the onset of epilepsy 1 year after TBI. We put together an interface that includes the most common calculations currently used and introduced new concepts such as our event detection and the normalization process.

Among the many conditions associated with TBI, disturbances in sleeping patterns and frequent periods of fatigue are extremely common. A recent review of TBI and sleep disturbance suggests the need for objective, measurable deficits in sleep [7, 8]. To study sleeping patterns, long term observation under standard living conditions is necessary. Subjects would not need to be kept under strict observation for the entire recording period, a few lines of code could total the number of delta events occurring by themselves and compare these findings to control groups. Fatigue could be studied in a similar manner by examining the number of events in the alpha band.

The use of power in the 100–200 Hz band for the normalization process has yet to be proven. There are at least two known events that can occur in this frequency range in awake, freely moving animals. The first is muscle contraction artifact and the second is high frequency oscillations known as “ripples” [9]. There is literature that suggests the development of ripples in the 100–200 Hz range are associated with epilepsy and thus might occur after TBI [1012]. Both are periodic and recognizable through spectra analysis. As a result, they can be easily avoided by selecting short records of ECoG devoid of these components for normalization. Normalization might not be necessary if the raw ECoG signals among subjects are identical. However, power is a function of the square of ECoG amplitude and thus small differences in raw signal amplitude among subjects will result in large differences in power.

As mentioned above, the FFT provides no time information from the original signal. This is improved by advancing through the signal in short epochs, taking the FFT of each individual time segment. This method is still not perfect as it is limited by the uncertainty principle. If the epochs are too short, slow frequencies cannot be resolved. If epochs are too long, the specific time in which frequencies are present is lost.

Another transform to be considered which may be an improvement to the STFT is the Wavelet Transform. The Wavelet Transform has a varying epoch, a long epoch to accommodate slow frequencies and a short epoch for fast frequencies. Because the wavelet transform does not require a fixed epoch, it has advantages over STFT for signals with features that vary significantly in duration and frequency, such as ECoG [13]. This transform could prove to be a substantial improvement for our event detection and entropy analysis. MATLAB also provides a Wavelet toolbox that could be implemented into ECoG analysis. We are currently working on introducing a wavelet analysis package into our GUI.

As our study progresses, we expect to use these calculations to automate behavior recognition through ECoG analysis. An automated process will allow us to efficiently quantify changes in sleep patterns and epileptic activity, as well as identify other changes which might not be currently known to be associated with TBI.

Abbreviations

EEG: 

Electroencephalography

ECoG: 

Electrocorticographic

TBI: 

Traumatic brain injury

PSD: 

Power Spectral Density

SE: 

Spectral Entropy

GUI: 

Graphical user interface

FFT: 

Fast Fourier Transform

STFT: 

Short Time Fourier Transform

2.0 ATM: 

Moderate

mTBI: 

Parasagittal fluid percussion brain injury

Fs: 

Sampling frequency

NaN: 

Not a number

PTZ: 

Pentylene tetrazol

EMG: 

Electromyographic.

Declarations

Acknowledgements

This works was supported by Veterans Affairs grant 1 I01 BX000521-01 and NIH NS030291.

We thank Jessie Truettner for analyzing the videotaped behavior of the animals and Alexandra Wick for observing the recordings and producing recording hardware.

Authors’ Affiliations

(1)
The Miami Project to Cure Paralysis, University of Miami, Miller School of Medicine
(2)
Department of Neurosurgery, Lois Pope LIFE Center
(3)
Department of Neurology, Lois Pope LIFE Center
(4)
Bruce W. Carter Department of Veterans Affairs

References

  1. Pitkänen A, Bolkvadze T:Head Trauma and Epilepsy. Jasper’s Basic Mechanisms of the Epilepsies. Edited by: Noebels JL, Massimo A, Rogawski MA, Olsen RW, Delgado-Escueta AV. 2012, USA: Oxford University Press, 1-14. 4.Google Scholar
  2. Atkins CM, Truettner JS, Lotocki G, Sanchez-Molano J, Kang Y, Alonso OF, Sick TJ, Dietrich WD, Bramlett HM: Post-traumatic seizure susceptibility is attenuated by hypothermia therapy. Eur J Neurosci. 2010, 32: 1912-1920. 10.1111/j.1460-9568.2010.07467.x.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Vespa PM, McArthur DL, Xu Y, Eliseo M, Etchepare M, Dinov I, Alger J, Glenn TP, Hovda D: Nonconvulsive seizures after traumatic brain injury are associated with hippocampal atrophy. Neurology. 2010, 75: 792-798. 10.1212/WNL.0b013e3181f07334.PubMed CentralView ArticlePubMedGoogle Scholar
  4. Racine RJ: Modification of seizure activity by electrical stimulation: II. Motor seizure. Electroencephalogr Clin Neurophysiol. 1972, 32: 281-294. 10.1016/0013-4694(72)90177-0.View ArticlePubMedGoogle Scholar
  5. Viertio-Oja H, Maja V, Sarkela M, Tenkanen N, Tolvanen-Laakso H, Paloheimo M, Vakkuri A, Yli-Hankala A, Merilainen P: Description of the Entropy algorithm as applied in the Datex-Ohmeda S/5 Entropy Module. Acta Anaesthesiol Scand. 2004, 48: 154-161. 10.1111/j.0001-5172.2004.00322.x.View ArticlePubMedGoogle Scholar
  6. Singer W: Synchronization of cortical activity and its putative role in information processing and learning. Annu Rev Physiol. 1993, 55: 349-374. 10.1146/annurev.ph.55.030193.002025.View ArticlePubMedGoogle Scholar
  7. Orff HJ, Ayalon L, Drummond SPA: Traumatic brain injury and sleep disturbance: a review of current research. J Head Trauma Rehabil. 2009, 24: 155-165. 10.1097/HTR.0b013e3181a0b281.View ArticlePubMedGoogle Scholar
  8. Beaulieu-Bonneau S, Morin CM: Sleepiness and fatigue following traumatic brain injury. Sleep Med. 2012, 13: 598-605. 10.1016/j.sleep.2012.02.010.View ArticlePubMedGoogle Scholar
  9. Jefferys JGR, Menendez de la Prida L, Wendling F, Bragin A, Avoli M, Timofeev I, Lopes da Silva FH: Mechanisms of physiological and epileptic HFO generation. Prog Neurobiol. 2012, 98: 250-264. 10.1016/j.pneurobio.2012.02.005.View ArticlePubMedGoogle Scholar
  10. Jiruska P, Bragin A: High-frequency activity in experimental and clinical epileptic foci. Epilepsy Res. 2011, 97: 300-307. 10.1016/j.eplepsyres.2011.09.023.PubMed CentralView ArticlePubMedGoogle Scholar
  11. Bragin A, Engel J, Staba RJ: High-frequency oscillations in epileptic brain. Curr Opin Neurol. 2010, 23: 151-156. 10.1097/WCO.0b013e3283373ac8.PubMed CentralView ArticlePubMedGoogle Scholar
  12. Engel J, Bragin A, Staba R, Mody I: High-frequency oscillations: what is normal and what is not?. Epilepsia. 2009, 50: 598-604. 10.1111/j.1528-1167.2008.01917.x.View ArticlePubMedGoogle Scholar
  13. Schiff SJ, Aldroubi A, Unswer M, Sato S: Fast wavelet transformation of EEG. Electroencephalogr Clin Neurophysiol. 1994, 91: 442-455. 10.1016/0013-4694(94)90165-1.View ArticlePubMedGoogle Scholar

Copyright

© Sick et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement