System identification of the nociceptive function

Research output: Chapter in Book/Report/Conference proceedingConference contributionAcademic

Abstract

Pain is essential for survival but persistent pain reduces quality of life. One clinical symptom is hyperalgesia which results in either a decrease in threshold or an increase in suprathreshold response\cite{Sandkühler}. Hyperalgesia can be a consequence of peripheral and/or central sensitization. This perturbs the nociceptive system that is then malfunctioning. Our goal is to identify the nociceptive pathway and to characterize which parts are (mal)functioning. To this end, we adopt a system identification approach by integrating input-output measurements for estimating the parameters in a mathematical model of the ascending nociceptive pathway. We consider data from a (binary) detection task with suitable variation in stimulus parameters. This is achieved by application of electrocutaneous stimulation using an intra-epidermal needle electrode. In each trial, current is delivered as a pulse train characterized by the amplitude and the temporal combination of the number of pulses, the inter-pulse interval and the pulse width. Then, the subject reports a binary response, perceived or not. From one trial to the next, the amplitude of the current is be adjusted using a staircase procedure \cite{Doll}. The set of stimulus-response pairs (SRP) is the input-output data for our identification method. In a previous study, we constructed a mathematical model of the ascending nociceptive pathway considering peripheral and central mechanisms. Analogous to the experimental paradigm, one can define the perception threshold (PT) for the model: the amplitude such that given the experimental temporal parameters and model parameters, half of the stimuli are detected \cite{Yang}. The model reproduces PT variations similar as in experiments and depends on the parameters of the model. Based on this we were able to estimate the model parameters given the PTs. We found, however, that the addition of (measurement) noise in the PTs led to large estimation errors. Indeed, a relative error of only 0.1\% in the PTs resulted in a 10\% relative error in the estimation. We suspect this is due to losing information when transforming the SRP into a set of PTs. Therefore, we aim to improve the estimation by using SRP directly. The staircase procedures can be formulated as Markov models \cite{Wetherill}. This describes the eventual distribution of the amplitude of applied currents. For the staircase procedures we consider, the Markov models are ergodic, hence, a single realization or measurement is informative if the number of trials is large enough. Based on the final asymptotic distribution and the SRP, we formulate a maximum likelihood estimator for identification of the model parameters. We show that on simulated data this direct approach with SRP outperforms the indirect method with PTs. We study the effect of the number of trials and the choice of stimulus parameters, i.e. the pulse width and the inter-pulse interval, on the estimation performance.
Original languageUndefined
Title of host publicationProceedings of the 29th Annual Meeting of the International Society for Psychophysics
EditorsJ. Wackermann, M. Wittman, W. Skrandies
Place of PublicationFreiburg
PublisherInternational Society for Psychophysics
Pages38
Number of pages1
ISBN (Print)not assigned
Publication statusPublished - Oct 2013
Event29th Annual Meeting of the International Society for Psychophysics: Fechner Day 2013 - Hotel Rheingold, Freiburg im Breisgau, Germany
Duration: 21 Oct 201325 Oct 2013

Publication series

Name
PublisherInternational Society for Psychophysics

Conference

Conference29th Annual Meeting of the International Society for Psychophysics
CountryGermany
CityFreiburg im Breisgau
Period21/10/1325/10/13

Keywords

  • EWI-23902
  • IR-87908
  • METIS-302552

Cite this

Yang, H., Buitenweg, J. R., & Meijer, H. G. E. (2013). System identification of the nociceptive function. In J. Wackermann, M. Wittman, & W. Skrandies (Eds.), Proceedings of the 29th Annual Meeting of the International Society for Psychophysics (pp. 38). Freiburg: International Society for Psychophysics.
Yang, H. ; Buitenweg, Jan R. ; Meijer, Hil Gaétan Ellart. / System identification of the nociceptive function. Proceedings of the 29th Annual Meeting of the International Society for Psychophysics. editor / J. Wackermann ; M. Wittman ; W. Skrandies. Freiburg : International Society for Psychophysics, 2013. pp. 38
@inproceedings{7f1fa8c2cb1345a8ae05b28dfe171d67,
title = "System identification of the nociceptive function",
abstract = "Pain is essential for survival but persistent pain reduces quality of life. One clinical symptom is hyperalgesia which results in either a decrease in threshold or an increase in suprathreshold response\cite{Sandk{\"u}hler}. Hyperalgesia can be a consequence of peripheral and/or central sensitization. This perturbs the nociceptive system that is then malfunctioning. Our goal is to identify the nociceptive pathway and to characterize which parts are (mal)functioning. To this end, we adopt a system identification approach by integrating input-output measurements for estimating the parameters in a mathematical model of the ascending nociceptive pathway. We consider data from a (binary) detection task with suitable variation in stimulus parameters. This is achieved by application of electrocutaneous stimulation using an intra-epidermal needle electrode. In each trial, current is delivered as a pulse train characterized by the amplitude and the temporal combination of the number of pulses, the inter-pulse interval and the pulse width. Then, the subject reports a binary response, perceived or not. From one trial to the next, the amplitude of the current is be adjusted using a staircase procedure \cite{Doll}. The set of stimulus-response pairs (SRP) is the input-output data for our identification method. In a previous study, we constructed a mathematical model of the ascending nociceptive pathway considering peripheral and central mechanisms. Analogous to the experimental paradigm, one can define the perception threshold (PT) for the model: the amplitude such that given the experimental temporal parameters and model parameters, half of the stimuli are detected \cite{Yang}. The model reproduces PT variations similar as in experiments and depends on the parameters of the model. Based on this we were able to estimate the model parameters given the PTs. We found, however, that the addition of (measurement) noise in the PTs led to large estimation errors. Indeed, a relative error of only 0.1\{\%} in the PTs resulted in a 10\{\%} relative error in the estimation. We suspect this is due to losing information when transforming the SRP into a set of PTs. Therefore, we aim to improve the estimation by using SRP directly. The staircase procedures can be formulated as Markov models \cite{Wetherill}. This describes the eventual distribution of the amplitude of applied currents. For the staircase procedures we consider, the Markov models are ergodic, hence, a single realization or measurement is informative if the number of trials is large enough. Based on the final asymptotic distribution and the SRP, we formulate a maximum likelihood estimator for identification of the model parameters. We show that on simulated data this direct approach with SRP outperforms the indirect method with PTs. We study the effect of the number of trials and the choice of stimulus parameters, i.e. the pulse width and the inter-pulse interval, on the estimation performance.",
keywords = "EWI-23902, IR-87908, METIS-302552",
author = "H. Yang and Buitenweg, {Jan R.} and Meijer, {Hil Ga{\'e}tan Ellart}",
year = "2013",
month = "10",
language = "Undefined",
isbn = "not assigned",
publisher = "International Society for Psychophysics",
pages = "38",
editor = "J. Wackermann and M. Wittman and W. Skrandies",
booktitle = "Proceedings of the 29th Annual Meeting of the International Society for Psychophysics",

}

Yang, H, Buitenweg, JR & Meijer, HGE 2013, System identification of the nociceptive function. in J Wackermann, M Wittman & W Skrandies (eds), Proceedings of the 29th Annual Meeting of the International Society for Psychophysics. International Society for Psychophysics, Freiburg, pp. 38, 29th Annual Meeting of the International Society for Psychophysics, Freiburg im Breisgau, Germany, 21/10/13.

System identification of the nociceptive function. / Yang, H.; Buitenweg, Jan R.; Meijer, Hil Gaétan Ellart.

Proceedings of the 29th Annual Meeting of the International Society for Psychophysics. ed. / J. Wackermann; M. Wittman; W. Skrandies. Freiburg : International Society for Psychophysics, 2013. p. 38.

Research output: Chapter in Book/Report/Conference proceedingConference contributionAcademic

TY - GEN

T1 - System identification of the nociceptive function

AU - Yang, H.

AU - Buitenweg, Jan R.

AU - Meijer, Hil Gaétan Ellart

PY - 2013/10

Y1 - 2013/10

N2 - Pain is essential for survival but persistent pain reduces quality of life. One clinical symptom is hyperalgesia which results in either a decrease in threshold or an increase in suprathreshold response\cite{Sandkühler}. Hyperalgesia can be a consequence of peripheral and/or central sensitization. This perturbs the nociceptive system that is then malfunctioning. Our goal is to identify the nociceptive pathway and to characterize which parts are (mal)functioning. To this end, we adopt a system identification approach by integrating input-output measurements for estimating the parameters in a mathematical model of the ascending nociceptive pathway. We consider data from a (binary) detection task with suitable variation in stimulus parameters. This is achieved by application of electrocutaneous stimulation using an intra-epidermal needle electrode. In each trial, current is delivered as a pulse train characterized by the amplitude and the temporal combination of the number of pulses, the inter-pulse interval and the pulse width. Then, the subject reports a binary response, perceived or not. From one trial to the next, the amplitude of the current is be adjusted using a staircase procedure \cite{Doll}. The set of stimulus-response pairs (SRP) is the input-output data for our identification method. In a previous study, we constructed a mathematical model of the ascending nociceptive pathway considering peripheral and central mechanisms. Analogous to the experimental paradigm, one can define the perception threshold (PT) for the model: the amplitude such that given the experimental temporal parameters and model parameters, half of the stimuli are detected \cite{Yang}. The model reproduces PT variations similar as in experiments and depends on the parameters of the model. Based on this we were able to estimate the model parameters given the PTs. We found, however, that the addition of (measurement) noise in the PTs led to large estimation errors. Indeed, a relative error of only 0.1\% in the PTs resulted in a 10\% relative error in the estimation. We suspect this is due to losing information when transforming the SRP into a set of PTs. Therefore, we aim to improve the estimation by using SRP directly. The staircase procedures can be formulated as Markov models \cite{Wetherill}. This describes the eventual distribution of the amplitude of applied currents. For the staircase procedures we consider, the Markov models are ergodic, hence, a single realization or measurement is informative if the number of trials is large enough. Based on the final asymptotic distribution and the SRP, we formulate a maximum likelihood estimator for identification of the model parameters. We show that on simulated data this direct approach with SRP outperforms the indirect method with PTs. We study the effect of the number of trials and the choice of stimulus parameters, i.e. the pulse width and the inter-pulse interval, on the estimation performance.

AB - Pain is essential for survival but persistent pain reduces quality of life. One clinical symptom is hyperalgesia which results in either a decrease in threshold or an increase in suprathreshold response\cite{Sandkühler}. Hyperalgesia can be a consequence of peripheral and/or central sensitization. This perturbs the nociceptive system that is then malfunctioning. Our goal is to identify the nociceptive pathway and to characterize which parts are (mal)functioning. To this end, we adopt a system identification approach by integrating input-output measurements for estimating the parameters in a mathematical model of the ascending nociceptive pathway. We consider data from a (binary) detection task with suitable variation in stimulus parameters. This is achieved by application of electrocutaneous stimulation using an intra-epidermal needle electrode. In each trial, current is delivered as a pulse train characterized by the amplitude and the temporal combination of the number of pulses, the inter-pulse interval and the pulse width. Then, the subject reports a binary response, perceived or not. From one trial to the next, the amplitude of the current is be adjusted using a staircase procedure \cite{Doll}. The set of stimulus-response pairs (SRP) is the input-output data for our identification method. In a previous study, we constructed a mathematical model of the ascending nociceptive pathway considering peripheral and central mechanisms. Analogous to the experimental paradigm, one can define the perception threshold (PT) for the model: the amplitude such that given the experimental temporal parameters and model parameters, half of the stimuli are detected \cite{Yang}. The model reproduces PT variations similar as in experiments and depends on the parameters of the model. Based on this we were able to estimate the model parameters given the PTs. We found, however, that the addition of (measurement) noise in the PTs led to large estimation errors. Indeed, a relative error of only 0.1\% in the PTs resulted in a 10\% relative error in the estimation. We suspect this is due to losing information when transforming the SRP into a set of PTs. Therefore, we aim to improve the estimation by using SRP directly. The staircase procedures can be formulated as Markov models \cite{Wetherill}. This describes the eventual distribution of the amplitude of applied currents. For the staircase procedures we consider, the Markov models are ergodic, hence, a single realization or measurement is informative if the number of trials is large enough. Based on the final asymptotic distribution and the SRP, we formulate a maximum likelihood estimator for identification of the model parameters. We show that on simulated data this direct approach with SRP outperforms the indirect method with PTs. We study the effect of the number of trials and the choice of stimulus parameters, i.e. the pulse width and the inter-pulse interval, on the estimation performance.

KW - EWI-23902

KW - IR-87908

KW - METIS-302552

M3 - Conference contribution

SN - not assigned

SP - 38

BT - Proceedings of the 29th Annual Meeting of the International Society for Psychophysics

A2 - Wackermann, J.

A2 - Wittman, M.

A2 - Skrandies, W.

PB - International Society for Psychophysics

CY - Freiburg

ER -

Yang H, Buitenweg JR, Meijer HGE. System identification of the nociceptive function. In Wackermann J, Wittman M, Skrandies W, editors, Proceedings of the 29th Annual Meeting of the International Society for Psychophysics. Freiburg: International Society for Psychophysics. 2013. p. 38