Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Software Tool Article

AR2, a novel automatic artifact reduction software method for ictal EEG interpretation: Validation and comparison of performance with commercially available software

[version 1; peer review: 2 approved with reservations]
PUBLISHED 10 Jan 2017
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

This article is included in the INCF gateway.

Abstract

Objective: To develop a novel software method (AR2) for reducing muscle contamination of ictal scalp electroencephalogram (EEG), and validate this method on the basis of its performance in comparison to a commercially available software method (AR1) to accurately depict seizure-onset location.
Methods: A blinded investigation used 23 EEG recordings of seizures from 8 patients. Each recording was uninterpretable with digital filtering because of muscle artifact and processed using AR1 and AR2 and reviewed by 26 EEG specialists. EEG readers assessed seizure-onset time, lateralization, and region, and specified confidence for each determination. The two methods were validated on the basis of the number of readers able to render assignments, confidence, the intra-class correlation (ICC), and agreement with other clinical findings.
Results: Among the 23 seizures, two-thirds of the readers were able to delineate seizure-onset time in 10 of 23 using AR1, and 15 of 23 using AR2 (p<0.01). Fewer readers could lateralize seizure-onset (p<0.05). The confidence measures of the assignments were low (probable-unlikely), but increased using AR2 (p<0.05). The ICC for identifying the time of seizure-onset was 0.15 (95% confidence interval (CI), 0.11-0.18) using AR1 and 0.26 (95% CI 0.21-0.30) using AR2.  The EEG interpretations were often consistent with behavioral, neurophysiological, and neuro-radiological findings, with left sided assignments correct in 95.9% (CI 85.7-98.9%, n=4) of cases using AR2.
Conclusions: EEG artifact reduction methods for localizing seizure-onset does not result in high rates of interpretability, reader confidence, and inter-reader agreement. However, the assignments by groups of readers are often congruent with other clinical data. Utilization of the AR2 software method may improve the validity of ictal EEG artifact reduction.

Keywords

scalp EEG, electroencephalogram, muscle artifact, independent component analysis, seizure

Introduction

The scalp electroencephalogram (EEG) is a critical diagnostic tool in the evaluation of seizures, but artifact from muscle contraction often limits its use because of the obscuring of the cerebrally generated potentials. This problem is present in 11% of ictal EEGs overall and up to 70% of frontal lobe seizures13. The inability to discern the seizure-onset zone from scalp EEG often necessitates additional testing, including (positron emission tomography) PET, magnetoencephalography, ictal Single-photon emission computed tomography (SPECT), and intracranial EEG4. Each of these tests adds undesired time and cost to the evaluation.

Digital filters are the common approach to maximizing the likelihood of identifying a seizure-onset zone from EEG with muscle artifact. This filtering reduces muscle artifact by attenuating all frequencies beyond a selected value5, but it may impair the integrity of the EEG recording since brain-generated potentials may be in the same frequency band6,7. Recently, new technologies to reduce muscle artifact based on independent component analysis (ICA)810 have become available. ICA removes artifacts based on source-related features instead of frequencies1114. Prior studies have demonstrated that ICA-based methods improve the interpretation of artifact-laden ictal EEG recordings; in these studies researchers manually performed the ICA analysis prior to performing the EEG interpretation15,16. Automatic artifact reduction using ICA8 has become commercially available and is included in the latest versions of popular EEG viewer software17.

Despite the utilization of these software products by neurologists around the globe, the clinical benefit has not been established. It is also unknown if the new approaches introduce confounding artifacts that may lead to erroneous interpretations.

The goal of this study was to assess the validity of a commercially available EEG artifact reduction tool (AR1)17 and compare its validity to a novel automatic artifact reduction tool (AR2), which was developed at the University of California Los Angeles on the basis of inter-reader agreement, confidence, and congruence with other clinical findings, and which we are describing here.

Methods

Implementation

The custom software algorithm involved importing EEG scalp recordings as European Data Format (EDF) files in Matlab 8.4 (Mathworks, Natick, MA). The imported EEG was band pass filtered (16–70 Hz) using a 500th order finite impulse response filter, i.e. FIR1 in referential montage. We then applied a power spectral density algorithm to find extended intervals of elevated high frequency power across electrodes. We next calculated the normalized mutual information (MI)18 adjacency matrix across all scalp electrode contacts during the (16–70 Hz) band-pass filtered artifact epoch of greatest duration and assigned each scalp EEG electrode a single MI value derived from the maximum pairwise MI values in the adjacency matrix. We then determined if this maximum mutual information value exceeded a threshold value, and if that electrode should be included in subsequent artifact reduction processing. If the recording lacked an artifact epoch, or all channels were excluded, artifact reduction was applied to the referential recordings from all recording electrodes.

The high pass filtered (>16 Hz) scalp EEG was then separated into consecutive 120-second trials and each trial was processed using CUDAICA19. The purpose of the ICA was to separate the (>16 Hz) seizure activity, from the (>16 Hz) muscle artifact. The 16 Hz cut-off for the filter was chosen to isolate the vast majority of the muscle artifact. Independent components that explained an amount of variance above a particular threshold were excluded from the signal. The threshold was selected on the basis of the values of the raw and normalized mixing matrix (i.e. inverse weight matrix) calculated in each of the ICA iterations. We assumed that the last myogenic component and first neurogenic component can be differentiated on the basis of the inverse weight matrix, which provides the spatial distribution of each component, and identifying the independent component of greatest order with a focal spatial topography defined on the basis of exceeding a normalized threshold in at least one electrode of the inverse weight matrix.

The pruned EEG calculated for each 120 second trial of EEG (i.e iteration of CUDAICA) was concatenated, and subsequently the entire raw ictal EEG was low pass filtered (<16 Hz) using a 500th order symmetric digital FIR filter, and the resulting low pass filtered EEG was reconstituted with the high pass (>16 Hz) filtered EEG, following the exclusion of the independent components suspected to represent muscle artifact. The reconstituted and modified ictal EEG was exported from Matlab format to EDF for subsequent visual analysis.

Operation

All computations were carried out using compiled Matlab 8.4 custom scripts on a cluster of HP SL230s Gen 8 ES-2670 nodes with dual-eight-core 2.6 GHz Intel ES-2670 central processing units, 4 GB of memory per core, and NVIDIA Tesla graphics processing units. Minimal system requirements for operating AR2 include Matlab v8.4 or above, an Intel Xeon CPU, 2 GB of memory, a NVIDIA GPU, which is CUDA compatible, and CUDAICA. For scalp EEG files exported from Neuroworkbench (Nihon-Kohden, Irvine, CA, USA), executing the AR2 software method requires only inputting the file name of the EDF file of interest at the command line. For EDF files exported from other equipment manufacturers, full automation of the AR2 software method can be easily accomplished with simple modifications of the input parameters.

Patients and sample selection

We tested AR2 retrospectively using 23 seizures from eight adult patients with suspected focal-onset seizures undergoing evaluation at the UCLA Seizure Disorder Center. The patients and seizures were selected by S.A.W, whom was not a reviewer, from a review of consecutive clinical neurophysiology case conference presentations between January 1, 2014 and December 1, 2015 and based on case conference consensus that the ictal EEG records were uninterpretable due to muscle artifact contamination when reviewed with conventional digital filtering. For each of these patients, between 1–4 uninterpretable seizures were selected for inclusion in the study on the basis of a lack of ECG, electrode, and salt bridge artifact by S.A.W. Clinical data for each patient included seizure semiology, inter-ictal epileptiform abnormality, unobscured findings and radiological reports from MRI, PET, SPECT scans. The EEG and clinical records were deidentified and research informed consent was not required. This study was approved under UCLA IRB#15-001481. The video EEGs were acquired using a EEG-1200 amplifier (Nihon-Kohden, Irvine, CA) at a sampling rate of 200 Hz. Electrodes were placed according to the 10–20 international system with the additional anterotemporal electrodes at T1/T2. The duration of the exported EEG recording included the entire seizure and a 2-3 minute peri-ictal epoch.

Muscle artifact removal

AR1 was the commercially available Persyst v12 artifact reduction software17 (Persyst Development, San Diego, CA). The methods are proprietary. AR2 was developed by S.A.W and involved a two-step procedure consisting of a custom algorithm. EEG processed by AR2 was also interpreted using the Persyst v12 artifact reduction software.

Performance measures of AR1 and AR2

The ictal recordings for AR1 and AR2 were reviewed in Persyst v12 without video by 26 neurologists with a specialization in EEG. The readers were blinded to which records received AR1 or AR2, and each reader reviewed the 46 seizures in random. Following review of each ictal record, the reader completed a multiple choice questionnaire (Supplementary File 1), which assessed ability to visualize seizure-onset (Y,N) lateralize seizure-onset (L,R,N), locate the region of ictal onset (anterior temporal, anterior frontal, mid-temporal, temporal-parietal-occipital, occipital, none), and self-identify confidence of interpretation on a 5 point scale [(5) entirely confident (4) somewhat sure (3) probable (2) not confident (1) unlikely i.e. slight probability] for each measure. When time of onset, laterality, or the seizure onset region was not assigned the confidence was taken as (0). Readers were not provided with a definition of seizure-onset.

EEG analysis

During the interpretation of the ictal EEG processed by AR1 or AR2, no restrictions were placed on the use of Persyst v12 built in EEG filters (low-pass, high-pass, band-pass), or changes to montage. A comment in each recording was used to demarcate the time prior to the clinical seizure but not the EEG onset. The assessment was not time limited.

Statistical analysis

Differences in EEG interpretation utilizing AR1 and AR2 were assessed using the Wilcoxon signed rank test and the McNemar test on paired nominal data. Agreement across readers (Y,N,L,R), using either AR1 or AR2, was calculated using the inter-class correlation coefficient (ICC). For these outcomes, missing values were imputed to be in between non-missing values, and were analyzed using cumulative logit mixed effects models, which capture this ordering in the values and accounts for the clustering of readings into patients, and seizures within patients. Agreement across readers for onset region was calculated using a Fleiss kappa and treating the missing values as a category of response. Errors are given as standard error of the mean (s.e.m), unless otherwise specified.

Results

Implementation of the AR2 method

We applied the AR2 method developed at UCLA to the 23 seizures in the dataset. The method was automatic and unsupervised and separated the high-pass filtered (> 16 Hz) scalp EEG recordings into putative neurogenic and myogenic components (Figure 1). After pruning the putative myogenic components, the putative neurogenic components were reconstituted with the low-pass filtered (< 16 Hz) scalp EEG (Figure 2). The AR2 and AR1 processed scalp EEG recordings were subsequently inspected by the 26 specialists (Figure 3).

be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure1.gif

Figure 1. The AR2 method automatically separates independent components containing myogenic from neurogenic potentials.

The AR2 method automatically separates independent components containing myogenic from neurogenic potentials in the beta and gamma band on the basis of spatial topography and explained variance. A. Unprocessed scalp ictal EEG recording that was deemed uninterpretable. B. The same epoch after applying a low pass (<16 Hz) filter demonstrating a lack of a convincing ictal rhythm. C. The ictal epoch after applying a high pass (> 16 Hz) filter demonstrating dense muscle artifact. D. An example of a mutual information adjacency matrix calculated during an epoch of artifact in the high pass (> 16 Hz) filtered scalp EEG recording. Three scalp electrode recordings exhibited relatively low mutual information with all other electrodes and were designated poor quality and excluded from further processing to optimize INFO-MAX ICA based artifact reduction. E. The inverse weight matrix, and normalized inverse weight matrix (panel F) of all independent components across scalp electrode recordings for the seizure in panel A. Independent components 1-13 exhibited strong focality and were designated as containing myogenic potentials, while independent components 14 and above were designated neurogenic.

be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure2.gif

Figure 2. Ictal onset is revealed with reconstitution of the low pass (<16 Hz) ictal scalp EEG with the high pass (>16 Hz) neurogenic independent components.

Reconstitution of the low pass (<16 Hz) ictal scalp EEG with the high pass (>16 Hz) neurogenic independent components reveals a clear ictal onset. A. The tentative neurogenic independent components (A1) and myogenic independent components (A2) derived from INFO-MAX ICA processing of the high pass (> 16 Hz) filtered ictal scalp EEG recording. The largest amplitude activity in the neurogenic components are evident frontally and in the left hemisphere. B. The low pass filtered ictal scalp EEG suggests a possible left frontal onset but a convincing ictal rhythm is lacking. C. Reconstitution of the low pass EEG with the neurogenic high pass (> 16 Hz) independent components results in an ictal EEG that demonstrates a more convincing left frontal onset consisting of beta-gamma oscillations with some clear phase reversals in F3 and F7.

be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure3.gif

Figure 3. A comparison of the results of artifact reduction methodologies.

Ictal scalp EEG recording from seizure 18 prior to artifact reduction processing (top), after processing with artifact reduction methodology 1 (AR1, middle), and after processing with artifact reduction methology 2 (AR2, bottom). Only processing with AR2 reveals a right hemispheric onset followed by clear spread to right frontal regions.

Identifying time of seizure-onset

Across the 23 seizures considered previously uninterpretable with digital filtering (Table 1) two-thirds of the readers were able to delineate the time of seizure-onset for 10 of the 23 using AR1, and 15 of the 23 using AR2 (Figure 4A, n=23, p<0.01). Across the 23 seizures, the mean confidence measure for the determination of seizure-onset was 2.68 +/- 0.19 (probable-not confident) when AR2 was utilized and 2.19 +/- 0.18 (not confident) with AR1 (Figure 5A, p<0.01). The inter-class coefficient (ICC) was 0.26 (95% Confidence Interval (CI) 0.21-0.30) with AR2, and 0.15 (95% CI 0.11-0.18) with AR1 (p=0.333).

Table 1. Clinical description of patients.

Clinical description of patients and ictal EEG laterality and focus assignments with AR1 and AR2. Abbreviations (L:left, R:right), PET findings refer to hypometabolism, SPECT findings to hyperperfusion. The focus was determined on a majority basis across all the assignments made by the readers for a subject’s seizure(s).

PatientAge
Gender
Aura/SemiologyIEDsUn-
Obscured
Seizure
Onset
Laterality
sMRIPET/SPECTSeizure Onset
or Spread
Laterality (AR2)
Seizure
Onset or
Spread
Laterality
(AR1)
AR2 focusAR1 focus
#1 46MSomato-sensory
(warmth)/arousal from
sleep, hyperkinetic,
noneleft frontal
ictal rhythm
nonlesionalnormal PET,
SPECT left insula
1. 14/21 L1. 17/19 Lant/mid
temporal
ant
temporal
#2 32MSomato-sensory
(discomfort)/right facial
grimacing, right leg
elevation, breath holding
nonenonenonlesionalPET right
temporal, SPECT
bilateral frontal
lobes
2. 6/10 L
3. 6/8 L
4. 13/16 L
2. 6/7 R
3. 6/6 R
4. 7/7 R
ant.Frontalmid.
Temporal
#3 23MTachycardia/arousal from
sleep, hyperkinetic, b/l
dystonic posturing
nonenonenonlesionalnormal PET5. 8/10 R
6. 14/16R
7. 7/11 R
5. 6/11 L
6. 6/6 L
7. 6/6 R
frontal/mid
temporal
mid
temporal
#4 53MVisual disturbance/
behavioral arrest, cursing,
right arm dystonic
posturing
L temporalL anterior
temporal
L MTS, L parietal
encephalomalacia
PET L parietal8. 18/19 L
9. 16/21L
10. 17/21 L
8. 14/17 L
9. 10/13 L
10. 21/23 L
ant/mid
temporal
ant
temporal
#5 20MVague/right head and eye
version, right arm clonic
movements,
L temporalL temporalL frontal
polymicrogyria
normal PET11. 11/18 R
12. 7/11 L
13. 9/10 L
11. 9/18 R
12. 10/13 L
13. 7/11 R
ant/mid
temporal
ant/mid
temporal
#6 27MNone/arousal from sleep,
dyscognitive, right head
and body version.
L frontalL frontalnormalPET L inferior
frontal
14. 20/22 L
15. 21/24 L
16. 22/24 L
14. 25/25 L
15. 24/24 L
16. 24/25 L
ant frontal/
ant temporal
ant frontal/
ant
temporal
#7 26FNone/nocturnal arousal
or daytime events,
hyperkinetic, right or left
dystonic posturing
L and R
temporal
NoneRight middle cranial
fossa arachnoid
cyst
PET R parietal
lobe
17. 21/23 R
18. 12/16 R
19. 21/23 R
20. 20/21 R
17. 20/23 R
18. 18/21 R
19. 14/16 R
20. 6/11 R
ant/mid
temporal
ant/mid
temporal
#8 19MLightheaded/loss
of consciousness,
right > left arm clonic
movements, and
posturing
L and R
temporal
NoneL mesial
temporal CD,
R>L gyrus rectus
encephalomalacia
PET L>R
temporal lobe
21. 12/16 R
22. 12/23 R
23. 23/26 L
21. 14/16 L
22. 23/24 L
23. 22/22 L
ant/mid
temporal
ant/mid
temporal
be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure4.gif

Figure 4. More readers could lateralize seizure onset utilizing AR2 as compared to AR1.

More readers could visualize the time of seizure onset, and assign laterality to seizure onset utilizing AR2 as compared to AR1, and the assigned laterality of seizure onset sometimes differed between the two methods. A. Bar plot of the number of readers whom visualized the time of onset for each seizure utilizing AR1 (blue) or AR2 (red). Across seizures more readers visualized seizure onset utilizing AR2 compared with AR1 (p<0.01). Asterisks indicate statistically significant differences between the two methods in individual seizures (McNemar, p<0.05). B. Stacked bar plot of the number of readers selecting a left- or right-sided seizure onset utilizing AR1 (light blue, left; light yellow, right) or AR2 (dark blue, left; yellow, right). Across seizures more readers lateralized seizure onset utilizing AR2 compared with AR1 (p<0.01). Asterisks indicate statistically significant differences in individual seizures (McNemar, p<0.05), number sign indicates a significant change in the determination of laterality utilizing AR2 compared to AR1 (McNemar, p<0.05).

be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure5.gif

Figure 5. Confidence in the interpretation of ictal EEG onset improves with utilization of AR2 as compared to AR1.

A. Bar plot of the mean confidence scale values for visualizing the time of seizure onset for the 23 seizures interpreted utilizing AR1 (blue), and AR2 (red). Across seizures, confidence scale values were greater when AR2 was utilized as compared with AR1 (p<0.01). Asterisks indicate differences in confidence values in individual seizures (p<0.05). Error bars are calculated as s.e.m. B. The respective mean confidence scale values for seizure onset lateralization. C. The respective mean confidence scale values for seizure focus localization. Across seizures, confidence scale values for lateralizing seizure onset, and identifying the seizure focus were greater when AR2 was utilized as compared with AR1 (p<0.05).

Lateralizing and localizing seizure-onset

Compared with identifying the time of seizure-onset, fewer readers could lateralize seizure-onset after either AR1 or AR2 (Figure 4B, p<0.05). However, more readers were able to lateralize seizure-onset using AR2 compared to AR1 (Figure 4B, p<0.05) and readers were more confident with AR2, although both methods did not produce high levels of confidence. The mean confidence measure for seizure-onset lateralization was 1.87+/- 0.198 (not confident-unlikely) for AR2 and 1.54+/- 0.176 (not confident-unlikely) for AR1 (Figure 5B, p<0.01). The ICC was equivalent (p=0.501) for AR1 (ICC=0.33 95% CI 0.30-0.37) and AR2 (ICC=0.28 95% CI 0.25-0.31). For localizing the region of seizure-onset reader confidence (Figure 5C), and agreement was very low (Figure 6, AR1 Fleiss’ kappa = 0.1199, 95% CI = 0.116-0.124, AR2 Fleiss’ kappa = 0.121, 95% CI =0.118-0.125). For two of the seizures, the laterality assignments were different when AR2 was used as compared to AR1 (Figure 4B, McNemar p<0.05).

be01c09c-d7e5-493f-bf78-4aac56d43e4f_figure6.gif

Figure 6. Differences in ictal onset region assignments using AR1 or AR2.

Stacked bar plot of the ictal onset region assignments using either AR1 (lighter colors) or AR2 (darker colors) for all 23 seizures. Overall, across seizures, more readers were able to render an assignment using AR2 as compared to AR1 (p<0.05). Inter-reader agreement using for assigning the ictal onset region was marginal using either AR1 or AR2.

Comparison of seizure-onset lateralization assignments with other clinical findings

We identified the patients with at least two consistent clinical findings that lateralized the suspected seizure-onset zone (SOZ). Compared to AR1, more readers were able to render seizure-onset laterality assignments using AR2, and these assignments were more often congruent with other clinical data (Table 2). These clinical findings included seizure semiology, onset of seizures without EEG obscuration, structural MRI, PET, or SPECT findings. If any of the clinical findings were contradictory with respects to the laterality of the suspected SOZ, the SOZ was designated unknown. Overall, 4 patients (#1,4,5,6) had clinical findings that supported a left-hemispheric SOZ, and 1 patient (#7) had clinical findings that supported a right-hemispheric SOZ (Table S1). Among the 8 patients, if the reader lateralized the seizure-onset to the left using AR2 they were correct in 95.9% (95% CI 85.7-98.9%) of cases, but using AR1 they were correct in 91.9% (95% CI 77.0-97.5%) of cases (Table 3, p<0.0607).

Table 2. Contingency table of agreement between assigned seizure onset laterality and other clinical findings.

Contingency table of the agreement between seizure-onset laterality using AR1 (left), and AR2 (right) and the laterality of seizure-onset assigned on the basis of other clinical data for all the study patients and seizures. Note that clinical seizure-onset lateralization was not available for all patients, and when readers rendered a laterality decision that matched the laterality based on other clinical data, the assignments “agreed”.

AR1AR2
EEG seizure-onset
lateralization
EEG seizure-onset
lateralization
YNYN
Clinical seizure-
onset
lateralization



Y
 
Agree

145
 
Disagree

32
 


187
 
Clinical seizure-
onset
lateralization



Y
 
Agree

171
 
Disagree

39
 


154
 
N
 
83
 
151
 
N
 
107
 
127
 

Table 3. Agreement between seizure-onset laterality and other clinical findings.

Agreement between seizure-onset laterality assignments using either AR1 or AR2 and the suspected laterality of the SOZ assigned on the basis of other clinical data. Parentheses indicate the 95% confidence interval. “n” refers to the number of subjects.

Artifact
Reduction
Method
Reader Assignment
of Seizure-Onset
Laterality
Percentage of reader assignments
in concordance with SOZ
laterality defined by other clinical
criteria.
AR1Right 59.3 (28.5-84.2) (n=1)
Unknown66.8 (38.1-86.9) (n=3)
Left91.9 (77.0-97.5) (n=4)
AR2Right61.8 (31.3-85.1) (n=1)
Unknown71.4 (42.8-89.3) (n=3)
Left95.9 (85.7-98.9) (n=4)

Discussion

In this study, we present a new artifact reduction software, AR2, and its application compared with a commercially available tool, AR1. 26 neurologists used the two methods to interpret 23 ictal EEG recordings that were uninterpretable due to muscle artifact when reviewed with conventional filtering. The major findings from this study include: 1) the utilization of artifact reduction software results in non-uniform interpretation of ictal EEG, with many readers not able to render assignments; 2) when readers did render seizure-onset laterality assignments it often agreed with other clinical findings; 3) although the study size was small, the AR2 software method increased the number of readers that rendered assignments, and reader confidence suggesting it aids in diagnosis.

Both AR1 and AR2 are digital signal processing software tools8,15,17 that may confound accurate ictal EEG interpretation by altering the appearance of the EEG. Digital filtering also can mislead5. One concern about AR1 and AR2 relates to the lack of understanding of the waveform alteration. Specifically, the readers were not confident in their interpretations, and the determination of seizure lateralization sometimes differed between the AR1 or AR2 methods. As such, the artifact reduction methods may introduce false positive findings. This demonstrates the limits of EEG artifact reduction approaches and puts the advantages into perspective.

Neurologists often disagree on the interpretation of ictal EEG processed with artifact reduction software, however the seizure-onset laterality assignments rendered by a quorum are often correct. Further refinement of this technology may successfully improve the efficiency of video-EEG monitoring and the utilization of epilepsy surgery; however, correlation with epilepsy resective surgery outcomes will be required for further validation.

With regard to AR2, the novel software method developed for this study, the slight improvement seen in ictal EEG interpretability after applying the method suggests that the algorithm can (1) reliably produce signals that are, exclusively or mainly, EEG or EMG, and (2) identify which signals are of brain origin and which are contaminant.

One explanation for AR2’s ability to isolate myogenic from neurogenic independent components may be that scalp EEG electrodes record weighted and summated far-field signals from all brain and muscle sources, as well as near-field electrode noise generated at the electrode/skin interface. The decomposition of scalp EEG data into components with maximally independent time courses using independent component analysis results in time series that may resemble single equivalent dipoles because of the bias towards increased local connectivity in neurons and myocytes as compared to long distance connectivity14.

Data and software availability

All software code for the new AR2 software developed by S.A.W. is openly and permanently available at https://github.com/shennanw/AR2.

Archived source code as at time of publication: doi, 10.5281/zenodo.22989321

License: GNU Public License 3.

The raw scalp ictal EEG files that were analyzed in this study using AR2, as well as the scalp ictal EEG files following processing using AR2 are available from Zenodo: Dataset 1. Validity of two automatic artifact reduction software methods in ictal EEG interpretation. Doi, 10.5281/zenodo.22109522 (https://www.zenodo.org/record/221095#.WF63m7YrLdR)

The raw data used for the comparative assessments are available from Zenodo: Dataset 2. Validity of two automatic artifact reduction software methods in ictal EEG interpretation. Doi. 10.5281/zenodo.223329 (https://zenodo.org/record/223329#.WHN-HLYrLdQ)

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 10 Jan 2017
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Weiss SA, Asadi-Pooya AA, Vangala S et al. AR2, a novel automatic artifact reduction software method for ictal EEG interpretation: Validation and comparison of performance with commercially available software [version 1; peer review: 2 approved with reservations] F1000Research 2017, 6:30 (https://doi.org/10.12688/f1000research.10569.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 10 Jan 2017
Views
19
Cite
Reviewer Report 16 Mar 2017
David M. Groppe, Department of Psychology, Toronto, ON, Canada 
Approved with Reservations
VIEWS 19
In this manuscript, Weiss and colleagues present a novel algorithm for removing electromyographic (EMG) artifacts from ictal EEG recordings, called AR2. Moreover, they evaluate the performance of the algorithm on data from 8 patients and compare it to a similar ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Groppe DM. Reviewer Report For: AR2, a novel automatic artifact reduction software method for ictal EEG interpretation: Validation and comparison of performance with commercially available software [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:30 (https://doi.org/10.5256/f1000research.11390.r20245)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 04 Apr 2017
    Shennan Weiss, Department of Neurology, Jefferson Comprehensive Epilepsy Center, Thomas Jefferson University, Philadelphia, USA
    04 Apr 2017
    Author Response
    Dear Dr. David Groppe,
     
    We are grateful for your insightful and thoughtful comments and suggestions. Appended below are answers to your inquiries, and changes we have made to the ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 04 Apr 2017
    Shennan Weiss, Department of Neurology, Jefferson Comprehensive Epilepsy Center, Thomas Jefferson University, Philadelphia, USA
    04 Apr 2017
    Author Response
    Dear Dr. David Groppe,
     
    We are grateful for your insightful and thoughtful comments and suggestions. Appended below are answers to your inquiries, and changes we have made to the ... Continue reading
Views
27
Cite
Reviewer Report 22 Feb 2017
Patrícia Figueiredo, Institute for Systems and Robotics and Department of Bioengineering, Instituto Superior Técnico University of Lisbon, Lisbon, Portugal 
Rodolfo Abreu, Institute for Systems and Robotics and Department of Bioengineering, Instituto Superior Técnico University of Lisbon, Lisbon, Portugal 
Approved with Reservations
VIEWS 27
Referee report
The manuscript needs careful revision by a native English speaker within the scientific community. Although I feel that the performance measures used by the authors are adequate, and that a substantial number of EEG specialists quantified them, ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Figueiredo P and Abreu R. Reviewer Report For: AR2, a novel automatic artifact reduction software method for ictal EEG interpretation: Validation and comparison of performance with commercially available software [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:30 (https://doi.org/10.5256/f1000research.11390.r20461)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 04 Apr 2017
    Shennan Weiss, Department of Neurology, Jefferson Comprehensive Epilepsy Center, Thomas Jefferson University, Philadelphia, USA
    04 Apr 2017
    Author Response
    Dear Dr. Figueiredo and Dr. Abreu,

    Thank you very much for your thoughtful and helpful comments and suggestions. We have substantially revised the manuscript according to your feedback as ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 04 Apr 2017
    Shennan Weiss, Department of Neurology, Jefferson Comprehensive Epilepsy Center, Thomas Jefferson University, Philadelphia, USA
    04 Apr 2017
    Author Response
    Dear Dr. Figueiredo and Dr. Abreu,

    Thank you very much for your thoughtful and helpful comments and suggestions. We have substantially revised the manuscript according to your feedback as ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 10 Jan 2017
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions