# A MISO-ARX-Based Method for Single-Trial Evoked Potential Extraction.

1. IntroductionEvoked potentials (EPs) are localized potential changes generated by the central nervous system when stimulated by well-defined external stimuli (such as electrical, light, sound, and other stimuli) [1]. Thus, EPs can be categorized as auditory evoked potentials (AEPs), visual evoked potentials (VEPs), somatosensory evoked potentials (SEPs), and motor evoked potentials (MEPs) according to the modality of stimulation. Depending on the experimental paradigm, EPs may include a complex of partially overlapping components [2], reflecting different processing stages along the neural pathways. The latency variations of specific components can objectively reflect changes in the underlying state of the neural pathways, which is very meaningful in cognitive science research and clinical applications [3]. Many single-trial EP extracting methods have been proposed in order to enhance the ability to track latency variations.

Parametric modeling using autoregression with exogenous inputs (ARX) is a commonly used method for extracting single-trial EPs over the conventional moving time average [4]. ARX modeling for single-trial EP estimation was first proposed by Cerutti et al. [5]. In ARX, the electroencephalogram (EEG) can be viewed as an autoregression (AR) model driven by white noise, and the EP can be accurately modeled by an autoregressive-moving-average (ARMA) filter with a known signal [6]. The known signal is typically the average of the reference EPs (AREP). The order and parameters of the AR and ARMA models can be estimated by utilizing various optimization techniques, such as the final prediction error (FPE) [5] and the least-squares (LS) method [7]. The EPs can then be reconstructed by ARMA filtering with the AREP. ARX modeling has been widely adopted by researchers to rapidly extract middle latency AEPs, VEPs, and SEPs. For example, Mainardi et al. [8] used the ARX model to quantify changes in auditory N100 for the monitoring of sedation in cardiac surgery patients. Rossi et al. [9] extracted single-trial SEPs with ARX filtering for monitoring the functional integrity of the spinal cord during surgery. Lange and Inbar [10] further extended the ARX estimator to make the single-trial estimation process resistant to noise present in the system using a robust evoked potential estimator (REPE). However, Cerutti et al. [5] recently found, by systemic experimentation, that EP extraction using ARX modeling is completely invalid when latency varies greatly compared with the AREP. We carried out a further study on the single-trial extracting experiments made by De Silva et al. [11]. We found that they always assumed temporal lag between the input and the output of the ARX model equaled zero, which causes significant error when latency varies greatly. In addition, they limited the ARX method to yield a waveform similar to the average response, differing only in global latency. Thus, their procedure cannot demonstrate the method's tracking capabilities of specific components of the EPs. An EP complex may contain components that originated from different functional brain sites [12]. The summation of these components results in component overlap, which may cause partial occlusion of the desired component's features. Because of this, the tracking of latency variations of specific components is very difficult.

In this paper, we present a novel single-trial evoked potential estimation method based on multiple-input single-output ARX (MISO-ARX). In MISO-ARX, each component of the EP is individually modeled by an ARMA filter with a reference signal to avoid different components interfering with each other as in ARX. In addition, all parameters are calculated synchronously to guarantee that the estimated EP is optimal overall. Moreover, as EPs have been proven (in our previous paper) to have strong sparsity over an appropriate dictionary, we first roughly estimate the temporal lag of specific components with sparse coding before calculating the parameters of MISO-ARX in order to improve robustness against great latency variations. A series of experiments carried out on simulated and human test responses confirmed the superior performance of our MISO-ARX method for tracking latency variations even in situations of extremely low SNR. The rest of this paper is organized as follows. Section 2 gives a detailed description of our single-trial estimation algorithm. Section 3 contains our experimental results obtained by using the MISO-ARX method and a comparison with ARX and REPE methods. Section 4 presents our conclusions.

2. Single-Trial Evoked Potential Extraction with MISO-ARX

EPs are always embedded in the ongoing spontaneous EEG background, and the SNR is extremely low (below 0dB). The main parts of our method consist of removing the EEG e(t) from the measurement y(t) and then reconstructing the single-trial EP s(t) [13]. The measurement y(t) is

y(t) = s(t)+e(t). (1)

2.1. The EP Signal. The single-trial EP is considered as a complex containing many components; these components may originate from different functional brain sites and can be distinguished according to their respective latencies and amplitudes. The EP waveform s(t) is assumed to be a superposition of Q components:

s(t) = [[Q.summation over (q=1)][k.sub.i][v.sub.q] (t-[[tau].sub.q]), (2)

where [v.sub.q](t) is the basic shape of the qth component, [[tau].sub.q] is the component's latency, and [k.sub.q] indicates the component's amplitude.

In MISO-ARX, each component [v.sub.q] is derived by filtering the reference [u.sub.q](t) using the ARMA model parameters, as shown in

[v.sub.q](t) = [B.sub.q] (z-1)/[A.sub.q]([z.sup.-1][u.sub.q], (3)

where [mathematical expression not reproducible]. Thus, (2) can be rewritten as

[mathematical expression not reproducible]. (4)

[mathematical expression not reproducible] can be parameterized as

[mathematical expression not reproducible]. (5)

Then, we assume [[bar.b].sup.b.sub.q], = [b.sup.q.sub.j][k.sub.q] and [[bar.d].sup.q] = [d.sup.q] + [[tau].sub.q], which yields

[mathematical expression not reproducible]. (6)

Thus, [mathematical expression not reproducible], and we obtain

[mathematical expression not reproducible]. (7)

2.2. Reference Signal for Each Component. The selection of reference signals, which directly affects the accuracy of EP extraction, is a very important process. In MISO-ARX, we extract the reference signal [u.sub.q](t) for each component from the AREP using a specific filtering window function, such as the Hamming window or the Blackman window. The central location and width of the window are determined by the location and width of the peak, respectively. An example of a simulated decomposition is provided in Figure 1. It can be seen that the simulated waveform consists of three underlying components, which were extracted using a Hamming window.

2.3. The EEG Signal. In this paper, the EEG signal e(t) is viewed as an AR model driven by white noise w(t), defined as

e(t) = 1/[A.sub.e]([z.sup.-1]) w(t), (8)

where [mathematical expression not reproducible].

2.4. Estimation of Temporal Lag. In ARX, most researchers assume that temporal lag d between the input and output of the model is equal to zero before estimating the model orders m and n. This assumption is not in accordance with practice. In this study, we used sparse coding to roughly estimate the value of d. Sparse coding has had significant success in signal denoising and separation. In addition, in our previous paper, EPs were proven to have strong sparsity over an appropriate dictionary [14]. Assuming D and [theta] are the dictionary and the sparse coefficients, respectively, s(t) can be expressed as s(t) = D[theta] [15]. Thus, y(t) is

y(t) = s (t) + e(t) = D[theta] + e (t). (9)

The estimator for [theta] is calculated by solving

[mathematical expression not reproducible], (10)

where [[epsilon].sub.0] is determined by the variance of the EEG. Equation (10) can be solved by using optimization methods, such as basis pursuit [16], orthonormal matching pursuit [17], and Lasso [18]. Since the atoms of D are constructed by left or right translation of the basic components of EPs, we use the location of nonzero values in d to estimate the temporal lag d.

2.5. Single-Trial Extraction. By replacing (7) and (8) in (1), we get

[mathematical expression not reproducible]. (11)

In order to simplify this model, we assume [A.sub.q]([z.sup.-1]) = [A.sub.e]([z.sup.-1]) = A([z.sup.-1]). This implies a partial loss of generality with respect to the AR model for the noise and the ARMA model for the signal, completely independent of each other. Nevertheless, the MISO-ARX model requires more complex algorithms to be characterized. In this paper, this MISO system is characterized with the global separable nonlinear multi-innovation recursive least-squares-identification method [19]. Then, parameters A([z.sup.-1]) and [[bar.B].sub.q]B([z.sup.-1]) are estimated. The EP can be reconstructed as

[mathematical expression not reproducible]. (12)

3. Experimental Results

3.1. Simulation Experiment. A computer simulation was conducted to verify the performance of our MISO-ARX method for EP signal estimation. Three single-trial-estimation methods for EP signals, namely, ARX, REPE, and SP, were compared in the following simulations. In order to measure the performance of each of the methods, we measured the SNR of the estimated EPs and the accuracy of the latencies for different SNRs to evaluate the quality of these methods. All experiments were implemented in Matlab R2006b on a Pentium 2.7 GHz PC with 4 GB RAM.

The reference signals [v.sub.i](k) were simulated by the superimposition of three basic components, which can be represent by the Gaussian distribution function [20]; thus

u (t) = 0.25 sinc [0.13[pi] (4t - 8)]

+ 0.5 sinc [0.13[pi](4t-16)] (13)

+ sinc [0.13[pi](4t-24 + l)].

The synthetic reference signals (l = 0 and l = 3) are shown in Figure 2.

From Figure 2, it can be seen that all three components have changed and the latency of the third wave varies significantly. In this study, we used an EP with l = 0 as the reference signal. The background EEG superimposed on the EP signal was simulated by an autoregressive process [21], as shown in the following equation:

q(t) = 1.5084q(t- 1) - 0.1587q(t -2) - 0.3109q (t -3) - 0.0510q (t-4) + w (t), (14)

where w(t) is Gaussian white noise. During the process of estimation, the SNR of the observations may change over time due to the nonstationary characteristics of the EEG. Therefore, in this experiment, the performance of the four different methods was examined under various SNR conditions. The SNRs of the observations were changed from 0 dB to -10 dB, and I was changed from -5 to 5. For each SNR value, 100 pairs of observations were generated. The average results for different SNR and I values in 100 independent runs are shown in Figure 3.

It is clear that, with the decrease of the value of the SNR, estimation performance declines. However, in MISOARX, the changes in I had hardly any impact on the estimation results. This illustrates that our method is apt for tracking latency variation of EPs. Since the SNR is defined by the complete signal and not by a specific feature of it, we measured the MISO-ARX method's ability of tracking latency variations. We compared ARX and MISO-ARX in the case of three different SNR values (0 dB, -5 dB, and -10 dB) and four different latencies (4.17 ms (l = 7), 5.29 ms (l = 3), 6.26 ms (l = -1), and 7.24 ms (l = -5)). Results are shown in Figure 4. We can see from Figures 4(a) and 4(b) that, for high SNR values (0 dB and -5 dB), our method had strong tracking capability for all latencies. Figure 4(c) indicates a deterioration of latency tracking performance when the SNR decreased, but the MISO-ARX method's accuracy rate still exceeded 70%. With ARX, even for high SNR values, for l = -5 and I = 7 the estimations were totally wrong, suggesting that tracking such variations is not possible.

3.2. Real Data. For further evaluation of the performance of our method, VEPs were collected from six eyes belonging to three human subjects during pattern reversal VEP experiments. This study was conducted with the approval of the local ethics committee, and all experiments with human participants were performed according to the National Institutes of Health Guidelines. VEP signals were recorded via a surface electrode placed on the occipital region of the scalp (Oz), using the right earlobe as a reference, with the forehead grounded. Subjects were required to gaze at a cross on the stimulus screen. The stimulus pattern was a conventional black-and-white checkerboard, which was reversed every half a second. Recordings were made using a digital EEG recording system (NuAmps EEG Amplifier, NeuroScan, USA) with a sampling rate of 1000 Hz and stored on a computer. Signals were bandpass filtered in the range of 0.05-450 Hz.

In the VEP signal, critical responses were located at approximately 100 ms after stimulation, where the positive peak (P100) occurred. With the stimulus being delivered at 0 ms, we processed each VEP trial from 0 ms to 200 ms. We used data from the first 50 pairs of trials to perform the experiment, from six eyes belonging to three subjects, as shown in Figure 5. As shown in the figure, it is clear that most estimated VEPs have a peak at around 100 ms.

For each eye, we used data from the first 50 pairs of trials to calculate the average of the estimated VEPs. The results are shown in Figure 6. The estimates using MISOARX are indicated by solid line, and the averages of the measurements are indicated by dashed-dotted line. Clearly, as shown in Figure 6, the average of the estimated results for each eye is very similar to the average of the measurements.

Then, we estimated the latency of the P100 for each trial, as shown in Figure 7. The estimates are indicated by dots, the averages of the estimates are indicated by solid line, and the P100 latencies of the average of the measurements are indicated by dashed-dotted line. We can see that the trial-to-trial variation in the latencies of P100 is large. However, the averages of the estimates are close to the results of the average of the measurements.

4. Conclusion

In this paper, we presented a novel single-trial EP extraction method based on MISO-ARX. This method considers the single-trial EP as a complex containing many components, and each component can be modeled by ARMA. In addition, in order to improve the accuracy of the model parameters estimation, we used sparse coding to roughly estimate temporal lag. Since each component is modeled individually, our method has greater tracking capabilities of specific components of the EP complex. We conducted a series of experiments on synthetic and real data, and the results were evaluated using waveform observations and several metrics. From point of view of the experimental results, our method achieved a better and more favorable estimation performance than other currently used state-of-the-art methods in single-trial EP estimations.

http://dx.doi.org/10.1155/2017/7395385

Competing Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work; there are no professional or other personal interests of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, this manuscript. And the funding in "Acknowledgments" would not lead to any conflict of interests regarding the publication of this manuscript.

Acknowledgments

This work was supported by the Nature Science Foundation of China (Grant no. 61401181) and Nature Science Foundation of China (Grant no. 61403174).

References

[1] Q. Wang, Y. Wu, W. Liu, and L. Gao, "Dominant eye and visual evoked potential of patients with myopic anisometropia," BioMed Research International, vol. 2016, Article ID 5064892, 6 pages, 2016.

[2] A. S. Gevins, "Analysis of the electromagnetic signals of the human brain: milestones, obstacles, and goals," IEEE Transactions on Biomedical Engineering, vol. 31, no. 12, pp. 833-850, 1984.

[3] M. H. Costa, "Estimation of the noise autocorrelation function in auditory evoked potential applications," Biomedical Signal Processing & Control, vol. 7, no. 5, pp. 542-548, 2012.

[4] S. F. Regacone, D. D. B. D. Lima, V. E. Valenti, and A. C. F. Frizzo, "Resting heart rate and auditory evoked potential," BioMed Research International, vol. 2015, Article ID 847506, 6 pages, 2015.

[5] S. Cerutti, G. Baselli, D. Liberati, and G. Pavesi, "Single sweep analysis of visual evoked potentials through a model of parametric identification," Biological Cybernetics, vol. 56, no. 2, pp. 111-120, 1987.

[6] B. Berger, T. Minarik, G. Liuzzi, F. C. Hummel, and P. Sauseng, "EEG oscillatory phase-dependent markers of corticospinal excitability in the resting brain," BioMed Research International, vol. 2014, Article ID 936096, 8 pages, 2014.

[7] H. Kumru, D. Soler, J. Vidal, J. M. Tormos, A. Pascual-Leone, and J. Valls-Sole, "Evoked potentials and quantitative thermal testing in spinal cord injury patients with chronic neuropathic pain," Clinical Neurophysiology, vol. 123, no. 3, pp. 598-604, 2012.

[8] L. T. Mainardi, J. Kupila, K. Nieminen et al., "Single sweep analysis of event related auditory potentials for the monitoring of sedation in cardiac surgery patients," Computer Methods & Programs in Biomedicine, vol. 63, no. 3, pp. 219-227, 2000.

[9] L. Rossi, A. M. Bianchi, A. Merzagora, A. Gaggiani, S. Cerutti, and F. Bracchi, "Single trial somatosensory evoked potential extraction with ARX filtering for a combined spinal cord intraoperative neuromonitoring technique," BioMedical Engineering Online, vol. 6, article no. 2, 2007.

[10] D. H. Lange and G. F. Inbar, "A robust parametric estimator for single-trial movement related brain potentials," IEEE Transactions on Biomedical Engineering, vol. 43, no. 4, pp. 341-347,1996.

[11] A. C. De Silva, N. C. Sinclair, and D. T. J. Liley, "Limitations in the rapid extraction of evoked potentials using parametric modeling," IEEE Transactions on Biomedical Engineering, vol. 59, no. 5, pp. 1462-1471, 2012.

[12] D. H. Lange, H. Pratt, and G. F. Inbar, "Modeling and estimation of single evoked brain potential components," IEEE Transactions on Biomedical Engineering, vol. 44, no. 9, pp. 791-799,1997.

[13] J. L. K. Kramer, J. Haefeli, A. Curt, and J. D. Steeves, "Increased baseline temperature improves the acquisition of contact heat evoked potentials after spinal cord injury," Clinical Neurophysiology, vol. 123, no. 3, pp. 582-589, 2012.

[14] N. Yu, H. Liu, X. Wang, and H. Lu, "A joint sparse representation-based method for double-trial evoked potentials estimation," Computers in Biology & Medicine, vol. 43, no. 12, pp. 2071-2078, 2013.

[15] N. Yu, F. Hu, D. Zou, Q. Ding, and H. Lu, "Single-trial sparse representation-based approach for VEP extraction," BioMed Research International, vol. 2016, Article ID 8569129, 9 pages, 2016.

[16] O. Shental, "Sparse representation of white Gaussian noise with application to L0-norm decoding in noisy compressed sensing," in Proceedings of the Qualcomm Technology Forum (QTECH '11), San Diego, Calif, USA, June 2011.

[17] S. S. Chen, D. L. Donoho, and M. A. Saunders, "Atomic decomposition by basis pursuit," SIAM Review, vol. 43, no. 1, pp. 129-159, 2001.

[18] Y. C. Pati, R. Rezaiifar, and P. S. Krishnaprasad, "Orthogonal matching pursuit: recursive function approximation with applications to wavelet decomposition," in Proceedings of the 27th Asilomar Conference on Signals, Systems and Computers, vol. 1, pp. 40-44, Pacific Grove, Calif, USA, November 1993.

[19] J. H. Wang, D. B. Wang, and Z. S. Wang, "Recursive identification of MISO systems with multiple unknown time delays," Control and Decision, vol. 25, no. 1, pp. 93-98, 2010.

[20] D. L. Jewett and J. S. Williston, "Auditory-evoked far fields averaged from the scalp of humans," Brain, vol. 94, no. 4, pp. 681-696, 1971.

[21] N. Kamel, M. Z. Yusoff, and A. F. M. Hani, "Single-trial subspace-based approach for VEP extraction," IEEE Transactions on Biomedical Engineering, vol. 58, no. 5, pp. 1383-1393, 2011.

Nannan Yu, (1) Lingling Wu, (1) Dexuan Zou, (1) Ying Chen, (1) and Hanbing Lu (2)

(1) School of Electrical Engineering and Automation, Jiangsu Normal University, Xuzhou 221116, China

(2) Department of Internal Neurology, Xuzhou Central Hospital, Xuzhou 221116, China

Correspondence should be addressed to Hanbing Lu; luhanbing111@126.com

Received 15 November 2016; Accepted 9 January 2017; Published 8 February 2017

Academic Editor: Victor H. C. de Albuquerque

Caption: FIGURE 1: The reference signal extracted from the AREP.

Caption: FIGURE 2: The waves of the reference signals in 1 = 0 and 1 = 3.

Caption: FIGURE 3: Performance evaluation of four methods with different SNR and I values.

Caption: FIGURE 4: Latency tracking of the third wave using ARX and MISO-ARX.

Caption: Figure 5: Estimation performance.

Caption: FIGURE 6: Average of the estimated VEPs and the average of the measurements.

Caption: FIGURE 7: Estimation of latencies of P100 by MISO-ARX.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Yu, Nannan; Wu, Lingling; Zou, Dexuan; Chen, Ying; Lu, Hanbing |

Publication: | BioMed Research International |

Date: | Jan 1, 2017 |

Words: | 3444 |

Previous Article: | GosB Inhibits Triacylglycerol Synthesis and Promotes Cell Survival in Mouse Mammary Epithelial Cells. |

Next Article: | Apoptotic Signaling Pathways in Glioblastoma and Therapeutic Implications. |

Topics: |