- Open Access
Single trial somatosensory evoked potential extraction with ARX filtering for a combined spinal cord intraoperative neuromonitoring technique
BioMedical Engineering OnLinevolume 6, Article number: 2 (2007)
When spinal cord functional integrity is at risk during surgery, intraoperative neuromonitoring is recommended. Tibial Single Trial Somatosensory Evoked Potentials (SEPs) and H-reflex are here used in a combined neuromonitoring method: both signals monitor the spinal cord status, though involving different nervous pathways. However, SEPs express a trial-to-trial variability that is difficult to track because of the intrinsic low signal-to-noise ratio. For this reason single trial techniques are needed to extract SEPs from the background EEG.
The analysis is performed off line on data recorded in eight scoliosis surgery sessions during which the spinal cord was simultaneously monitored through classical SEPs and H-reflex responses elicited by the same tibial nerve electrical stimulation. The single trial extraction of SEPs from the background EEG is here performed through AutoRegressive filter with eXogenous input (ARX). The electroencephalographic recording can be modeled as the sum of the background EEG, which can be described as an autoregressive process not related to the stimulus, and the evoked potential (EP), which can be viewed as a filtered version of a reference signal related to the stimulus. The choice of the filter optimal orders is based on the Akaike Information Criterion (AIC). The reference signal used as exogenous input in the ARX model is a weighted average of the previous SEPs trials with exponential forgetting behavior.
The moving average exponentially weighted, used as reference signal for the ARX model, shows a better sensibility than the standard moving average in tracking SEPs fast inter-trial changes. The ability to promptly detect changes allows highlighting relations between waveform changes and surgical maneuvers. It also allows a comparative study with H-reflex trends: in particular, the two signals show different fall and recovery dynamics following stressful conditions for the spinal cord.
The ARX filter showed good performances in single trial SEP extraction, enhancing the available information concerning the current spinal cord status. Moreover, the comparison between SEPs and H-reflex showed that the two signals are affected by the same surgical maneuvers, even if they monitor the spinal cord through anatomically different pathways.
The monitoring of the functionality of vital parameters is known to be a basic aspect in up-to-date surgical techniques. This is particularly true in the case of surgery performed on the vertebral column (i.e.: correction of serious scoliosis), when a continuous monitoring of the spinal cord functionality is usually required [1–4].
In fact, during surgery the integrity and functionality of the spinal cord may be compromised by surgical maneuvers to such an extent as to cause prolonged insufficient blood supply to the cord, or mechanical compression and, as a consequence, its temporary malfunctioning. Sometimes this malfunctioning, if undetected at an early stage, can lead to irreversible spinal cord damage [5–7].
The aim of this study is to present a signal processing method used for extracting a single trial potential in a new combined neuromonitoring technique , based on the simultaneous monitoring of parameters extracted from two neurophysiologic signals: the single trial somatosensory evoked potentials (SEPs) [9, 10] and the soleus muscle H-reflex [11, 12], both elicited by the same electrical stimulus to the posterior tibial nerve (Fig. 1). The advantage of observing two different signals both elicited by the same stimulus, instead of observing the variation of only one of them, is that the variability of the monitoring is decreased. Moreover the SEP provides information on the functionality of the ascending pathways, up to the somatosensory cortex, while the H-reflex involves mainly the descending motor pathways. Thus the pieces of information coming from the two systems are complementary, and the use of both can provide a more reliable and robust monitoring, especially when real-time processing and immediate decisions are required as in a surgical theatre. Indeed the SEP and the H-reflex involve different nervous fibers which however share the same peripheral nerve and both depend on the integrity of spinal cord pathways. It is therefore reasonable to expect that these two signals should show an independent variability but a damage to the peripheral or medullar function would modify the two signals in the same direction hence, it is hypothesized that the real time comparison of the SEP and the H-reflex, both elicited by the same stimulus, should significantly increase the monitoring reliability and sensitivity, as suggested in [13–15].
However, whereas the H-reflex analysis does not require refined processing methods, the SEPs due to their poor signal to noise ratio (SNR), require proper tools for the analysis of the EEG in order to extract the evoked potential waveform.
The current clinical SEP extraction techniques usually average several trials in order to enhance the SNR. However, the SEP has a trial-to-trial variability which is lost if it is extracted using the averaging technique. This would lead to the loss of information about the trial-to-trial changes of the SEP signal.
When combining SEPs and H-reflex, a limit is introduced by the long recovery time of the H-reflex, so that the inter-stimulus interval (ISI) must be at least 10 s -. In case of a SEPs' sudden changes the averaging technique provides intolerably delayed information, for this reason, a single trial SEP analysis is needed in order to compare it with the H-reflex obtained from the same electrical stimulus.
In this study, recordings from 8 scoliosis correction surgery cases are utilized. The correction was performed at the same spinal level with similar surgical procedures: 6 of them were corrected with Cotrel-Dubousset technique; the other 2 were corrected with the Harrington technique (Table 1). In all cases the same i.v. Propofol anesthesia was used. All analyzed cases did not present any post operative complication after recovery from anesthesia.
EEG trials for the evoked potentials and EMG for the H-reflex were acquired during surgery form each of eight patients. The signals were elicited stimulating the posterior tibial nerve with a 1 ms rectangular voltage pulse up to 150 V at 0.1 Hz repetition rate.
The EEG signal was recorded through needle electrodes placed on the scalp in Cz' with an auricular reference. The H-reflex was acquired using a surface electrode (RedDot, 3 M, USA) on the soleus muscle and a reference one over the Achilles tendon.
The EEG and EMG acquisition were performed by means of two analog ampli-filters ICP 511-AP by Grass Instruments, with different gains and filter settings (EEG: gain = 80dB, bandpass 1 Hz – 300 Hz; EMG: gain = 40 dB, bandpass 30 Hz – 3 kHz).
The stimulus triggered EMG recording lasts 62.5 ms, sampled at 10 kHz; the EEG recording lasts 125 ms, sampled at 2.5 kHz by PCI board (type PCI-MIO-16E-4, National Instruments, USA). The acquisition board range is +/-1.25 V with a 12-bit resolution.
The analysis software used for off-line processing the recorded data has been implemented by means of Virtual Instruments developed in LabView with integrated Matlab processing scripts. The software performs the single trial extraction of SEP waveform and the detection of amplitude and latency both of the SEP and of the H-Reflex simulating an on-line neuromonitoring system.
EEG signal processing
An Autoregressive Model with an eXogenous Input (ARX) has been chosen, because it has been successfully applied and tested in the analysis of different kinds of evoked potentials (visual, acoustic, somatosensory [23–25]).
According to the ARX model, the signal y(k), recorded after a somatosensory stimulation, can be modeled as the sum of two different contributions: a former not related to the stimulus (the background EEG) that, according to the literature , can be described as an autoregressive process driven by a white noise e(k); a latter related to the stimulus (the EP), that can be viewed as a filtered version of a reference signal u(k). The general equation of an ARX process is:
The model orders (n, m, d) are defined as follows:
n and m, that are respectively the model orders of the autoregressive and of the moving average sections;
d, takes into account the possible temporal anticipation between the reference input and the output of the filter (the filter is not causal).
In the z-transform domain, where z-1 is the unit delay operator, the ARX model is expressed as:
A (z)·Y = B (z)·U + E (2)
This last expression (3) is visualized by the block diagram shown in Fig. 2.
If u(k) and y(k) are known, then the a i and b j coefficients can be estimated together with the variance of the input white noise, for every pre-fixed model order set and delay (n, m, d), by means of a least squares approach, that minimizes the following figure of merit:
where err(k) is the prediction error of the model and N is the number of samples in the trial. The prediction error is the difference between the recorded signal and its prediction according the model expressed in eq. 1.
Choice of the model order set (n, m, d)
The method used for selecting the optimal model order set has been implemented in the following steps:
Choice of an initial set of model orders (n, m, d);
Minimization of a figure of merit J in order to find the best ARX coefficients for the chosen order.
Check the whiteness of prediction error err(k): the set of identified parameters is accepted for the following calculations when the err(k) is a white sequence with a confidence level of 95%; the whiteness of the residual is tested using the cumulative Anderson test.
If the whiteness test is not verified, the model orders are increased until the prediction error is a white noise.
The starting triplet (n, m, d), selected in the first step, is calculated for each subject before starting the monitoring phase of the risky surgical procedures.
The first 50 sweeps are recorded and stored, and the optimal model orders, among the ones that satisfy the whiteness criterion, are selected by minimizing the Akaike Information Criterion (AIC) function; though there are several methods for finding the optimal orders for model (Final Prediction Error, Minimum Description Length, etc.), for the specific application the best one is the AIC because where the model orders are low the AIC shows better performances .
The AIC function is expressed as:
where σ2 is the variance of the prediction error err(k).
The d parameter is chosen as m/2 considering a previous investigation .
On the basis of the AIC function, evaluated on the first 50 recorded sweeps, the triplet (n, m, d) is selected for the specific patient.
The reference signal
In many applications of the ARX filter for the extraction of the single trial EP, the reference u(k) is the average of a sufficient number of trials, representing a standard EP shape. The choice of the reference signal is fundamental for the ARX filter capability in tracking the inter-trials EP variations: in fact the filter tries to identify a waveform, similar to the input reference, in the recorded signal trial y(t). Thus, a fixed reference could be misleading in the case of a dynamic process. For this reason a fast adapting reference could be a better choice in order to be able to track also sudden changes in the waveform.
However, the ability of this reference to adapt to morphological variations in the EP wave depends on the number of sweeps included: the moving average introduces a delay in the updating of the reference that depends on the number of averaged trials. In order to enhance the adaptation ability of the reference signal, a moving average with an exponential decay is here proposed and tested.
The reference is computed as follows:
u i = μ·u i - 1+ (1 - μ)·y i (6)
The μ coefficient (forgetting coefficient) provides an exponential windowing on the averaged sweeps, so that only the most recent trials contribute significantly to the reference u i , while the oldest ones are progressively forgotten, with an exponential decay. Different values of μ have been tested in order to have a reference signal that adaptively varies according to the possible dynamical changes in the evoked potentials. Off-line tests led to the value μ = 0.95 .
In order to compare the ARX performances with the reference built on a static average and with reference built on a moving average driven by a forgetting coefficient, as defined by the eq. (6), signal processing simulations have been performed.
Three known input signals, representing the x i (k) target signals, have been built from one generic EP that has been multiplied for 3 different coefficients: k 1 = 20, k 2 = 15, k 3 = 1 in order to simulate different amplitudes.
An AR model has been identified from an EEG trial of 125 ms not related to any stimulus. A white Gaussian noise has been utilized as the AR model input in order to generate EEG signals n i (k). With this method it was possible to collect, as the sum of random n i (k) and known targets signals kj*x i (k), 50 trials using the k 1 coefficient, 1 using the k 2 and 1 using the k 3 . These input signals have been utilized for simulating a possible variation of the EP due to a stressful maneuver applied at the 51st trial that produces a strong EP decrease at the 52nd trial.
The ARX filter using the two different reference signals, i.e. the moving average and the moving average with forgetting coefficient, have been compared according to their ability to extract the input signals kj*x i (k).
Correlation with the H-reflex amplitude
The degree of correlation between the amplitude of SEP and H-reflex has been computed through the Pearson's correlation coefficient in the 8 surgery sessions studied.
The Pearson's coefficient is defined as:
where As and Ah are the amplitudes of SEP and H-reflex, A and A are the averages of the respective amplitudes, Ss and Sh are the standard deviations and n is the number of the considered trials. Pearson's coefficient value indicates the correlation level for a specific freedom degree given by the number of trials analyzed; a coefficient value of 1 indicates a perfect correlation whereas a value of 0 indicates no correlation.
Data analysis protocol
The sequence of the off-line data analysis steps have been structured in order to comply with the requirements for a data processing during a surgery session.
The data analysis steps proposed are:
Cleaning corrupted recorded trials
This procedure automatically deletes the trials corrupted by environmental noise by checking 2 different parameters:
presence of saturated samples: if there is at least one sample with a value of +/-1.25 V (the maximum/minimum acquisition board voltage input), the trial is rejected;
single-trial amplitude range: given the maximum sample value (Max) and the minimum sample value (Min) of the ith-trial, if |Max-Min|>0.7 V, the trial is rejected.
All the parameter checking and its respective threshold values have been experimentally found.
Building the initial reference signal
The first reference signal is obtained from 50 trials recorded before the beginning of risky surgical maneuvers. Then, the reference signal is updated at each trial throughout the surgery according to the weighted moving average procedure described above.
Choosing the optimal order for the model
The trial set, used to build the initial reference signal, is also used to evaluate the best model order. The following range of model orders have been used in order to calculate the AIC function:
[1<n<21], [2<m<20] and d = m/2.
Single sweep analysis.
The trials acquired during the surgical procedure are evaluated in combination with the H-reflex signal.
The 3 input test trials computed by the AR model built from a generic EEG trial have been utilized as target input signals for single trials extracting simulation (Fig. 3). It is very difficult to valuate without any signal processing that the 52nd unfiltered trial has a known input EP signal one tenth smaller than the 50th trial.
The 52nd EP input is correctly extracted from the raw signal by the ARX filter both using, as input model reference, the average of the previous 50 trials or using the moving average with a forgetting coefficient μ = 0.95 (Fig. 4).
However the difference due to the reference is clearly visible in the extraction of the 51st trial. The 51st raw trial, as described above, is built from an EP reduced by 20% in respect to the EP of the 50th trial. Utilizing the moving average with a forgetting coefficient, it is possible to notice a variation of extracted EP of about 50% on P30 peak amplitude whereas using as reference the average of the previous 50 trails the peak shows a decrease of only 5%.
Three significant trials have been chosen from one of the eight cases studied (described in further detail in the next section) in order to show how the reference signal, chosen for the ARX filter, affects the tracking of the changes in the SEP waveform after a sequence of hammer blows at the vertebral column (Fig. 5). Using the exponential average with the forgetting factor as ARX reference signal (Fig. 5A) the tracking of the sudden waveform change is faster than using the classical moving average (Fig. 5B). This increased sensibility to the EP variation is also observed in the signal processing simulation previously presented (Fig. 4).
Considering the correlation of the P30 EP amplitude variation with respect to that of the H-reflex one, the Pearson's coefficient shows that in all of the 8 cases analyzed the two sets of data are correlated with a confidence level p < 0.01 (Table 2).
A surgery case
The procedure followed in the signal processing method is here outlined using data obtained from one of the eight cases presented in Table 1.
The initial ARX reference signal has been built using 50 trials recorded before the beginning of the surgical procedures (Fig. 6). Reference values, specific to the subject, are obtained from the evaluation of the SEP main peak amplitude (P30) and its latency. The optimum model order is then chosen analyzing the AIC function: Fig. 7 shows the trend of AIC when n varies between 1 and 10 and for two different values of m, m = 2 and m = 4. The optimum values in this case are: n = 2, m = 4 and d = 2.
Fig. 8a shows a sequence of single trial SEPs during different surgery phases, while in Fig. 8b the corresponding H-reflex responses are plotted. The surgery history summary can be observed in Table 3.
In step 1 and step 2 there is no detectable changes in the monitored waveforms as they are recorded during the initial surgery phase and no stressful maneuvers have been performed yet.
At step 3 the vertebral column begins to be hammered and both signal amplitudes decrease considerably. The two signals, however, show different behaviors. The H-reflex recovers very quickly, whereas the SEP shows a slower recovery dynamic. The same behavior can be observed after other stressful maneuvers, as decortications and tractions.
Fig. 9 shows in better detail the amplitude trends of the two signals. As observed in Fig. 8 and 9, both signals are affected by surgical maneuvers such as hammer strokes and decortications. However, the amplitude of the H-reflex restores quickly to the initial values, while the SEP amplitude requires a longer recovery time and does not restore completely, showing a globally decreasing trend.
Conclusions and discussion
The application target was to implement a signal processing method able to monitor single trial SEPs during surgery and to be very sensitive to the evoked potential variations. The ARX filter, using a reference generated by a recursive exponential averaging, fulfills the requirements and provides better performances than the classical moving average reference. In order to test the method we used a lower high pass filter (1 Hz) than the one utilized in usual neuromonitoring applications (30 Hz) in order to test the signal processing method in less favorable conditions and demonstrating its capability to detect also small single potential variations.
The paper presents off-line results, but the method is designed for online application. During surgery, the time available to complete a single trial extraction is 10 sec, i.e the interstimuli interval. The proposed algorithm completes the analysis in less than 1 sec (running on Laptop with Intel P4 512MB) so it is suitable for online application. Another critical point is the choice of the model orders (n, m, d). The inter-subject variability imposes the calculation of the optimal orders for each new patient as described in the method section and this procedure has to be completed between the beginning of the surgery, when the patient falls asleep, and the starting of the surgical risky maneuvers (about 15 minutes). Actually this procedure is completed in less than 2 minute because, as described in the method section, it needs only 50 trials obtained at higher stimulation frequency (5 Hz). There is also some intra-individual variability, even if less critical than inter-subject variability. For this reason the test of the model is carried out each trial through the whiteness test of the prediction error err(k). In the event of a not whiteness of the err(k) the model orders are increased.
Other signal processing methods [29, 30] have been proposed in literature for SEP monitoring during surgery, but the presented application compared the single trial SEP with the H-reflex using the same stimulus delivered at very low frequency. For this reason we develop a signal processing method, based on previous works [22–26], which gives information at each single trial recorded. The method improves the sensibility to the SEP variation as presented in the result section. Though 8 surgery sessions are not enough to have clinical and medical results, some comments are still possible, even at this early stage of research. The analysis of these 8 clinical cases recorded engenders a general consideration about the correlation between SEP and H-reflex amplitudes that are statistically correlated as shown by the Pearson's coefficient. As expected, the Pearson's coefficient shows a fair correlation between the amplitude of the single trial SEP and the H-reflex. As described previously, the behavior of the two signals is quite different regarding the response of the spinal cord to a stressful maneuver. However, a correlation exists between the EP and H-reflex. The example above, in fact, shows that the same maneuver affects both signals. The presence of correlation between the changes that occur in both signals means that the changes are related to the same cause: in this case, this could mean that modification of the spinal cord functionality is induced by the same surgical maneuvers. Thus, the two signals, anatomically and functionally independent, generated by a single transcutaneous stimulus, are both related to the functionality of the spinal cord and are both affected by the same cause. Regarding the neuromonitoring technique where the method has been applied, the combined neuromonitoring has shown the expected behavior: even if SEPs and H-reflex monitor anatomically different pathways, when they are affected by the same surgical maneuvers, they show significantly correlated amplitude changes.
A larger set of monitored surgeries will be collected in order to obtain a strong and significant clinical result that confirms a higher reliability of this new combined monitoring technique.
Guérit JM: Neuromonitoring in the operating room: why, when, and how to monitor? Electroencephalography & Clinical Neurophysiology 1998, 106: 1–21. 10.1016/S0013-4694(97)00077-1
Laureau E, Marciniak B, Hebrad A, Herbaux B, Guieu JD: Neuromonitoring and anesthesia in surgery of the spine. Clin Neurophysiol 1998, 28: 299–320. 10.1016/S0987-7053(98)80002-6
Minahan RE: Intraoperative neuromonitoring. Neurologist 2002, 8: 209–226. 10.1097/00127893-200207000-00001
Guérit JM, Dion RA: State-of-art of neuromonitoring for prevention of immediate and delayed paraplegia in thoracic and thoracoabdominal aorta surgery. Ann Thorac Surg 2002, 74: 1867–1869. 10.1016/S0003-4975(02)04130-9
Lin J, Zhang J, Ye Q, Sheng J, Qiu G: Failure and complication following surgical treatment of scoliosis – an analysis of 101 cases. Chin Med Sci J 1999, 14: 174–179.
Cervellati S, Bettini N, Bianco T, Parisini P: Neurological complications in segmental spinal instrumentation: analysis of 750 patients. Eur Spine J 1996, 5: 161–166. 10.1007/BF00395507
Wiedemayer H, Fauser B, Sandalcioglu IE, Schafer H, Stolke D: The impact of neurophysiological intraoperative monitoring on surgical decisions: a critical analysis of 423 cases. J Neurosurg Spine 2000, 96: 255–262.
Bracchi F, Rossi L, Merzagora AC: Apparatus and Methods for Monitoring the Functionality of Spinal Cord. International Patent Pending: PCT/IT2004/000440
Nuwer MR: Spinal Cord Monitoring With Somatosensory Techniques. J Clin Neurophysiol 1998, 15: 183–193. 10.1097/00004691-199805000-00002
Lille F, Petit B, Margules S, Mazel C, Roy-Camille R: Somatosensory evoked potentials during spinal orthopedic surgery in adult patients. Neurophysiol Clin 1993, 23: 179–92. 10.1016/S0987-7053(05)80229-1
Bracchi F, Grossi P, Trovati L, Viganò P: H-Reflex spinal cord monitoring during vertebral column stabilization surgery. In Handbook of Spinal Cord Monitoring. Edited by: Boyd J. London: Kluwer; 1992:253–258.
Bracchi F, Guerrieri M, Frassi P, Vigano P, Robinson MA: H-Reflex intraoperative spinal cord monitoring. In Neurological Complications of Spinal Surgery Edited by: G.I.C.D. Paris. 1994, 113–119.
Pelosi L, Lamb J, Grevitt M, Mehdian SM, Web JK, Blumhardt LD: Combined monitoring of motor and somatosensory evoked potential in orthopaedic spinal surgery. Clin Neurophysiol 2002, 113: 1028–1091. 10.1016/S1388-2457(02)00027-5
Slimp JC: Electrophysiologic intraoperative monitoring for spine procedures. Phys Med Rehabil Clin N Am 2004, 15: 85–105.
Balzer JR, Rose RD, Welch WC, Sclabassi RJ, Robert JMD: Simultaneous somatosensory evoked potential and electromyographic recording during lumbo-sacral decompression and instrumentation. Neurosurgery 1998, 42: 1318–1324. 10.1097/00006123-199806000-00074
Schieppati M: The Hoffmann Reflex: a means of assessing spinal reflex excitability and its descending control in man. Prog Neurobiol 1987, 28: 345–376. 10.1016/0301-0082(87)90007-4
Fung KSM, Chan FHY, Lam FK, Poon PWF: A tracing evoked potential estimator. Med Biol Eng Comput 1999, 37: 218–227. 10.1007/BF02513290
Cerutti S, Baselli G, Liberati D, Pavesi G: Single sweep analysis of visual evoked potentials through a model of parametric identification. Biol Cybern 1987, 56: 111–120. 10.1007/BF00317986
Laguna P, Meste O, Poon PW, Caminal P, Rix H, Thakor NT: Adaptive filter for event-related bioelectric signals using an impulse correlated reference input: comparison with signal averaging technique. IEEE Trans Biomed Eng 1992, 39: 1032–1044. 10.1109/10.161335
Walter DO: A posteriori "Wiener filtering" of average evoked response. Electroencephalography & Clinical Neurophysiology 1968, (Suppl 27):61–70.
Heinze HJ, Kunkel H: ARMA-filtering of evoked potentials. Methods Inf Med 1984, 33: 29–36.
Spreckelsen MV, Bromm B: Estimation of single-evoked cerebral potentials by means of parametric modeling and Kalman filtering. IEEE Trans Biomed Eng 1988, 35: 691–700. 10.1109/10.7270
Mainardi LT: Single sweep analysis of event related auditory potential for the monitoring of sedation in cardiac surgery patients. Comput Methods Programs Biomed 2000, 63: 219–227. 10.1016/S0169-2607(00)00112-7
Magni R, Giunti S, Bianchi A, Reni G, Bandello F, Durante A, Cerutti S, Brancato R: Single-sweep analysis using an autoregressive with exogenous input (ARX) model. Doc Ophthalmol 1994, 86: 95–104. 10.1007/BF01224631
Liberati D, Dicorrado S, Mandelli S: Topographic mapping of single sweep evoked potentials in the brain. IEEE Trans Biomed Eng 1992, 39: 943–51. 10.1109/10.256428
Cerutti S, Liberati D, Mascellani P: Autoregressive modelling and filtering of EEG signal generating mechanism in patients under surgical interventions. Model Simulat Contr 1985, 2: 11–23.
Haykin S: Adaptive filter theory. Prentice Hall; 2002.
Bracchi F, Perale G, Rossi L, Gaggiani A, Bianchi A: A PC-Based System for H-Reflex and Single Sweep SEP Coupled Monitoring of Spinal Cord Function in Vertebral Column Surgery. Proceedings of the 25th Annual International Conference of IEEE Engineering in Medicine and Biology Society: 17–21 September 2003; Cancun
Lam BSC, Hu Y, Lu WW, Luk KDK, Chang CQ, Qui W, Chan FHY: Multi-adaptive filtering technique for surface somatosensory evoked potentials processing. Med Eng Phys 2007, 27: 257–266. 10.1016/j.medengphy.2004.09.007
Qui W, Fung KS, Chan FHY, Lam FK, Poon PW, Hamernik RP: Adaptive filtering of evoked potential with radial-basis-function neural network prefilter. IEEE Trans Biomed Eng 2002, 49: 255–232.
FB proposed the neuromonitoring technique and was in charge of the neuromonitoring in the operating theater. SC proposed the signal processing method. AMB and AM performed the signal processing and statistical analysis. LR worked on the simulation of the off-line data processing. AG helped to record data in surgery sessions.
All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.