title: MCMC can detect nonidentifiable models.
abstract: Continuous-time Markov models have been considered the best representation for the stochastic dynamics of ion channels for more than thirty years. For most single-channel data sets, several open and closed states are required for accurately representing the dynamics. However, each data point only shows if the channel is open or closed but not in which state it is. Consequently, some model structures are inherently overparameterized and therefore, in principle, unsuitable for representing any data--those models are called "nonidentifiable". As of this writing, it seems to be poorly understood which continuous-time Markov models are identifiable and which are not, therefore the unconscious use of a nonidentifiable model is a considerable concern. To address this problem, an improved variant of a recently published Markov-chain Monte Carlo method is presented. The algorithm is tested using test data as well as experimental data. We demonstrate that, opposed to a widely used maximum-likelihood estimator, it gives clear warning signs when a nonidentifiable model is used for fitting. Furthermore, for test data that was generated from a nonidentifiable model, the Markov-chain Monte Carlo results recover much more information from the data than maximum-likelihood estimation.
authors.affiliation: Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand. ivo.siekmann@auckland.ac.nz
authors.forename: Ivo
authors.initials: I
authors.lastname: Siekmann
authors.forename: James
authors.initials: J
authors.lastname: Sneyd
authors.forename: Edmund J
authors.initials: EJ
authors.lastname: Crampin
mesh_terms: D000465:Algorithms; D000818:Animals; D003198:Computer Simulation; D006801:Humans; D015640:Ion Channel Gating / Q000502:physiology*; D007473:Ion Channels / Q000502:physiology*; D008390:Markov Chains; D008954:Models, Biological; D015233:Models, Statistical; D009010:Monte Carlo Method
chemical_list: D007473:Ion Channels