Abstract

This paper describes a Bayesian approach to determining the order of a finite state Markov chain whose transition probabilities are themselves governed by a homogeneous finite state Markov chain. It extends previous work on homogeneous Markov chains to more general and applicable hidden Markov models. The method we describe uses a Markov chain Monte Carlo algorithm to obtain samples from the (posterior) distribution for both the order of Markov dependence in the observed sequence and the other governing model parameters. These samples allow coherent inferences to be made straightforwardly in contrast to those which use information criteria. The methods are illustrated by their application to both simulated and real data sets.

Highlights

  • Markov chains are commonly used as models for data which are observed in discrete time and have a discrete and finite state space

  • The Markov chain Monte Carlo (MCMC) algorithm was run for 110000 iterations with the first 10000 being discarded as burn-in

  • The MCMC algorithm was again run for 110000 iterations with the first 10000 being discarded as burnin

Read more

Summary

Introduction

Markov chains are commonly used as models for data which are observed in discrete time and have a discrete and finite state space. The ease of their approach is a direct consequence of their choice of prior distribution for the transition probabilities governing the evolution of the underlying process; we will return to this point later Their method can be extended to the more general HMM context (r > 1) if the configuration of hidden states s is known. A fundamental drawback of using HMMs is that the configuration s is unknown and has to be determined from the observed data y This complication precludes a fully analytic treatment of the model and so we resort to computer intensive Markov chain Monte Carlo (MCMC) methods.

Model description
Bayesian inference
Posterior inference
Posterior summaries
Simulated data
DNA sequence data
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call