1 / 19

Evoked Potentials Noise Reduction by Ensemble Averaging

This chapter explores the techniques and methods for noise reduction in evoked potentials through ensemble averaging. Topics covered include spike artifacts and robust averaging, estimation of latency shifts, and weighting of averaged EPs using ensemble correlation.

acary
Download Presentation

Evoked Potentials Noise Reduction by Ensemble Averaging

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. EVOKED POTENTIALSNOISE REDUCTION BY ENSEMBLE AVERAGINGChapters 4.3.4 (case 2) - 4.3.8Seppo Mattila (BRU)

  2. Overview • Averaging of Inhomogenious Ensembles • Spike Artifacts and Robust Averaging • The effect of Latency Shifts • Estimation of Latency Shifts • Weighting of Averaged EPs Using Ensemble Correlation

  3. Averaging of Inhomogenious Ensembles • Varying noise variance (case 1) • Varying signal amplitude but constant noise variance from potential to potential (case 2)

  4. Varying signal amplitude assume that signal amplitude differs from potential to potential: optimal weights from the eigenvalue problem: where all eigenvalues equal to zero, except: and optimla weight vector proportional to corresponding eigenvector: normalise: wi = ai

  5. Varying signal amplitude II ensemble weighted

  6. Gaussian noise with varying variance Weights from maximum likelyhood estimation. The joint PDF of the potentials xi(n) at time n: We maximise its logarithm: by setting its derivative wrt s(n) to zero:

  7. Gaussian noise with varying variance II weighted average of xi(n) i.e. each potential weighted by identical to the result from SNR maximisation moreover, (ensemble average)

  8. Spike artifacts & robust averaging • ensemble & exponential averaging perform well when Gaussian noise • spike (outlier) artifacts degrade performance • need more robust methods: • ensemble averaging with outlier rejection • recursive, robust averaging with outlier rejection

  9. Ensemble averaging with outlier rejection Consider the generalised Gaussian PDF: where the Gamma function The Laplacian PDF (v = 1) for the noise sample:

  10. Ensemble averaging with outlier rejection II ML estimate from again, by setting its derivative wrt s(n) to zero choose s(n) such that exactly half of the sample values greater and half smaller Gaussian noise Laplacian noise ML estimator of s(n) MEDIAN when Laplacian noise

  11. Trimmed means • ensemble average and ensemble median special cases of where K is the largest integer less than or equal to vM v = 0 for enemble average v = 0.5 for ensemble median

  12. Recursive, robust averaging with outlier rejection closely related to exponential average but has an updated part modified by the influence function influence functions sign limiter hard limiter

  13. Effect of latency shifts • variations in latency distortion in ensemble average • shifts in continuous-time signals • caused by biological mechanisms • not constrained to sampling time grid • shifts in discrete-time signals • variations taking place in sampled signal ensemble average

  14. Shifts in continuous time signals for the expected value of ensemble average: The convolution integral can be expressed as a product in the frequency domain: Zero-mean Gaussian PDF is an example of the characteristic function:

  15. Shifts in continuous time signals Latency shifts can act as a lowpass filter on s(t) -3 dB cut-off frequency associated with the low-pass filtering effect due to latency shifts Gaussian sigma sampling interval

  16. Estimation of latency shifts • need to find the shift in each individual potential • compute latency corrected ensamble average • Woody method most well-known • estimates individual shifts • estimates latency corrected ensemble average • interative procedure for improving the estimates

  17. Woody method EP affected by latency shift: PDF of observed signal: ML estimate from: best cross-correlation between s(n) and xi(n) update iteratively the ensemble average:

  18. Weighting of averaged EPs using Ensemble correlation Weight the individual samples such that the difference between s(n) and the average minimised Sample-by-sample weighted ensemble average:

More Related