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