A Bayesian change-point analysis of electromyographic data: detecting muscle activation patterns and associated applications.

Many facets of neuromuscular activation patterns and control can be assessed via electromyography and are important for understanding the control of locomotion. After spinal cord injury, muscle activation patterns can affect locomotor recovery. We present a novel application of reversible jump Markov chain Monte Carlo simulation to estimate activation patterns from electromyographic data. We assume the data to be a zero-mean, heteroscedastic process. The variance is explicitly modeled using a step function. The number and location of points of discontinuity, or change-points, in the step function, the inter-change-point variances, and the overall mean are jointly modeled along with the mean and variance from baseline data. The number of change-points is considered a nuisance parameter and is integrated out of the posterior distribution. Whereas current methods of detecting activation patterns are deterministic or provide only point estimates, ours provides distributional estimates of muscle activation. These estimates, in turn, are used to estimate physiologically relevant quantities such as muscle coactivity, total integrated energy, and average burst duration and to draw valid statistical inferences about these quantities.

[1]  W. Snodgrass.  Physiology , 1897, Nature.

[2]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[3]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[4]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[5]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[6]  R. Engle.  Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation , 1982 .

[7]  T. Bollerslev,et al.  Generalized autoregressive conditional heteroskedasticity , 1986 .

[8]  G. Loeb,et al.  Electromyography for Experimentalists , 1986 .

[9]  P. Thompson.  Electromyography for Experimentalists , 1987 .

[10]  J. Todd.  Book Review: Digital image processing (second edition). By R. C. Gonzalez and P. Wintz, Addison-Wesley, 1987. 503 pp. Price: £29.95. (ISBN 0-201-11026-1) , 1988 .

[11]  G. Wahba.  Spline models for observational data , 1990 .

[12]  William K. Pratt,et al.  Digital image processing, 2nd Edition , 1991, A Wiley-Interscience publication.

[13]  J. Perry,et al.  Gait Analysis , 2024 .

[14]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[15]  R. Chou,et al.  ARCH modeling in finance: A review of the theory and empirical evidence , 1992 .

[16]  S. Simon.  Gait Analysis, Normal and Pathological Function. , 1993 .

[17]  A comparison of methods that characterize pulses in a time series. , 1995, Statistics in medicine.

[18]  P. Green.  Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[19]  P. Hodges,et al.  A comparison of computer-based methods for the determination of onset of muscle contraction using electromyography. , 1996, Electroencephalography and clinical neurophysiology.

[20]  R. Neptune,et al.  The effect of pedaling rate on coordination in cycling. , 1997, Journal of biomechanics.

[21]  R. Eubank.  Nonparametric Regression and Spline Smoothing , 1999 .

[22]  G Staude,et al.  Objective motor response onset detection in surface myoelectric signals. , 1999, Medical engineering & physics.

[23]  Saad T. Bakir,et al.  Nonparametric Regression and Spline Smoothing , 2000, Technometrics.

[24]  Susan J. Harkema,et al.  Clonus after human spinal cord injury cannot be attributed solely to recurrent muscle-tendon stretch , 2003, Experimental Brain Research.