Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
173 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
46 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Inference for the stochastic FitzHugh-Nagumo model from real action potential data via approximate Bayesian computation (2405.17972v2)

Published 28 May 2024 in stat.CO and math.DS

Abstract: The stochastic FitzHugh-Nagumo (FHN) model is a two-dimensional nonlinear stochastic differential equation with additive degenerate noise, whose first component, the only one observed, describes the membrane voltage evolution of a single neuron. Due to its low-dimensionality, its analytical and numerical tractability and its neuronal interpretation, it has been used as a case study to test the performance of different statistical methods in estimating the underlying model parameters. Existing methods, however, often require complete observations, non-degeneracy of the noise or a complex architecture (e.g., to estimate the transition density of the process, "recovering" the unobserved second component) and they may not (satisfactorily) estimate all model parameters simultaneously. Moreover, these studies lack real data applications for the stochastic FHN model. The proposed method tackles all challenges (non-globally Lipschitz drift, non-explicit solution, lack of available transition density, degeneracy of the noise and partial observations). It is an intuitive and easy-to-implement sequential Monte Carlo approximate Bayesian computation algorithm, which relies on a recent computationally efficient and structure-preserving numerical splitting scheme for synthetic data generation and on summary statistics exploiting the structural properties of the process. All model parameters are successfully estimated from simulated data and, more remarkably, real action potential data of rats. The presented novel real-data fit may broaden the scope and credibility of this classic and widely used neuronal model.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (57)
  1. A stochastic version of the Jansen and Rit neural mass model: analysis and numerics. J. Math. Neurosci., 7(8), 2017.
  2. Approximate Bayesian Computation with the Wasserstein Distance. J. Royal Stat. Soc. B, 81(2):235–269, 2019.
  3. S. Bonaccorsi and E. Mastrogiacomo. Analysis of the stochastic FitzHugh-Nagumo system. Infin. Dimens. Anal. Quantum Probab. Relat. Top., 11(03):427–446, 2008.
  4. Splitting schemes for FitzHugh–Nagumo stochastic partial differential equations. Discrete Cont Dyn-B, 29(1):214–244, 2024.
  5. A splitting method for SDEs with locally Lipschitz drift: Illustration on the FitzHugh-Nagumo model. Appl. Numer. Math., 179:191–220, 2022.
  6. Spectral density-based and measure-preserving ABC for partially observed diffusion processes. An illustration on Hamiltonian SDEs. Stat. Comput., 30(3):627–648, 2020.
  7. Q. Clairon and A. Samson. Optimal control for parameter estimation in partially observed hypoelliptic stochastic differential equations. Comput. Stat., 37:2471–2491, 2022.
  8. An adaptive sequential Monte Carlo method for approximate Bayesian computation. Stat. Comput., 22:1009–1020, 2012.
  9. S. Ditlevsen and A. Samson. Hypoelliptic diffusions: filtering and inference from complete and partial observations. J. Royal Stat. Soc. B, 81(2):361–384, 2019.
  10. Network inference in a stochastic multi-population neural mass model via approximate Bayesian computation. arXiv, 2306.15787v2, 2023.
  11. R. O. Doruk and L. Abosharb. Estimating the Parameters of Fitzhugh-Nagumo Neurons from Neural Spiking Data. Brain Sci., 9:364, 2019.
  12. Exact and approximate Bayesian inference for low integer-valued time series models with intractable liklihoods. Bayesian Anal., 11:325–352, 2016.
  13. D. Eddelbuettel and R. François. Rcpp: seamless R and C++ integration. J. Stat. Soft., 40(8):1–18, 2011.
  14. M. Fasiolo. An introduction to mvnfast. R package version 0.2.8. University of Bristol, 2014.
  15. On optimality of kernels for approximate Bayesian computation using sequential Monte Carlo. Stat. Appl. Genet. Mol. Biol., 12(1):87–107, 2013.
  16. R. FitzHugh. Impulses and physiological states in theoretical models of nerve membrane. Biophys. J., 1(6):445–466, 1961.
  17. A. Gloter and N. Yoshida. Adaptive estimation for degenerate diffusion processes. Elec. J. Stat., 15:1424–1472, 2021.
  18. Manifold markov chain monte carlo methods for bayesian inference in diffusion models. J. Roy. Stat. Soc. B, 84:1229–1256, 2022.
  19. Bayesian inference for SDE models: a case study for an excitable stochastic-dynamical model. NextJournal, 2020.
  20. The partially truncated Euler–Maruyama method and its stability and boundedness. Appl. Numer. Math., 115:235–251, 2017.
  21. A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol., 117(4):500–544, 1952.
  22. M. Hutzenthaler and A. Jentzen. Numerical approximations of stochastic differential equations with non-globally Lipschitz continuous coefficients. Mem. Am. Math. Soc., 236, 2012.
  23. Strong and weak divergence in finite time of Euler’s method for stochastic differential equations with non-globally Lipschitz continuous coefficients. Proc. R. Soc. A, 467:1563–1576, 2011.
  24. Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients. Ann. Appl. Probab., 22(4):1611–1641, 2012.
  25. Y. Iguchi and A. Beskos. Parameter inference for hypo-elliptic diffusions under a weak design condition. submitted, 2023.
  26. Parameter estimation with increased precision for elliptic and hypo-elliptic diffusions. submitted, 2022.
  27. A. Jasra. Approximate Bayesian computation for a class of time series models. Int Stat. Rev., 83:405–435, 2015.
  28. Towards data-conditional simulation for abc inference in stochastic differential equations. arXiv, preprint, 2023.
  29. R. Khasminskii and F. C. Klebaner. Long term behavior of solutions of the Lotka-Volterra system under small random perturbations. Ann. Appl. Probab., 11(3):952–963, 2001.
  30. A tutorial introduction to Bayesian inference for stochastic epidemic models using approximate Bayesian computation. Math Biosci., 287:42–53, 2017.
  31. J.R. León and A. Samson. Hypoelliptic stochastic FitzHugh-Nagumo neuronal model: mixing, up-crossing and estimation of the spike rate. Ann. Appl. Probab., 28(4):2243–2274, 2018.
  32. Stochastic resonance in an underdamped system with FitzHug-Nagumo potential for weak signal detection. Journal of Sound and Vibration, 411:34–46, 2017.
  33. X. Mao. The truncated Euler–Maruyama method for stochastic differential equations. J. Comput. Appl. Math., 290:370–384, 2015.
  34. X. Mao. Convergence rates of the truncated Euler–Maruyama method for stochastic differential equations. J. Comput. Appl. Math., 296:362–375, 2016.
  35. Approximate Bayesian computational methods. Stat Comput, 22:1167–1180, 2012.
  36. Auxiliary likelihood-based approximate Bayesian computation in state space models. J Comput. Graph. Stat., 28:1–31, 2019.
  37. Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise. Stoch. Process. Their Appl., 101(2):185–232, 2002.
  38. Fast approximate Bayesian inference for stable differential equation models. submitted, 2017.
  39. A. Melnykova. Parametric inference for hypoelliptic ergodic diffusions with full observations. Stat. Inf. Stoc. Proc., 23:595–635, 2020.
  40. A dataset of action potentials recorded from the L5 dorsal rootlet of rat using a multiple electrode array. Data in Brief, 33:106561, 2020.
  41. C. Morris and H. Lecar. Voltage oscillations in the barnacle giant muscle fiber. Biophys. J., 35:193–213, 1981.
  42. An active pulse transmission line simulating nerve axon. Proc. IRE, 50(10):2061–2070, 1962.
  43. U. Picchini. Inference for SDE models via Approximate Bayesian Computation. J Comput. Graph. Stat., 23:1080–1100, 2014.
  44. U. Picchini and Forman J. L. Accelerating inference for diffusions observed with measurement error and large sample sizes using approximate Bayesian computation. J. Stat. Comput. Simul., 86(1):195–213, 2016.
  45. U. Picchini and A. Samson. Coupling stochastic EM and approximate Bayesian computation for parameter inference in state-space models. Comput. Stat., 33(1):179–212, 2018.
  46. U. Picchini and M. Tamborrino. Guided sequential ABC schemes for intractable Bayesian models. Accepted for publication in Bayesian Analysis. Preprint at ArXiv, 2206.12235, 2022.
  47. Parameter estimation for partially observed hypoelliptic diffusions. J. Royal Stat. Soc. B, 71(1):49–73, 2009.
  48. O. Pons. Functional Estimation for Density, Regression Models and Processes. World Scientific Publishing, Singapore, 2011.
  49. On the periodogram estimators of periods from interleaved sparse, noisy timing data. IEEE Stat. Signal Processing Workshop, pages 232–235, 2014.
  50. R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria, 2011.
  51. S. Sabanis. Euler approximations with varying coefficients: the case of superlinearly growing diffusion coefficients. Ann. Appl. Probab., 26(4):2083–2105, 2016.
  52. A. Samson and M. Thieullen. A contrast estimator for completely or partially observed hypoelliptic diffusion. Stoc. Proc. Appl., 122:2521–2552, 2012.
  53. Handbook of Approximate Bayesian Computation. Chapman & Hall/CRC Handbooks of Modern Statistical Methods. CRC Press, Taylor & Francis Group, 2018.
  54. Parameter inference and model selection in deterministic and stochastic dynamical models via approximate Bayesian computation: Modeling a wildlife epidemic. Environmetrics, 26(7):451–462, 2015.
  55. A. Tancredi. Approximate Bayesian inference for discretely observed continuous-time multi-state models. Biometrics, 75:966–977, 2019.
  56. M. V. Tretyakov and Z. Zhang. A fundamental mean-square convergence theorem for SDEs with locally Lipschitz coefficients and its applications. SIAM J. Numer. Anal., 51(6):3135–3162, 2013.
  57. Z. Zhang and H. Ma. Order–preserving strong schemes for SDEs with locally Lipschitz coefficients. Appl. Numer. Math., 112:1–16, 2017.

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now
X Twitter Logo Streamline Icon: https://streamlinehq.com

Tweets