Given a sequence of states arising from a stationary state, it fits the underlying Markov chain distribution with higher order.

1 2 3 |

`sequence` |
A character list. |

`order` |
Markov chain order |

A list containing lambda, Q, and X.

This function is written in Rcpp.

Giorgio Spedicato, Tae Seung Kang

Ching, W. K., Huang, X., Ng, M. K., & Siu, T. K. (2013). Higher-order markov chains. In Markov Chains (pp. 141-176). Springer US.

Ching, W. K., Ng, M. K., & Fung, E. S. (2008). Higher-order multivariate Markov chains and their applications. Linear Algebra and its Applications, 428(2), 492-507.

package version 0.2.5

`markovchain`

1 2 3 | ```
sequence<-c("a", "a", "b", "b", "a", "c", "b", "a", "b", "c", "a", "b", "c", "a"
, "b", "c", "a", "b", "a", "b")
fitHigherOrder(sequence)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.