Abstract
Acupuncturing the ST36 acupoint can evoke a responding activity in the spinal dorsal root ganglia and generate spikes. In order to identify the responding mechanism of different acupuncture manipulations, in this paper the spike history of neurons is taken as the starting point and the coupling generalized linear model is adopted to encode the neuronal spiking activity evoked by different acupuncture manipulations. Then, maximum likelihood estimation is used to fit the model parameters and estimate the coupling parameters of stimulus, the self-coupling parameters of the neuron's own spike history and the cross-coupling parameters of other neurons' spike history. We use simulation data to test the estimation algorithm's effectiveness and analyze the main factors that evoke neuronal responding activity. Finally, we use the coupling generalized linear model to encode neuronal spiking activity evoked by two acupuncture manipulations, and estimate the coupling parameters of stimulus, the self-coupling parameters and the cross-coupling parameters. The results show that in acupuncture experiments, acupuncture stimulus is the inducing factor of neuronal spiking activity, and the cross-coupling of other neurons' spike history is the main factor of neuronal spiking activity. Additionally, the higher the amplitude of the neuronal spiking waveform, the greater the cross-coupling parameter. This lays a theoretical foundation for the scientific application of acupuncture therapy.