Maximum likelihood estimation of motor unit firing pattern statistics

Exact inter-discharge interval distribution of motor unit firing patterns with gamma model

Inter-discharge interval distribution modeling of the motor unit firing pattern plays an important role in electromyographic decomposition and the statistical analysis of firing patterns. When modeling firing patterns obtained from automatic procedures, false positives and false negatives can be taken into account to enhance performance in estimating firing pattern statistics. Available models of this type, however, are only approximate and use Gaussian distributions, which are not strictly suitable for modeling renewal point processes. In this paper, the theory of point processes is used to derive an exact solution to the distribution when a gamma distribution is used to model the physiological firing pattern. Besides being exact, the solution provides a way to model the skewness of the inter-discharge distribution, and this may make it possible to obtain a better fit with available experimental data. In order to demonstrate potential applications of the model, we use it to obtain a maximum likelihood estimator of firing pattern statistics. Our tests found this estimator to be reliable over a wide range of firing conditions, whether dealing with real or simulated firing patterns, the proposed solution had better agreement than other models.

Graphical Abstract

Comparison of methods for estimating motor unit firing rate time series from firing times

The central nervous system regulates recruitment and firing of motor units to modulate muscle tension. Estimation of the firing rate time series is typically performed by decomposing the electromyogram (EMG) into its constituent firing times, then lowpass filtering a constituent train of impulses. Little research has examined the performance of different estimation methods, particularly in the inevitable presence of decomposition errors. The study of electrocardiogram (ECG) and electroneurogram (ENG) firing rate time series presents a similar problem, and has applied novel simulation models and firing rate estimators. Herein, we adapted an ENG/ECG simulation model to generate realistic EMG firing times derived from known rates, and assessed various firing rate time series estimation methods. ENG/ECG-inspired rate estimation worked exceptionally well when EMG decomposition errors were absent, but degraded unacceptably with decomposition error rates of ⩾1%. Typical EMG decomposition error rates-even after expert manual review-are 3-5%. At realistic decomposition error rates, more traditional EMG smoothing approaches performed best, when optimal smoothing window durations were selected. This optimal window was often longer than the 400ms duration that is commonly used in the literature. The optimal duration decreased as the modulation frequency of firing rate increased, average firing rate increased and decomposition errors decreased. Examples of these rate estimation methods on physiologic data are also provided, demonstrating their influence on measures computed from the firing rate estimate.

Inter-discharge interval distribution of motor unit firing patterns with detection errors

Inter-discharge interval (IDI) distribution analysis of motor unit firing patterns is a valuable tool in EMG decomposition and analysis. However, the firing pattern obtained by EMG decomposition may have detection errors: false positives (incorrectly classified firings) and false negatives (missed firings). In this paper, the mathematical derivation of an IDI distribution model that accommodates false positives and false negatives of the detection process is presented. An approximation of the general model to adapt to specific EMG decomposition conditions is also presented. To illustrate the usefulness of the model, the obtained distribution is used to derive the maximum likelihood estimates of the statistics of motor unit firing patterns, the IDI mean and standard deviation, and estimates of the false negative and false positive ratios. Results obtained from simulation experiments and tests with real motor unit firing patterns show an enhanced estimation performance when compared to previously available algorithms. Goodness-of-fit tests applied to estimations for real data corrupted with false positives showed that the model-driven estimations fitted the uncorrupted data better than EFE estimations: 82% versus 52% not rejectable, respectively, when false positives were about 10% of IDIs. With about 5% false positives, the not rejectable estimations were 85% versus 70%.