Deconvolution of the vestibular evoked myogenic potential using the power spectrum of the electromyogram
- Bernd Lütkenhöner^{1}Email author
https://doi.org/10.1186/s12976-015-0018-x
© Lütkenhöner. 2015
Received: 16 April 2015
Accepted: 21 September 2015
Published: 6 October 2015
Abstract
Background
The vestibular evoked myogenic potential (VEMP) can be modelled reasonably well by convolving two functions: one representing an average motor unit action potential (MUAP), the other representing the temporal modulation of the MUAP rate (rate modulation). It is the latter which contains the information of interest, and so it would be desirable to be able to estimate this function from a combination of the VEMP with some other data. As the VEMP is simply a stimulus-triggered average of the electromyogram (EMG), a supplementary, easily accessible source of information is the EMG power spectrum, which can be shown to be roughly proportional to the squared modulus of the Fourier transform of the MUAP. But no phase information is available for the MUAP so that a straightforward deconvolution is not possible.
Methods
To get around the problem of incomplete information, the rate modulation is described by a thoughtfully chosen function with just a few adjustable parameters. The convolution model is then used to make predictions as to the energy spectral density of the VEMP, and the parameters are optimized using a cost function that quantifies the difference between model prediction and data.
Results
The workability of the proposed approach is demonstrated by analysing Monte Carlo simulated data and exemplary data from patients who underwent VEMP testing as part of a clinical evaluation of their dizziness symptoms.
Conclusions
The approach is suited, for example, to estimate the duration of the inhibition causing the VEMP or to disentangle a VEMP consisting of more than one component.
Keywords
EMG MUAP VEMPIntroduction
The otolith organs in the inner ear, saccule and utricle, are sensors for linear accelerations and, as such, are important for controlling posture and eye movement. To assess their function in clinical settings it can be exploited that they also respond to brief high-intensity sounds. Muscle reflexes elicited that way give rise to short-latency myogenic responses, which can be recorded using surface electrodes placed over the respective muscles. This type of signal is called vestibular evoked myogenic potential (VEMP). Two varieties of VEMP are to be distinguished: The cVEMP represents a vestibulo-collic reflex and is recorded from a cervical muscle such as the sternocleidomastoid, whereas the oVEMP represents a vestibulo-ocular reflex and arises from the extraocular muscles [1, 2]. The two signals provide complementary diagnostic information: While the sound-evoked cVEMP is assumed to originate primarily in the ipsilateral saccule, the oVEMP is attributed mainly to the contralateral utricle [3]. Only the cVEMP is considered in this article, but the theoretical framework developed is applicable to the oVEMP as well.
where V, R, and H are the Fourier transforms of v, r, and h, respectively. The quantity of interest is r(t), and therefore the question arises as to how to estimate this quantity from the data. A method presented in a previous article [6] exploited the fact that r(t) effects not only the mean of the EMG, i.e., the VEMP, but causes also a modulation of the variance. Unfortunately, this VEMP-associated variance modulation is not generally workable in standard VEMP investigations, because it typically has a much lower signal to noise ratio than the VEMP itself (as yet, the method was applied only to grand-averaged data). In the present article we therefore exploit a different source of information which promises to be much more robust: the power spectrum of the EMG. A simulation conducted by Wit and Kingma [4] suggested that the overall shape of this spectrum corresponds to the energy density spectrum of the MUAP. Thus, by implication, the power spectrum of the electromyogram allows us to estimate |H(f)|, except for a factor that is essentially independent of frequency. Moreover, from Eq. (2) it follows that dividing |V(f)| by this estimate provides an estimate of |R(f)|, at least conceptually.
The problem remains that this procedure provides no phase information. A conceivable solution is to resort to higher-order spectra [7]. However, caution is advised because such spectra are highly susceptible to noise [8]. One might also consider reconstructing H(f) from its modulus by making certain assumptions about the phase. It could be assumed, for example, that H(f) represents, in good approximation, a minimum phase filter. However, such assumptions would not be well-founded. Therefore yet another approach is proposed here: Assumptions are made about R(f) rather than H(f), which appears to be a much easier task. It will be shown that this approach allows us, indeed, to interpret the VEMP using information from the EMG. Depending on the quality of the data, it becomes possible, for example, to estimate the duration of an inhibitory or excitatory modulation and to get hints as to a possible multi-component structure of the VEMP.
Theory
Power spectral density of the electromyogram
A simulation by Wit and Kingma [4] suggests that |S(f)| gets its overall shape from the factor |H(f)| and its noisy character from the factor |A _{ N }(f)|.
with sinc(x)= sin(x)/x.
is the L ^{2} norm of h(t).
In practice, the power spectrum has to be estimated from a finite sample of the EMG, and so the question arises whether it is justified to neglect the second term on the right-hand side of Eq. (10). As \(\bar {a}\) and \( {\bar {\bar {{a}}}}\) typically have about the same value, the condition to be fulfilled is N/(π f T)^{2}≪1. With N=ρ _{0} T, as already assumed above, the condition can be rewritten as f ^{2}≫ρ _{0}/(π ^{2} T). Making the more or less realistic assumptions ρ _{0}=1000 s ^{−1} (see, e.g., [4]) and T=10 s, the requirement can be approximated as f ^{2}≫10 s ^{−2}. This simple estimation suggests that proportionality between |H(f)|^{2} and power spectrum of the EMG can be expected only at sufficiently high frequencies. In this article we generally choose a lower frequency limit of 10 Hz.
Optimization problem
For symmetry reasons there is no need to consider negative frequencies. A lower frequency limit, f _{min}, is defined to avoid the problem considered at the end of the previous subsection. Defining an upper frequency limit, f _{max}, may be useful to exclude noise outside the frequency range of interest. The normalization overcomes the problem that \(|\hat {V}(f)|\) is arbitrarily scaled. The nature of the cost function can be varied by means of the parameters κ and λ, which are assumed to be integers. Extensive simulations (presented below) suggest that, all in all, the cost function Q _{2,1} is a good default choice.
Modelling the rate modulation
This function has more parameters than can be determined by minimizing the right-hand side of Eq. (16). The reason is that r(t) is in this context indistinguishable from c·r(t−t _{0}), assuming c≠0. We fix this problem by setting r _{1}=−1 and t _{1}=0. Thus, in the simplest case (K=1) there is only a single parameter to be optimized: τ _{1}.
The three rate modulation functions presented in Fig. 1 a agree regarding the first peak, for which we assumed r _{1}=−1, t _{1}=0, and τ _{1}=2 ms (the apparent minor differences between the curves go back to the normalization). While the grey curve shows only this one peak, the other two curves show a second peak of opposite polarity, centered at t _{2}=20 ms. In the case of the solid black curve we assumed τ _{2}=τ _{1}, whereas in the case of the dashed curve we assumed τ _{2}=2τ _{1}. Moreover, as to the peak amplitude we assumed r _{2} τ _{2}=0.8 ms, which ensures that the area under the second peak is the same for the two cases.
The exponential term clearly explains the bell shape of the grey curve, whereas the cosine term explains the oscillations of the other two curves. Equation (21) predicts a spectral minimum at f=1/(t _{2}−t _{1}), which is 50 Hz in the present example. This minimum is a potentially important fingerprint when analysing and interpreting real data. Hence, it is noteworthy that the minimum is found also in the dashed curve, although the assumption τ _{1}=τ _{2}, on which Eq. (21) is based, is not fulfilled.
Figure 1 c illustrates how the magnitude of the Fourier transform depends on the parameter τ, which was varied between 1 and 8 ms. All curves refer to a single-peaked rate modulation; the grey curve (τ=2 ms) is identical with the grey curve in Fig. 1 b. A comparison of the curves for τ=2 ms and τ=4 ms gives a hint as to why the oscillation of the dashed curve in Fig. 1 b fades more rapidly with increasing frequency than does the oscillation of the solid black curve: The broader the second peak is in the time domain, the less does it contribute to higher frequencies.
In Fig. 1 d, the parameter r _{2} was systematically varied between -1 and 1, whereas τ _{2} was kept constant at 2 ms. The thick grey curve, representing the case r _{2}=0, is again identical with the bell-shaped grey curve in Fig. 1 b, whereas the black and the thin grey curves were obtained for r _{2}>0 and r _{2}<0, respectively. As to be expected in view of Eq. (21), the parameter r _{2} determines magnitude and polarity of the deviation from the bell-shaped “baseline”.
VEMP deconvolution
where R _{max} is the maximum value of |R(f)| and ε≪1 is a positive real number. The first factor has basically no effect for \(|R(f)|^{2}\gg \epsilon R_{\max }^{2}\), whereas for \(|R(f)|^{2}\ll \epsilon R_{\max }^{2}\) it causes \(|\hat {H}(f)|\) to have a value close to zero. Lütkenhöner and Basel [6] showed that this kind of regularization can be understood as optimal (Wiener) filtering. In this article we defined ε as 10^{−2}.
Non-uniqueness
The setup of the optimization procedure already accounted for the fact that the rate modulation can only be determined up to an unknown scale factor and that any time shift of the rate modulation can be compensated for by a reverse time shift of the MUAP. But there is yet another kind of non-uniqueness inherent to the inverse problem considered here: Time reversal of the rate modulation function, which in the frequency domain corresponds to substituting R(f) by its complex conjugate, \({R}^{*}(f)=|R(f)|e^{-\mathfrak {i}\arg (R(f))}\), has no effect on the cost function defined in Eq. (16). Using R ^{∗}(f) rather than R(f) in Eq. (22) yields \(\hat {H}_{{R}^{*}}(f)=\hat {H}(f)e^{2\mathfrak {i}\arg (R(f))}\). The fact that \(\hat {H}(f)\) and \(\hat {H}_{{R}^{*}}(f)\) usually differ fundamentally could help to resolve the ambiguity related to time reversal, provided that a rough idea about the phase of H(f) is available. For example, it might be reasonable to assume that H(f) more likely resembles a minimum rather than a maximum phase filter.
Simulation results
Simple Gaussian model
In the present example, τ and θ have the values 3 and 4 ms, respectively, so that θ _{ v } is 5 ms. The extrema of h(t) and v(t) are found at the times ±θ and ±θ _{ v }, respectively (dotted vertical lines).
Calculating the modulus of the Fourier transform yields the curves shown on the right of Fig. 2. In view of the above time-domain considerations it is not surprising that the curves for MUAP (black) and VEMP (gray) are similar. The differences between these curves can be exploited to draw conclusions as to the rate modulation. This is most easily done by comparing the locations of the spectral maxima, which, for the situation considered here, are found at the frequencies 1/(2π θ) and 1/(2π θ _{ v }), respectively [10]. Numerical evaluation yields 39.8 and 31.8 Hz, respectively (dotted vertical lines). By reversing the line of thought, the frequencies of experimentally determined spectral maxima of MUAP and VEMP can be converted into estimates of θ and θ _{ v }, and τ can subsequently be calculated using Eq. (23).
Estimation of the parameter τ of the simple Gaussian model
Q _{1,1} | Q _{1,2} | Q _{1,4} | Q _{2,1} | Q _{2,2} | Q _{2,4} | |
---|---|---|---|---|---|---|
Median | 2.47 | 2.61 | 2.70 | 2.80 | 2.70 | 2.61 |
Mean | 2.49 | 2.65 | 2.71 | 2.83 | 2.74 | 2.62 |
SD | 0.62 | 0.59 | 0.68 | 0.72 | 0.83 | 0.97 |
RMSD | 0.80 | 0.69 | 0.74 | 0.74 | 0.87 | 1.04 |
More sophisticated model
The grey curves in Fig. 3 visualize theoretical functions which are, of course, not available when dealing with real data. Instead, the analysis has to begin with two curves that are generally quite noisy: the estimated VEMP and the estimated power spectral density of the EMG. Here we consider the square root of the latter, which serves as a surrogate for the modulus of the Fourier transform of the MUAP. The examples shown as black curves in Fig. 3 c and 3 e, respectively, were obtained by Monte Carlo simulation. As in the simulations before, the MUAP rate ρ _{0} was set to 1000 s ^{−1}, the data was sampled at a rate of 2000 s ^{−1}, and stimuli were presented at a rate of 4 s ^{−1}. The number of stimulus presentations was raised to 200, though. In accordance with the Theory section of the article, where it was explained that the EMG power spectrum is basically proportional to the energy spectral density of the MUAP, the grey curve in the background of Fig. 3e appears, indeed, as a smoothed version of the black curve.
What remains to be done now is to choose a rate modulation so that, in the frequency-domain, multiplication with the estimate obtained for the MUAP (square root of the EMG power) optimally fits the observed VEMP. The frequency-domain representation of the VEMP is shown as a black curve in Fig. 3 f, whereas the best fitting model (the one that minimizes the cost function Q _{2,1}) is represented by the dotted (red) curve. The corresponding rate modulation, represented by the black curve in Fig. 3 a, is in good agreement with the theoretical counterpart (grey curve).
Estimation of the MUAP function h(t) is optional, but can help to check the plausibility of the rate modulation obtained. The dashed curve in Fig. 3 b shows the Fourier transform of \(\hat {H}(f)\) as calculated using Eq. (22). The curve apparently agrees well with the theoretical counterpart (grey curve).
Parameter estimation for the two-component model
Cost function | |||||||
---|---|---|---|---|---|---|---|
Q _{1,1} | Q _{1,2} | Q _{1,4} | Q _{2,1} | Q _{2,2} | Q _{2,4} | ||
τ _{1} | median | 1.77 | 1.82 | 1.85 | 1.92 | 1.81 | 1.32 |
(ms) | mean | 1.78 | 1.84 | 1.86 | 1.94 | 1.86 | 1.59 |
SD | 0.20 | 0.24 | 0.41 | 0.28 | 0.51 | 0.66 | |
RMSD | 0.29 | 0.29 | 0.43 | 0.28 | 0.53 | 0.78 | |
τ _{2} | median | 4.23 | 3.86 | 3.61 | 4.04 | 3.76 | 3.14 |
(ms) | mean | 4.38 | 4.17 | 3.88 | 4.15 | 3.92 | 3.54 |
SD | 1.21 | 1.76 | 2.09 | 1.14 | 2.24 | 2.62 | |
RMSD | 1.27 | 1.77 | 2.09 | 1.15 | 2.24 | 2.66 | |
r _{2} | median | 0.25 | 0.26 | 0.26 | 0.25 | 0.26 | 0.21 |
mean | 0.25 | 0.27 | 0.29 | 0.25 | 0.30 | 0.28 | |
SD | 0.08 | 0.09 | 0.13 | 0.09 | 0.16 | 0.18 | |
RMSD | 0.08 | 0.09 | 0.14 | 0.09 | 0.17 | 0.18 | |
t _{2} | median | 19.87 | 19.90 | 19.86 | 19.98 | 19.88 | 19.85 |
(ms) | mean | 19.82 | 19.80 | 19.76 | 19.99 | 19.82 | 19.83 |
SD | 1.62 | 1.65 | 1.61 | 1.52 | 1.64 | 1.85 | |
RMSD | 1.63 | 1.66 | 1.63 | 1.52 | 1.65 | 1.86 |
Analysis of exemplary real data
Discussion and conclusions
This study started from the conjecture that the differences in the spectral densities of VEMP and EMG represent a signature which, if interpreted appropriately, can tell something about the generation of the VEMP. Model simulations building on previous theoretical work [4, 5, 10] allowed us to make predictions about this signature. Our current understanding is that the VEMP represents a brief reduction of the firing probability of the motor units. To implement this idea in the convolution model given by Eq. (1), the rate modulation is defined as a negative pulse. The width of the pulse can principally be estimated by comparing the locations of the spectral maxima of VEMP and EMG: the more they differ, the broader the peak (Fig. 1 c). But when working with real data such as those presented in Fig. 4, noise typically prevents a precise estimation of spectral maxima. The problem can be overcome by resorting to a parameter optimization approach which accounts for the whole (or nearly the whole) spectrum. This way we derived the rate modulation shown in Fig. 4 a. For a Gaussian function with the standard deviation τ, the full width at half maximum is \(2\sqrt {2\ln 2} \tau \approx 2.355 \tau \). Inserting the estimated τ we get 5.9 ms, which is consistent with Colebatch and Rothwell [12]: In a study of single motor units they found inhibition windows with a duration between 2 and 8 ms (mean 3.6 ms).
A second VEMP component, which follows roughly 20 ms after the first one and has opposite polarity, is found in 55 % [13] to 76 % [14] of healthy subjects. Because this component does not depend on the integrity of the vestibular nerve, it is generally assumed to be of cochlear rather than vestibular origin [11]. But the component was identified also in deaf ears, which could be indicative of a dual origin: cochlear as well as vestibular [13]. No matter what the origin of the second component is, the analysis of a VEMP consisting of two components requires to choose a rate modulation with two peaks. In that case, the spectral density of the VEMP differs from the spectral density of the EMG by a more or less pronounced spectral minimum. This kind of signature is so salient in Fig. 5 f that the evaluation can be done just by visual inspection. The location of the spectral minimum is easy to interpret, because it corresponds to the reciprocal of the latency difference between the two VEMP components. Other frequency-domain features are more difficult to understand. Thus, a parameter optimization has to be performed to explain given data in terms of a model. In the example presented in Fig. 5, the assumption of a two-component rate modulation was indeed suitable to explain the data. It was even possible to reduce the number of parameters by assuming that the two peaks of the rate modulation have the same width.
A key feature of the proposed method is that the rate modulation is characterized in terms of just a few parameters, whereas no assumptions whatsoever are made about the MUAP. The latter is considered as an arbitrary function to be estimated from the data (more precisely, only the modulus of the Fourier transform is relevant). Thus, the finding that the estimated MUAP (Figs. 4 b and 5 b) is roughly biphasic, in good agreement with the theoretical functions used in the model simulations, does not reflect assumptions inherent to the model, but is based on the measured data. The reason why the rate modulation rather than the MUAP was parametrized (in theory, the latter possibility would work equally well) is that parameters closely related to the questions of interest can be chosen. If the VEMP has two components, the questions concern the latency difference as well as the amplitude ratio of the components, and the half width of each component. Inhibitory and excitatory rate modulations are distinguished by the sign, but apart from that they are handled in the same way. Although the method cannot be used to confirm the generally accepted opinion that the first VEMP component (p13-n23) is of an inhibitory nature, it is possible to relate the “polarity” (excitatory versus inhibitory) of any other component to that of the first component.
A drawback of the normalization in Eq. (16) is that the estimated functions are unscaled. In this respect the present approach is clearly inferior to a method developed by Lütkenhöner and Basel [6], where the deconvolution algorithm exploits the relationship between VEMP and associated variance modulation. But that method is applicable only in well-chosen cases, because a variance modulation with a sufficient signal-to-noise ratio is rather the exception. The present approach, by contrast, promises to be workable whenever an acceptable VEMP was recorded, because this normally means that the EMG power spectrum can be estimated reasonably well, too. Considering the fact that both deconvolution approaches have advantages and disadvantages, the question arises as to whether they can be combined. This should indeed be possible. The general idea is that the unscaled functions obtained with the method proposed here are scaled using the previously developed algorithm. Conceptually, given unscaled estimates of rate modulation and MUAP, a simulation as presented in Fig. 3 could be run, provided that suitable assumptions are made about what is unknown. Detailed knowledge of the statistical distribution of the MUAP amplitudes is not essential, because a normalization of the EMG (such that it has unit variance) prior to the estimation of the VEMP basically eliminates the dependence of the VEMP on the MUAP amplitudes [5]. In essence, only two unknown parameters remain: the mean MUAP rate and the scaling factor for the rate modulation. These two parameters could be optimized by comparing the variance modulation estimated from simulated EMG data with the variance modulation derived from the experimental data. Even if the signal-to-noise ratio of the latter is too bad for a stand-alone deconvolution as proposed in our previous article [6], there may be enough information for optimizing only two parameters. This consideration shows that the deconvolution method developed in this study could eventually become a key component of a more comprehensive deconvolution approach that simultaneously works on VEMP, associated variance modulation, and power spectrum of the EMG.
Methods
Monte Carlo experiments
The application of Eq. (3) requires that random numbers a _{ n } and t _{ n } (1≤n≤N) are generated. Strictly speaking, the total number of contributing MUAPs, N, is a random number as well. As to the MUAP occurrence times t _{ n }, the task was accomplished by considering the MUAP generation as a time-dependent Poisson process (see, e.g., [15]). More specifically, the mean number of MUAPs occurring between t and t+Δ t was assumed to be ρ(t)Δ t, where ρ(t) is a rate function; the Matlab function POISSRND converted the mean number into a random number of MUAPs (mostly 0 or 1, provided that Δ t is sufficiently small; the simulations presented here were performed with Δ t=0.5 ms). To generate all the required MUAP occurrence times, the time range of interest (with some extra time on both sides, to avoid edge effects) was sampled at intervals of Δ t, and the idea of a Poisson process was successively applied to each sampling time. The amplitude values associated with the occurrence times were drawn from a gamma distribution with shape parameter 2, as already suggested by Wit and Kingma [4]. For the sake of convenience, the value 1/2 was chosen for the scale parameter of the distribution, which has the consequence that the mean amplitude is \(\bar {a}=1\). Gamma distributed random numbers were generated using the Matlab function GAMRND.
where ρ _{0} denotes, as in the context of Eq. (12), the mean MUAP rate of the undisturbed EMG (no VEMP eliciting stimulus). To simulate a typical VEMP experiment, the stimulus was assumed to be presented at intervals of 250 ms, and a corresponding periodicity was implemented also into the rate ρ(t).
under the conditions considered here [5]. The convolution integral was calculated using the Matlab function CONV.
Exemplary real data
Two exemplary sets of real data were taken from an archive compiled by Lütkenhöner et al. [16], which comprises data from patients who underwent VEMP testing as part of a clinical evaluation of their dizziness symptoms. In brief, VEMPs were elicited by 500-Hz Gaussian tone pulses with a full width at half maximum of 4 ms, which were presented at a rate of 4/s (peak-equivalent sound pressure level of 107 dB). The electromyogram, recorded from the sternocleidomastoid muscle, was continuously digitized at a rate of 10 kHz. Each data set comprises 200 stimulus presentations.
Estimation of the rate modulation from given data
The data required by the algorithm consist of the measured VEMP, \(\bar {v}(t)\), and an estimate of the magnitude of the Fourier transformed MUAP, \(|\bar {H}(f)|\), which is basically the square root of the EMG power spectrum. The latter was estimated by means of the Matlab function PWELCH. Figure 7 illustrates that the algorithm uses the two types of data in a quite different way: \(|\bar {H}(f)|\) is multiplied by |R(f)| to predict the magnitude of the Fourier transformed VEMP, whereas \(\bar {v}(t)\) directly enters the calculation of a figure of merit (after transformation into the frequency domain).The figure of merit quantifies the difference between predicted and measured VEMP using the cost function defined in Eq. (16). This measure eventually determines the next action taken by the nonlinear optimization routine (i.e., either continuation with a new set of parameters or exit).
Constraints were implemented by means of parameter transformations. More specifically, to ensure that a certain parameter a is always greater than a lower limit a _{0}, the optimization procedure internally worked with a parameter α, assuming that a=a _{0}+α ^{2}. This way the model parameters in Eq. (18) were constrained as follows: r _{ k }≥−1 and τ _{ k }≥1 ms (1≤k≤K), and t _{ k }≥t _{ k−1} (2≤k≤K).
Declarations
Acknowledgements
This work was supported by the Deutsche Forschungsgemeinschaft (grant LU 342/12-1) and the Open Access Publication Fund of the University of Munster.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Rosengren SM, Welgampola MS, Colebatch JG. Vestibular evoked myogenic potentials: past, present and future. Clin Neurophysiol. 2010; 121:636–51.View ArticlePubMedGoogle Scholar
- Curthoys IS. A critical review of the neurophysiological evidence underlying clinical vestibular testing using sound, vibration and galvanic stimuli. Clin Neurophysiol. 2010; 121:132–44.View ArticlePubMedGoogle Scholar
- Rosengren SM, Kingma H. New perspectives on vestibular evoked myogenic potentials. Curr Opin Neurol. 2013; 26:74–80.View ArticlePubMedGoogle Scholar
- Wit HP, Kingma CM. A simple model for the generation of the vestibular evoked myogenic potential (VEMP). Clin Neurophysiol. 2006; 117:1354–8.View ArticlePubMedGoogle Scholar
- Lütkenhöner B, Stoll W, Basel T. Modeling the vestibular evoked myogenic potential. J Theor Biol. 2010; 263:70–8.View ArticlePubMedGoogle Scholar
- Lütkenhöner B, Basel T. Deconvolution of the vestibular evoked myogenic potential. J Theor Biol. 2012; 294:87–97.View ArticlePubMedGoogle Scholar
- Shahid S, Walker J, Lyons GM, Byrne CA, Nene AV. Application of higher order statistics techniques to EMG signals to characterize the motor unit action potential. IEEE Trans Biomed Eng. 2005; 52(7):1195–209.View ArticlePubMedGoogle Scholar
- Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical recipes. The art of scientific computing, 3rd edn. New York: Cambridge University Press; 2007.Google Scholar
- Brigham EO. The Fast Fourier Transform. Englewood Cliffs, NJ: Prentice-Hall; 1974.Google Scholar
- Lütkenhöner B, Basel T. An analytical model of the vestibular evoked myogenic potential. J Theor Biol. 2011; 286:41–9.View ArticlePubMedGoogle Scholar
- Colebatch JG, Halmagyi GM, Skuse NF. Myogenic potentials generated by a click-evoked vestibulocollic reflex. J Neurol Neurosurg Psychiatry. 1994; 57:190–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Colebatch JG, Rothwell JC. Motor unit excitability changes mediating vestibulocollic reflexes in the sternocleidomastoid muscle. Clin Neurophysiol. 2004; 115:2567–73.View ArticlePubMedGoogle Scholar
- Wu CC, Young YH. Vestibular evoked myogenic potentials are intact after sudden deafness. Ear Hear. 2002; 23:235–8.View ArticlePubMedGoogle Scholar
- Eleftheriadou A, Deftereos S, Zarikas V, Panagopoulos G, Sfetsos S, Karageorgiou K, et al. Test-retest reliability of p13n23 and n34p44 components of vestibular evoked myogenic potentials in a large healthy population. J Otolaryngol Head Neck Surg. 2009; 38:462–7.PubMedGoogle Scholar
- Cox DR, Miller HD. The theory of stochastic processes. London: Chapman and Hall; 1965.Google Scholar
- Lütkenhöner B, Rudack C, Basel T. The variance modulation associated with the vestibular evoked myogenic potential. Clin Neurophysiol. 2011; 122:1448–56.View ArticlePubMedGoogle Scholar
- Lagarias J, Reeds J, Wright M, Wright P. Convergence properties of the Nelder–Mead simplex method in low dimensions. SIAM J Optim. 1998; 9:112–47.View ArticleGoogle Scholar