Statistical Identification of Markov Chain on Trees

The theoretical study of continuous-time homogeneous Markov chains is usually based on a natural assumption of a known transition rate matrix (TRM). However, the TRM of a Markov chain in realistic systems might be unknown and might even need to be identified by partially observable data. Thus, an issue on how to identify the TRM of the underlying Markov chain by partially observable information is derived from the great significance in applications. That is what we call the statistical identification of Markov chain. The Markov chain inversion approach has been derived for basic Markov chains by partial observation at few states. In the current letter, a more extensive class of Markov chain on trees is investigated. Firstly, a type of a more operable derivative constraint is developed. Then, it is shown that all Markov chains on trees can be identified only by such derivative constraints of univariate distributions of sojourn time and/or hitting time at a few states. A numerical example is included to demonstrate the correctness of the proposed algorithms.

[1]  R Horn,et al.  Estimating kinetic constants from single channel data. , 1983, Biophysical journal.

[2]  Frank G. Ball,et al.  Marked Continuous-Time Markov Chain Modelling of Burst Behaviour for Single Ion Channels , 2007, Adv. Decis. Sci..

[3]  R J French,et al.  Software for acquisition and analysis of ion channel data: choices, tasks, and strategies. , 1992, Methods in enzymology.

[4]  Christopher Shelley,et al.  Linking Exponential Components to Kinetic States in Markov Models for Single-Channel Gating , 2008, The Journal of general physiology.

[5]  O Flomenbom,et al.  Toolbox for analyzing finite two-state trajectories. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Ophir Flomenbom,et al.  Utilizing the information content in two-state trajectories. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Bret Larget,et al.  A canonical representation for aggregated Markov processes , 1998, Journal of Applied Probability.

[8]  U. Seifert,et al.  Communications: Can one identify nonequilibrium in a three-state system by analyzing two-state trajectories? , 2010, The Journal of chemical physics.

[9]  Xiang Xu-yan,et al.  Markov chain inversion approach to identify the transition rates of ion channels , 2012 .

[10]  A. Hawkes,et al.  On the stochastic properties of bursts of single ion channel openings and of clusters of bursts. , 1982, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[11]  Jian-Qiang Feng,et al.  Identifying transition rates of ionic channels of star-graph branch type , 2006 .

[12]  J. Rice,et al.  On aggregated Markov processes , 1986 .

[13]  Nicholas P. Jewell,et al.  Mixtures of Exponential Distributions , 1982 .

[14]  Jens Timmer,et al.  Estimating transition rates in aggregated Markov models of ion channel gating with loops and with nearly equal dwell times , 1999, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[15]  Clive L. Keatinge MODELING LOSSES WITH THE MIXED EXPONENTIAL DISTRIBUTION , 2000 .

[16]  Minping Qian,et al.  Identifying transition rates of ionic channels via observations at a single state , 2003 .

[17]  P. Moran,et al.  Reversibility and Stochastic Networks , 1980 .

[18]  Chen Jia,et al.  A second perspective on the Amann–Schmiedl–Seifert criterion for non-equilibrium in a three-state system , 2015 .

[19]  A G Hawkes,et al.  The quality of maximum likelihood estimates of ion channel rate constants , 2003, The Journal of physiology.

[20]  Tzu-kʿun Wang,et al.  Birth and death processes and Markov chains , 1992 .